The annotation of repetitive elements in the genome of channel catfish (Ictalurus punctatus).

Channel catfish (Ictalurus punctatus) is a highly adaptive species and has been used as a research model for comparative immunology, physiology, and toxicology among ectothermic vertebrates. It is also economically important for aquaculture. As such, its reference genome was generated and annotated with protein coding genes. However, the repetitive elements in the catfish genome are less well understood. In this study, over 417.8 Megabase (MB) of repetitive elements were identified and characterized in the channel catfish genome. Among them, the DNA/TcMar-Tc1 transposons are the most abundant type, making up ~20% of the total repetitive elements, followed by the microsatellites (14%). The prevalence of repetitive elements, especially the mobile elements, may have provided a driving force for the evolution of the catfish genome. A number of catfish-specific repetitive elements were identified including the previously reported Xba elements whose divergence rate was relatively low, slower than that in untranslated regions of genes but faster than the protein coding sequences, suggesting its evolutionary restrictions.

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.

Construction of a high-density, high-resolution genetic map and its integration with BAC-based physical map in channel catfish.

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.

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

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

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.

SNP discovery in wild and domesticated populations of blue catfish, Ictalurus furcatus, using genotyping-by-sequencing and subsequent SNP validation.

Blue catfish, Ictalurus furcatus, are valued in the United States as a trophy fishery for their capacity to reach large sizes, sometimes exceeding 45 kg. Additionally, blue catfish × channel catfish (I. punctatus) hybrid food fish production has recently increased the demand for blue catfish broodstock. However, there has been little study of the genetic impacts and interaction of farmed, introduced and stocked populations of blue catfish. We utilized genotyping-by-sequencing (GBS) to capture and genotype SNP markers on 190 individuals from five wild and domesticated populations (Mississippi River, Missouri, D&B, Rio Grande and Texas). Stringent filtering of SNP-calling parameters resulted in 4275 SNP loci represented across all five populations. Population genetics and structure analyses revealed potential shared ancestry and admixture between populations. We utilized the Sequenom MassARRAY to validate two multiplex panels of SNPs selected from the GBS data. Selection criteria included SNPs shared between populations, SNPs specific to populations, number of reads per individual and number of individuals genotyped by GBS. Putative SNPs were validated in the discovery population and in two additional populations not used in the GBS analysis. A total of 64 SNPs were genotyped successfully in 191 individuals from nine populations. Our results should guide the development of highly informative, flexible genotyping multiplexes for blue catfish from the larger GBS SNP set as well as provide an example of a rapid, low-cost approach to generate and genotype informative marker loci in aquatic species with minimal previous genetic information.

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.

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.

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.

Characterization of a mannose-binding lectin from channel catfish (Ictalurus punctatus).

Mannose-binding lectin (MBL) is an important component of innate immunity capable of activating the lectin pathway of the complement system. A MBL gene was isolated from channel catfish (Ictalurus punctatus). The deduced protein contains a canonical collagen-like domain, a carbohydrate recognition domain (CRD), and a neck region similar to mammalian mannose-binding lectin. The catfish mannose-binding lectin CRD contains the EPN motif shown previously to mediate mannose specificity. The catfish mannose-binding lectin showed 30-43% identity with MBL protein sequences of rainbow trout, zebrafish, common carp, and goldfish, and 33-35% identity with sequences of mammalian species. In this study, while liver was the predominant source of mannose-binding lectin gene expression in healthy tissues, mannose-binding lectin expression in spleen rose sharply following challenge with a Gram-negative bacterium.

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

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

Microfibrillar-associated protein 4 (MFAP4) genes in catfish play a novel role in innate immune responses.

The lectin pathway of the complement system is characterized by two groups of soluble pattern recognition molecules, mannose-binding lectins (MBLs) and ficolins. These molecules recognize and bind carbohydrates in pathogens and activate complement leading to opsonization, leukocyte activation, and direct pathogen killing. While MBLs have been reported in many fish species, ficolins do not appear to be present in the teleost lineage, despite their importance in invertebrate and higher vertebrate innate immunity. A protein with a similar fibrinogen-like domain, microfibrillar-associated protein 4, MFAP4, is present in fish, albeit with no described immune function. We examined whether MFAP4 genes in fish may potentially act as pathogen receptors in the absence of ficolin. We isolated and characterized five MFAP4 genes from channel catfish. Linkage mapping and phylogenetic analysis indicated that at least three of the catfish MFAP4 genes are tightly clustered on a single chromosome, suggesting that they may have arisen through tandem duplication. Divergent, duplicated families of MFAP4 genes are also present in other teleost species. Expression analysis of the catfish MFAP4 transcripts revealed unique patterns of homeostatic expression among the genes in gill, spleen, skin, liver, and muscle. Expression of the five MFAP4 transcripts showed significant changes in expression as soon as 4h after infection with either Edwardsiella ictaluri or Flavobacterium columnare with modulation of expression continuing up to 7 d following pathogen exposure. Several different tissues and gene-specific patterns were captured and transcript expression changes of >30-fold were observed over the course of the bacterial challenges. Our results suggest a novel role for MFAP4 in teleost immune responses.

Molecular responses of calreticulin genes to iron overload and bacterial challenge in channel catfish (Ictalurus punctatus).

Infection and inflammation are often accompanied by oxidative stress caused by the accumulation of reactive oxygen species which can be deleterious to the health of the host. Antioxidant defense mechanisms and components are crucial in limiting cellular and tissue-level damage and restoring homeostasis. In mammals, calreticulin is a 46-kDa multifunctional calcium binding protein of the endoplasmic reticulum that has many critical functions in the eukaryotic cell including regulation of intracellular calcium homoeostasis, lectin binding and chaperoning, and oxidative stress responses. In previous studies from our lab, the calreticulin gene was observed to be strongly upregulated in catfish during challenge with infectious Gram-negative bacteria. However, little is known about the function of this gene in teleost fish. The objective of this study, therefore, was to characterize the calreticulin gene from channel catfish, to determine its genomic organization, to profile its patterns of tissue expression, and to establish its potential for physiological antioxidant and immune responses in catfish after bacterial infection with Edwardsiella ictaluri and iron treatment. Our results indicate that there are at least three calreticulin related genes in the catfish genome. The three calreticulin genes are widely expressed in various tissues under homeostatic conditions and their expression showed significant upregulation following infection and/or iron level changes.

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

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

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.

Channel catfish, Ictalurus punctatus, CD4-like molecules.

Two CD4-like (CD4L) molecules have been identified in channel catfish, Ictalurus punctatus. Although phylogenetically related to other vertebrate CD4 molecules, they exhibit only 19% amino acid identity to each other. IpCD4L-1 encodes a predicted protein containing four immunoglobulin domains, a transmembrane region and a cytoplasmic tail containing a p56(lck) binding site. In contrast, IpCD4L-2 encodes for a similar protein with three immunoglobulin domains. The gene organization of IpCD4L-1 is very similar to that of other vertebrate CD4 genes, while the genomic organization of IpCD4L-2 is different. Southern blots indicate both catfish molecules are likely single copy genes and mapping studies show that both are found on a single Bacterial Artificial Chromosome suggesting close linkage. At the message level, IpCD4L-1 and -2 are expressed in various catfish lymphoid tissues and in non-B-cell peripheral blood leukocytes (PBL). Both messages are upregulated in concanavalin A (ConA) and alloantigen stimulated PBL, but not in lipopolysaccharide (LPS)-stimulated cultures. In contrast, they are differentially expressed among the catfish clonal T cell lines. While both molecules appear to be T cell specific, their functional significance in catfish is unknown.

In vivo and in vitro CYP1B mRNA expression in channel catfish.

Our goal was to study the induction of CYP1B mRNA expression in channel catfish (Ictalurus punctatus). CYP1B belongs to the cytochrome P450 superfamily of genes, is involved in the oxidation of endogenous and exogenous compounds, and could potentially be a useful biomarker in fish for exposure to AhR ligands. The full-length catfish CYP1B cDNA is 2417 nt to the polyA tail and encodes a putative protein of 536 amino acids. It has 67% amino acid similarity to carp and zebrafish CYP1B and 68% similarity to carp CYP1B2. Male channel catfish were collected from three Mississippi Delta sites: Lake Roebuck, Itta Bena; Bee Lake, Thornton; and Sunflower River, Indianola. Total RNA was isolated from wild-caught catfish gill, blood, gonad and liver tissues. Quantitative real-time reverse transcriptase PCR was used to determine relative induction of CYP1B in wild catfish compared to laboratory control and BaP-exposed catfish (20mg/kg i.p. after 4 days). BaP exposure significantly induced CYP1B message in blood, gonad, and liver of laboratory catfish. In these same tissues of wild catfish from sites with relatively low sediment contaminants, CYP1B message was not statistically increased relative to laboratory control catfish. CYP1B transcript abundance was higher in gills compared to other tissues in both laboratory and wild catfish. When primary cultured gill cells were treated with increasing concentrations of BaP, TCDD, and PCBs 77, 126 and 169, CYP1B mRNA was induced more than 10-fold while PCB153 and 4,4'DDT did not cause significant CYP1B induction. Our results suggest that catfish CYP1B is induced by the classic AhR ligands.

Analysis of a catfish gene resembling interleukin-8: cDNA cloning, gene structure, and expression after infection with Edwardsiella ictaluri.

Chemokines are important mediators for innate immunity involved in recruitment, activation and adhesion of a variety of leukocyte types to inflammatory foci. While almost all chemokines have been identified from mammals, only a handful of fish chemokines have been identified. Here we report molecular cloning, sequence analysis, and expression of a channel catfish gene resembling interleukin-8 (IL-8). The gene has two alternatively spliced transcripts encoding 114 and 111 amino acids, respectively. The gene has four exons and three introns, typical of the CXC chemokine gene organization. In spite of the structural conservation through evolution, the piscine IL-8 genes showed a much greater sequence divergence than their counterparts among mammals. RT-PCR indicated that both spliced forms were expressed. Expression of the IL-8 like gene was up-regulated 3-5-fold in channel catfish and blue catfish after infection with pathogenic bacteria Edwardsiella ictaluri.