Gene duplication in the genome of parasitic Giardia lamblia

Background: Giardia are a group of widespread intestinal protozoan parasites in a number of vertebrates. Much evidence from G. lamblia indicated they might be the most primitive extant eukaryotes. When and how such a group of the earliest branching unicellular eukaryotes developed the ability to successfully parasitize the latest branching higher eukaryotes (vertebrates) is an intriguing question. Gene duplication has long been thought to be the most common mechanism in the production of primary resources for the origin of evolutionary novelties. In order to parse the evolutionary trajectory of Giardia parasitic lifestyle, here we carried out a genome-wide analysis about gene duplication patterns in G. lamblia. Results: Although genomic comparison showed that in G. lamblia the contents of many fundamental biologic pathways are simplified and the whole genome is very compact, in our study 40% of its genes were identified as duplicated genes. Evolutionary distance analyses of these duplicated genes indicated two rounds of large scale duplication events had occurred in G. lamblia genome. Functional annotation of them further showed that the majority of recent duplicated genes are VSPs (Variant-specific Surface Proteins), which are essential for the successful parasitic life of Giardia in hosts. Based on evolutionary comparison with their hosts, it was found that the rapid expansion of VSPs in G. lamblia is consistent with the evolutionary radiation of placental mammals. Conclusions: Based on the genome-wide analysis of duplicated genes in G. lamblia, we found that gene duplication was essential for the origin and evolution of Giardia parasitic lifestyle. The recent expansion of VSPs uniquely occurring in G. lamblia is consistent with the increment of its hosts. Therefore we proposed a hypothesis that the increment of Giradia hosts might be the driving force for the rapid expansion of VSPs.

BarraCUDA – a Fast Sequence Mapping Software using Graphics Processing Units

High-throughput DNA sequencing (HTS) instruments today are capable of generating millions of sequencing reads in a short period of time, and this represents a serious challenge to current bioinformatics pipeline in processing such an enormous amount of data in a fast and economical fashion. Modern graphics cards are powerful processing units that consist of hundreds of scalar processors in parallel in order to handle the rendering of high-definition graphics in real-time. It is this computational capability that we propose to harness in order to accelerate some of the time-consuming steps in analyzing data generated by the HTS instruments. We have developed BarraCUDA, a novel sequence mapping software that utilizes the parallelism of NVIDIA CUDA graphics cards to map sequencing reads to a particular location on a reference genome. While delivering a similar mapping fidelity as other mainstream programs , BarraCUDA is a magnitude faster in mapping throughput compared to its CPU counterparts. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the mapping throughput. BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the mapping of millions of sequencing reads generated by HTS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology. BarraCUDA is currently available at http://seqbarracuda.sf.net

Discovery of a male-biased mutant family and identification of a male-specific SCAR marker in gynogenetic gibel carp Carassius auratus gibelio

Gibel carp ( Carassius auratus gibelio) is a uniquely gynogenetic species with a minor ratio of males in natural habitats, but its male origin and sex determination mechanisms have been unknown. In this study, a male-biased mutant family was discovered from the gynogenetic gibel carp, and a male-specific SCAR marker was identified from the mutant family. Normal spermatogenesis was observed in the male testes by immuno. fluorescence histochemistry. Nearly identical AFLP profiles were observed between males and females, but a male-specific 86 bp AFLP fragment was screened by sex-pool bulked segregant analysis and individual screening. Based on the male-specific AFLP fragment, a total of 579 bp sequences were cloned by genome walking. Subsequently, a male-specific SCAR marker was designed, and the male-specific DNA fragment was confirmed to be steadily transmitted to the next generation and consistently detected only in males. (C) 2009 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Genome-wide identification and characterization of cytochrome P450 monooxygenase genes in the ciliate Tetrahymena thermophila

Background: Cytochrome P450 monooxygenases play key roles in the metabolism of a wide variety of substrates and they are closely associated with endocellular physiological processes or detoxification metabolism under environmental exposure. To date, however, none has been systematically characterized in the phylum Ciliophora. T. thermophila possess many advantages as a eukaryotic model organism and it exhibits rapid and sensitive responses to xenobiotics, making it an ideal model system to study the evolutionary and functional diversity of the P450 monooxygenase gene family. Results: A total of 44 putative functional cytochrome P450 genes were identified and could be classified into 13 families and 21 sub-families according to standard nomenclature. The characteristics of both the conserved intron-exon organization and scaffold localization of tandem repeats within each P450 family clade suggested that the enlargement of T. thermophila P450 families probably resulted from recent separate small duplication events. Gene expression patterns of all T. thermophila P450s during three important cell physiological stages (vegetative growth, starvation and conjugation) were analyzed based on EST and microarray data, and three main categories of expression patterns were postulated. Evolutionary analysis including codon usage preference, sit-especific selection and gene-expression evolution patterns were investigated and the results indicated remarkable divergences among the T. thermophila P450 genes. Conclusion: The characterization, expression and evolutionary analysis of T. thermophila P450 monooxygenase genes in the current study provides useful information for understanding the characteristics and diversities of the P450 genes in the Ciliophora, and provides the baseline for functional analyses of individual P450 isoforms in this model ciliate species.

Antimicrobial activity-specific to Gram-negative bacteria and immune modulation-mediated NF-kappa B and Sp1 of a medaka beta-defensin

Defensins are a group of cationic antimicrobial peptides which play an important role in the innate immune system by exerting their antimicrobial activity against pathogens. In this study, we cloned a novel beta-defensin cDNA from medaka (Oryzias latipes) by rapid amplification of cDNA ends (RACE) technique. The full-length cDNA consists of 480 bp, and the open reading frame (CRF) of 189 bp encodes a polypeptide of 63 amino acids (aa) with a predicted molecular weight of 7.44 kDa. Its genomic organization was analyzed, and Southern blot detection confirmed that only one copy of beta-defensin exists in the medaka HNI strain. RT-PCR, Western blot and immunohistochemistry detections showed that the beta-defensin transcript and protein could be detected in eyes, liver, kidney, blood, spleen and gill, and obviously prevalent expression was found in eyes. Antimicrobial activity of the medaka beta-defensin was evaluated, and the antibacterial activity-specific to Gram-negative bacteria was revealed. Furthermore, the lipopolysaccharide (LPS), a major component of the outer membrane of Gram-negative bacteria, was demonstrated to be able to induce about 13-fol up-regulation of the beta-defensin within first 12 h. In addition, promoter and promoter mutagenesis analysis were performed in the medaka beta-defensin. A proximal 100 base pair(bp) sequence (+26 to -73)and the next 1700 bp sequence (-73 to -1755) were demonstrated to be responsible for the basal promoter activity and for the transcription regulation. Three nuclear factor kappa B (NF-kappa B) cis-elements and a Sp1 cis-element were revealed by mutagenesis analysis to exist in the 5' flanking sequence, and they were confirmed to be responsible for the up-regulation of medaka beta-defensin stimulated by LPS. And, the Sp1 cis-element was further revealed to be related to the basal promoter activity, and transcriptional factor II D (TFIID) was found to be in charge of the gene transcription initiation. All the obtained data suggested that the novel medaka beta-defensin should have antimicrobial activity-specific to Gram-negative bacteria, and the antibacterial immune function should be modulated by NF-kappa B and Sp1. (C) 2008 Elsevier Ltd. All rights reserved.

Genomic sequence of mandarin fish rhabdovirus with an unusual small non-transcriptional ORF

The complete genome of mandarin fish Siniperca chuatsi rhabdovirus (SCRV) was cloned and sequenced. It comprises 11,545 nucleotides and contains five genes encoding the nucleoprotein N, the phosphoprotein P, the matrix protein M, the glycoprotein G, and the RNA-dependent RNA polymerase protein L. At the 3' and 5' termini of SCRV genome, leader and trailer sequences show inverse complementarity. The N, P, M and G proteins share the highest sequence identities (ranging from 14.8 to 41.5%) with the respective proteins of rhabdovirus 903/87, the L protein has the highest identity with those of vesiculoviruses, especially with Chandipura virus (44.7%). Phylogenetic analysis of L proteins showed that SCRV clustered with spring vireamia of carp virus (SVCV) and was most closely related to viruses in the genus Vesiculovirus. In addition, an overlapping open reading frame (ORF) predicted to encode a protein similar to vesicular stomatitis virus C protein is present within the P gene of SCRV. Furthermore, an unoverlapping small ORF downstream of M ORF within M gene is predicted (tentatively called orf4). Therefore, the genomic organization of SCRV can be proposed as 3' leader-N-P/C-M-(orf4)-G-L-trailer 5'. Orf4 transcription or translation products could not be detected by northern or Western blot, respectively, though one similar mRNA band to M mRNA was found. This is the first report on one small unoverlapping ORF in M gene of a fish rhabdovirus. (c) 2007 Elsevier B.V. All rights reserved.

Site-directed gene integration in transgenic zebrafish mediated by cre recombinase using a combination of mutant Lox sites

With current gene-transfer techniques in fish, insertion of DNA into the genome occurs randomly and in many instances at multiple sites. Associated position effects, copy number differences, and multiple gene interactions make gene expression experiments difficult to interpret and fish phenotype less predictable. To meet different fish engineering needs, we describe here a gene targeting model in zebrafish. At first, four target zebrafish lines, each harboring a single genomic lox71 target site, were generated by zebrafish transgenesis. The zygotes of transgenic zebrafish lines were coinjected with capped Cre mRNA and a knockin vector pZklox66RFP. Site-specific integration event happened from one target zebrafish line. In this line two integrant zebrafish were obtained from more than 80,000 targeted embryos (integrating efficiency about 10(-4) to 10(-5)) and confirmed to have a sole copy of the integrating DNA at the target genome site. Genomic polymerase chain reaction analysis and DNA sequencing verified the correct gene target events where lox71 and lox66 have accurately recombined into double mutant lox72 and wild-type loxP. Each integrant zebrafish chosen for analysis harbored the transgene rfp at the designated egfp concatenates. Although the Cre-mediated recombination is site specific, it is dependent on a randomly placed target site. That is, a genomic target cannot be preselected for integration based solely on its sequence. Conclusively, an rfp reporter gene was successfully inserted into the egfp target locus of zebrafish genome by Cre-lox-mediated recombination. This site-directed knockin system using the lox71/lox66 combination should be a promising gene-targeting platform serving various purposes in fish genetic engineering.

The complete mitochondrial genome of the Chinese hook snout carp Opsariichthys bidens (Actinopterygii : Cyprinifonnes) and an altemative pattem of mitogenomic evolution in vertebrate

The complete mitochondrial genome sequence of the Chinese hook snout carp, Opsariichthys bidens, was newly determined using the long and accurate polymerase chain reaction method. The 16,611-nucleotide mitogenome contains 13 protein-coding genes, two rRNA genes (12S, 16S) 22 tRNA genes, and a noncoding control region. We use these data and homologous sequence data from multiple other ostariophysan fishes in a phylogenetic evaluation to test hypothesis pertaining to codon usage pattern of O. bidens mitochondrial protein genes as well as to re-examine the ostariophysan phylogeny. The mitochondrial genome of O. bidens reveals an alternative pattern of vertebrate mitochondrial evolution. For the mitochondrial protein genes of O. bidens, the most frequently used codon generally ends with either A or C, with C preferred over A for most fourfold degenerate codon families; the relative synonymous codon usage of G-ending codons is greatly elevated in all categories. The codon usage pattern of O. bidens mitochondrial protein genes is remarkably different from the general pattern found previously in the relatively closely 9 related zebrafish and most other vertebrate mitochondria. Nucleotide bias at third codon positions is the main cause of codon bias in the mitochondrial protein genes of O. bidens, as it is biased particularly in favor of C over A. Bayesian analysis of 12 concatenated mitochondrial protein sequences for O. bidens and 46 other teleostean taxa supports the monophyly of Cypriniformes and Otophysi and results in a robust estimate of the otophysan phylogeny. (C) 2007 Published by Elsevier B.V.

Identification of genome organization in the unusual allotetraploid form of Carassius auratus gibelio

The unusual allotetraploid form with unequal contribution of chromosome sets was discovered from the gynogenetic offspring of Carassius auratus gibelio stimulated by red common carp sperm. In this study, genomic in situ hybridization (GISH) and fluorescent in situ hybridization (FISH) with 45S rDNA probe are used. The GISH results lead to the identification of species-specific chromosomes, which permits to demonstrate the origin and genome organization in the allotetraploid form. Moreover, chromosome localization of 45S rDNA and co-localizations of 45S rDNA and Cyprinus carpio genomic DNA further confirm that one extra 45S rDNA positive chromosome in the allotetraploid form originates from the paternal haploid genome of C carpio, and other 5 45S rDNA-containing chromosomes are from the maternal genome of Carassius auratus gibelio. And, the correlation between 45 rDNA and the nucleolar organizer regions (NORs) is confirmed by silver nitrate staining. The data provide direct experiment evidence that the allotetraploid actually contains three chromosome sets of Carassius auratus gibelio and one chromosome set of C carpio, and will be a useful genetic material for both basic research and breeding practice. (c) 2006 Elsevier B.V. All rights reserved.

Identification of a putative oocyte-specific small nuclear ribonucleoprotein polypeptide C in gibel carp

Previous studies have demonstrated that germinal vesicle of amphibian oocyte contains small nuclear ribonucleoprotein polypeptide C (SNRPC). In this study, a putative member of SNRPC was identified from Carassius auratus gibelio oocyte cDNA library. Its full-length cDNA has an open reading frame of 201 nt for encoding a peptide of 66 an, a short 5'-UTR of 19 nt and a long 3'-UTR of 347 nt including a polyadenylation signal and poly- (A) tail, and the deduced amino acid sequence has 47% identity with the C-terminal of the zebrafish small nuclear ribonucleoprotein polypeptide C. Western blot analysis revealed its oocyte-specific expression. Immunofluorescence localization indicated that its gene product localized to numerous nucleoli within the oocytes and showed dynamic changes with the nucleoli during oocyte maturation. RT-PCR and Western blot analysis further revealed its constant presence in the oocytes and in the embryos until hatching. The data suggested that the newly identified CagOSNRPC might be a nucleolar protein. (c) 2006 Elsevier Inc. All rights reserved.

Sox genes in grass carp (Ctenopharyngodon idella) with their implications for genome duplication and evolution

The Sox gene family is found in a broad range of animal taxa and encodes important gene regulatory proteins involved in a variety of developmental processes. We have obtained clones representing the HMG boxes of twelve Sox genes from grass carp (Ctenopharyngodon idella), one of the four major domestic carps in China. The cloned Sox genes belong to group B1, B2 and C. Our analyses show that whereas the human genome contains a single copy of Sox4, Sox11 and Sox14, each of these genes has two co-orthologs in grass carp, and the duplication of Sox4 and Sox11 occurred before the divergence of grass carp and zebrafish, which support the "fish-specific whole-genome duplication" theory. An estimation for the origin of grass carp based on the molecular clock using Sox1, Sox3 and Sox11 genes as markers indicates that grass carp (subfamily Leuciscinae) and zebrafish (subfamily Danioninae) diverged approximately 60 million years ago. The potential uses of Sox genes as markers in revealing the evolutionary history of grass carp are discussed.

Upstream regulatory region of zebrafish lunatic fringe: Isolation and promoter analysis

Lunatic fringe (Lfng), one modulator of Notch signaling, plays an essential part in demarcation of tissues boundaries during animal early development, especially somitogenesis. To characterize the promoter of zebrafish 1fng and generate somite-specific transgenic zebrafish, we isolated the upstream regulatory region of zebrafish 1fng by blast search at the Ensembl genome database (http://www. ensembl.org) and analyzed the promoter activity using green fluorescent protein (GFP) as a reporter. Promoter activity assay in zebrafish shows that the 0.2-kb fragment containing GC-box, CAAT-box, and TATA-box can direct tissue-specific GFP expression, while the 0.4-kb and 1.2-kb fragments with further upstream sequence included drive GFP expression more efficiently. We produced 1fngEGFP-transgenic founders showing somite-specific expression of GFP and consequently generated a hemizygous 1fngEGFP-transgenic line. The eggs from 1fngEGFP-transgenic female zebrafish show strong GFP expression, which is consistent to the reverse-transcription polymerase chain reaction PCR (RT-PCR) detection of 1fng transcripts in the fertilized eggs. This reveals that zebrafish 1fng is a maternal factor existing in matured eggs, suggesting that fish somitogcnesis may be influenced by maternal factors.

The complete mitochondrial genome of the helmet catfish Cranoglanis bouderius (Silurifonnes : Cranoglanididae) and the phylogeny of otophysan fishes

The complete sequence of the 16,539 nucleotide mitochondrial genome from the single species of the catfish family Cranoglanididae, the helmet catfish Cranoglanis bouderius, was determined using the long and accurate polymerase chain reaction (LA PCR) method. The nucleotide sequences of C. bouderius mitochondrial DNA have been compared with those of three other catfish species in the same order. The contents of the C. bouderius mitochondrial genome are 13 protein-coding genes, two ribosomal RNA and 22 transfer RNA genes, and a non-coding control region, the gene order of which is identical to that observed in most other vertebrates. Phylogenetic analyses for 13 otophysan fishes were performed using Bayesian method based on the concatenated mtDNA protein-coding gene sequence and the individual protein-coding gene sequence data set. The competing otophysan topologies were then tested by using the approximately unbiased test, the Kishino-Hasegawa test, and the Shimodaira-Hasegawa test. The results show that the grouping ((((Characifonnes, Gymnotiformes), Siluriformes), Cyprinifionnes), outgroup) is the most likely but there is no significant difference between this one and the other alternative hypotheses. In addition, the phylogenetic placement of the family Cranoglanididae among siluriform families was also discussed. (c) 2006 Elsevier B.V. All rights reserved.

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.