The skeletal muscle alpha-actin gene of channel catfish (Ictalurus punctatus) and its association with piscine specific SINE elements

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

The Regulatory Region of Muscle-Specific Alpha Actin 1 Drives Fluorescent Protein Expression in Olive Flounder Paralichthys olivaceus

To develop a promoter capable of driving transgene expression in non-model fish, we identified and characterized the muscle-specific alpha-actin gene in olive flounder, Paralichthys olivaceus (PoACTC1). The regulatory region of PoACTC1 includes putative regulatory elements such as a TATA box, two MyoD binding sites, three CArG boxes, and a CCAAT box. Microinjection experiments demonstrated that the regulatory region of PoACTC1, covering from -2,126 bp to +751 bp, just prior to the start codon, drove the expression of red fluorescent protein in developing zebrafish embryos and hatching olive flounder. These results suggest that the regulatory region of PoACTC1 may be useful in developing a promoter for biotechnological applications such as transgene expression in olive flounder.

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.

Generation of physical map contig-specific sequences useful for whole genome sequence scaffolding

Along with the rapid advances of the nextgen sequencing technologies, more and more species are added to the list of organisms whose whole genomes are sequenced. However, the assembled draft genome of many organisms consists of numerous small contigs, due to the short length of the reads generated by nextgen sequencing platforms. In order to improve the assembly and bring the genome contigs together, more genome resources are needed. In this study, we developed a strategy to generate a valuable genome resource, physical map contig-specific sequences, which are randomly distributed genome sequences in each physical contig. Two-dimensional tagging method was used to create specific tags for 1,824 physical contigs, in which the cost was dramatically reduced. A total of 94,111,841 100-bp reads and 315,277 assembled contigs are identified containing physical map contig-specific tags. The physical map contig-specific sequences along with the currently available BAC end sequences were then used to anchor the catfish draft genome contigs. A total of 156,457 genome contigs (~79% of whole genome sequencing assembly) were anchored and grouped into 1,824 pools, in which 16,680 unique genes were annotated. The physical map contig-specific sequences are valuable resources to link physical map, genetic linkage map and draft whole genome sequences, consequently have the capability to improve the whole genome sequences assembly and scaffolding, and improve the genome-wide comparative analysis as well. The strategy developed in this study could also be adopted in other species whose whole genome assembly is still facing a challenge.

Determination and quantification of γ-glutamyl-valyl-glycine in commercial fish sauces

It was recently reported that kokumi substances such as glutathione are perceived through the calcium-sensing receptor (CaSR). In addition, screening by the CaSR assay and sensory evaluation revealed that γ-glutamyl-valyl-glycine (γ-Glu-Val-Gly) was a potent kokumi peptide. In this study, the quantities of γ-Glu-Val-Gly in various commercial fish sauces originating from Vietnam (Nuoc Mum), Thailand (Nampra), China (Yu-lu), Korea, Japan (Shottsuru and Ikanago-shoyu), and Italy (Garum) were investigated using a LC/MS/MS method followed by derivatization with 6-aminoquinoyl-N-hydroxysuccinimidyl-carbamate (AQC). The analyses revealed γ-Glu-Val-Gly at concentrations ranging from 0.04 to 1.26 mg/dL, indicating that γ-Glu-Val-Gly is widely distributed among various commercial fish sauces.

A pilot study for channel catfish whole genome sequencing and de novo assembly

Background

Recent advances in next-generation sequencing technologies have drastically increased throughput and significantly reduced sequencing costs. However, the average read lengths in next-generation sequencing technologies are short as compared with that of traditional Sanger sequencing. The short sequence reads pose great challenges for de novo sequence assembly. As a pilot project for whole genome sequencing of the catfish genome, here we attempt to determine the proper sequence coverage, the proper software for assembly, and various parameters used for the assembly of a BAC physical map contig spanning approximately a million of base pairs.

Results

A combination of low sequence coverage of 454 and Illumina sequencing appeared to provide effective assembly as reflected by a high N50 value. Using 454 sequencing alone, a sequencing depth of 18 X was sufficient to obtain the good quality assembly, whereas a 70 X Illumina appeared to be sufficient for a good quality assembly. Additional sequencing coverage after 18 X of 454 or after 70 X of Illumina sequencing does not provide significant improvement of the assembly. Considering the cost of sequencing, a 2 X 454 sequencing, when coupled to 70 X Illumina sequencing, provided an assembly of reasonably good quality. With several software tested, Newbler with a seed length of 16 and ABySS with a K-value of 60 appear to be appropriate for the assembly of 454 reads alone and Illumina paired-end reads alone, respectively. Using both 454 and Illumina paired-end reads, a hybrid assembly strategy using Newbler for initial 454 sequence assembly, Velvet for initial Illumina sequence assembly, followed by a second step assembly using MIRA provided the best assembly of the physical map contig, resulting in 193 contigs with a N50 value of 13,123 bp.

Conclusions

A hybrid sequencing strategy using low sequencing depth of 454 and high sequencing depth of Illumina provided the good quality assembly with high N50 value and relatively low cost. A combination of Newbler, Velvet, and MIRA can be used to assemble the 454 sequence reads and the Illumina reads effectively. The assembled sequence can serve as a resource for comparative genome analysis. Additional long reads using the third generation sequencing platforms are needed to sequence through repetitive genome regions that should further enhance the sequence assembly.

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.

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.

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.

The catfish liver-expressed antimicrobial peptide 2 (LEAP-2) gene is expressed in a wide range of tissues and developmentally regulated

Antimicrobial peptides (AMPs) are important components of the host's innate immune response against microbial invasion. The cysteine-rich AMPs such as defensin and hepcidin have been extensively studied, but the recently identified cysteine-rich liver-expressed antimicrobial peptide 2 (LEAP-2) has been characterized from only a few organisms. Here we cloned and sequenced the LEAP-2 cDNAs from both Channel catfish and Blue catfish. The LEAP-2 gene from Channel catfish was also sequenced and characterized. Channel catfish LEAP-2 gene consists of two introns and three exons that encode a peptide of 94 amino acids with a 28 amino acid signal peptide and a mature peptide of 41 amino acids. The amino acid sequences and gene organization were conserved between catfish and other organisms. The Channel catfish LEAP-2 gene is expressed in a wide range of tissues except brain and stomach. Its expression is developmentally regulated with no detection of mature mRNA in early stages of development. It appears that the catfish LEAP-2 gene is regulated at the level of splicing; it is constitutively transcribed, but remains unspliced until 6 days after hatching. The expression of LEAP-2 was induced in a tissue-specific manner. Its expression was upregulated in the spleen, but not in the liver and head kidney, after challenge with Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish (ESC).

SINEs and LINEs: symbionts of eukaryotic genomes with a common tail

Many SINEs and LINEs have been characterized to date, and examples of the SINE and LINE pair that have the same 3' end sequence have also increased. We report the phylogenetic relationships of nearly all known LINEs from which SINEs are derived, including a new example of a SINE/LINE pair identified in the salmon genome. We also use several biological examples to discuss the impact and significance of SINEs and LINEs in the evolution of vertebrate genomes.

Catfish hepcidin gene is expressed in a wide range of tissues and exhibits tissue-specific upregulation after bacterial infection

Antimicrobial peptides (AMPs) are important components of the host innate immune response against microbial invasion. The cysteine-rich AMPs such as defensin and hepcidin have been extensively studied from various organisms, but their role in disease defense in catfish is unknown. As a first step, we sequenced a hepcidin cDNA from both channel catfish and blue catfish, and characterized the channel catfish hepcidin gene. The channel catfish hepcidin gene consists of two introns and three exons that encode a peptide of 96 amino acids. The amino acid sequences and gene organization were conserved between catfish and other organisms. In contrast to its almost exclusive expression in the liver in humans, the channel catfish hepcidin gene was expressed in a wide range of tissues except brain. Its expression was detected early during embryonic and larval development, and induced after bacterial infection with Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish (ESC) in a tissue-specific manner. The upregulation was observed in the spleen and head kidney, but not in the liver. The expression of hepcidin was upregulated 1--3 days after challenge, but returned to normal levels at 7 days after challenge. The expression profile of the catfish hepcidin gene during the course of bacterial infection mirrors those of inflammatory proteins such as chemokines, suggesting an important role for hepcidin during inflammatory responses.

Tissue engineering of fish skin: behavior of fish cells on poly(ethylene glycol terephthalate)/poly(butylene terephthalate) copolymers in relation to the composition of the polymer substrate as an initial step in constructing a robotic/living tissue hybrid

This study presents the development of a biosynthetic fish skin to be used on aquatic robots that can emulate fish. Smoothness of the external surface is desired in improving high propulsive efficiency and maneuvering agility of autonomous underwater vehicles such as the RoboTuna (Triantafyllou, M., and Triantafyllou, G. Sci. Am. 272, 64, 1995). An initial step was to determine the seeding density and select a polymer for the scaffolds. The attachment and proliferation of chinook salmon embryo (CHSE-214) and brown bullhead (BB) cells were studied on different compositions of a poly(ethylene glycol terephthalate) (PEGT) and poly(butylene terephthalate) (PBT) copolymer (Polyactive). Polymer films were used, cast of three different compositions of PEGT/PBT (weight ratios of 55/45, 60/40, and 70/30) and two different molecular masses of PEGT (300 and 1000 Da). When a 55 wt% and a 300-Da molecular mass form of PEGT was used, maximum attachment and proliferation of CHSE-214 and BB cells were achieved. Histological studies and immunostaining indicate the presence of collagen and cytokeratins in the extracellular matrix formed after 14 days of culture. Porous scaffolds of PEGT/PBT copolymers were also used for three-dimensional tissue engineering of fish skin, using BB cells. Overall, our results indicate that fish cells can attach, proliferate, and express fish skin components on dense and porous Polyactive scaffolds.

Comparing genomic and histologic correlations to radiographic changes in tumors: a murine SCC VII model study

RATIONALE AND OBJECTIVES

To investigate the correlation between the temporal changes in T1- and T2-weighted contrast-enhanced magnetic resonance imaging (MRI), histologic evaluation, and genomic analysis using oligonucleotide microarrays in a murine squamous cell carcinoma tumor models.

MATERIALS AND METHODS

The squamous cell carcinoma (SCC VII) cell line was used to initiate subcutaneous tumors in mice. This mouse model has been used as a model for human head and neck carcinomas. Animals were imaged using contrast enhanced MRI (CE-MRI). Different stages of tumor growth were defined based on changes in the T1- and T2-weighted MRI patterns. The contrast enhancing (CE) and nonenhancing (NE) regions of the tumors were marked and biopsied for oligonucleotide microarray and histologic analysis. Tumors with no differential contrast enhancement were used as controls.

RESULTS

Distinct temporal stages of tumor progression can be defined using both T1- and T2-weighted CE-MRI and microarray analysis. The early stage tumors show a homogeneous contrast enhancement pattern in the T1- and T2-weighted images with no significant differential gene expression from the center and periphery of the tumor. The more advanced tumors that show discrete regions of contrast enhancement in the post-contrast T1-weighted MRIs and tissues from the CE and NE regions show distinctly differential gene expression profiles. Histologic analysis (hematoxylin-eosin stain) showed that the samples obtained from the periphery and center of the early stage tumors and the CE and NE regions from these more advanced tumors were similar. The gene expression profiles of late-stage tumors that showed changes in T2-weighted MRI signal intensity were consistent with tissue degradation in the NE region, which also showed characteristic signs of tissue necrosis in histologic analysis.

CONCLUSION

These results show that temporal changes in T1- and T2-weighted CE-MRI are related to distinct gene expression profiles, and histologic analysis may not be sufficient to detect these detailed changes. As tumors progress, discrete regions of post-contrast T1 enhancement are identified; these regions have distinct gene expression patterns despite similar histologic features. In late-stage tumors, regions of T2 signal changes are observed which correspond with tissue necrosis.

An AFLP-Based Genetic Linkage Map of Channel Catfish (Ictalurus punctatus) Constructed by Using an Interspecific Hybrid Resource Family

Catfish is the major aquaculture species in the United States. The hybrid catfish produced by crossing channel catfish females with blue catfish males exhibit a number of desirable production traits, but their mass production has been difficult. To introduce desirable genes from blue catfish into channel catfish through introgression, a genetic linkage map is helpful. In this project, a genetic linkage map was constructed using amplified fragment length polymorphism (AFLP). A total of 607 AFLP markers were analyzed using 65 primer combinations and an interspecific backcross resource family. A total of 418 AFLP markers were assigned to 44 linkage groups. Among the remaining 189 markers, 101 were not used because of significant segregation distortion, 29 were unlinked, and 59 were eliminated because they span very large distances. The 418 AFLP markers covered 1593 cM Kosambi. The AFLP markers showed a high level of clustering that appears to be related to certain primer combinations. This linkage map will serve as the basis for mapping a greater number of markers to provide a map with high enough resolution for it to be useful for selective breeding programs using introgression.

Molecular characterization and differential expression of the myostatin gene in channel catfish (Ictalurus punctatus)

Myostatin is a recently discovered gene that inhibits muscle growth. In the present study, we characterized the myostatin locus and its expression in channel catfish (Ictalurus punctatus). The genomic DNA and cDNA encoding the channel catfish myostatin were cloned and sequenced. The myostatin gene has three exons encoding a protein of 389 amino acids. Comparison of the genomic sequences with those of the cDNA revealed that the myostatin cDNA was 1673 base pair (bp) long with a 5'-untranslated region (UTR) and 3'-UTR of 180 and 323 bp, respectively. The deduced amino acid sequences of the catfish myostatin is highly conserved with those of other organisms. The myostatin locus is highly polymorphic in channel catfish because of the presence of several microsatellites and single nucleotide polymorphic sites. The myostatin gene was expressed in various tissues and developmental stages at differential levels, suggesting complex regulation of this gene and perhaps roles for myostatin in addition to those originally suggested.

V-SINEs: a new superfamily of vertebrate SINEs that are widespread in vertebrate genomes and retain a strongly conserved segment within each repetitive unit

We have identified a new superfamily of vertebrate short interspersed repetitive elements (SINEs), designated V-SINEs, that are widespread in fishes and frogs. Each V-SINE includes a central conserved domain preceded by a 5'-end tRNA-related region and followed by a potentially recombinogenic (TG)(n) tract, with a 3' tail derived from the 3' untranslated region (UTR) of the corresponding partner long interspersed repetitive element (LINE) that encodes a functional reverse transcriptase. The central domain is strongly conserved and is even found in SINEs in the lamprey genome, suggesting that V-SINEs might be approximately 550 Myr old or older in view of the timing of divergence of the lamprey lineage from the bony fish lineage. The central conserved domain might have been subject to some form of positive selection. Although the contemporary 3' tails of V-SINEs differ from one another, it is possible that the original 3' tail might have been replaced, via recombination, by the 3' tails of more active partner LINEs, thereby retaining retropositional activity and the ability to survive for long periods on the evolutionary time scale. It seems plausible that V-SINEs may have some function(s) that have been maintained by the coevolution of SINEs and LINEs during the evolution of vertebrates.

A microsatellite-based genetic linkage map for channel catfish, Ictalurus punctatus

Microsatellite loci were identified in channel catfish gene sequences or random clones from a small insert genomic DNA library. Outbred populations of channel catfish contained an average of eight alleles per locus and an average heterozygosity of 0.70. A genetic linkage map of the channel catfish genome (N = 29) was constructed from two reference families. A total of 293 microsatellite loci were polymorphic in one or both families, with an average of 171 informative meioses per locus. Nineteen type I loci, 243 type II loci, and one EST were placed in 32 multipoint linkage groups covering 1958 cM. Nine more type II loci were contained in three two-point linkage groups covering 24.5 cM. Twenty-two type II loci remained unlinked. Multipoint linkage groups ranged in size from 11.9 to 110.5 cM with an average intermarker distance of 8.7 cM. Seven microsatellite loci were closely linked with the sex-determining locus. The microsatellite loci and genetic linkage map will increase the efficiency of selective breeding programs for channel catfish.