Zhenxia Sha

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.

Quality assessment parameters for EST-derived SNPs from catfish

BackgroundSNPs are abundant, codominantly inherited, and sequence-tagged markers. They are highly adaptable to large-scale automated genotyping, and therefore, are most suitable for association studies and applicable to comparative genome analysis. However, discovery of SNPs requires genome sequencing efforts through whole genome sequencing or deep sequencing of reduced representation libraries. Such genome resources are not yet available for many species including catfish. A large resource of ESTs is to become available in catfish allowing identification of large number of SNPs, but reliability of EST-derived SNPs are relatively low because of sequencing errors. This project was designed to answer some of the questions relevant to quality assessment of EST-derived SNPs.Resultswo factors were found to be most significant for validation of EST-derived SNPs: the contig size (number of sequences in the contig) and the minor allele sequence frequency. The larger the contigs were, the greater the validation rate although the validation rate was reasonably high when the contigs contain four or more EST sequences with the minor allele sequence being represented at least twice in the contigs. Sequence quality surrounding the SNP under test is also crucially important. PCR extension appeared to be limited to a very short distance, prohibiting successful genotyping when an intron was present, a surprising finding.ConclusionStringent quality assessment measures should be used when working with EST-derived SNPs. In particular, contigs containing four or more ESTs should be used and the minor allele sequence should be represented at least twice. Genotyping primers should be designed from a single exon, completely avoiding introns. Application of such quality assessment measures, along with large resources of ESTs, should provide effective means for SNP identification in species where genome sequence resources are lacking.

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

BackgroundThrough the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energys Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification.ResultsA 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.ConclusionsThis 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.

Nod-like subfamily of the nucleotide-binding domain and leucine-rich repeat containing family receptors and their expression in channel catfish.

The NLRs (nucleotide-binding domain and leucine-rich repeat containing family receptors) are a recently identified family of pattern recognition receptors in vertebrates. Several subfamilies of NLRs have been characterized in human, mouse, and zebrafish, but studies of NLRs in other species, especially teleost species, have been lacking. Here we report characterization of five NLRs from channel catfish: NOD1, NOD2, NLRC3, NLRC5, and NLRX1. Structural analysis indicated that the genes were organized in a similar fashion as in the mammals and in zebrafish. Phylogenetic analysis suggested that they were orthologous to the NOD-like subfamily of NLRs. All five NOD-like genes exist as a single copy gene in the catfish genome. Hybridization of gene-specific probes allowed mapping of three NLR genes to the catfish physical map, laying a foundation for genome characterization and for establishing orthologies with NLR genes from other species. These genes are widely expressed in various tissues and leukocyte cell lines. While the majority of the NLR genes appeared to be constitutively expressed, NOD1 was induced after infection with a bacterial pathogen, Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish (ESC), suggesting its involvement in immunity against the intracellular pathogen.

Establishment of a pluripotent embryonic cell line from sea perch (Lateolabrax japonicus) embryos

Abstract Pluripotent embryonic stem (ES) cells provide an efficient approach for genome manipulation with many applications in marine biotechnology and development studies. To develop this technology, we have worked to derive fish ES cells for in vitro studies of embryo cell growth and differentiation and for the generation of transgenic fish. A pluripotent cell line, LJES1, was established from blastula-stage embryos of a cultured marine fish, Lateolabrax japonicus . The LJES1 cells were cultured in Dulbeccos modified Eagles medium (DMEM) supplemented with fetal bovine serum (FBS), marine fish serum, sea perch embryo extract, selenium, basic fibroblast growth factor (bFGF) and leukemia inhibitory factor. The ES-like cells were small and round or polygonal, and grew actively and stably in culture. The cells exhibited a positive alkaline phosphatase activity upon histochemical staining. When the cells were treated with all- trans retinoic acid, differentiation into various types of cells, including neuron-like cells, muscle cells and some unidentified cells were observed, suggesting that the LJES1 cells remained pluripotent in culture. Chromosome analysis revealed that LJES1 cells have a normal diploid karyotype with 2 n =48. Up to now, LJES1 cells have been continuously cultured for more than 150 days with more than 40 passages. High survival rate has been obtained after cryopreservation of cell cultures. GFP reporter gene were transferred into LJES1 cells and successfully expressed.

Pathogen recognition receptors in channel catfish: I. Identification, phylogeny and expression of NOD-like receptors.

Innate immune system plays a significant role in all multicellular organisms. The key feature of the system is its ability to recognize and respond to invading microorganisms. Vertebrates including teleost fish have evolved an array of pathogen recognition receptors (PRRs) for detecting and responding to various pathogen-associated molecular patterns (PAMPs), including Toll-like receptors (TLRs), nucleotide-binding domain, leucine-rich repeat containing receptors (NLRs), and the retinoic acid inducible gene I (RIG-I) like receptors (RLRs). In this study, we identified 22 NLRs including six members of the NLR-A subfamily (NODs), two members of the NLR-B subfamily, 11 members of the NLR-C subfamily, and three genes that do not belong to any of these three subfamilies: Apaf1, CIITA, and NACHT-P1. Phylogenetic analysis indicated that orthologs of the mammalian NOD1, NOD2, NOD3, NOD4, and NOD5 were all identified in catfish. In addition, an additional truncated NOD3-like gene was also identified in catfish. While the identities of subfamily A NLRs could be established, the identities of the NLR-B and NLR-C subfamilies were inconclusive at present. Expression of representative NLR genes was analyzed using RT-PCR and qRT-PCR. In healthy catfish tissues, all the tested NLR genes were found to be ubiquitously expressed in all 11 tested catfish tissues. Analysis of expression of these representative NLR genes after bacterial infection with Edwardsiella ictaluri revealed a significant up-regulation of all tested genes in the spleen and liver, but a significant down-regulation in the intestine and head kidney, suggesting their involvement in the immune responses of catfish against the intracellular bacterial pathogen in a tissue-specific manner. The up-regulation and down-regulation of the tested genes exhibited an amazing similarity of expression profiles after infection, suggesting the co-regulation of these genes.

Transcriptome analysis revealed changes of multiple genes involved in immunity in Cynoglossus semilaevis during Vibrio anguillarum infection

Half-smooth tongue sole (Cynoglossus semilaevis) is one of the most valuable marine aquatic species in Northern China. Given to the rapid development of aquaculture industry, the C. semilaevis was subjected to disease-causing bacteria Vibrio anguillarum. It therefore is indispensable and urgent to understand the mechanism of C. semilaevis host defense against V. anguillarum infection. In the present study, the extensively analysis at the transcriptome level for V. Anguillarum disease in tongue sole was carried out. In total, 94,716 high quality contigs were generated from 75,884,572 clean reads in three libraries (HOSG, NOSG, and CG). 22,746 unigenes were identified when compared with SwissProt, an NR protein database and NT nucleotide database. 954 genes exhibiting the differentially expression at least one pair of comparison in all three libraries were identified. GO enrichment for these genes revealed gene response to biotic stimulus, immune system regulation, and immune response and cytokine production. Further, the pathways such as complement and coagulation cascades and Vibrio cholerae infection pathways were enriched in defensing of pathogen. Besides, 13,428 SSRs and 118,239 SNPs were detected in tongue sole, providing further support for genetic variation and marker-assisted selection in future. In summary, this study identifies several putative immune pathways and candidate genes deserving further investigation in the context of development of therapeutic regimens and lays the foundation for selecting resistant lines of C. semilaevis against V. anguillarum.

Structure and expression of transferrin gene of channel catfish, Ictalurus punctatus.

Transferrin is important in iron metabolism and has been reported to be involved in disease defence responses after bacterial infection. In this study, we identified, sequenced, and characterized the transferrin gene from channel catfish, Ictalurus punctatus. The catfish transferrin gene was similar to those of other vertebrate species with 17 exons and 16 introns. Sequence analysis indicated the presence of the two duplicated lobes, each containing two sub-domains separated by a cleft harboring the iron-binding site, suggesting their structural conservation. The channel catfish transferrin cDNA encodes 679 amino acids with 42-56% similarity to known transferrin genes from various species. Southern blot analysis suggested the presence of two copies of the transferrin gene in the catfish genome, perhaps arranged in a tandem fashion. The catfish transferrin gene was mapped to a catfish BAC-based physical map. The catfish transferrin gene was highly expressed in the liver, but expression was low in most other tested tissues. Transferrin expression was significantly up-regulated after infection with Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish. Such induction was also found with co-injection of iron-dextran and E. ictaluri, while transferrin expression was not significantly induced with the injection of iron-dextran alone.

Molecular cloning, genomic structure, polymorphism and expression analysis of major histocompatibility complex class IIA and IIB genes of half-smooth tongue sole (Cynoglossus semilaevis).

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Its function is to present foreign peptide to the T-cell. In order to study the function and molecular polymorphism of class II genes in teleost, the full lengths of MHC class IIA and IIB cDNA were cloned from half-smooth tongue sole by homology cloning and rapid amplification of cDNA ends polymerase chain reaction (RACE-PCR). Genomic organizations, molecular polymorphism, and expression profiles of class IIA and IIB were examined to study the function in fish. As in other teleost, four exons and three introns were identified in half-smooth tongue sole class IIA gene, five exons and four introns were identified in class IIB gene. The deduced amino acid sequence of class IIA had 27.3-69.8% identity with those of mammal and teleost. Nine class IIA alleles were identified from four individuals. Four different alleles observed in a single individual may infer the existence of two loci at least. The deduced amino acid sequence of class IIB had 7.9-71.9% identity with those of other species. Fifteen class IIB alleles were identified. Six different alleles observed in a single individual may suggest that there are at least three loci in class IIB genes. Real-time quantitative RT-PCR demonstrated that the MHC class IIA and IIB were ubiquitously expressed in twelve normal tissues. Challenge of half-smooth tongue sole with the pathogenic bacteria, Vibrio anguillarum, resulted in significant changes in the expression of MHC IIA and IIB mRNA in three tissues.

Second-Generation Genetic Linkage Map of Catfish and Its Integration with the BAC-Based Physical Map

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