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.

Molecular cloning, promoter analysis and induced expression of the complement component C9 gene in the grass carp Ctenopharyngodon idella

Complement-mediated killing of pathogens through lytic pathway is an important effector mechanism of innate immune response. C9 is the ninth member of complement components, creating the membrane attack complex (MAC). In the present study, a putative cDNA sequence encoding the 650 amino acids of C9 and its genomic organization were identified in grass carp Ctenopharyngodon idella. The deduced amino acid sequence of grass carp C9 (gcC9) showed 48% and 38.5% identity to Japanese flounder and human C9, respectively. Domain search revealed that gcC9 contains a LDL receptor domain, an EGF precursor domain, a MACPF domain and two TSP domain located in the N-terminal and C-terminal, respectively. Phylogenetic analysis demonstrated that gcC9 is clustered in a same clade with Japanese flounder, pufferfish and rainbow trout C9. The gcC9 gene consists of 11 exons with 10 introns, spacing over approximately 7 kb of genomic sequence. Analysis of gcC9 promoter region revealed the presence of a TATA box and some putative transcription factor such as C/EBP, HSF, NF-AT, CHOP-C, HNF-3B, GATA-2, IK-2, EVI- 1, AP-1, CP2 and OCT-1 binding sites. The first intron region contains C/EBPb, HFH-1 and Oct-1 binding sites. RT-PCR and Western blotting analysis demonstrated that the mRNA and protein of gcC9 gene have similar expression patterns, being constitutively expressed in all organs examined of healthy fish, with the highest level in hepatopancreas. By real-time quantitative RT-PCR analysis, gcC9 transcripts were significantly up-regulated in head kidney, spleen, hepatopancreas and down-regulated in intestine from inactivated fish bacterial pathogen Flavobacterium columnare-stimulated fish, demonstrating the role of C9 in immune response. (c) 2007 Elsevier B.V. All rights reserved.

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.

Molecular evidence for the monophyly of East Asian groups of Cyprinidae (Teleostei : Cypriniformes) derived from the nuclear recombination activating gene 2 sequences

The family Cyprinidae is one of the largest families of fishes in the world and a well-known component of the East Asian freshwater fish fauna. However, the phylogenetic relationships among cyprinids are still poorly understood despite much effort paid on the cyprinid molecular phylogenetics. Original nucleotide sequence data of the nuclear recombination activating gene 2 were collected from 109 cyprinid species and four non-cyprinid cypriniform outgroup taxa and used to infer the cyprinid phylogenetic relationships and to estimate node divergence times. Phylogenetic reconstructions using maximum parsimony, maximum likelihood, and Bayesian analysis retrieved the same clades, only branching order within these clades varied slightly between trees. Although the morphological diversity is remarkable, the endemic cyprinid taxa in East Asia emerged as a monophyletic clade referred to as Xenocypridini. The monophyly for the subfamilies including Cyprininae and Leuciscinae, as well as the tribes including Labeonini, Gobionini, Acheilognathini, and Leuciscini, was also well resolved with high nodal support. Analysis of the RAG2 gene supported the following cyprinid molecular phylogeny: the Danioninae is the most basal subfamily within the family Cyprinidae and the Cyprininae is the sister group of the Leuciscinae. The divergence times were estimated for the nodes corresponding to the principal clades within the Cyprinidae. The family Cyprinidae appears to have originated in the mid-Eocene in Asia, with the cladogenic event of the key basal group Danioninae occurring in the early Oligocene (about 31-30 MYA), and the origins of the two subfamilies, Cyprininae and Leuciscinae, occurring in the mid-Oligocene (around 26 MYA). (c) 2006 Elsevier Inc. 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.

NK-lysin of channel catfish: Gene triplication, sequence variation, and expression analysis

Antimicrobial peptides (AMPs) are important components of the host innate immune response against microbial invasion. In addition to the previously known four classes of antimicrobial peptides, a fifth class of antimicrobial peptides has been recently identified to include NK-lysins that have a globular three-dimensional structure and are larger with 74-78 amino acid residues. NK-lysin has been shown to harbor antimicrobial activities against a wide spectrum of microorganisms including bacteria, fungi, protozoa, and parasites. To date, NK-lysin genes have been reported from only a limited number of organisms. We previously identified a NK-lysin cDNA in channel catfish. Here we report the identification of two noveltypes of NK-lysin transcripts in channel catfish. Altogether, three distinct NK-lysin transcripts exist in channel catfish. In this work, their encoding genes were identified, sequenced, and characterized. We provide strong evidence that the catfish NK-lysin gene is tripled in the same genomic neighborhood. All three catfish NK-lysin genes are present in the same genomic region and are tightly linked on the same chromosome, as the same BAC clones harbor all three copies of the NK-lysin genes. All three NK-lysin genes are expressed, but exhibit distinct expression profiles in various tissues. In spite of the existence of a single copy of NK-lysin gene in the human genome, and only a single hit from the pufferfish,genome, there are two tripled clusters of NK-lysin genes on chromosome 17 of zebrafish in addition to one more copy on its chromosome 5. The similarity in the genomic arrangement of the tripled NK-lysin genes in channel catfish and zebrafish suggest similar evolution of NK-lysin genes. (c) 2005 Elsevier Ltd. All rights reserved.

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.

Molecular identification and expression analysis of natural resistance associated macrophage protein (Nramp) cDNA from Japanese flounder (Paralichthys olivaceus)

Natural resistance associated macrophage protein (Nramp) controls partially innate resistance to intracellular parasites. Its function is to enhance the ability of macrophages to kill pathogens. However, little is known about the structure and function of Nramp in lower vertebrates such as teleosts. We have recently isolated a cDNA encoding Nramp from Japanese flounder (Paratichthys olivaceus). The full-length cDNA of the Nramp is 3066 bp in length, including 224 bp 5' terminal UTR, 1662 bp encoding region and 1180 bp 3' terminal UTR. The 1662-nt open reading frame was found to code for a protein with 554 amino acid residues. Comparison of amino acid sequence indicated that Japanese flounder Nramp consists of 12 transmembrane (TM) domains. A consensus transport motif (CTM) containing 20 residues was observed between transmembrane domains 8 and 9. The deduced amino acid sequence of Japanese flounder had 77.30%, 82.71%, 82.67%, 79.64%, 80.72%, 90.97%, 91.16%, 60.14%, 71.48%, 61.69%, 72.37% identity with that of rainbow trout Nramp alpha and beta, channel catfish Nramp, fathead minnow Nramp, common carp Nramp, striped sea bass Nramp, red sea bream Nramp, mouse Nramp 1 and 2, human Nramp 1 and 2, respectively. RT-PCR indicated that Nramp transcripts were highly abundant in spleen, head kidney, abundant in intestine, liver and gill, and less abundant in heart. The level of Nramp mRNA in embryos gradually increases during embryogenesis from 4 h (8 cell stage) to 80 h (hatched stage) after fertilization. (c) 2005 Elsevier Ltd. All rights reserved.

Molecular phylogenetics of Gymnocypris (Teleostei : Cyprinidae) in Lake Qinghai and adjacent drainages

149 complete mitochondrial DNA (mtDNA) cytochrome b (Cyt b) genes (1140 bp) of Gymnocypris przewalskii, Gymnocypris eckloni and Gymnocyptis scolistomus from the Lake Qinghai, Yellow River and Qaidam Basin were sequenced and analyzed. Consistent dendrogram indicated that the samples collected from the same species do not constitute a separate monophyletic group and all the samples were grouped into three highly divergent lineages (A, B and C). Among them, Lineage A contained all samples of G. przewalskii from the Lake Qinghai and partial samples of the G. eckloni from the Yellow River. Lineage B contained the remaining samples of G. eckloni from the Yellow River. Lineage C was composed of a monophyletic group by G. eckloni from the Qaidam Basin. Analysis of molecular variance (AMOVA) indicated that most of genetic variations were detected within these three mtDNA lineages (93.12%), suggesting that there are three different lineages of Gymnocypris in this region. Our Cyt b sequence data showed that G. przewalskii was not a polytypic species, and G. scolistomus was neither an independent species nor a subspecies of G. eckloni. The divergent mtDNA lineages of G. eckloni from the Yellow River suggested that gene flow between the different populations was restricted to a certain extent by several gorges on the upper reach of the Yellow River. Lineage B of G. eckloni might be the genetic effect from the ancestor which was incorporated with the endemic schizothoracine fishes when the headward erosion of the Yellow River reached to its current headwaters of late. The G. eckloni from Basin Qaidam was a monophyletic group (lineage C) and F-st values within G. eckloni from the Yellow River were higher than 0.98, suggesting that the gene flow has been interrupted for a long time and the G. eckloni from Basin Qaidam might have been evolved into different species by ecology segregation. The correlation between the rakers number of Gymnocypris and population genetic variation was not significant. All Gymnocypris populations exhibited a low nucleotide diversity (pi = 0.00096-0.00485). Therefore the Gymnocyptis populations from Basin Qaidam could have experienced severe bottleneck effect in history. Our result suggested Gym-nocypris populations of Basin Qaidam should give a high priority in conservation programs.

Molecular cloning and characterization of the copper/zinc and manganese superoxide dismutase genes from the human parasite Clonorchis sinensis

A copper/zinc superoxide dismutase (Cu/ZnSOD) gene and a manganese superoxide dismutase (MnSOD) gene of the human parasite Clonorchis sinensis have been cloned and their gene products functionally characterized. Genes Cu/ZnSOD and MnSOD encode proteins of 16 kDa and 25.4 kDa, respectively. The deduced amino acid sequences of the two genes contained highly conserved residues required for activity and secondary structure formation of Cu/ZnSOD and MnSOD, respectively, and show up to 73.7% and 75.4% identities with their counterparts in other animals. The genomic DNA sequence analysis of Cu/ZnSOD gene revealed this as an intronless gene. Inhibitor studies with purified recombinant Cu/ ZnSOD and MnSOD, both of which were functionally expressed in Escherichia coli, confirmed that they are copper/zinc and manganese-containing SOD, respectively. Immunoblots showed that both C. sinensis Cu/ZnSOD and MnSOD should be antigenic for humans, and both, especially the C. sinensis MnSOD, exhibit extensive cross-reactions with sera of patients infected by other trematodes or cestodes. RT-PCR and SOD activity staining of parasite lysates indicate that there are no significant differences in mRNA level or SOD activity for both species of SOD, indicating cytosolic Cu/ZnSOD and MnSOD might play a comparatively important role in the C. sinensis antioxidant system.

Molecular cloning of TRAF2 binding protein gene and its promoter region from the grass carp Ctenopharyngodon idellus

A tumor necrosis factor receptor-associated factor 2 binding protein (T2BP) gene was isolated from the grass carp (Ctenopharyngodon idellus) by utilizing suppression subtractive hybridization (SSH) and rapid amplification of cDNA ends (RACE). The grass carp T2BP (GT2BP) gene contains an open reading frame of 579 nucleotide(s) (nt), encoding 193 amino acids, with 23 nt 5'-untranslated region and a long 3'-untranslated region of 434 nt including poly (A), 1 AUUUA motif and 4 AUUUUA motifs. No signal peptide has been detected in the predicted GT2BP, but a characteristic forkhead associated domain is present. The GT2BP mRNA shares 83% identity with the zebrafish DNA sequence, and they both have no introns in the genomic DNA. The putative transcription factor binding sites of GT2BP include two C/EBP alpha binding sites, and one c-Jun binding, one AP-1 binding, and one nuclear factor kappa B (NF kappa B) binding sites. Southern blot analysis revealed that the GT2BP was a single-copy gene. Individual difference was observed in GT2BP expression in examined organs of healthy grass carp. However, the expression of GT2BP in all examined organs in a fish with the highest copepod infection level and the significantly higher expression level in spleen and liver in infected fish may indicate its up-regulation with the parasite infection. (c) 2005 Elsevier B.V. All rights reserved.

Molecular and expression characterization of three gonadotropin subunits common alpha, FSH beta and LH beta in groupers

A SMART cDNA plasmid library was constructed from protogyous greasy grouper (Epinephelus coioides) pituitary, and the full-length cDNAs of three gonadotropin (GTH) subunits common alpha, FSH beta and LH beta were cloned and sequenced from the library. The nucleotide sequences of common alpha, FSH beta and LH beta subunit cDNAs are 647, 594 and 574 bp in length, and encode for mature peptides of 94, 99 and 115 aa, respectively. High homology was observed by amino acid sequence alignment and identity comparison of the grouper mature peptides of common alpha, FSH beta and LH beta with that of other fishes. Phylogenetic tree analyses of the three GTH mature subunits revealed similar phylogeny relationships among the studied fish species. Three polyclonal antibodies were prepared from the in vitro expressed common alpha, FSH beta and LH beta mature proteins, respectively. Western blot analysis and immunofluoresence localization were performed on two typical stages of ovarian development stages in red-spotted grouper. Significant differences in protein expression levels of three gonadotropin subunits were revealed between the two ovarian development stages. In the individuals with resting ovary, common alpha was almost not detected in pituitaries, and FSH beta and LH beta expression levels were very low. While in the individuals with developing ovary, the expression of all three gonadotropin subunits reached to a high level. Immunofluoresence localization indicated that the grouper FSH beta cells mainly distributed in the middle area of PPD, while the LH beta cells distributed more widely, including in the area similar to the FSH beta cells and at the external periphery of pituitary near to the PI side. The common alpha might be expressed in both FSH beta and LH beta cells. Double immunofluoresence localization further demonstrated FSH beta and LH beta expression in distinct cells in the PPD area, although the FSH beta and LH beta cells were detected in the identical area of PPD. (c) 2005 Elsevier Ireland Ltd. All rights reserved.

Mitochondrial 16S rRNA sequence variations and phylogeny of the Chinese sisorid catfishes

Partial sequences of mitochondrial 16S rRNA gene were obtained by PCR amplification for comparisons among nine species of glyptosternoid fishes and six species of non-glyptosternoids representing 10 sisorid genera. There are compositional biases in the A-rich impaired regions and G-rich paired regions. A-G transitions are primarily responsible for the Ts/Tv bias in impaired regions. The overall substitution rate in impaired regions is almost two times higher than that in the paired regions. Saturation plots at comparable levels of sequence divergence demonstrate no saturation effects. Phylogenetic analyses using both maximum likelihood and Bayesian methods support the monophyly of Sisoridae. Chinese sisorid catfishes are composed of two major lineages, one represented by (Gagata (Bagarius, Glyptothorax)) and the other by "glyptosternoids + Pseudecheneis". The glyptosternoids may not be a monophyletic group. A previous hypothesis referring to Pseudecheneis as the sister group of monophyletic glyptosternoids, based on morphological evidence, is not supported by the molecular data. Pseudecheneis is shown to be a sister taxon of Glaridoglanis. Pareuchiloglanis might be paraphyletic with Pseudexostoma and Euchiloglanis. Our results also support the hypothesis that Pareuchiloglanis anteanalis might be considered as the synonyms of Pareuchiloglanis sinensis, and genus Euchiloglanis might have only one valid species, Euchiloglanis davidi.

Molecular cloning of the viperin gene and its promoter region from the mandarin fish Siniperca chuatsi

A viperin gene has been cloned from the mandarin fish (Siniperca chuatsi). From the first transcription initiation site, the mandarin fish viperin gene extends 3163 nucleotides to the end of the 3' untranslated region, and it contains six exons and five introns. The open reading frame of the viperin transcript has 1062 nucleotides which encode a 354 amino acid peptide. The amino acid sequence of mandarin fish viperin shows high identities with its homologues in teleosts and mammals except for the first 70 amino acids. A characteristic feature in the viperin promoter region was the presence of five putative ICSBP (IRF8) binding sites and one IRFI binding site. The viperin gene expressed mainly in lymphoid tissues before stimulation, but its expression can be examined in almost all the organs investigated after stimulation with virus or Poly I:C. The expression pattern and promoter sequence may be considered as the indirect evidence that the transcription of viperin is regulated by interferons or interferon induced genes. (C) 2004 Elsevier B.V. All rights reserved.

Molecular phylogeny of the specialized schizothoracine fishes (Teleostei : Cyprinidae), with their implications for the uplift of the Qinghai-Tibetan Plateau

Molecular phylogeny of three genera containing nine species and subspecies of the specialized schizothoracine fishes are investigated based on the complete nucleotide sequence of mitochondrial cytochrome b gene. Meantime relationships between the main cladogenetic events of the specialized schizothoracine fishes and the stepwise uplift of the Qinghai-Tibetan Plateau are also conducted using the molecular clock, which is calibrated by geological isolated events between the upper reaches of the Yellow River and the Qinghai Lake. Results indicated that the specialized schizothoracine fishes are not a monophyly. Five species and subspecies of Ptychobarbus form a monophyly. But three species of Gymnodiptychus do not form a monophyly. Gd. integrigymnatus is a sister taxon of the highly specialized schizothoracine fishes while Gd. pachycheilus has a close relation with Gd. dybowskii, and both of them are as a sister group of Diptychus maculatus. The specialized schizothoracines fishes might have originated during the Miocene (about 10 MaBP), and then the divergence of three genera happened during late Miocene (about 8 MaBP). Their main specialization occurred during the late Pliocene and Pleistocene (3.54-0.42 MaBP). The main cladogenetic events of the specialized schizothoracine fishes are mostly correlated with the geological tectonic events and intensive climate shift happened at 8, 3.6, 2.5 and 1.7 MaBP of the late Cenozoic. Molecular clock data do not support the hypothesis that the Qinghai-Tibetan Plateau uplifted to near present or even higher elevations during the Oligocene or Miocene, and neither in agreement with the view that the plateau uplifting reached only to an altitude of 2000 in during the late Pliocene (about 2.6 MaBP).