Molecular evolution of connective tissue growth factor in Cyprinidae (Teleostei : Cypriniformes)

Connective tissue growth factor (CTGF) plays an important role in regulation of cell growth, differentiation, apoptosis and individual development in animals. The study of sequences variation and molecular evolution of CTGF gene across various species of the cyprinid could be helpful for understanding of speciation and gene divergence in this kind of fish. In this study, 19 novel sequences of CTGF gene were obtained from the representative species of the family Cyprinidae using PCR amplification, cloning and sequencing. Phylogenetic relationships of Cyprinidae were reconstructed by neighbor-joining (NJ), maximum parsimony (MP), maximum likelihood (ML) and Bayesian method. Oryzias latipes from the family Cyprinodontidae was assigned to be the outgroup taxon. Leuciscini and Barbini were clustered into the monophyletic lineages, respectively, with the high nodal supports. The estimation of the ratio of non-synonymous to synonymous substitution (dN/dS) for the various branches indicated that there stood the different evolution rates between the Leuciscini and the Barbini. With the ratio of dN/dS of the Leuciscini being lower than that of the Barbini, species within the Barbini were demonstrated to be subjected to the relatively less selection pressure and under the relaxable evolution background. A 6 by indel (insertion/deletion) was found at the 5' end of CTGF gene of Cyprinidae, and this 6 by deletion only appeared in the Leuciscini, which is a typical characteristic of the Leuciscini and provides evidence for the monophylogeny of the Leuciscini. For the amino acid sequences of CTGF protein, the most variations and indels were distributed in the signal region and IGFBP region of this protein, implying that these variations were correlated with the regulation of the CTGF gene expression and protein activity. (c) 2007 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Cloning and characterization of interferon stimulated genes Viperin and ISG15, and their promoters from snakehead Channa argus

By suppression subtractive hybridization, rapid amplification of cDNA ends and gene walking methods, interferon stimulated genes (ISGs), Viperin and ISG15, and their promoters have been cloned and characterized from snakehead Channa argus. The Viperin cDNA was found to be 1474 nt and contain an open reading frame (ORF) of 1059 nt that translates into a putative peptide of 352 amino acid (aa). The putative peptide of Viperin shows high identity to that in teleosts and mammals except for the N-terminal 70 aa. The ISG15 cDNA was found to be 758 nt and contain an ORF of 468 nt that translates into a putative peptide of 155 aa. The putative peptide of ISG15 is composed of two tandem repeats of ubiquitin-like (UBL) domains, and a canonical conjugation motif (LRGG) at C-terminal. Viperin and ISG15 promoter regions were characterized by the presence of interferon stimulating response elements (ISRE) and gamma-IFN activation sites (GAS). ISRE is a feature of IFN-induced gene promoter and partially overlaps interferon regulatory factor (IRF) 1 and IRF2 recognition sites. GAS is responsible for the gamma-IFN mediated transcription. One conserved site for NF-kappa B was found in the promoter region of Viperin. This is the first report of conservative binding motif for NF-kappa B in accordance with the consensus sequence (GGGRN-NYYCC) among teleost ISG promoters. Moreover, there were also TATA, CAAT and Sp1 transcription factor sites in Viperin and ISG15 promoters. In 5' untranslated region (UTR), snakehead ISG15 gene contains a single intron, which differs from Viperin gene. The transcripts of Vipeirn and ISG15 mRNA were mainly expressed in head kidney, posterior kidney, spleen and gill. The expression levels in liver were found to increase obviously in response to induction by IFN-inducer poly I : C.

Characterization and expression analysis of TNFR-associated factor 1 (TRAF1) in grass carp Ctenopharyngodon idella

TNF receptor associated factor 1 (TRAF1) plays an important role in regulating the TNF signaling and protecting cells from apoptosis. In the present study, a TRAF1 gene has been cloned from grass carp (Ctenopharyngodon idella) by reverse transcription (RT)-PCR and rapid amplification of cDNA ends (RACE). The full-length cDNA is 2235 bp, including a 250 bp 5' UTR (untranslated region), a 1659 bp open reading frame, and a 326 bp 3'UTR. The polyadenylation signal (AATAAA, AATAA) and one mRNA instability motif (AUUUA) were found followed by a poly (A) tail in the 3'UTR. No signal peptide or transmembrane region has been found in the putative amino acids of grass carp TRAF1 (gcTRAF1). The putative amino acids of gcTRAF1 share 72% identity with the homologue in zebrafish. It is characterized by a zinc finger at the N-terminus and a TRAF domain (contains one TRAF-C and one TRAF-N) at the C-terminus. The identity of the TRAF domain among all the TRAF1 homologues in vertebrates varies from 52% to 58%, while the identities of TRAF-C were almost the same as 70%. The recombinant gcTRAF1 has been constructed successfully and expressed in Escherichia coli by using pET-32a expression vector. The polyclonal antibody for rabbit has been successfully obtained. The expression of gcTRAF1 in different organs was examined by real-time quantitative PCR and Western blotting, respectively. It was widely distributed in heart, head kidney, thymus, brain, gill, liver, spleen, and trunk kidney. This is the first report of TRAF1 homologue molecule found in fish. (c) 2007 Elsevier B.V. All rights reserved.

Increasing the genetic uniformity of bighead carp [Aristichthys nobilis (Richardson)] by means of spontaneous diploidization of gynogenetically activated eggs

Three groups of gynogenetic diploid bighead carp were successfully obtained by means of artificial gynogenesis. The activation rates of gynogenesis varied from 75.9% to 98.8%, and the frequency of spontaneous diploidization was around 0.4%. Over 2000 normally gynogenetic diploid fry were obtained in three gynogenetic groups. The haploid karyotype consisted of nine metacentric, 12 submetacentric, three subtelocentric chromosomes and 45 arms. The chromosome number was 48 from gynogenetic diploid. The results showed that the genetic material of offspring was maternal. The aneuploid hybrid embryos of bighead carp and Xingguo red common carp with chromosome numbers ranging from 28 to 73 did not survive post hatch, likely the result of incompatibility between the nucleus and the cytoplasm of two parents. Sixty RAPD primers from three groups were used for total DNA amplification of gynogenetic offspring, maternal and 'paternal' fish. A total of 451 bands were amplified from three kinds of samples above. From maternal bighead carp, 256 bands were amplified; however, there were 251 shared bands between maternal and gynogenetic bighead carp. From artificial gynogenetic offspring, two 'paternal' DNA segments without an expression function were found. An UPGMA tree showed that gynogenetic offspring were closely clustered and the genetic identity among them was very high (0.956).

Molecular identification and expression analysis of tumor necrosis factor receptor-associated factor 2 in grass carp Ctenopharyngodon idella

Tumor necrosis factor receptor-associated factor 2 (TRAF2) is a crucial component of almost the entire tumor necrosis factor receptor superfamily signaling pathway. In the present study, a TRAF2 gene has been cloned from grass carp (Ctenopharyngodon idella) by reverse transcription-polymerase chain reaction and rapid amplification of cDNA ends. The full-length cDNA is 3162 bp, including a 60 bp 5' untranslated region (UTR), a 1611 bp open reading frame, and a 1491 bp 3' UTR. The polyadenylation signal (AATAAA) and the mRNA instability motifs (ATTTTA, ATTTA) were followed by a poly(A) tail in the 3' UTR. No signal peptide or transmembrane region has been found in the putative amino acids of grass carp TRAF2 (gcTRAF2). Phylogenetic tree analysis clearly showed that gcTRAF2 is nearest to the TRAF2 gene of goldfish. The identity of gcTRAF2 with its homologs in other vertebrates ranges from 56% to 97%. It is characterized by one RING-type signature at the N-terminus, one zinc finger in the middle part, and one conserved TRAF domain consisting of a C-proximal (TRAF-C) subdomain and a N-proximal (TRAF-N) subdomain. The identity of TRAF-C among all TRAF2 homologs in vertebrates varies from 78% to 97%, whereas the identity of TRAF-N ranges from 56% to 100%. The recombinant gcTRAF2 has been expressed in Escherichia coli using pET-32a expression vector. The rabbit anti-gcTRAF2 polyclonal antibody was obtained. The expression of gcTRAF2 in different organs was examined by real-time quantitative polymerase chain reaction and Western blot analysis. It was widely distributed in heart, head kidney, thymus, brain, gill, liver, spleen, and trunk kidney. This is the first report of a TRAF2 homolog molecule in fish.

Phylogenetic relationships of Cyprinidae (Teleostei : Cypriniformes) inferred from the partial S6K1 gene sequences and implication of indel sites in intron 1

The family Cyprinidae is widely distributed in East Asia, and has the important phylogenetic significance in the fish evolution. In this study, the 5' end partial sequences (containing exon 1, exon 2 and indel 1) of S6K1 gene were obtained from 30 representative species in Cyprinidae and outgroup using PCR amplification and sequencing. The phylogenetic relationships of Cyprinidae were reconstructed with neighbor joining (NJ), maximum parsimony (MP), maximum likelihood (ML), and Bayesian methods. Myxocyprinus asiaticus (Catostomidae) was assigned to the outgroup taxon. Similar phylogenetic relationships within the family Cyprinidae were achieved with the four analyses. Leuciscini and Barbini were monophyletic lineages respectively with the high nodal supports. Leuciscini comprises Hypophthalmichthyinae, Xenocyprinae, Cultrinae, Gobioninae, Acheilognathinae and East Asian species of Leuciscinae and Danioninae. Monophyly of East Asian clade was supported with high nodal support. Barbini comprises Schizothoracinae, Barbinae, Cyprininae and Labeoninae. The monophyletic lineage consisting of Danio rerio, D. myersi, and Rasbora trilineata was basal in the tree. In addition, the large fragment indels in intron 1 were analyzed to improve the understanding of Cyprinidae relationships. The results showed that the large fragment indels were correlated with the relations among species. Some conserved regions in intron 1 were thought to be involved in the functional regulation. However, no correlation was found between sequence variations and species characteristic size.

Polymorphic microsatellites in largemouth bronze gudgeon (Coreius guichenoti) developed from repeat-enriched libraries and cross-species amplifications

Largemouth bronze gudgeon (Coreius guichenoti) is a medium-sized fish endemic from the upper Yangtze River of China and its survival is threatened by the construction of the Three Gorges Dam. This study reports 20 new polymorphic microsatellites from a repeat-enriched genomic library with a mean number allele of 5.2, and observed and expected heterozygosities ranging from 0.035 to 1, and from 0.13 to 0.917, respectively. In a cross-species amplification test, nine of the 37 tested loci were found to be also polymorphic in a congeneric species, brass gudgeon (C. heterodon). In addition, other four loci from common carp (Cyprinus carpio) were also polymorphic in C. guichenoti. Out of these 24 polymorphic microsatellites, only three loci significantly deviated from Hardy-Weinberg equilibrium in the sampled population (P < 0.0025), and all pairwise tests for linkage disequilibrium among loci were nonsignificant after applying sequential Bonferroni correction (P > 0.0026). These novel microsatellites provide sufficient levels of polymorphism for studies on population genetics and conservation in C. guichenoti and its related species.

Development of novel EST-SSR markers in common carp by data mining from public EST sequences

Expressed sequence tags (ESTs) are a source for microsatellite development. In the present study, EST-derived microsatelltes (EST-SSRs) were generated and characterized in the common carp (Cyprinus carpio) by data mining from updated public EST databases and by subsequent testing for polymorphism. About 5.5% (555) of 10,088 ESTs contain repeat motifs of various types and lengths with CA being the most abundant dinucleotide one. Out of the 60 EST-SSRs for which PCR primers were designed, 25 loci showed polymorphism in a common carp population with the alleles per locus ranging from 3 to 17 (mean 7). The observed (H-O) and expected (HE) heterozygosities of these EST-SSRs were 0.13-1.00 and 0.12-0.91, respectively. Six EST-SSR loci significantly deviated from the Hardy-Weinberg equilibrium (HWE) expectation, and the remaining 19 loci were in HWE. Of the 60 primer sets, the rates of polymorphic EST-SSRs were 42% in common carp, 17% in crucian carp (Carassius auratus), and 5% in silver carp (Hypophthalmichthys molitrix), respectively. These new EST-SSR markers would provide sufficient polymorphism for population genetic studies and genome mapping of the common carp and its closely related fishes. (c) 2007 Published by Elsevier B.V.

Isolation and characterization of polymorphic microsatellite loci in Wuchang bream (Megalobrama amblycephala)

Wuchang bream (Megalobrama amblycephala) is an economically important fish in China. From a (GT)(13)-enriched genomic library, 20 microsatellites were developed. Nine of these 20 loci were polymorphic in a test population with allele numbers ranging from two to four, and the observed and expected heterozygosities ranging from 0.2609 to 0.7826 and from 0.3739 to 0.7546, respectively. In the cross-species amplifications, six of these nine loci were also polymorphic in white amur bream (Parabramis pekinensis). These polymorphic microsatellite loci are potentially useful for population genetics of Wuchang bream and its closely related species.

The c-myc coding DNA sequences of cyprinids (Teleostei : Cypriniformes): Implications for phylogeny

The family Cyprinidae is one of the largest fish families in the world, which is widely distributed in East Asian, with obvious difference in characteristic size among species. The phylogenetic analysis of cyprinid taxa based on the functionally important genes can help to understand the speciation and functional divergench of the Cyprinidae. The c-myc gene is an important gene regulating individual growth. In the present study, the sequence variations of the cyprinid c-myc gene and their phylogenetic significance were analyzed. The 41 complete sequences of the c-myc gene were obtained from cyprinids and outgroups through PCR amplification and clone. The coding DNA sequences of the c-myc gene were used to infer molecular phylogenetic relationships within the Cyprinidae. Myxocyprinus asiaticus (Catostomidae), Misgurnus anguillicaudatus (Cobitidae) and Hemimyzon sinensis (Homalopteridae) were assigned to the outgroup taxa. Phylogenetic analyses using maximum parsimony (MP), maximum likelihood (ML), and Bayesian retrieved similar topology. Within the Cyprinidae, Leuciscini and Barbini formed the monophyletic lineage respectively with high nodal supports. Leuciscini comprises Xeno-cyprinae, Cultrinae, East Asian species of Leuciscinae and Danioninae, Gobioninae and Acheilognathinae, and Barbini contains Schizothoracinae, Barbinae, Cyprininae and Labeoninae. Danio rerio, D. myersi and Rasbora trilineata were supposed to separate from Leuciscinae and Barbini and to form another lineage, The positions of some Danioninae species were still unresolved. Analyses of both amino acid variation with parsimony information and two high variation regions indicated that there is no correlation between variations of single amino acid or high variation regions and characteristic size of cyprinids. In,addition, the species with smaller size were usually found to be basal within clades in the tree, which might be the results of the adaptation to the primitive ecology and survival pressure.

Construction and characterization of a Lipotes vexillifer genomic DNA BAC library

We constructed a genomic DNA library for Lipotes vexillifer (L. vexillifer), the Baiji or Yangtze River dolphin, one of the most endangered mammals in the world. The library consists of 149,000 BAC clones, with an average insert size of 83 kb, representing approximately 3.4 haploid genome equivalents. PCR amplification of four known L. vexillifer genes yielded two to four positive clones each. To demonstrate the utility of this library, we isolated and sequenced the L. vexillifer alpha lactalbumin gene, which is a gene specific to mammals and one which has been widely used as molecular tool in phylogenetic analysis. We also end-sequenced 20 randomly selected clones, resulting in the identification of at least five new L. vexilliter genes, five SSR loci, and one SINE locus. These results suggest that this library is a valuable resource for candidate gene cloning, physical mapping, and genome sequencing of this important and threatened species.

Molecular characterization of common carp (Cyprinus carpio) Sonic Hedgehog and discovery of its maternal expression

Sonic hedgehog (Shh), one of important homologous members of the hedgehog (Hh) family in vertebrates, encodes a signaling molecule that is involved in short- or long-range patterning processes during embryogenesis. In zebrafish, maternal activity of Hh was found to be contributing to the formation of primary motoneurons. However, we found that all of the known Hh members were not maternally expressed in zebrafish. In the present study, full-length cDNA of common carp (Cyprinus carpio) Shh (cShh) was gained by degenerate reverse-transcription PCR (RT-PCR) and rapid amplification of cDNA ends. Sequence comparison shows that cShh coding sequence shares 93.4% identity with zebrafish Shh coding sequence, and their corresponding protein sequences have 91.9% similarity. Comparative analysis of Shh genomic sequences and Hh protein sequences from different species revealed that the genomic structures of Hh are conserved from invertebrate to vertebrate. In contrast to zebrafish Shh, cShh transcripts were detectable from one-cell stage by RT-PCR analysis. Whole mount in situ hybridization verified the maternal expression of Shh in common carp, which is, to our knowledge, the first report of that in vertebrates, suggesting that Shh might be responsible for the maternal Hh activity in common carp.

Isolation and characterization of a rhabdovirus from co-infection of two viruses in mandarin fish

Co-infection of two viruses has been observed in mandarin fish (Siniperca chuatsi), but the two viruses have not been characterized. In this study, a rhabdovirus has been isolated from the co-infected two viruses extracted from the diseased mandarin fish, and its morphological structure and partial biochemical and biophysical characteristics have been observed and analyzed. The isolated rhabdovirus has a typical bullet shape, and is therefore called S. chttatsi rhabdovirus (SCRV). And, the isolated rhabdovirus produced a higher titer (10(8.5) TCID50 ml(-1)) than did the co-infecting viruses (10(6.5) TCID50 ml(-1)). Subsequently, the viral genome RNA was extracted, and used as template to clone the complete nucleoprotein (N) gene by RT-PCR amplification. Cloning and sequencing of the SCRV N protein revealed 42%-31% amino acid identities to that of trout rhabdovirus 903/87 and the rhabdoviruses in genus Vesiculovirus. SDS-PAGE separation of the isolated SCRV and other two rhabdoviruses also revealed obvious polypeptide profile difference. Moreover, the anti-SCRV N protein antibody was prepared, and the anti-SCRV N protein antibody only could recognize the SCRV N protein, whereas no antigenicity was detected in other two rhabdoviruses. The data suggested that the SCRV should be a rhabdovirus member related to the genus Vesiculovirus in the Rhabdoviridae. (c) 2006 Elsevier B.V. All rights reserved.

Molecular cloning of growth hormone receptor (GHR) from common carp (Cyprinus carpio L.) and identification of its two forms of mRNA transcripts

The cDNA of growth hormone receptor (GHR) was cloned from the liver of 2-year common carp (Cyprinus carpio L.) by reverse transcription-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA end (RACE). Its open reading frame (ORF) of 1806 nucleotides is translated into a putative peptide of 602 amino acids, including an extracellular ligand-binding domain of 244 amino acids (aa), a single transmembrane domain of 24 aa and an intracellular signal-transduction domain of 334 aa. Sequence analysis indicated that common carp GHR is highly homologous to goldfish (Carassius auratus) GHR at both gene and protein levels. Using a pair of gene-specific primers, a GHR fragment was amplified from the cDNA of 2-year common carp, a 224 bp product was identified in liver and a 321 bp product in other tissues. The sequencing of the products and the partial genomic DNA indicated that the difference in product size was the result of a 97 bp intron that alternatively spliced. In addition, the 321 bp fragment could be amplified from all the tissues of 4-month common carp including liver, demonstrating the occurrence of the alternative splicing of this intron during the development of common carp. Moreover, a semi-quantitative RT-PCR was performed to analyze the expression level of GHR in tissues of 2-year common carp and 4-month common carp. The result revealed that in the tissues of gill, thymus and brain, the expression level of GHR in 2-year common carp was significantly tower than that of 4-month common carp.

Molecular characterization of three Rana grylio virus (RGV) isolates and Paralichthys olivaceus lymphocystis disease virus (LCDV-C) in iridoviruses

Three Rana grylio virus (RGV) isolates and lymphocystis disease virus (LCDV-C) were molecularly characterized by antigenicity comparison, Western blot detection of viral polypeptides, restriction fragment length polymorphism analysis of viral genomes, and MCP sequence analysis. Significant antigenicity differences existed among the three RGV isolates and LCDV-C. Western blot detection indicated that the viral polypeptides of three RGV isolates could be recognized by the anti-RGV9807 serum, whereas no bands were observed in the LCDV-C, and significant differences exist among the band patterns of three RGV isolates. Restriction fragment length polymorphism (RFLP) analysis was performed by digesting genomic DNA of the four iridovirus isolates with restriction endonucleases HindIII, KpnI, XbaI and BamHI. On the whole, obvious discrepancies existed between LCDV-C and RGV isolates, and some significant band pattern differences were also revealed between RGV9808 and RGV9506 (or RGV9807) in the profiles of restriction endonucleases Xbal, Kpn I and BamHI. PCR amplification and sequence analysis of MCP gene sequence further revealed their phylogenetic relationship among the three RGV isolates, LCDV-C and other iridoviruses. RGV9506, RGV9807 and RGV9808 are clustered together with other ranaviruses, such as FV3, BIV, TFV and ENHV, although the RGV9808 is more close to EHNV than to other ranaviruses. Additionally, LCDV-C is clustered with LCDV-1, the type species of genus Lymphocystisvirus. The current study provides clear evidence that significant genetic difference exists among the three RGV isolates. Therefore, further work on comparative genomic studies will contribute significantly to understanding of their taxonomic position and pathological mechanism. (C) 2005 Elsevier B.V. All rights reserved.