Transcriptome of channel catfish (Ictalurus punctatus): initial analysis of genes and expression profiles of the head kidney

Molecular Identification and Dual Functions of Two Different CXC Chemokines in Nile Tilapia (Oreochromis niloticus) against Streptococcus agalactiae and Flavobacterium columnare

Two CXC chemokines in Nile tilapia (On-CXC1 and On-CXC2) were identified at both the genomic and proteomic levels. A southern blot analysis and comparison searching in Ensembl confirmed the typical structure of the CXC chemokine genes and provided evidence for unusual mechanisms used to generate the two different CXC chemokine transcripts that have not been reported in other vertebrate species so far. The expression levels of On-CXC1 and On-CXC2 were analyzed by quantitative real-time PCR. These two mRNAs were detected in various tissues of normal Nile tilapia, especially in the spleen, heart, and head kidney, indicating a homeostatic function in immunosurveillance. A time-course experiment clearly demonstrated that these two transcripts were effectively enhanced in the head kidney, spleen and trunk kidney of Nile tilapia 6, 12 and 24 h after injection with Streptococcus agalactiae but were down-regulated in all tested tissues at 48 h, reflecting the fact that they have short half-lives during the crucial response to pathogens that is characteristic of CXC chemokine genes in other vertebrates. Functional analyses obviously exhibited that these two CXC chemokines at concentrations of 1–10 μg strongly inactivated S. agalactiae and Flavobacterium columnare and effectively induced phagocytosis of leukocytes in vitro.

Comparative analysis of the transcriptome of the Amazonian fish species Colossoma macropomum (tambaqui) and hybrid tambacu by next generation sequencing

Background

The C. macropomum is a characiform fish from the Amazon basin that has been hybridized with other pacu species to produce commercial hybrids, such as the tambacu. However, little is known of the functional genomics of the parental species or these hybrid forms. The transcriptome of C. macropomum and tambacu were sequenced using 454 Roche platform (pyrosequencing) techniques to characterize the domains of Gene Ontology (GO) and to evaluate the levels of gene expression in the two organisms.

Results

The 8,188,945 reads were assembled into 400,845 contigs. A total of 58,322 contigs were annotated with a predominance of biological processes for both organisms, as determined by Gene Ontology (GO). Similar numbers of metabolic pathways were identified in both the C. macropomum and the tambacu, with the metabolism category presenting the largest number of transcripts. The BUSCO analysis indicated that our assembly was more than 40% complete. We identified 21,986 genes for the two fishes. The P and Log2FC values indicated significant differences in the levels of gene expression, with a total of 600 up-regulated genes.

Conclusion

In spite of the lack of a reference genome, the functional annotation was successful, and confirmed a considerable difference in the specificity and levels of gene expression between the two organisms. This report provides a comprehensive baseline for the genetic management of these commercially important fishes, in particular for the identification of specific genes that may represent markers involved in the immunity, growth, and fertility of these organisms, with potential practical applications in aquaculture management.

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.

Molecular characterization, functional analysis, and defense mechanisms of two CC chemokines in Nile tilapia (Oreochromis niloticus) in response to severely pathogenic bacteria

Two full-length cDNAs encoding CC chemokine genes in Nile tilapia (Oreochromis niloticus) (On-CC1 and On-CC2) were cloned and characterized. On-CC1 and On-CC2 showed signature cysteine motifs consisting of four cysteines. The expression levels of On-CC1 and On-CC2 were analyzed by RT-PCR, which showed that low expression of these two genes was only observed in the peripheral blood leukocytes (PBLs) and spleen of normal fish. Expression levels of these two molecules were quantified in 13 tissues of fish infected with virulent strains of Streptococcus agalactiae and Flavobacterium columnare. Most tissues, especially PBLs, the spleen and the liver, expressed significantly higher mRNA levels than the controls, particularly at 12 and 24 h after infection (P < 0.05). The current study strongly indicates that CC chemokine genes in Nile tilapia are crucially involved in the early immune responses to pathogens. Functional analyses clearly demonstrated that 10 and 100 μg/ml of recombinant rOn-CC1 and rOn-CC2 proteins efficiently enhanced the phagocytic activity (in vitro) of Nile tilapia phagocytes. Finally, Southern blot analysis and searching in Ensembl databases demonstrated that two different functional CC chemokine genes and other pseudogene fragments were discovered in the Nile tilapia genome.

Transcriptome and expression profiling analysis of the hemocytes reveals a large number of immune-related genes in mud crab Scylla paramamosain during Vibrio parahaemolyticus infection

BACKGROUND

Mud crab Scylla paramamosain is an economically important marine species in China. However, frequent outbreaks of infectious diseases caused by marine bacteria, such as Vibrio parahaemolyticus, result in great economic losses.

METHODOLOGY/PRINCIPAL FINDINGS

Comparative de novo transcriptome analysis of S. paramamosain infected with V. parahaemolyticus was carried out to investigate the molecular mechanisms underlying the immune response to pathogenic bacteria by using the Illumina paired-end sequencing platform. A total of 52,934,042 clean reads from the hemocytes of V. parahaemolyticus-infected mud crabs and controls were obtained and assembled into 186,193 contigs. 59,120 unigenes were identified from 81,709 consensus sequences of mud crabs and 48,934 unigenes were matched proteins in the Nr or Swissprot databases. Among these, 10,566 unigenes belong to 3 categories of Gene Ontology, 25,349 to 30 categories of KEGG, and 15,191 to 25 categories of COG database, covering almost all functional categories. By using the Solexa/Illumina's DGE platform, 1213 differentially expressed genes (P<0.05), including 538 significantly up-regulated and 675 down-regulated, were detected in V. parahaemolyticus-infected crabs as compared to that in the controls. Transcript levels of randomly-chosen genes were further measured by quantitative real-time PCR to confirm the expression profiles. Many differentially expressed genes are involved in various immune processes, including stimulation of the Toll pathway, Immune Deficiency (IMD) pathway, Ras-regulated endocytosis, and proPO-activating system.

CONCLUSIONS/SIGNIFICANCE

Analysis of the expression profile of crabs under infection provides invaluable new data for biological research in S. paramamosain, such as the identification of novel genes in the hemocytes during V. parahaemolyticus infection. These results will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection and will be helpful for diseases prevention in crab aquaculture.

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.

Transcriptome and expression profiling analysis of Leuciscus waleckii: an exploration of the alkali-adapted mechanisms of a freshwater teleost

The strategies by which freshwater teleosts maintain acid-base homeostasis under alkaline stress are attractive and have been explored for a long time. In this study, a cyprinid fish that tolerates extremely alkaline environments (pH 9.6), Leuciscus waleckii, was used as a model to explore the molecular mechanisms of acid-base regulation. Using a lab-controlled alkaline challenge test and 454 sequencing, the transcriptomes of their gills and kidney were profiled and compared. mRNA profiling produced 1 826 022 reads, generated 30 606 contigs with an average length of 1022 bp, of which 19 196 were annotated successfully. Comparative analysis of the expression profiles between alkaline and freshwater L. waleckii habitats revealed approximately 4647 and 7184 genes that were differentially expressed (p < 0.05) in gills and kidney, respectively, of which 2398 and 5127 had more than twofold changes in expression. Gene ontology analysis and KEGG enrichment analysis were conducted. Comprehensive analysis found that genes involved in ion transportation, ammonia transportation, and arachidonic acid metabolism pathways changed dramatically and played important roles in acid-base homeostasis in fish under alkaline stress. These results support the existing hypotheses about candidate genes involved in acid-base regulation under alkaline stress and prompt several new hypotheses. The large transcriptome dataset collected in this study is a useful resource for the exploration of homeostasis modulation in other fish species.

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.

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.

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.

Transcriptome analysis of head kidney in grass carp and discovery of immune-related genes

BACKGROUND

Grass carp (Ctenopharyngodon idella) is one of the most economically important freshwater fish, but its production is often affected by diseases that cause serious economic losses. To date, no good breeding varieties have been obtained using the oriented cultivation technique. The ability to identify disease resistance genes in grass carp is important to cultivate disease-resistant varieties of grass carp.

RESULTS

In this study, we constructed a non-normalized cDNA library of head kidney in grass carp, and, after clustering and assembly, we obtained 3,027 high-quality unigenes. Solexa sequencing was used to generate sequence tags from the transcriptomes of the head kidney in grass carp before and after grass carp reovirus (GCRV) infection. After processing, we obtained 22,144 tags that were differentially expressed by more than 2-fold between the uninfected and infected groups. 679 of the differentially expressed tags (3.1%) mapped to 483 of the unigenes (16.0%). The up-regulated and down-regulated unigenes were annotated using gene ontology terms; 16 were annotated as immune-related and 42 were of unknown function having no matches to any of the sequences in the databases that were used in the similarity searches. Semi-quantitative RT-PCR revealed four unknown unigenes that showed significant responses to the viral infection. Based on domain structure predictions, one of these sequences was found to encode a protein that contained two transmembrane domains and, therefore, may be a transmembrane protein. Here, we proposed that this novel unigene may encode a virus receptor or a protein that mediates the immune signalling pathway at the cell surface.

CONCLUSION

This study enriches the molecular basis data of grass carp and further confirms that, based on fish tissue-specific EST databases, transcriptome analysis is an effective route to discover novel functional genes.

A review of catfish genomics: progress and perspectives

Catfish is one of the lower teleosts whose genome research is important for evolutionary genomics. As the major aquaculture species in the USA, its genome research also has practical and economical implications. Much progress has been made in recent years, including the development of large numbers of molecular markers, the construction of framework genetic linkage maps, the identification of putative markers involved in performance traits, and the development of genomic resources. Repetitive elements have been identified and characterized in the catfish genome that should facilitate physical analysis of the catfish genome. A large number of genes or full-length cDNAs have been analysed using genomic approaches, providing information on gene structure, gene evolution and gene expression in relation to functions. Catfish genome research has come to a stage when physical mapping through BAC contig construction is greatly demanded, in order to develop regional markers for QTL analysis and for large-scale comparative mapping. The current effort in large-scale EST analysis and type I marker mapping should further enhance research efficiency through comparative mapping. Candidate gene identification is being accelerated through the use of cDNA microarrays.

Genomics and mapping of teleostei (bony fish)

Until recently, the Human Genome Project held centre stage in the press releases concerning sequencing programmes. However, in October 2001, it was announced that the Japanese puffer fish (Takifugu rubripes, Fugu) was the second vertebrate organism to be sequenced to draft quality. Briefly, the spotlight was on fish genomes. There are currently two other fish species undergoing intensive sequencing, the green spotted puffer fish (Tetraodon nigroviridis) and the zebrafish (Danio rerio). But this trio are, in many ways, atypical representations of the current state of fish genomic research. The aim of this brief review is to demonstrate the complexity of fish as a group of vertebrates and to publicize the ‘lesser-known’ species, all of which have something to offer.

Molecular responses of ceruloplasmin to Edwardsiella ictaluri infection and iron overload in channel catfish (Ictalurus punctatus)

Ceruloplasmin is a serum ferroxidase that carries more than 90% of the copper in plasma and has documented roles in iron homeostasis as well as antioxidative functions. In our previous studies, it has been shown that the ceruloplasmin gene is strongly up-regulated in catfish during challenge with Edwardsiella ictaluri. However, little is known about the function of this gene in teleost fish. The objective of this study, therefore, was to characterize the ceruloplasmin gene from channel catfish, determine its genomic organization, profile its patterns of tissue expression, and establish its potential for physiological antioxidant responses in catfish after bacterial infection with E. ictaluri and iron treatment. The genomic organization suggested that the catfish ceruloplasmin gene had 20 exons and 19 introns, encoding 1074 amino acids. Exon sizes of the catfish ceruloplasmin gene were close to or identical with mammalian and zebrafish homologs. Further phylogenetic analyses suggested that the gene was highly conserved through evolution. The catfish ceruloplasmin gene was mapped to both the catfish physical map and linkage map. The catfish ceruloplasmin gene was mainly expressed in liver with limited expression in other tissues, and it was significantly up-regulated in the liver after bacterial infection alone or after co-injection with bacteria and iron-dextran, while expression was not significantly induced with iron-dextran treatment alone.

The catfish genome database cBARBEL: an informatic platform for genome biology of ictalurid catfish

The catfish genome database, cBARBEL (abbreviated from catfish Breeder And Researcher Bioinformatics Entry Location) is an online open-access database for genome biology of ictalurid catfish (Ictalurus spp.). It serves as a comprehensive, integrative platform for all aspects of catfish genetics, genomics and related data resources. cBARBEL provides BLAST-based, fuzzy and specific search functions, visualization of catfish linkage, physical and integrated maps, a catfish EST contig viewer with SNP information overlay, and GBrowse-based organization of catfish genomic data based on sequence similarity with zebrafish chromosomes. Subsections of the database are tightly related, allowing a user with a sequence or search string of interest to navigate seamlessly from one area to another. As catfish genome sequencing proceeds and ongoing quantitative trait loci (QTL) projects bear fruit, cBARBEL will allow rapid data integration and dissemination within the catfish research community and to interested stakeholders. cBARBEL can be accessed at http://catfishgenome.org.

Transcriptome and expression profiling analysis revealed changes of multiple signaling pathways involved in immunity in the large yellow croaker during Aeromonas hydrophila infection

BACKGROUND

The large yellow croaker (Pseudosciaena crocea) is an economically important marine fish in China suffering from severe outbreaks of infectious disease caused by marine bacteria such as Aeromonas hydrophila (A. hydrophila), resulting in great economic losses. However, the mechanisms involved in the immune response of this fish to bacterial infection are not fully understood. To understand the molecular mechanisms underlying the immune response to such pathogenic bacteria, we used high-throughput deep sequencing technology to investigate the transcriptome and comparative expression profiles of the large yellow croaker infected with A. hydrophila.

RESULTS

A total of 13,611,340 reads were obtained and assembled into 26,313 scaffolds in transcriptional responses of the A. hydrophila-infected large yellow croaker. Via annotation to the NCBI database, we obtained 8216 identified unigenes. In total, 5590 (68%) unigenes were classified into Gene Ontology, and 3094 unigenes were found in 20 KEGG categories. These genes included representatives from almost all functional categories. By using Solexa/Illumina's DeepSAGE, 1996 differentially expressed genes (P value < 0.05) were detected in comparative analysis of the expression profiles between A. hydrophila-infected fish and control fish, including 727 remarkably upregulated genes and 489 remarkably downregulated genes. Dramatic differences were observed in genes involved in the inflammatory response. Bacterial infection affected the gene expression of many components of signaling cascades, including the Toll-like receptor, JAK-STAT, and MAPK pathways. Genes encoding factors involved in T cell receptor (TCR) signaling were also revealed to be regulated by infection in these fish.

CONCLUSION

Based on our results, we conclude that the inflammatory response may play an important role in the early stages of infection. The signaling cascades such as the Toll-like receptor, JAK-STAT, and MAPK pathways are regulated by A. hydrophila infection. Interestingly, genes encoding factors involved in TCR signaling were revealed to be downregulated by infection, indicating that TCR signaling was suppressed at this early period. These results revealed changes of multiple signaling pathways involved in immunity during A. hydrophila infection, which will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection in the large yellow croaker.

Molecular characterization and gene expression of the channel catfish ferritin H subunit after bacterial infection and iron treatment

Ferritins are the major iron storage protein in the cytoplasm of cells, responsible for regulating levels of intracellular iron. Ferritin genes are widely distributed in both prokaryotes and eukaryotes. In mammals, ferritin molecules are composed of heavy- (H) and light- (L) chain subunits; amphibian genomes contain three ferritin-type genes (H; middle, M; and L subunits); and teleost genomes to date contain H and M subunits. The objective of this study was to characterize the ferritin H gene in channel catfish (Ictalurus punctatus) to determine its genomic organization and copy numbers, to determine its patterns of tissue expression, and to establish if it is involved in defense responses of catfish after bacterial infection. The catfish ferritin H gene was completely sequenced and characterized, using both mRNA and genomic DNA. Catfish ferritin H gene has a full-length mRNA sequence of 999 bp, an open reading frame of 534 bp, and 4,704 bp genomic DNA sequence. Catfish ferritin H has a 5 exon and 4 intron genetic organization, containing a long 5'-untranslated region, which shares high similarity with mammalian and zebrafish genes. Based on phylogenetic analyses, the catfish ferritin H gene is highly conserved throughout evolution. Southern blot analysis suggested that the ferritin H gene has only one copy in the catfish genome. The catfish ferritin H gene was widely expressed in various healthy tissues. The catfish ferritin H gene was significantly up-regulated in the liver after intraperitoneal injection of iron dextran and coinjection of Edwardsiella ictaluri and iron dextran treatment, suggesting its role in iron metabolism and immunity.

Development of genomic resources in support of sequencing, assembly, and annotation of the catfish genome

Major progress has been made in catfish genomics including construction of high-density genetic linkage maps, BAC-based physical maps, and integration of genetic linkage and physical maps. Large numbers of ESTs have been generated from both channel catfish and blue catfish. Microarray platforms have been developed for the analysis of genome expression. Genome repeat structures are studied, laying grounds for whole genome sequencing. USDA recently approved funding of the whole genome sequencing project of catfish using the next generation sequencing technologies. Generation of the whole genome sequence is a historical landmark of catfish research as it opens the real first step of the long march toward genetic enhancement. The research community needs to be focused on aquaculture performance and production traits, take advantage of the unprecedented genome information and technology, and make real progress toward genetic improvements of aquaculture brood stocks.

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

Twelve cDNA libraries from two species of catfish have been sequenced, resulting in the generation of nearly 500,000 ESTs.

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.

Comparative analysis of catfish BAC end sequences with the zebrafish genome

Background

Comparative mapping is a powerful tool to transfer genomic information from sequenced genomes to closely related species for which whole genome sequence data are not yet available. However, such an approach is still very limited in catfish, the most important aquaculture species in the United States. This project was initiated to generate additional BAC end sequences and demonstrate their applications in comparative mapping in catfish.

Results

We reported the generation of 43,000 BAC end sequences and their applications for comparative genome analysis in catfish. Using these and the additional 20,000 existing BAC end sequences as a resource along with linkage mapping and existing physical map, conserved syntenic regions were identified between the catfish and zebrafish genomes. A total of 10,943 catfish BAC end sequences (17.3%) had significant BLAST hits to the zebrafish genome (cutoff value ≤ e^-5), of which 3,221 were unique gene hits, providing a platform for comparative mapping based on locations of these genes in catfish and zebrafish. Genetic linkage mapping of microsatellites associated with contigs allowed identification of large conserved genomic segments and construction of super scaffolds.

Conclusion

BAC end sequences and their associated polymorphic markers are great resources for comparative genome analysis in catfish. Highly conserved chromosomal regions were identified to exist between catfish and zebrafish. However, it appears that the level of conservation at local genomic regions are high while a high level of chromosomal shuffling and rearrangements exist between catfish and zebrafish genomes. Orthologous regions established through comparative analysis should facilitate both structural and functional genome analysis in catfish.

Construction of genetic linkage maps and comparative genome analysis of catfish using gene-associated markers

A genetic linkage map of the channel catfish genome (N=29) was constructed using EST-based microsatellite and single nucleotide polymorphism (SNP) markers in an interspecific reference family. A total of 413 microsatellites and 125 SNP markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 allowed mapping of 331 markers (259 microsatellites and 72 SNPs) to 29 linkage groups. Each linkage group contained 3-18 markers. The largest linkage group contained 18 markers and spanned 131.2 cM, while the smallest linkage group contained 14 markers and spanned only 7.9 cM. The linkage map covered a genetic distance of 1811 cM with an average marker interval of 6.0 cM. Sex-specific maps were also constructed; the recombination rate for females was 1.6 times higher than that for males. Putative conserved syntenies between catfish and zebrafish, medaka, and Tetraodon were established, but the overall levels of genome rearrangements were high among the teleost genomes. This study represents a first-generation linkage map constructed by using EST-derived microsatellites and SNPs, laying a framework for large-scale comparative genome analysis in catfish. The conserved syntenies identified here between the catfish and the three model fish species should facilitate structural genome analysis and evolutionary studies, but more importantly should facilitate functional inference of catfish genes. Given that determination of gene functions is difficult in nonmodel species such as catfish, functional genome analysis will have to rely heavily on the establishment of orthologies from model species.

Systematic sequencing of mRNA from the Antarctic krill (Euphausia superba) and first tissue specific transcriptional signature

Background

Little is known about the genome sequences of Euphausiacea (krill) although these crustaceans are abundant components of the pelagic ecosystems in all oceans and used for aquaculture and pharmaceutical industry. This study reports the results of an expressed sequence tag (EST) sequencing project from different tissues of Euphausia superba (the Antarctic krill).

Results

We have constructed and sequenced five cDNA libraries from different Antarctic krill tissues: head, abdomen, thoracopods and photophores. We have identified 1.770 high-quality ESTs which were assembled into 216 overlapping clusters and 801 singletons resulting in a total of 1.017 non-redundant sequences. Quantitative RT-PCR analysis was performed to quantify and validate the expression levels of ten genes presenting different EST countings in krill tissues. In addition, bioinformatic screening of the non-redundant E. superba sequences identified 69 microsatellite containing ESTs. Clusters, consensuses and related similarity and gene ontology searches were organized in a dedicated E. superba database .

Conclusion

We defined the first tissue transcriptional signatures of E. superba based on functional categorization among the examined tissues. The analyses of annotated transcripts showed a higher similarity with genes from insects with respect to Malacostraca possibly as an effect of the limited number of Malacostraca sequences in the public databases. Our catalogue provides for the first time a genomic tool to investigate the biology of the Antarctic krill.

Divergent Toll-like receptors in catfish (Ictalurus punctatus): TLR5S, TLR20, TLR21

Toll-like receptors (TLR) mediate pathogen recognition in vertebrate species through detection of conserved microbial ligands. Families of TLR molecules have been described from the genomes of the teleost fish model species zebrafish and Takifugu, but much research remains to characterize the full length sequences and pathogen specificities of individual TLR members in fish. While the majority of these pathogen receptors are conserved among vertebrate species with clear orthologues present in fish for most mammalian TLRs, several interesting differences are present in the TLR repertoire of teleost fish when compared to that of mammals. A soluble form of TLR5 has been reported from salmonid fish and Takifugu rubripes which is not present in mammals, and a large group of TLRs (arbitrarily numbered 19-23) was identified from teleost genomes with no easily discernible orthologues in mammals. To better understand these teleost adaptations to the TLR family, we have isolated, sequenced, and characterized the full-length cDNA and gene sequences of TLR5S, TLR20, and TLR21 from catfish as well as studied their expression pattern in tissues. We also mapped these genes to bacterial artificial chromosome (BAC) clones for genome analysis. While TLR5S appeared to be common in teleost fish, and TLR21 is common to birds, amphibians and fish, TLR20 has only been identified in zebrafish and catfish. Phylogenetic analysis of catfish TLR20 indicated that it is closely related to murine TLR11 and TLR12, two divergent TLRs about which little is known. All three genes appear to exist in catfish as single copy genes.

Characterization of a BAC library from channel catfish Ictalurus punctatus: indications of high levels of chromosomal reshuffling among teleost genomes

The CHORI-212 bacterial artificial chromosome (BAC) library was constructed by cloning EcoRI/EcoRI partially digested DNA into the pTARBAC2.1 vector. The library has an average insert size of 161 kb, and provides 10.6-fold coverage of the channel catfish haploid genome. Screening of 32 genes using overgo or cDNA probes indicated that this library had a good representation of the genome as all tested genes existed in the library. We previously reported sequencing of approximately 25,000 BAC ends that generated 20,366 high-quality BAC end sequences (BES) and identified a large number of sequences similar to known genes using BLASTX searches. In this work, particular attention was given to identification of BAC mate pairs with known genes from both ends. When identified, comparative genome analysis was conducted to determine syntenic regions of the catfish genome with the genomes of zebrafish and Tetraodon. Of the 141 mate pairs with known genes from channel catfish, conserved syntenies were identified in 34 (24.1%), with 30 conserved in the zebrafish genome and 14 conserved in the Tetraodon genome. Additional analysis of three of the 34 conserved syntenic groups by direct sequencing indicated conserved gene contents in all three species. This indicates that comparative genome analysis may provide shortcuts to genome analysis in catfish, especially for short genomic regions once the conserved syntenies are identified.

A BAC-based physical map of the channel catfish genome

Catfish is the major aquaculture species in the United States. To enhance its genome studies involving genetic linkage and comparative mapping, a bacterial artificial chromosome (BAC) contig-based physical map of the channel catfish (Ictalurus punctatus) genome was generated using four-color fluorescence-based fingerprints. Fingerprints of 34,580 BAC clones (5.6x genome coverage) were generated for the FPC assembly of the BAC contigs. A total of 3307 contigs were assembled using a cutoff value of 1x10(-20). Each contig contains an average of 9.25 clones with an average size of 292 kb. The combined contig size for all contigs was 0.965 Gb, approximately the genome size of the channel catfish. The reliability of the contig assembly was assessed by both hybridization of gene probes to BAC clones contained in the fingerprinted assembly and validation of randomly selected contigs using overgo probes designed from BAC end sequences. The presented physical map should greatly enhance genome research in the catfish, particularly aiding in the identification of genomic regions containing genes underlying important performance traits.

Generation and analysis of expressed sequence tags from the ciliate protozoan parasite Ichthyophthirius multifiliis

BACKGROUND

The ciliate protozoan Ichthyophthirius multifiliis (Ich) is an important parasite of freshwater fish that causes 'white spot disease' leading to significant losses. A genomic resource for large-scale studies of this parasite has been lacking. To study gene expression involved in Ich pathogenesis and virulence, our goal was to generate expressed sequence tags (ESTs) for the development of a powerful microarray platform for the analysis of global gene expression in this species. Here, we initiated a project to sequence and analyze over 10,000 ESTs.

RESULTS

We sequenced 10,368 EST clones using a normalized cDNA library made from pooled samples of the trophont, tomont, and theront life-cycle stages, and generated 9,769 sequences (94.2% success rate). Post-sequencing processing led to 8,432 high quality sequences. Clustering analysis of these ESTs allowed identification of 4,706 unique sequences containing 976 contigs and 3,730 singletons. These unique sequences represent over two million base pairs (~10% of Plasmodium falciparum genome, a phylogenetically related protozoan). BLASTX searches produced 2,518 significant (E-value < 10-5) hits and further Gene Ontology (GO) analysis annotated 1,008 of these genes. The ESTs were analyzed comparatively against the genomes of the related protozoa Tetrahymena thermophila and P. falciparum, allowing putative identification of additional genes. All the EST sequences were deposited by dbEST in GenBank (GenBank: EG957858-EG966289). Gene discovery and annotations are presented and discussed.

CONCLUSION

This set of ESTs represents a significant proportion of the Ich transcriptome, and provides a material basis for the development of microarrays useful for gene expression studies concerning Ich development, pathogenesis, and virulence.

Towards the ictalurid catfish transcriptome: generation and analysis of 31,215 catfish ESTs

Background

EST sequencing is one of the most efficient means for gene discovery and molecular marker development, and can be additionally utilized in both comparative genome analysis and evaluation of gene duplications. While much progress has been made in catfish genomics, large-scale EST resources have been lacking. The objectives of this project were to construct primary cDNA libraries, to conduct initial EST sequencing to generate catfish EST resources, and to obtain baseline information about highly expressed genes in various catfish organs to provide a guide for the production of normalized and subtracted cDNA libraries for large-scale transcriptome analysis in catfish.

Results

A total of 17 cDNA libraries were constructed including 12 from channel catfish (Ictalurus punctatus) and 5 from blue catfish (I. furcatus). A total of 31,215 ESTs, with average length of 778 bp, were generated including 20,451 from the channel catfish and 10,764 from blue catfish. Cluster analysis indicated that 73% of channel catfish and 67% of blue catfish ESTs were unique within the project. Over 53% and 50% of the channel catfish and blue catfish ESTs, respectively, had significant similarities to known genes. All ESTs have been deposited in GenBank. Evaluation of the catfish EST resources demonstrated their potential for molecular marker development, comparative genome analysis, and evaluation of ancient and recent gene duplications. Subtraction of abundantly expressed genes in a variety of catfish tissues, identified here, will allow the production of low-redundancy libraries for in-depth sequencing.

Conclusion

The sequencing of 31,215 ESTs from channel catfish and blue catfish has significantly increased the EST resources in catfish. The EST resources should provide the potential for microarray development, polymorphic marker identification, mapping, and comparative genome analysis.

Structural characterisation and expression analysis of toll-like receptor 2 gene from catfish

Toll-like receptors (TLRs) are important components of innate immunity. They were found to recognise specific structures on pathogens termed pathogen-associated molecular patterns (PAMPs) and utilise conserved signaling pathways to activate pro-inflammatory cytokines and type-1 interferons. In spite of much understanding gained from the mammalian systems, many fish TLRs are unknown. Recent studies in Japanese flounder as well as in zebrafish suggested that the ligand binding and activation of inflammatory responses in fish may be different from and more complex than those found in mammals. In channel catfish, the major aquaculture species in the United States, only partial sequences of TLR3 and TLR5 were reported. As a part of efforts to characterise the innate immune components in channel catfish, here we cloned and sequenced both the cDNA and the gene for TLR2, a receptor believed mostly responsible for recognition of lipopeptides on the surface of most Gram-positive bacteria. However, expression analysis after infection with a Gram-negative bacterium, Edwardsiella ictaluri indicated that TLR2 was modestly down-regulated in the head kidney tissue of blue catfish, and with a similar pattern in the head kidney of channel catfish though the down-regulation in channel catfish was not statistically significant. In the spleen, an insignificant down-regulation was initially observed early after infection, with an increase of TLR expression later after infection. These results suggest the involvement of TLR2 in the responses after the bacterial infection. As LPS is believed to be the major PAMP for Gram-negative bacteria, additional research is warranted to determine the functions and mechanisms of TLR2 in infections of Gram-negative bacteria.

Repeat structure of the catfish genome: a genomic and transcriptomic assessment of Tc1-like transposon elements in channel catfish (Ictalurus punctatus)

We have assessed the distribution and diversity of members of the Tc1/mariner superfamily of transposable elements in the channel catfish (Ictalurus punctatus) genome as well as evaluating the extent of transcription of Tc1 transposases in the species. Through use of PCR amplification and sequencing, assessment of random BAC end sequences (BES) equivalent to 1.2% genome coverage, and screening of over 45,000 catfish ESTs, a significant proportion of Tc1-like elements and their associated transcripts were captured. Up to 4.2% of the catfish genome in base pairs appears to be composed of Tc1-like transposon-related sequences and a significant fraction of the catfish cellular mRNA, approximately 0.6%, was transcribed from transposon-related sequences in both sense and antisense orientations. Based on results of repeat-masking, as much as 10% of BAC end sequences from catfish, which is a random survey of the genome, contain some remnant of Tc1 elements, suggesting that these elements are present in the catfish genome as numerous, small remnants of the transposons. Phylogenetic analysis allowed comparison of catfish Tc1 transposase types with those found in other vertebrate and invertebrate species. In spite of the existence of many types of Tc1-like sequences that are not yet able to be placed in clades with strong statistical support, it is clear that multiple families of Tc1-like elements exist in channel catfish.

Channel catfish BAC-end sequences for marker development and assessment of syntenic conservation with other fish species

In the present study, 25 195 BAC ends for channel catfish (Ictalurus punctatus) were sequenced, generating 20 366 clean BAC-end sequences (BES), with an average read length of 557 bp after trimming. A total of 11 414 601 bp were generated, representing approximately 1.2% of the catfish genome. Based on this survey, the catfish genome was found to be highly AT-rich, with 60.7% A+T and 39.3% G+C. Approximately 12% of the catfish genome consisted of dispersed repetitive elements, with the Tc1/mariner transposons making up the largest percentage by base pair (4.57%). Microsatellites were detected in 17.5% of BES. Catfish BACs were anchored to the zebrafish and Tetraodon genome sequences by BLASTN, generating 16% and 8.2% significant hits (E < e(-5)) respectively. A total of 1074 and 773 significant hits were unique to the zebrafish and Tetraodon genomes, respectively, of which 417 and 406, respectively, were identified as known genes in other species, providing a major genome resource for comparative genomic mapping.

Toll-like receptor 3 and TICAM genes in catfish: species-specific expression profiles following infection with Edwardsiella ictaluri

Toll-like receptors (TLRs) are a family of transmembrane proteins that recognize specific pathogen-associated molecular patterns and use conserved signaling pathways to activate proinflammatory cytokines and type-1 interferons to fight infection. TLR3 in mammals is best known for its recognition of dsRNA as ligand and its MyD88-independent signaling. TLR3, upon recognition of dsRNA, recruits and binds its adaptor protein TIR domain-containing adapter molecule (TICAM) 1. Here we report the genomic sequences and structures of TLR3 and a TICAM adaptor from channel catfish (Ictalurus punctatus). Whereas a partial TLR3 cDNA sequence has been reported from channel catfish, and complete TLR3 genes are known from other teleost fish species, a complete TICAM sequence has not been previously reported from a nonmammalian species. Analysis of catfish TLR3 and TICAM expression after infection with Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish (ESC), suggested a conserved TLR3-TICAM receptor-adaptor relation in catfish. Comparison of TLR3 and TICAM expression profiles in channel catfish with those from the closely related blue catfish species (Ictalurus furcatus), which exhibits strong resistance to ESC, revealed a striking pattern of species-specific expression. A dramatic downregulation of TLR3 and TICAM gene expression was observed in blue catfish head kidney and spleen, which we speculate may be the result of maturation and migration of different cell types to and from the lymphoid tissues following infection.

In silico identification and expression analysis of 12 novel CC chemokines in catfish

Chemokines, a superfamily of chemotactic cytokines involved in recruitment, activation, and adhesion of a variety of leukocyte types to inflammatory foci, are a crucial component of the immune system of Sarcopterygiian vertebrates. Although all mammalian chemokines are believed to have been found, the status of these molecules in Actinopterygii was unknown until recently. The identification of chemokines in fish species has been complicated by low sequence conservation and confusion over expected numbers. Earlier discoveries of single fish chemokines coupled with rapidly expanding genetic resources in these species have recently provided a foundation for large-scale in silico discoveries of these important immune regulators. We report here the identification and expression analysis of 12 new CC chemokine sequences from catfish. When added to our previous report of 14 catfish CC chemokines, the number of CC chemokines in catfish now stands at 26, two more than known from humans. Establishing orthologous relationships among the majority of catfish CC chemokines, a newly available set of chicken CC chemokines, and their mammalian counterparts remain difficult, suggesting high levels of duplication and divergence within individual species.

Cloning and characterization of myogenic regulatory genes in three Ictalurid species

We report sequence, tissue expression and map-position data for myogenin, MYOD1, myostatin and follistatin in three Ictalurid catfish species: channel catfish (Ictalurus punctatus), blue catfish (I. furcatus) and white catfish (Ameiurus catus). These genes are involved in muscle growth and development in mammals and may play similar roles in catfish. Amino acid sequences were highly conserved among the three Ictalurid species (>95% identity), moderately conserved among catfish and zebrafish (approximately 80% identity), and less conserved among catfish and humans (approximately 40-60% identity) for all four genes. Gene structure (number of exons and introns and exon-intron boundaries) was conserved between catfish and other species for all genes. Myogenin and MYOD1 expression was limited to skeletal muscle in juvenile channel catfish, similar to expression patterns for these genes in other fish and mammalian species. Myostatin was expressed in a variety of tissues in juvenile channel catfish, a pattern common in other fish species but contrasting with data from mammals where myostatin is primarily expressed in skeletal muscle. Follistatin was expressed in juvenile catfish heart, testes and spleen. All four genes contained polymorphic microsatellite repeats in non-coding regions and linkage analysis based on inheritance of these microsatellite loci was used to place the genes on the channel catfish linkage map. Information provided in this study will be useful in further studies to determine the role these genes play in muscle growth and development in catfish.

Constitutive expression of three novel catfish CXC chemokines: homeostatic chemokines in teleost fish

Chemokines are best known for their vital role in leukocyte chemotaxis, as part of the larger inflammatory response. Expression analysis and functional characterization of chemokines in mammalian species have often overlooked the role of these proteins under homeostatic conditions. Recent investigations of chemokine diversity in teleost fish have also centered on the immune-related functions of chemotactic cytokines, such as CXCL8 and CXCL10. While a disease-based approach to chemokines is essential to the development of remediative therapies for both human and animal infections, it may be a poor measure of the overall complexity of chemokine functions. As part of a larger effort to assess the conservation of chemokine diversity in teleost fish, we report here the identification of three novel, constitutively expressed CXC chemokines from channel catfish (Ictalurus punctatus). Phylogenetic analyses indicated that two of the three CXC chemokines were orthologues for mammalian CXCL12 and CXCL14, respectively. Whereas a clear orthology could not yet be established for the third CXC chemokine, it shared highest amino acid identity with mammalian CXCL2. All three CXC chemokines show expression in a wide range of tissues, and early expression during development was observed for CXCL12. The expression of this new set of catfish CXC chemokines was not induced during challenge by infection of Edwardsiella ictaluri, the causative agent of the fish pathogen enteric septicemia of catfish. In contrast to the gene duplication of CXCL12 in carp and zebrafish, Southern blot analysis indicated that all three catfish CXC chemokines exist as single copy genes in the catfish genome suggesting that gene duplication of CXC chemokines in specific teleost fish was a recent evolutionary event.

Bioinformatic mining of type I microsatellites from expressed sequence tags of channel catfish (Ictalurus punctatus)

Gene-derived markers are pivotal to the analysis of genome structure, organization, and evolution and necessary for comparative genomics. However, gene-derived markers are relatively difficult to develop. This project utilized the genomic resources of channel catfish expressed sequence tags (ESTs) to identify simple sequence repeats (SSRs), or microsatellites. It took the advantage of ESTs for the establishment of gene identities, and of microsatellites for the acquisition of high polymorphism. When microsatellites are tagged to genes, the microsatellites can then be used as gene markers. A bioinformatic analysis of 43,033 ESTs identified 4855 ESTs containing microsatellites. Cluster analysis indicated that 1312 of these ESTs fell into 569 contigs, and the remaining 3534 ESTs were singletons. A total of 4103 unique microsatellite-containing genes were identified. The dinucleotide CA/TG and GA/TC pairs were the most abundant microsatellites. AT-rich microsatellite types were predominant among trinucleotide and tetranucleotide microsatellites, consistent with our earlier estimation that the catfish genome is highly AT-rich. Our preliminary results indicated that the majority of the identified microsatellites were polymorphic and, therefore, useful for genetic linkage mapping of catfish. Mapping of these gene-derived markers is under way, which will set the foundation for comparative genome analysis in catfish.

Putative SNP discovery in interspecific hybrids of catfish by comparative EST analysis

In this study, we identified putative SNP markers within genes by comparative analysis of expressed sequence tags (ESTs). Comparison of 849 ESTs from blue catfish (Ictalurus furcatus) with >11,000 ESTs from channel catfish (I. punctatus) deposited in GenBank resulted in the identification of 1020 putative SNPs within 161 genes, of which 145 were nuclear genes of known function. The observed frequency of SNPs within ESTs of the two closely related catfish species was 1.32 SNP per 100 bp. The majority of identified SNPs differed between the two species and, therefore, these SNPs are useful for mapping genes in channel catfish x blue catfish interspecific resource families. The SNPs that differed within species were also observed; these can be applied to genome scans in channel catfish resource families.

Analysis of genes expressed in head kidney of common carp Cyprinus carpio L. treated with cortisol

We analyzed genes expressed from head kidney of common carp Cyprinus carpio L. treated with cortisol. The results of single-pass sequencing of expressed sequence tags (ESTs) from 188 clones (AU240288-AU240367, AU301120-AU301227) from kidney cDNA are presented. One-hundred-twenty-seven clones (67.6%) were completely unknown and are likely to represent newly described genes, whereas 61 clones (32.4%) were identified based on matches to sequences in the database. The putative genes contain several ribosomal proteins, cytochrome oxidase subunits. Immune related cDNA clones identified from kidney were immunoglobulin light chain (n=4), FK506/rapamycin-binding protein (FKBP), CXC chemokine receptor type 4, complement factor B/C2-A3, peptidylprolyl isomerase (cyclophilin; Cyp)-like1, cyclophilin S1, heat shock-70 kDa protein-4, stress-activated protein kinase-3 (n=2). FKBP and cyclophilin genes expressed in normal tissues (head kidney, spleen, liver, brain and heart). Expression of FKBP and cyclophilin genes were not detected in liver, brain and heart when treated with cortisol for 16 h.

Towards a catalogue of genes transcribed in multiple tissues of Mytilus galloprovincialis

Mussels are among the most studied marine organisms because of their ecological role, economic importance and long-standing use in coastal pollution biomonitoring. Despite the bulk of data derived from traditional investigation areas, only limited knowledge is available on mussel genes and their expression in ordinary or stressing conditions. We present here the first data obtained by production and sequencing of 3' end-specific expressed sequence tags (3'ESTs) from multiple tissues of putatively unstressed mussels (Mytilus galloprovincialis). A total of 524 clusters (98 virtual consensuses and 426 singletons) derived from 829 reliable 3'ESTs were searched in the non-redundant (nr) sequence databases (National Centre for Biotechnology Information (Bethesda, USA), NCBI). Most contigs (59%) showed poor or no similarity, thus revealing unknown mussel genes, other contigs recognized some of the few available records of M. galloprovincialis (e.g. actin, cytochrome oxidase III and twitchin) and a number of significant similarities with mitochondrial (mt) or nuclear genes from other organisms were found. Actually, a variety of ribosomal as well as motility- and adhesion-related genes and some genes potentially involved in stress responses (e.g. myticin A, heat shock proteins, methallothionein) were putatively identified. Uniquely, mitochondrial transcripts were primarily represented by a cluster of 123 3'ESTs (1296 base pairs (bp) of the mussel 16S rRNA). In our knowledge, these results provide the first systematic production and annotation of ESTs in M. galloprovincialis. They also represent the first report of a wider project, based on the strategy of 3'EST identification from normal and stressed mussels and intended to define a mussel complementary DNA (cDNA) microarray for genome-wide transcription studies.

Translational machinery of channel catfish: II. Complementary DNA and expression of the complete set of 47 60S ribosomal proteins

Ribosomal protein genes have become widely used as markers for phylogenetic studies and comparative genomics, but they have not been available in fish. We have cloned and sequenced a complete set of all 47 60S ribosomal protein cDNAs from channel catfish (Ictalurus punctatus), of which 43 included the complete protein encoding regions. Most ribosomal protein mRNAs in channel catfish are highly similar to their mammalian counterparts. However, L4, L14, and L29 are significantly shorter in channel catfish than in mammals due to deletions in the 3' end of the gene. Two distantly related L5 cDNAs, L5a and L5b, were found in channel catfish. L5a is more similar to L5 in other vertebrates, while L5b showed significant levels of divergence, suggesting independent evolution of the two L5-encoding genes. The 47 ribosomal protein genes are generally highly expressed and together account for 11-14% of overall gene expression, depending on the tissues. Expression levels were highly variable both within a single tissue among different ribosomal protein genes, and among tissues with regard to a single ribosomal protein gene. Strong tissue preference expression was also observed for some ribosomal proteins. This set of ribosomal protein gene sequences represents one of the most complete sets from any single organism and will aid in fish phylogenetic and comparative genomic studies.

Translational machinery of channel catfish: I. A transcriptomic approach to the analysis of 32 40S ribosomal protein genes and their expression

Ribosomal protein (RP) genes have become widely used as markers for phylogenetic studies and comparative genomics. However, they have not been available for evolutionary studies in fish although teleosts are the largest group of vertebrates with more than 23,000 species. Using a transcriptomic approach, we have cloned and sequenced 32 40S RP complementary DNAs (cDNAs) from channel catfish (Ictalurus punctatus), making them one of the most complete sets of 40S RP gene sequences from a single organism. Most 40S RPs in channel catfish are highly similar to their orthologues in mammalian species, but S19, S21, and S25 are highly divergent. Only one type of cDNA was found for all RP genes except S26 and S27, for which two cDNAs were found in channel catfish. Alternatively spliced transcripts for the S3 and alternatively polyadenylated transcripts for S19 and S21 were found. The 32 40S RP genes are generally highly expressed and together they account for 5.33-11.42% of expression depending on the tissues. Expression levels of the RP genes were highly variable both within a single tissue among different RP genes and among tissues with regard to a single RP gene. Taken together, these data strongly suggest post-transcriptional regulation of RP gene expression, particularly in consideration of the stoichiometry of their representation in ribosomes.

Transcriptome analysis of channel catfish (Ictalurus punctatus): initial analysis of gene expression and microsatellite-containing cDNAs in the skin

Previous molecular genetic studies on channel catfish (Ictalurus punctatus) have focused on limited number of genes and gene products. Recent advancement of molecular techniques made high throughput analysis of transcriptomes possible. As part of our transcriptome analysis of channel catfish, we have analyzed 1909 expressed sequence tags (ESTs) derived from a skin library. Of the 1909 ESTs analyzed, 1376 (72.1%) ESTs representing 496 unique genes had homologies with other organisms while 478 (25.0%) ESTs had no significant homologies and were designated as unknown. The remaining 55 (2.9%) EST clones were eliminated because of their low quality or short sequences. Of the 496 unique genes, 327 (65.9%) genes were singletons while 169 (34.1%) genes represented by two or more ESTs. A total of 1007 (52.8%) ESTs representing 235 unique genes matched previously reported channel catfish ESTs while 847 (44.4%) ESTs representing 261 unique genes were newly identified from this research. Functional categorization of the channel catfish genes indicated that the largest group was ribosomal proteins with 65 unique genes represented by 500 clones. The most abundantly expressed gene, the calcium binding protein ictacalcin, accounted for almost 5% of overall expression, indicating its important function in the skin. Sequence analysis of ESTs revealed the presence of 89 microsatellite-containing genes that may be valuable for future mapping studies.