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

Chromosome-level assembly and annotation of the blue catfish Ictalurus furcatus, an aquaculture species for hybrid catfish reproduction, epigenetics, and heterosis studies

Background

The blue catfish is of great value in aquaculture and recreational fisheries. The F₁ hybrids of female channel catfish (Ictalurus punctatus) × male blue catfish (Ictalurusfurcatus) have been the primary driver of US catfish production in recent years because of superior growth, survival, and carcass yield. The channel–blue hybrid also provides an excellent model to investigate molecular mechanisms of environment-dependent heterosis. However, transcriptome and methylome studies suffered from low alignment rates to the channel catfish genome due to divergence, and the genome resources for blue catfish are not publicly available.

Results

The blue catfish genome assembly is 841.86 Mbp in length with excellent continuity (8.6 Mbp contig N50, 28.2 Mbp scaffold N50) and completeness (98.6% Eukaryota and 97.0% Actinopterygii BUSCO). A total of 30,971 protein-coding genes were predicted, of which 21,781 were supported by RNA sequencing evidence. Phylogenomic analyses revealed that it diverged from channel catfish approximately 9 million years ago with 15.7 million fixed nucleotide differences. The within-species single-nucleotide polymorphism (SNP) density is 0.32% between the most aquaculturally important blue catfish strains (D&B and Rio Grande). Gene family analysis discovered significant expansion of immune-related families in the blue catfish lineage, which may contribute to disease resistance in blue catfish.

Conclusions

We reported the first high-quality, chromosome-level assembly of the blue catfish genome, which provides the necessary genomic tool kit for transcriptome and methylome analysis, SNP discovery and marker-assisted selection, gene editing and genome engineering, and reproductive enhancement of the blue catfish and hybrid catfish.

Environment-Dependent Heterosis and Transgressive Gene Expression in Reciprocal Hybrids between the Channel Catfish Ictalurus punctatus and the Blue Catfish Ictalurus furcatus

The hybrid between female channel catfish (Ictalurus punctatus) and male blue catfish (Ictalurus furcatus) is superior in feed conversion, disease resistance, carcass yield, and harvestability compared to both parental species. However, heterosis and heterobeltiosis only occur in pond culture, and channel catfish grow much faster than the other genetic types in small culture units. This environment-dependent heterosis is intriguing, but the underlying genetic mechanisms are not well understood. In this study, phenotypic characterization and transcriptomic analyses were performed in the channel catfish, blue catfish, and their reciprocal F₁s reared in tanks. The results showed that the channel catfish is superior in growth-related morphometrics, presumably due to significantly lower innate immune function, as investigated by reduced lysozyme activity and alternative complement activity. RNA-seq analysis revealed that genes involved in fatty acid metabolism/transport are significantly upregulated in channel catfish compared to blue catfish and hybrids, which also contributes to the growth phenotype. Interestingly, hybrids have a 40-80% elevation in blood glucose than the parental species, which can be explained by a phenomenon called transgressive expression (overexpression/underexpression in F₁s than the parental species). A total of 1140 transgressive genes were identified in F₁ hybrids, indicating that 8.5% of the transcriptome displayed transgressive expression. Transgressive genes upregulated in F₁s are enriched for glycan degradation function, directly related to the increase in blood glucose level. This study is the first to explore molecular mechanisms of environment-dependent heterosis/heterobeltiosis in a vertebrate species and sheds light on the regulation and evolution of heterosis vs. hybrid incompatibility.

Exposure to environmentally relevant cadmium concentrations negatively impacts early life stages of channel catfish (Ictalurus punctatus)

Cadmium is a persistent contaminant of surface waters. The effects of cadmium on early life stages of fish are not well understood, although they are often disproportionately affected by contaminants. The objectives of this study were to examine effects of chronic exposure to environmentally relevant concentrations on growth, development, cellular stress, and glucose metabolism of channel catfish, Ictalurus punctatus. Eggs were wet-fertilized in treatment water at concentrations of 0.4 (control), 2.2 (low), or 8.5 (high) μg L^-1 and monitored through swim-up, black fry stage. Eggs and fry accumulated cadmium dose-dependently. Fertilization rates were unaffected, yet hatch rate was significantly reduced in the high treatment. Survival to black fry and overall size and condition factor were not affected; however, differences in yolk sac size, and presumably energetics of yolk fry, was detected. Physiological pathways were also affected, demonstrated by altered gene expression, most notably in genes related to carbohydrate metabolism. Elevated expression of HK and G6PD, rather than G6P and GADPH, suggests glucose may be shunted towards the pentose-phosphate pathway. Overall, observations indicate cadmium negatively affects development in early life stages of channel catfish, which could lead to shifts in population structure and life history patterns in exposed populations of wild fish.

Comparison of Channel Catfish and Blue Catfish Gut Microbiota Assemblages Shows Minimal Effects of Host Genetics on Microbial Structure and Inferred Function

The microbiota of teleost fish has gained a great deal of research attention within the past decade, with experiments suggesting that both host-genetics and environment are strong ecological forces shaping the bacterial assemblages of fish microbiomes. Despite representing great commercial and scientific importance, the catfish within the family Ictaluridae, specifically the blue and channel catfish, have received very little research attention directed toward their gut-associated microbiota using 16S rRNA gene sequencing. Within this study we utilize multiple genetically distinct strains of blue and channel catfish, verified via microsatellite genotyping, to further quantify the role of host-genetics in shaping the bacterial communities in the fish gut, while maintaining environmental and husbandry parameters constant. Comparisons of the gut microbiota among the two catfish species showed no differences in bacterial species richness (observed and Chao1) or overall composition (weighted and unweighted UniFrac) and UniFrac distances showed no correlation with host genetic distances (Rst) according to Mantel tests. The microbiota of environmental samples (diet and water) were found to be significantly more diverse than that of the catfish gut associated samples, suggesting that factors within the host were further regulating the bacterial communities, despite the lack of a clear connection between microbiota composition and host genotype. The catfish gut communities were dominated by the phyla Fusobacteria, Proteobacteria, and Firmicutes; however, differential abundance analysis between the two catfish species using analysis of composition of microbiomes detected two differential genera, Cetobacterium and Clostridium XI. The metagenomic pathway features inferred from our dataset suggests the catfish gut bacterial communities possess pathways beneficial to their host such as those involved in nutrient metabolism and antimicrobial biosynthesis, while also containing pathways involved in virulence factors of pathogens. Testing of the inferred KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways by DESeq2 revealed minor difference in microbiota function, with only two metagenomic pathways detected as differentially abundant between the two catfish species. As the first study to characterize the gut microbiota of blue catfish, our study results have direct implications on future ictalurid catfish research. Additionally, our insight into the intrinsic factors driving microbiota structure has basic implications for the future study of fish gut microbiota.

Characterization of two channel catfish, Ictalurus punctatus, glucocorticoid receptors and expression following an acute stressor

Two channel catfish glucocorticoid receptor genes, ipGR1 (NR3C1_1) and ipGR2 (NR3C1_2) were partially characterized. Identification and analysis of the genomic organization of two channel catfish glucocorticoid (GC) receptors (GRs) revealed differences in the lengths of exons 1 and 2 and the addition of an extra 27-bp exon inserted after exon 2 in the GR1 gene, yielding a 9-aa insert in the receptor protein. Sequence of the 9-aa insert in ipGR1 (WRARQNTHG) is unique compared to other teleost fish GRs. Amino acid sequence alignment of the two channel catfish GRs, revealed 55% sequence identity between them, with a high degree of sequence conservation (82%) in the DNA binding and ligand binding domains. Real-time PCR indicated that ipGR1 and ipGR2 were expressed in all tissues evaluated. Channel catfish GR1 was predominantly expressed in the gills, nearly 25-fold higher than in the liver. GR1 expression was higher than GR2 expression in gills, intestine, head kidney and heart (P<0.05). Channel catfish hepatic GR1 mRNA expression was significantly (P<0.05) increased from pre-stress expression 30min following removal of the acute stressor. After 30min of stress and during the 2h recovery period, ipGR1 mRNA expression was higher relative to ipGR2 expression. Hepatic ipGR2 expression was not affected (P>0.05) by the acute stress event. The present study adds to the growing body of information on GR evolution and function and further demonstrates the unique regulation of the GC/GR system in teleost fish.

Development of a 690 K SNP array in catfish and its application for genetic mapping and validation of the reference genome sequence

Single nucleotide polymorphisms (SNPs) are capable of providing the highest level of genome coverage for genomic and genetic analysis because of their abundance and relatively even distribution in the genome. Such a capacity, however, cannot be achieved without an efficient genotyping platform such as SNP arrays. In this work, we developed a high-density SNP array with 690,662 unique SNPs (herein 690 K array) that were relatively evenly distributed across the entire genome, and covered 98.6% of the reference genome sequence. Here we also report linkage mapping using the 690 K array, which allowed mapping of over 250,000 SNPs on the linkage map, the highest marker density among all the constructed linkage maps. These markers were mapped to 29 linkage groups (LGs) with 30,591 unique marker positions. This linkage map anchored 1,602 scaffolds of the reference genome sequence to LGs, accounting for over 97% of the total genome assembly. A total of 1,007 previously unmapped scaffolds were placed to LGs, allowing validation and in few instances correction of the reference genome sequence assembly. This linkage map should serve as a valuable resource for various genetic and genomic analyses, especially for GWAS and QTL mapping for genes associated with economically important traits.

Identification of differentially expressed genes in the spleens of polyriboinosinic polyribocytidylic acid (poly I:C)-stimulated yellow catfish Pelteobagrus fulvidraco

The yellow catfish, Pelteobagrus fulvidraco (Siluriformes: Bagridae) is an economically important fish in China. However, genomic research and resources on this species are largely unavailable and still in infancy. In the present study, we constructed a cDNA library following poly I:C injection to screen for immune response genes in the spleens of P. fulvidraco using suppression subtractive hybridization (SSH). A total of 420 putative expressed sequence tag (EST) clones were identified at 24 h post-injection, which contain 103 genes consisting of 25 immune response genes, 12 cytoskeleton genes, 7 cell cycle and apoptosis genes, 7 respiration and energy metabolism genes, 7 transport genes, 26 metabolism genes, 10 stress response genes, 9 translational regulation genes, and 71 unknown genes. Real-time quantitative reverse transcription-PCR (qRT-PCR) results revealed that a set of randomly selected immune response genes were identified to be up-regulated after 24 h of poly I:C stimulation compared to controls. Our study provides an annotation of immune genes in detail and insight into fish immunity.

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.

Genome-wide identification of Hsp70 genes in channel catfish and their regulated expression after bacterial infection

Heat shock proteins 70/110 (Hsp70/110) are a family of conserved ubiquitously expressed heat shock proteins which are produced by cells in response to exposure to stressful conditions. Besides the chaperone and housekeeping functions, they are also known to be involved in immune response during infection. In this study, we identified 16 Hsp70/110 geness in channel catfish (Ictalurus punctatus) through in silico analysis using RNA-Seq and genome databases. Among them 12 members of Hsp70 (Hspa) family and 4 members of Hsp110 (Hsph) family were identified. Phylogenetic and syntenic analyses provided strong evidence in supporting the orthologies of these HSPs. In addition, we also determined the expression patterns of Hsp70/110 genes after Flavobacterium columnare and Edwardsiella ictaluri infections by meta-analyses, for the first time in channel catfish. Ten out of sixteen genes were significantly up/down-regulated after bacterial challenges. Specifically, nine genes were found significantly expressed in gill after F. columnare infection. Two genes were found significantly expressed in intestine after E. ictaluri infection. Pathogen-specific pattern and tissue-specific pattern were found in the two infections. The significantly regulated expressions of catfish Hsp70 genes after bacterial infections suggested their involvement in immune response in catfish.

High-density interspecific genetic linkage mapping provides insights into genomic incompatibility between channel catfish and blue catfish

Catfish is the leading aquaculture species in the United States. The interspecific hybrid catfish produced by mating female channel catfish with male blue catfish outperform both of their parent species in a number of traits. However, mass production of the hybrids has been difficult because of reproductive isolation. Investigations of genome structure and organization of the hybrids provide insights into the genetic basis for maintenance of species divergence in the face of gene flow, thereby helping develop strategies for introgression and efficient production of the hybrids for aquaculture. In this study, we constructed a high-density genetic linkage map using the hybrid catfish system with the catfish 250K SNP array. A total of 26,238 SNPs were mapped to 29 linkage groups, with 12,776 unique marker positions. The linkage map spans approximately 3240 cM with an average intermarker distance of 0.25 cM. A fraction of markers (986 of 12,776) exhibited significant deviation from the expected Mendelian ratio of segregation, and they were clustered in major genomic blocks across 15 LGs, most notably LG9 and LG15. The distorted markers exhibited significant bias for maternal alleles among the backcross progenies, suggesting strong selection against the blue catfish alleles. The clustering of distorted markers within genomic blocks should lend insights into speciation as marked by incompatibilities between the two species. Such findings should also have profound implications for understanding the genomic evolution of closely related species as well as the introgression of hybrid production programs in aquaculture.

Phenoscape: Identifying Candidate Genes for Evolutionary Phenotypes

Phenotypes resulting from mutations in genetic model organisms can help reveal candidate genes for evolutionarily important phenotypic changes in related taxa. Although testing candidate gene hypotheses experimentally in nonmodel organisms is typically difficult, ontology-driven information systems can help generate testable hypotheses about developmental processes in experimentally tractable organisms. Here, we tested candidate gene hypotheses suggested by expert use of the Phenoscape Knowledgebase, specifically looking for genes that are candidates responsible for evolutionarily interesting phenotypes in the ostariophysan fishes that bear resemblance to mutant phenotypes in zebrafish. For this, we searched ZFIN for genetic perturbations that result in either loss of basihyal element or loss of scales phenotypes, because these are the ancestral phenotypes observed in catfishes (Siluriformes). We tested the identified candidate genes by examining their endogenous expression patterns in the channel catfish, Ictalurus punctatus. The experimental results were consistent with the hypotheses that these features evolved through disruption in developmental pathways at, or upstream of, brpf1 and eda/edar for the ancestral losses of basihyal element and scales, respectively. These results demonstrate that ontological annotations of the phenotypic effects of genetic alterations in model organisms, when aggregated within a knowledgebase, can be used effectively to generate testable, and useful, hypotheses about evolutionary changes in morphology.

Testing the evolutionary conservation of vocal motoneurons in vertebrates

Medullary motoneurons drive vocalization in many vertebrate lineages including fish, amphibians, birds, and mammals. The developmental history of vocal motoneuron populations in each of these lineages remains largely unknown. The highly conserved transcription factor Paired-like Homeobox 2b (Phox2b) is presumed to be expressed in all vertebrate hindbrain branchial motoneurons, including laryngeal motoneurons essential for vocalization in humans. We used immunohistochemistry and in situ hybridization to examine Phox2b protein and mRNA expression in caudal hindbrain and rostral spinal cord motoneuron populations in seven species across five chordate classes. Phox2b was present in motoneurons dedicated to sound production in mice and frogs (bullfrog, African clawed frog), but not those in bird (zebra finch) or bony fish (midshipman, channel catfish). Overall, the pattern of caudal medullary motoneuron Phox2b expression was conserved across vertebrates and similar to expression in sea lamprey. These observations suggest that motoneurons dedicated to sound production in vertebrates are not derived from a single developmentally or evolutionarily conserved progenitor pool.

Second generation physical and linkage maps of yellowtail (Seriola quinqueradiata) and comparison of synteny with four model fish

Background

Physical and linkage maps are important aids for the assembly of genome sequences, comparative analyses of synteny, and to search for candidate genes by quantitative trait locus analysis. Yellowtail, Seriola quinqueradiata, is an economically important species in Japanese aquaculture, and genetic information will be useful for DNA-assisted breeding. We report the construction of a second generation radiation hybrid map, its synteny analysis, and a second generation linkage map containing SNPs (single nucleotide polymorphisms) in yellowtail.

Results

Approximately 1.4 million reads were obtained from transcriptome sequence analysis derived from 11 tissues of one individual. To identify SNPs, cDNA libraries were generated from a pool of 500 whole juveniles, and the gills and kidneys of 100 adults. 9,356 putative SNPs were detected in 6,025 contigs, with a minor allele frequency ≥25%. The linkage and radiation hybrid maps were constructed based on these contig sequences. 2,081 markers, including 601 SNPs markers, were mapped onto the linkage map, and 1,532 markers were mapped in the radiation hybrid map.

Conclusions

The second generation linkage and physical maps were constructed using 6,025 contigs having SNP markers. These maps will aid the de novo assembly of sequencing reads, linkage studies and the identification of candidate genes related to important traits. The comparison of marker contigs in the radiation hybrid map indicated that yellowtail is evolutionarily closer to medaka than to green-spotted pufferfish, three-spined stickleback or zebrafish. The synteny analysis may aid studies of chromosomal evolution in yellowtail compared with model fish.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1600-7) contains supplementary material, which is available to authorized users.

Transcriptome assembly, gene annotation and tissue gene expression atlas of the rainbow trout

Efforts to obtain a comprehensive genome sequence for rainbow trout are ongoing and will be complemented by transcriptome information that will enhance genome assembly and annotation. Previously, transcriptome reference sequences were reported using data from different sources. Although the previous work added a great wealth of sequences, a complete and well-annotated transcriptome is still needed. In addition, gene expression in different tissues was not completely addressed in the previous studies. In this study, non-normalized cDNA libraries were sequenced from 13 different tissues of a single doubled haploid rainbow trout from the same source used for the rainbow trout genome sequence. A total of ~1.167 billion paired-end reads were de novo assembled using the Trinity RNA-Seq assembler yielding 474,524 contigs > 500 base-pairs. Of them, 287,593 had homologies to the NCBI non-redundant protein database. The longest contig of each cluster was selected as a reference, yielding 44,990 representative contigs. A total of 4,146 contigs (9.2%), including 710 full-length sequences, did not match any mRNA sequences in the current rainbow trout genome reference. Mapping reads to the reference genome identified an additional 11,843 transcripts not annotated in the genome. A digital gene expression atlas revealed 7,678 housekeeping and 4,021 tissue-specific genes. Expression of about 16,000–32,000 genes (35–71% of the identified genes) accounted for basic and specialized functions of each tissue. White muscle and stomach had the least complex transcriptomes, with high percentages of their total mRNA contributed by a small number of genes. Brain, testis and intestine, in contrast, had complex transcriptomes, with a large numbers of genes involved in their expression patterns. This study provides comprehensive de novo transcriptome information that is suitable for functional and comparative genomics studies in rainbow trout, including annotation of the genome.

A transcriptomic analysis of striped catfish (Pangasianodon hypophthalmus) in response to salinity adaptation: De novo assembly, gene annotation and marker discovery

The striped catfish (Pangasianodon hypophthalmus) culture industry in the Mekong Delta in Vietnam has developed rapidly over the past decade. The culture industry now however, faces some significant challenges, especially related to climate change impacts notably from predicted extensive saltwater intrusion into many low topographical coastal provinces across the Mekong Delta. This problem highlights a need for development of culture stocks that can tolerate more saline culture environments as a response to expansion of saline water-intruded land. While a traditional artificial selection program can potentially address this need, understanding the genomic basis of salinity tolerance can assist development of more productive culture lines. The current study applied a transcriptomic approach using Ion PGM technology to generate expressed sequence tag (EST) resources from the intestine and swim bladder from striped catfish reared at a salinity level of 9ppt which showed best growth performance. Total sequence data generated was 467.8Mbp, consisting of 4,116,424 reads with an average length of 112bp. De novo assembly was employed that generated 51,188 contigs, and allowed identification of 16,116 putative genes based on the GenBank non-redundant database. GO annotation, KEGG pathway mapping, and functional annotation of the EST sequences recovered with a wide diversity of biological functions and processes. In addition, more than 11,600 simple sequence repeats were also detected. This is the first comprehensive analysis of a striped catfish transcriptome, and provides a valuable genomic resource for future selective breeding programs and functional or evolutionary studies of genes that influence salinity tolerance in this important culture species.

The cytochrome P450 genes of channel catfish: their involvement in disease defense responses as revealed by meta-analysis of RNA-Seq data sets

BACKGROUND

Cytochrome P450s (CYPs) encode one of the most diverse enzyme superfamily in nature. They catalyze oxidative reactions of endogenous molecules and exogenous chemicals.

METHODS

We identified CYPs genes through in silico analysis using EST, RNA-Seq and genome databases of channel catfish. Phylogenetic analyses and conserved syntenic analyses were conducted to determine their identities and orthologies. Meta-analysis of RNA-Seq databases was conducted to analyze expression profile of CYP genes following bacterial infection.

RESULTS

A full set of 61 CYP genes was identified and characterized in channel catfish. Phylogenetic tree and conserved synteny provided strong evidence of their identities and orthorlogy. Lineage-specific gene duplication was evident in a number of clans in channel catfish. CYP46A1 is missing in the catfish genome as observed with syntenic analysis and RT-PCR analysis. Thirty CYPs were found up- or down-regulated in liver, while seven and eight CYPs were observed regulated in intestine and gill following bacterial infection.

CONCLUSION

We systematically identified and characterized a full set of 61 CYP genes in channel catfish and studied their expression profiles after bacterial infection. While bacterial challenge altered the expression of large numbers of CYP genes, the mechanisms and significance of these changes are not known.

GENERAL SIGNIFICANCE

This work provides an example to systematically study CYP genes in non-model species. Moreover, it provides a basis for further toxicological and physiological studies in channel catfish.

A candidate gene association study for growth performance in an improved giant freshwater prawn (Macrobrachium rosenbergii ) culture line

A candidate gene approach using type I single nucleotide polymorphism (SNP) markers can provide an effective method for detecting genes and gene regions that underlie phenotypic variation in adaptively significant traits. In the absence of available genomic data resources, transcriptomes were recently generated in Macrobrachium rosenbergii to identify candidate genes and markers potentially associated with growth. The characterisation of 47 candidate loci by ABI re-sequencing of four cultured and eight wild samples revealed 342 putative SNPs. Among these, 28 SNPs were selected in 23 growth-related candidate genes to genotype in 200 animals selected for improved growth performance in an experimental GFP culture line in Vietnam. The associations between SNP markers and individual growth performance were then examined. For additive and dominant effects, a total of three exonic SNPs in glycogen phosphorylase (additive), heat shock protein 90 (additive and dominant) and peroxidasin (additive), and a total of six intronic SNPs in ankyrin repeats-like protein (additive and dominant), rolling pebbles (dominant), transforming growth factor-β induced precursor (dominant), and UTP-glucose-1-phosphate uridylyltransferase 2 (dominant) genes showed significant associations with the estimated breeding values in the experimental animals (P =0.001-0.031). Individually, they explained 2.6-4.8 % of the genetic variance (R²=0.026-0.048). This is the first large set of SNP markers reported for M. rosenbergii and will be useful for confirmation of associations in other samples or culture lines as well as having applications in marker-assisted selection in future breeding programs.

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.

Whole genome comparative analysis of channel catfish (Ictalurus punctatus) with four model fish species

BACKGROUND

Comparative mapping is a powerful tool to study evolution of genomes. It allows transfer of genome information from the well-studied model species to non-model species. Catfish is an economically important aquaculture species in United States. A large amount of genome resources have been developed from catfish including genetic linkage maps, physical maps, BAC end sequences (BES), integrated linkage and physical maps using BES-derived markers, physical map contig-specific sequences, and draft genome sequences. Application of such genome resources should allow comparative analysis at the genome scale with several other model fish species.

RESULTS

In this study, we conducted whole genome comparative analysis between channel catfish and four model fish species with fully sequenced genomes, zebrafish, medaka, stickleback and Tetraodon. A total of 517 Mb draft genome sequences of catfish were anchored to its genetic linkage map, which accounted for 62% of the total draft genome sequences. Based on the location of homologous genes, homologous chromosomes were determined among catfish and the four model fish species. A large number of conserved syntenic blocks were identified. Analysis of the syntenic relationships between catfish and the four model fishes supported that the catfish genome is most similar to the genome of zebrafish.

CONCLUSION

The organization of the catfish genome is similar to that of the four teleost species, zebrafish, medaka, stickleback, and Tetraodon such that homologous chromosomes can be identified. Within each chromosome, extended syntenic blocks were evident, but the conserved syntenies at the chromosome level involve extensive inter-chromosomal and intra-chromosomal rearrangements. This whole genome comparative map should facilitate the whole genome assembly and annotation in catfish, and will be useful for genomic studies of various other fish species.

Transcriptome-wide single nucleotide polymorphisms (SNPs) for abalone (Haliotis midae): validation and application using GoldenGate medium-throughput genotyping assays

Haliotis midae is one of the most valuable commercial abalone species in the world, but is highly vulnerable, due to exploitation, habitat destruction and predation. In order to preserve wild and cultured stocks, genetic management and improvement of the species has become crucial. Fundamental to this is the availability and employment of molecular markers, such as microsatellites and single nucleotide (SNPs). Transcriptome sequences generated through sequencing-by-synthesis technology were utilized for the in vitro and in silico identification of 505 putative SNPs from a total of 316 selected contigs. A subset of 234 SNPs were further validated and characterized in wild and cultured abalone using two Illumina GoldenGate genotyping assays. Combined with VeraCode technology, this genotyping platform yielded a 65%–69% conversion rate (percentage polymorphic markers) with a global genotyping success rate of 76%–85% and provided a viable means for validating SNP markers in a non-model species. The utility of 31 of the validated SNPs in population structure analysis was confirmed, while a large number of SNPs (174) were shown to be informative and are, thus, good candidates for linkage map construction. The non-synonymous SNPs (50) located in coding regions of genes that showed similarities with known proteins will also be useful for genetic applications, such as the marker-assisted selection of genes of relevance to abalone aquaculture.

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.

Comparative genomic analysis of catfish linkage group 8 reveals two homologous chromosomes in zebrafish and other teleosts with extensive inter-chromosomal rearrangements

Background

Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study.

Results

With existing genome resources, 287 unique genes were identified in LG8. Comparative genome analysis indicated that most of these 287 genes on catfish LG8 are located on two homologous chromosomes of zebrafish, medaka, stickleback, and three chromosomes of green-spotted pufferfish. Large numbers of conserved syntenies were identified. Detailed analysis of the conserved syntenies in relation to chromosome level similarities revealed extensive inter-chromosomal and intra-chromosomal rearrangements during evolution. Of the 287 genes, 35 genes were found to be duplicated in the catfish genome, with the vast majority of the duplications being interchromosomal.

Conclusions

Comparative genome analysis is a powerful tool even in the absence of a well-assembled whole genome sequence. In spite of sequence stacking due to low resolution of the linkage and physical maps, conserved syntenies can be identified although the exact gene order and orientation are unknown at present. Through chromosome-level comparative analysis, homologous chromosomes among teleosts can be identified. Syntenic analysis should facilitate annotation of the catfish genome, which in turn, should facilitate functional inference of genes based on their orthology.

Genome-wide identification, characterization and phylogenetic analysis of 50 catfish ATP-binding cassette (ABC) transporter genes

Background

Although a large set of full-length transcripts was recently assembled in catfish, annotation of large gene families, especially those with duplications, is still a great challenge. Most often, complexities in annotation cause mis-identification and thereby much confusion in the scientific literature. As such, detailed phylogenetic analysis and/or orthology analysis are required for annotation of genes involved in gene families. The ATP-binding cassette (ABC) transporter gene superfamily is a large gene family that encodes membrane proteins that transport a diverse set of substrates across membranes, playing important roles in protecting organisms from diverse environment.

Methodology/Principal Findings

In this work, we identified a set of 50 ABC transporters in catfish genome. Phylogenetic analysis allowed their identification and annotation into seven subfamilies, including 9 ABCA genes, 12 ABCB genes, 12 ABCC genes, 5 ABCD genes, 2 ABCE genes, 4 ABCF genes and 6 ABCG genes. Most ABC transporters are conserved among vertebrates, though cases of recent gene duplications and gene losses do exist. Gene duplications in catfish were found for ABCA1, ABCB3, ABCB6, ABCC5, ABCD3, ABCE1, ABCF2 and ABCG2.

Conclusion/Significance

The whole set of catfish ABC transporters provide the essential genomic resources for future biochemical, toxicological and physiological studies of ABC drug efflux transporters. The establishment of orthologies should allow functional inferences with the information from model species, though the function of lineage-specific genes can be distinct because of specific living environment with different selection pressure.

Transcriptome sequencing and analysis of wild Amur Ide (Leuciscus waleckii) inhabiting an extreme alkaline-saline lake reveals insights into stress adaptation

Background

Amur ide (Leuciscus waleckii) is an economically and ecologically important species in Northern Asia. The Dali Nor population inhabiting Dali Nor Lake, a typical saline-alkaline lake in Inner Mongolia, is well-known for its adaptation to extremely high alkalinity. Genome information is needed for conservation and aquaculture purposes, as well as to gain further understanding into the genetics of stress tolerance. The objective of the study is to sequence the transcriptome and obtain a well-assembled transcriptome of Amur ide.

Results

The transcriptome of Amur ide was sequenced using the Illumina platform and assembled into 53,632 cDNA contigs, with an average length of 647 bp and a N50 length of 1,094 bp. A total of 19,338 unique proteins were identified, and gene ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses classified all contigs into functional categories. Open Reading Frames (ORFs) were detected from 34,888 (65.1%) of contigs with an average length of 577 bp, while 9,638 full-length cDNAs were identified. Comparative analyses revealed that 31,790 (59.3%) contigs have a significant similarity to zebrafish proteins, and 27,096 (50.5%), 27,524 (51.3%) and 27,996 (52.2%) to teraodon, medaka and three-spined stickleback proteins, respectively. A total of 10,395 microsatellites and 34,299 SNPs were identified and classified. A dN/dS analysis on unigenes was performed, which identified that 61 of the genes were under strong positive selection. Most of the genes are associated with stress adaptation and immunity, suggesting that the extreme alkaline-saline environment resulted in fast evolution of certain genes.

Conclusions

The transcriptome of Amur ide had been deeply sequenced, assembled and characterized, providing a valuable resource for a better understanding of the Amur ide genome. The transcriptome data will facilitate future functional studies on the Amur ide genome, as well as provide insight into potential mechanisms for adaptation to an extreme alkaline-saline environment.

Characterization of the pearl oyster (Pinctada martensii) mantle transcriptome unravels biomineralization genes

Pearl oyster, Pinctada martensii, is a marine bivalve species widely distributed in tropic and subtropic marine coasts. Mantle is the special tissue of P. martensii that secretes biomineralization proteins inducing shell deposition as well as iridescent nacre both in the inner shell and artificial nucleus. The pearl oyster is very efficient for artificial pearl production and is therefore an ideal organism for studies into the processes of biomineralization. However, deficiency of transcriptome information limits the insight into biomineralization mechanisms and pearl formation. In this study, we sequenced and characterized the P. martensii mantle transcriptome using 454 pyrosequencing. A total of 25,723 unique transcripts were assembled from 220,824 quality reads, followed by annotation and Gene Ontology classification analysis. A total of 146 unique transcript segments homologous to 49 reference biomineralization genes were identified, including calcineurin-binding protein, amorphous calcium carbonate binding protein 1, calmodulin, calponin-like protein, carbonic anhydrase 1, glycine-rich shell matrix protein, lysine-rich matrix protein, mantle gene or protein, nacrein, pearlin, PIF, regucalcin, and shematrin. The sequence data enabled the identification of 10,285 potential single nucleotide polymorphism loci and 7,836 putative indels, providing a resource for molecular biomarker, population genetics, and functional genomic studies. A large number of candidate genes for biomineralization were identified, considerably enriching resources for the study of shell formation. These sequence data will notably advance biomineralization and transcriptome study in pearl oyster and other Pinctada species.

Evasion of mucosal defenses during Aeromonas hydrophila infection of channel catfish (Ictalurus punctatus) skin

The mucosal surfaces of fish serve as the first line of defense against the myriad of aquatic pathogens present in the aquatic environment. The immune repertoire functioning at these interfaces is still poorly understood. The skin, in particular, must process signals from several fronts, sensing and integrating environmental, nutritional, social, and health cues. Pathogen invasion can disrupt this delicate homeostasis with profound impacts on signaling throughout the organism. Here, we investigated the transcriptional effects of virulent Aeromonas hydrophila infection in channel catfish skin, Ictalurus punctatus. We utilized a new 8 × 60 K Agilent microarray for catfish to examine gene expression profiles at critical early timepoints following challenge--2 h, 8 h, and 12 h. Expression of a total of 2,168 unique genes was significantly perturbed during at least one timepoint. We observed dysregulation of genes involved in antioxidant, cytoskeletal, immune, junctional, and nervous system pathways. In particular, A. hydrophila infection rapidly altered a number of potentially critical lectins, chemokines, interleukins, and other mucosal factors in a manner predicted to enhance its ability to adhere to and invade the catfish host.

Early mucosal responses in blue catfish (Ictalurus furcatus) skin to Aeromonas hydrophila infection

Bacterial pathogens are well-equipped to detect, adhere to, and initiate infection in their finfish hosts. The mucosal surfaces of fish, such as the skin, function as the front line of defense against such bacterial insults that are routinely encountered in the aquatic environment. While recent progress has been made, and despite the obvious importance of mucosal surfaces, the precise molecular events that occur soon after encountering bacterial pathogens remain unclear. Indeed, these early events are critical in mounting appropriate responses that ultimately determine host survival or death. In the present study, we investigated the transcriptional consequences of a virulent Aeromonas hydrophila challenge in the skin of blue catfish, Ictalurus furcatus. We utilized an 8×60K Agilent microarray to examine gene expression profiles at key early timepoints following challenge (2 h, 12 h, and 24 h). A total of 1155 unique genes were significantly altered during at least one timepoint. We observed dysregulation in a number of genes involved in diverse pathways including those involved in antioxidant responses, apoptosis, cytoskeletal rearrangement, immunity, and extracellular matrix protein diversity and regulation. Taken together, A. hydrophila coordinately modulates mucosal factors across numerous cellular pathways in a manner predicted to enhance its ability to adhere to and infect the blue catfish host.

Comparative analysis of the transcriptome in tissues secreting purple and white nacre in the pearl mussel Hyriopsis cumingii

The triangle sail mussel Hyriopsis cumingii (Lea) is the most important mussel species used for commercial freshwater pearl production in China. Mussel color is an important indicator of pearl quality. To identify genes involved in the nacre coloring, we conducted RNA-seq and obtained 541,268 sequences (298 bp average size) and 440,034 sequences (293 bp average size) in secreting purple and white nacre libraries (P- and W-libraries), respectively. The 981,302 Expressed Sequence Tags (ESTs) were assembled into 47,812 contigs and 289,386 singletons. In BLASTP searches of the deduced protein, 22,495 were proteins with functional annotations. Thirty-three genes involved in pearl or shell formation were identified. Digital expression analysis identified a total of 358 differentially expressed genes, and 137 genes in the P-library and 221 genes in the W-library showed significantly higher expression. Furthermore, a set of SSR motifs and SNPs between the two samples was identified from the ESTs, which provided the markers for genetic linkage, QTL analysis and future breeding. These EST sequences provided valuable information to further understand the molecular mechanisms involved in the formation, color determination and evolution of the pearl or shell.

Phylogenetic analyses uncover a novel clade of transferrin in nonmammalian vertebrates

Transferrin is a protein super-family involved in iron transport, a central process in cellular homeostasis. Throughout the evolution of vertebrates, transferrin members have diversified into distinct subfamilies including serotransferrin, ovotransferrin, lactoferrin, melanotransferrin, the inhibitor of carbonic anhydrase, pacifastin, and the major yolk protein in sea urchin. Previous phylogenetic analyses have established the branching order of the diverse transferrin subfamilies but were mostly focused on the transferrin repertoire present in mammals. Here, we conduct a comprehensive phylogenetic analysis of transferrin protein sequences in sequenced vertebrates, placing a special focus on the less-studied nonmammalian vertebrates. Our analyses uncover a novel transferrin clade present across fish, sauropsid, and amphibian genomes but strikingly absent from mammals. Our reconstructed scenario implies that this novel class emerged through a duplication event at the vertebrate ancestor, and that it was subsequently lost in the lineage leading to mammals. We detect footprints of accelerated evolution following the duplication event, which suggest positive selection and early functional divergence of this novel clade. Interestingly, the loss of this novel class of transferrin in mammals coincided with the divergence by duplication of lactoferrin and serotransferrin in this lineage. Altogether, our results provide novel insights on the evolution of iron-binding proteins in the various vertebrate groups.

A standardized microsatellite marker panel for parentage and kinship analyses in channel catfish, Ictalurus punctatus

This research was designed to produce a standardized set of microsatellite loci for parentage and kinship analyses in channel catfish, the leading species of US aquaculture. Three panels of five to six markers each were developed that contained a total of two dinucleotide-, eight trinucleotide- and seven tetranucleotide-microsatellite loci respectively. The loci had a range of nine to 31 alleles per locus in an outbred population. Based on the allele frequencies measured in commercial randomly bred broodstock, the combined probability of non-exclusion of an unrelated candidate parent pair was 5.36e-18. The combined probability of non-exclusion of unrelated identical genotypes was 2.58e-08. The microsatellite panels were validated by parentage and kinship evaluation in three populations. A total of 697 spawns were collected from matings of outbred broodstock over three spawning seasons, and parents were determined unambiguously for all but three spawns. Genotype analysis also enabled the identification of half-sibling and full-sibling families produced by pond spawning. In a second experiment, parentage was unambiguously determined in nine spawns from a population consisting of broodstock derived from only four families. A third experiment demonstrated that all but one of 374 individuals from 10 full-sibling families could be assigned to a family after coculture in an earthen pond for 1 year. The standardized microsatellite panels enable the development of pedigreed catfish populations and large-scale performance evaluations in common environments to support the genetic improvement of cultured catfish through selective breeding.

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

Background

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.

Results

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.

Conclusions

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.

AmphiEST: Enabling comparative analysis of ESTs from five developmental stages of amphioxus

Amphioxus has been an important model for understanding the evolution of chordates and origin of vertebrates. Comparative transcriptome analysis can facilitate delineation of gene expression patterns of amphioxus at different developmental stages. So far, however, few such analyses have been performed. Here we have systematically compared amphioxus ESTs from five developmental stages. For the egg, gastrula, neurula, larva and adult stages, amphioxus ESTs were assembled, respectively, into 3364, 3230, 10,299, 4052 and 3866 contigs, and 193, 183, 933, 178 and 151 singlets. 25,796 ORFs were identified, of which 6,529 predicted ORFs (egg, 922; gastrula, 997; neurula, 2,159; larva, 1,014; adult, 1,437) have found matches from the UniProt database. For all five stages, more GO terms were assigned to stage-specific contigs than to repertoire ones. Respectively, 2, 5, 11, 4 and 2 miRNAs were matched with ESTs from the five stages, and 12,232, 10,147, 43,170, 8,049 and 12,811 miRNA target sites were predicted from the 3' UTR sequences. A web interface was developed for visualization of the contigs, miRNAs and annotations (http://www.rich.yunda.org/test/amphioxusest/). The constructed AmphiEST database may serve as an essential resource for future studies of the development and evolution of amphioxus.

Molecular characterization, phylogenetic analysis and expression patterns of five protein arginine methyltransferase genes of channel catfish, Ictalurus punctatus (Rafinesque)

Protein arginine methylation, catalyzed by protein arginine methyltransferases (PRMT), has recently emerged as an important modification in the regulation of gene expression. In this communication, we identified and characterized the channel catfish orthologs to human PRMT 1, 3, 4 and 5, and PRMT4 like. Each PRMT nucleic acid sequence has an open reading frame (ORF) and 3'-untranslated regions. Each ORF appears to encode 361, 587 and 458 amino acid residues for PRMT1, PRMT4 and variant, respectively. The partial ORF of PRMT3 and PRMT5 encode 292 and 563 amino acids, respectively. By comparison with the human counterparts, each channel catfish PRMT also has conserved domains. For expression profile, the channel catfish PRMT1 transcript was detected by RT-PCR in spleens, anterior kidneys, livers, intestines, skin and gills of fish examined. Except in liver, the PRMT3 transcript was detected in all catfish tissues examined. However, the PRMT4 cDNA was detected in livers from all three catfish and gills from two fish, but not other tissues. This information will enable us to further elucidate PRMT functions in channel catfish.

Channel catfish (Ictalurus punctatus) protein disulphide isomerase, PDIA6: molecular characterization and expression regulated by bacteria and virus inoculation

Protein disulfide isomerases (PDIs) are thought to aid protein folding and assembly by catalyzing formation and shuffling of cysteine disulfide bonds in the endoplasmic reticulum (ER). Currently, increasing evidence suggests PDIs play an important role in host cell invasion and they are relevant targets for the host immune response. However the roles of specific PDIs in teleosts are little known. Here, we characterized the Protein disulfide isomerase family A, member 6 (PDIA6) from channel catfish, Ictalurus punctatus (named as ccPDIA6). The catfish ccPDIA6 gene was homologous to those of other vertebrate species with 13 exons and 12 introns. The consensus full-length ccPDIA6 cDNA contained an ORF of 1320 bp encoding a putative protein of 439 amino acids. It had a 19 amino acid signal peptide and two active thioredoxin-like domains. Sequence of phylogenic analysis and multiple alignments showed that ccPDIA6 was conserved throughout vertebrate evolution. Southern blot analysis suggested the presence of one copy of the ccPDIA6 gene in the catfish genome. Tissue distribution shows that ccPDIA6 was expressed in all examined tissues at the mRNA level. When using the aquatic zoonotic pathogens such as Edwardsiella tara, Streptococcus iniae, and channel catfish reovirus （CCRV） to challenge channel catfish, ccPDIA6 expression was significant changed in immune-related tissues such as head kidney, intestine, liver and spleen. The results suggested that ccPDIA6 might play an important role in the immunity of channel catfish. This is the first report that the PDI gene may be involved in fish host defense against pathogen infection.

Transcriptomic characterization of the larval stage in gilthead seabream (Sparus aurata) by 454 pyrosequencing

Gilthead seabream (Sparus aurata) is a teleost belonging to the family Sparidae with a high economical relevance in the Mediterranean countries. Although genomic tools have been developed in this species in order to investigate its physiology at the molecular level and consequently its culture, genomic information on post-embryonic development is still scarce. In this study, we have investigated the transcriptome of a marine teleost during the larval stage (from hatching to 60 days after hatching) by the use of 454 pyrosequencing technology. We obtained a total of 68,289 assembled contigs, representing putative transcripts, belonging to 54,606 different clusters. Comparison against all S. aurata expressed sequenced tags (ESTs) from the NCBI database revealed that up to 34,722 contigs, belonging to about 61% of gene clusters, are sequences previously not described. Contigs were annotated through an iterative Blast pipeline by comparison against databases such as NCBI RefSeq from Danio rerio, SwissProt or NCBI teleost ESTs. Our results indicate that we have enriched the number of annotated sequences for this species by more than 50% compared with previously existing databases for the gilthead seabream. Gene Ontology analysis of these novel sequences revealed that there is a statistically significant number of transcripts with key roles in larval development, differentiation, morphology, and growth. Finally, all information has been made available online through user-friendly interfaces such as GBrowse and a Blast server with a graphical frontend.

BIRC7 gene in channel catfish (Ictalurus punctatus): identification and expression analysis in response to Edwardsiella tarda, Streptococcus iniae and Channel catfish Hemorrhage Reovirus

A family member of inhibitor of apoptosis protein (IAP) termed baculoviral IAP repeat-containing 7 (BIRC7) from channel catfish (Ictalurus punctatus) was identified, the full length cDNA sequence of channel catfish BIRC7 (CcBIRC7) was 1686 bp, containing a 5'UTR of 93 bp, a 3'UTR of 399 bp with a poly (A) tail and an ORF of 1194 bp encoding a putative protein of 398 amino acids. The putative CcBIRC7 protein contains two BIR super-family conservative domains and a C-terminal RING finger motif. Phylogenetic analysis showed that catfish CcBIRC7 was moderately conserved with other BIRC7. Quantitative real-time PCR was conducted to examine the expression profiles of CcBIRC7 in healthy tissues and responding to different pathogens (Edwardsiella tarda, Streptococcus iniae and Channel catfish Hemorrhage Reovirus (CCRV)). CcBIRC7 was widely expressed in healthy tissues of channel catfish and with the highest 37.28-fold expression in blood. E. tarda and S. iniae could induce CcBIRC7 gene expression drastically in head kidney, liver and spleen, which the peak value reached 31.6-fold, 613.9-fold and 34.4-fold increase by E. tarda infection, and 248.3-fold, 1540.3-fold and 120.4-fold increase post S. iniae challenge, respectively. While, CCRV virus could slightly induce CcBIRC7 expression in head kidney and liver but reduce it in spleen. The result suggested BIRC7 may play a potential role in channel catfish innate immune system against bacterial and virus infections, especially as the anti-bacteria immune gene. This is the first report of BIRC7 gene identification and its expression in fish.

Development and characterisation of an expressed sequence tags (EST)-derived single nucleotide polymorphisms (SNPs) resource in rainbow trout

Background

There is considerable interest in developing high-throughput genotyping with single nucleotide polymorphisms (SNPs) for the identification of genes affecting important ecological or economical traits. SNPs are evenly distributed throughout the genome and are likely to be functionally relevant. In rainbow trout, in silico screening of EST databases represents an attractive approach for de novo SNP identification. Nevertheless, EST sequencing errors and assembly of EST paralogous sequences can lead to the identification of false positive SNPs which renders the reliability of EST-derived SNPs relatively low. Further validation of EST-derived SNPs is therefore required. The objective of this work was to assess the quality of and to validate a large number of rainbow trout EST-derived SNPs.

Results

A panel of 1,152 EST-derived SNPs was selected from the INRA Sigenae SNP database and was genotyped in standard and double haploid individuals from several populations using the Illumina GoldenGate BeadXpress assay. High-quality genotyping data were obtained for 958 SNPs representing a genotyping success rate of 83.2 %, out of which, 350 SNPs (36.5 %) were polymorphic in at least one population and were designated as true SNPs. They also proved to be a potential tool to investigate genetic diversity of the species, as the set of SNP successfully sorted individuals into three main groups using STRUCTURE software. Functional annotations revealed 28 non-synonymous SNPs, out of which four substitutions were predicted to affect protein functions. A subset of 223 true SNPs were polymorphic in the two INRA mapping reference families and were integrated into the INRA microsatellite-based linkage map.

Conclusions

Our results represent the first study of EST-derived SNPs validation in rainbow trout, a species whose genome sequences is not yet available. We designed several specific filters in order to improve the genotyping yield. Nevertheless, our selection criteria should be further improved in order to reduce the observed high rate of false positive SNPs which results from the occurrence of whole genome duplications.

RNA-seq analysis of mucosal immune responses reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus

The mucosal surfaces of fish (gill, skin, gastrointestinal tract) are important sites of bacterial exposure and host defense mechanisms. In mammalian systems, the intestinal epithelium is well characterized as both a selectively permeable barrier regulated by junctional proteins and as a primary site of infection for a number of enteric pathogens including viruses, bacteria, and parasites. The causative bacterium of enteric septicemia of catfish, Edwardsiella ictaluri, is believed to gain entry through the intestinal epithelium, with previous research using a rat intestinal epithelial cell line (IEC-6) indicating actin polymerization and receptor-mediated endocytosis as potential mechanisms of uptake. Here, we utilized high-throughput RNA-seq to characterize the role of the intestinal epithelial barrier following E. ictaluri challenge. A total of 197.6 million reads were obtained and assembled into 176,481 contigs with an average length of 893.7 bp and N50 of 1676 bp. The assembled contigs contained 14,457 known unigenes, including 2719 genes not previously identified in other catfish transcriptome studies. Comparison of digital gene expression between challenged and control samples revealed 1633 differentially expressed genes at 3 h, 24 h, and 3 day following exposure. Gene pathway analysis of the differentially expressed gene set indicated the centrality of actin cytoskeletal polymerization/remodelling and junctional regulation in pathogen entry and subsequent inflammatory responses. The expression patterns of fifteen differentially expressed genes related to intestinal epithelial barrier dysfunction were validated by quantitative real-time RT-PCR (average correlation coeff. 0.92, p < 0.001). Our results set a foundation for future studies comparing mechanisms of pathogen entry and mucosal immunity across several important catfish pathogens including E. ictaluri, Edwardsiellatarda, Flavobacterium columnare, and virulent atypical Aeromonas hydrophila. Understanding of molecular mechanisms of pathogen entry during infection will provide insight into strategies for selection of resistant catfish brood stocks against various diseases.

Characterization of common carp transcriptome: sequencing, de novo assembly, annotation and comparative genomics

BACKGROUND

Common carp (Cyprinus carpio) is one of the most important aquaculture species of Cyprinidae with an annual global production of 3.4 million tons, accounting for nearly 14% of the freshwater aquaculture production in the world. Due to the economical and ecological importance of common carp, genomic data are eagerly needed for genetic improvement purpose. However, there is still no sufficient transcriptome data available. The objective of the project is to sequence transcriptome deeply and provide well-assembled transcriptome sequences to common carp research community.

RESULT

Transcriptome sequencing of common carp was performed using Roche 454 platform. A total of 1,418,591 clean ESTs were collected and assembled into 36,811 cDNA contigs, with average length of 888 bp and N50 length of 1,002 bp. Annotation was performed and a total of 19,165 unique proteins were identified from assembled contigs. Gene ontology and KEGG analysis were performed and classified all contigs into functional categories for understanding gene functions and regulation pathways. Open Reading Frames (ORFs) were detected from 29,869 (81.1%) contigs with an average ORF length of 763 bp. From these contigs, 9,625 full-length cDNAs were identified with sequence length from 201 bp to 9,956 bp. Comparative analysis revealed that 27,693(75.2%) contigs have significant similarity to zebrafish Refseq proteins, and 24,371(66.2%), 24,501(66.5%) and 25,025(70.0%) to teraodon, medaka and three-spined stickleback refseq proteins. A total of 2,064 microsatellites were initially identified from 1,730 contigs, and 1,639 unique sequences had sufficient flanking sequences on both sides for primer design.

CONCLUSION

The transcriptome of common carp had been deep sequenced, de novo assembled and characterized, providing the valuable resource for better understanding of common carp genome. The transcriptome data will facilitate future functional studies on common carp genome, and gradually apply in breeding programs of common carp, as well as closely related other Cyprinids.

Molecular characterization of three L-type lectin genes from channel catfish, Ictalurus punctatus and their responses to Edwardsiella ictaluri challenge

L-type lectins have a leguminous lectin domain and can bind to high-mannose type oligosaccharides. In the secretory pathway, L-type lectins play crucial roles in selective protein trafficking, sorting and targeting. Three L-type lectins were cloned in the channel catfish, Ictalurus punctatus, the 53 kDa endoplasmic reticulum ER-Golgi intermediate compartment protein (ERGIC-53), the vesicular integral protein of 36 kDa (VIP36) and VIP36-like. Phylogenetic analysis indicated that the catfish genes are orthologous to their counterparts in other species. Southern blot analysis demonstrated that all three L-type lectin genes are likely single-copy genes in the catfish genome. Analysis of expression in healthy tissues using quantitative real time RT-PCR indicated that all three genes are expressed widely in all tested tissues, but with strong tissue preference of expression: ERGIC-53 was found to be abundantly expressed in the liver, VIP36 was found to be abundantly expressed in the head-kidney, whereas VIP36-like was found to be abundantly expressed in the brain. Upon infection with Edwardsiella ictaluri, expressions of the three genes all had significant up-regulation in the head-kidney, but had distinct expression patterns: ERGIC-53 was gradually induced with the highest expression 7 days after challenge in the head-kidney, but was down-regulated in the liver, spleen, and brain. VIP36 was highly induced in the head-kidney, and 3 days after challenge in the brain, but was not up-regulated in any other tissues or timepoints after challenge. Expression levels of the catfish VIP36-like gene appeared to also respond to infection, albeit with differing patterns among the tested tissues. Taken together, our results indicate that all three L-type lectin genes may be involved in the immune responses of catfish after infection with E. ictaluri.

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.

Analysis of the Asian seabass transcriptome based on expressed sequence tags

Analysis of transcriptomes is of great importance in genomic studies. Asian seabass is an important fish species. A number of genomic tools in it were developed, while large expressed sequence tag (EST) data are lacking. We sequenced ESTs from nine normalized cDNA libraries and obtained 11 431 high-quality ESTs. We retrieved 8524 ESTs from dbEST database and analyzed all 19 975 ESTs using bioinformatics tools. After clustering, we obtained 8837 unique sequences (2838 contigs and 5999 singletons). The average contig length was 574 bp. Annotation of these unique sequences revealed that 48.9% of them showed significant homology to RNA sequences in GenBank. Functional classification of the unique ESTs identified a broad range of genes involved in different functions. We identified 6114 putative single-nucleotide polymorphisms and 634 microsatellites in ESTs. We discovered different temporal and spatial expression patterns of some immune-related genes in the Asian seabass after challenging with a pathogen Vibrio harveyi. The unique EST sequences are being used in developing a cDNA microarray to examine global gene expression and will also facilitate future whole-genome sequence assembly and annotation of Asian seabass and comparative genomics.

Construction and characterization of the BAC library for common carp Cyprinus carpio L. and establishment of microsynteny with zebrafish Danio rerio

A bacterial artificial chromosome (BAC) library of common carp Cyprinus carpio L. was constructed as a part of ongoing common carp genome project, which is aiming assembly of common carp genome. The library, containing a total of 92,160 BAC clones with an average insert size of 141 kb, was constructed into the restriction site of Hind III on BAC vector CopyControl pCC1BAC, covering 7.7 X haploid genome equivalents. Three dimension pools and superpools of the BAC library were established and 23 positive clones of 14 targets were identified from one-fifth of the BAC library. Pilot project of BAC end sequencing was conducted on 2,688 BAC ends from 1,344 clones and harvested 2,522 high-quality Q20 sequences with average length of 677 bp. The sequencing success rate was 93.8% and pair-end success rate was 92.3%. A total of 212 microsyntenies had been established between common carp and zebrafish genomes as a trial for genome-wide comparative genomics in these two closely related species.

Transcriptome sequencing, and rapid development and application of SNP markers for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae)

The legume pod borer, Maruca vitrata (Lepidoptera: Crambidae), is an insect pest species of crops grown by subsistence farmers in tropical regions of Africa. We present the de novo assembly of 3729 contigs from 454- and Sanger-derived sequencing reads for midgut, salivary, and whole adult tissues of this non-model species. Functional annotation predicted that 1320 M. vitrata protein coding genes are present, of which 631 have orthologs within the Bombyx mori gene model. A homology-based analysis assigned M. vitrata genes into a group of paralogs, but these were subsequently partitioned into putative orthologs following phylogenetic analyses. Following sequence quality filtering, a total of 1542 putative single nucleotide polymorphisms (SNPs) were predicted within M. vitrata contig assemblies. Seventy one of 1078 designed molecular genetic markers were used to screen M. vitrata samples from five collection sites in West Africa. Population substructure may be present with significant implications in the insect resistance management recommendations pertaining to the release of biological control agents or transgenic cowpea that express Bacillus thuringiensis crystal toxins. Mutation data derived from transcriptome sequencing is an expeditious and economical source for genetic markers that allow evaluation of ecological differentiation.

Transcriptomics for understanding marine fish larval development1This review is part of a virtual symposium on current topics in aquaculture of marine fish and shellfish

The larval phase is a crucial period in the life of marine fish. During this phase, the organism will acquire the phenotype of an adult fish through the development of tissues and organs and the maturation of some of the principal physiological functions. Many biological processes (differentiation, cellular proliferation, growth, etc.) are regulated during this period. These regulations take place at different biological levels and particularly concern the expression of genes involved in larval ontogenesis processes. The development of bioinformatic resources (DNA or cDNA sequences) and molecular tools enabling high throughput gene expression analysis (microarrays) have allowed the transcriptome of marine fish species to be studied. In the present review, we summarize the main findings from transcriptomic investigations of development of marine fish larvae. Special attention is paid to investigations of transcriptomic patterns during postembryonic development and to the impact of environmental or nutritional factors on the transcriptome of marine fish larvae. Transcriptomic approaches will be especially useful in the future for investigating the effect of temperature and water acidification (or pH) on the development of different fish species in the context of global climate change.

Molecular characterization and expression analysis of the channel catfish cathepsin D genes

Cathepsin D is a lysosomal aspartic proteinase that participates in various degradation functions of the cell. In this study, we characterized the cathepsin D genes in channel catfish and found two genes encoding catfish cathepsin D, referred to as cathepsin D1 and D2 genes. These two genes are highly similar in genomic structure and organization, sharing a moderate level of amino acid sequence similarity (56%). Genomic Southern analysis suggested the presence of a single copy of each of the cathepsin D1 and D2 genes. Phylogenetic analysis provided strong evidence that two cathepsin D genes are present in most of the teleost lineage, with cathepsin D2 likely having been lost in some higher vertebrate lineages. The catfish cathepsin D1 and D2 genes are expressed in virtually all the 11 tested tissues (brain, gill, heart, head kidney, trunk kidney, intestine, liver, muscle, skin, spleen, and stomach) on the transcript level, but appear to exhibit greater levels of expression in immune-related tissues and organs. Upon infection with Edwardsiella ictaluri, the expression of the catfish cathepsin D genes showed the most significant changes in liver and head kidney, with time points and magnitude of transcript changes varying between the two genes. We additionally examined bacterially-mediated changes of expression in gill, intestine, and trunk kidney. The fact that bacterial infection can induce expression of the cathepsin D genes and that they appeared to be expressed naturally at higher levels in immune-related organs may suggest that they are an important component of the innate immune response of catfish against bacterial infections.

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

Dicentrarchus labrax is one of the major marine aquaculture species in the European Union. In this study, we have developed a directed-sequencing strategy to sequence three sea bass chromosomes and compared results with other teleosts. Three BAC DNA pools were created from sea bass BAC clones that mapped to stickleback chromosomes/groups V, XVII and XXI. The pools were sequenced to 17-39x coverage by pyrosequencing. Data assembly was supported by Sanger reads and mate pair data and resulted in superscaffolds of 13.2 Mb, 17.5 Mb and 13.7 Mb respectively. Annotation features of the superscaffolds include 1477 genes. We analyzed size change of exon, intron and intergenic sequence between teleost species and deduced a simple model for the evolution of genome composition in teleost lineage. Combination of second generation sequencing technologies, Sanger sequencing and genome partitioning strategies allows "high-quality draft assemblies" of chromosome-sized superscaffolds, which are crucial for the prediction and annotation of complete genes.

mRNA/microRNA Profile at the Metamorphic Stage of Olive Flounder (Paralichthys olivaceus)

Flatfish is famous for the asymmetric transformation during metamorphosis. The molecular mechanism behind the asymmetric development has been speculated over a century and is still not well understood. To date, none of the metamorphosis-related genes has been identified in flatfish. As the first step to screen metamorphosis-related gene, we constructed a whole-body cDNA library and a whole-body miRNA library in this study and identified 1051 unique ESTs, 23 unique miRNAs, and 4 snoRNAs in premetamorphosing and prometamorphosing Paralichthys olivaceus. 1005 of the ESTs were novel, suggesting that there was a special gene expression profile at metamorphic stage. Four miRNAs (pol-miR-20c, pol-miR-23c, pol-miR-130d, and pol-miR-181e) were novel to P. olivaceus; they were characterized as highly preserved homologies of published miRNAs but with at least one nucleotide differed. Representative 24 mRNAs and 23 miRNAs were quantified during metamorphosis of P. olivaceus by using quantitative RT PCR or stem-loop qRT PCR. Our results showed that 20 of mRNAs might be associated with early metamorphic events, 10 of mRNAs might be related with later metamorphic events, and 16 of miRNAs might be involved in the regulation of metamorphosis. The data provided in this study would be helpful for further identifying metamorphosis-related gene in P. olivaceus.