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.

Sox genes evolution in closely related young tetraploid cyprinid fishes and their diploid relative

Previous study and analysis of cytochrome b suggested that polyploidization event in the genus Tor occurred about 10 Mya ago. In order to understand evolutionary fates of Sox gene in the early stage of genome duplication at the nucleotide level, PCR surveys for Sox genes in three closely related cyprinid fishes T douronensis (2n = 100), T qiaojiensis (2n = ?), T sinensis (2n = 100) and their relative T brevifilis (2n = 50) were performed. Totally, 52 distinct Sox genes were obtained in these four species, representing SoxB, SoxC, and SoxE group. As expected, isoforms of some Sox genes correspond with the ploidy of species, such as two copies of Sox9a exist in tetraploid species. Analysis indicated that duplicated Sox gene pairs caused by polyploidization evolved independently of each other within polyploid species. Results of substitution rate showed nearly equal rate of nonsynonymous substitution of duplicated Sox orthologs among different polyploid species and their diploid relative orthologs, suggesting at the early stage of genome duplicated Sox orthologs are under similar selective constraints in different polyploidy species and their diploid relative at the amino acid level. All PCR fragments of Sox genes obtained in this study are not accompanied by obvious increase in mutations and pseudogene formation which means that they are under strong purifying selection, suggesting that they are functional at the DNA level. Cenealogical analysis revealed that T qiaojiensis was tetraploid, and T douronensis, T qiaojiensis as well as T sinensis had an allotetraploid ancestor. (C) 2009 Elsevier B.V. All rights reserved.

Origin and evolution of the RIG-I like RNA helicase gene family

Background: The DExD/H domain containing RNA helicases such as retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5) are key cytosolic pattern recognition receptors (PRRs) for detecting nucleotide pathogen associated molecular patterns (PAMPs) of invading viruses. The RIG-I and MDA5 proteins differentially recognise conserved PAMPs in double stranded or single stranded viral RNA molecules, leading to activation of the interferon system in vertebrates. They share three core protein domains including a RNA helicase domain near the C terminus (HELICc), one or more caspase activation and recruitment domains (CARDs) and an ATP dependent DExD/H domain. The RIG-I/MDA5 directed interferon response is negatively regulated by laboratory of genetics and physiology 2 (LGP2) and is believed to be controlled by the mitochondria antiviral signalling protein (MAVS), a CARD containing protein associated with mitochondria. Results: The DExD/H containing RNA helicases including RIG-I, MDA5 and LGP2 were analysed in silico in a wide spectrum of invertebrate and vertebrate genomes. The gene synteny of MDA5 and LGP2 is well conserved among vertebrates whilst conservation of the gene synteny of RIG-I is less apparent. Invertebrate homologues had a closer phylogenetic relationship with the vertebrate RIG-Is than the MDA5/LGP2 molecules, suggesting the RIG-I homologues may have emerged earlier in evolution, possibly prior to the appearance of vertebrates. Our data suggest that the RIG-I like helicases possibly originated from three distinct genes coding for the core domains including the HELICc, CARD and ATP dependent DExD/H domains through gene fusion and gene/domain duplication. Furthermore, presence of domains similar to a prokaryotic DNA restriction enzyme III domain (Res III), and a zinc finger domain of transcription factor (TF) IIS have been detected by bioinformatic analysis. Conclusion: The RIG-I/MDA5 viral surveillance system is conserved in vertebrates. The RIG-I like helicase family appears to have evolved from a common ancestor that originated from genes encoding different core functional domains. Diversification of core functional domains might be fundamental to their functional divergence in terms of recognition of different viral PAMPs.

Population genetic structure and geographical differentiation of the Chinese catfish Hemibagrus macropterus (Siluriformes, Bagridae): Evidence for altered drainage patterns

National Natural Science Foundation of China (NSFC) [2007CB411600, 30530120]

Identification of an additional two-cysteine containing type I interferon in rainbow trout Oncorhynchus mykiss provides evidence of a major gene duplication event within this gene family in teleosts

Multiple type I interferons (IFNs) have recently been identified in salmonids, containing two or four conserved cysteines. In this work, a novel two-cysteine containing (2C) IFN gene was identified in rainbow trout. This novel trout IFN gene (termed IFN5) formed a phylogenetic group that is distinct from the other three salmonid IFN groups sequenced to date and had a close evolutionary relationship with IFNs from advanced fish species. Our data demonstrate that two subgroups are apparent within each of the 2C and 4C type I IFNs, an evolutionary outcome possibly due to two rounds of genome duplication events that have occurred within teleosts. We have examined gene expression of the trout 2C type I IFN in cultured cells following stimulation with lipopolysaccharide, phytohaemagglutinin, polyI:C or recombinant IFN, or after transfection with polyI:C. The kinetics of gene expression was also studied after viral infection. Analysis of the regulatory elements in the IFN promoter region predicted several binding sites for key transcription factors that potentially play an important role in mediating IFN5 gene expression.

Bead-probe complex capture a couple of SINE and LINE family from genomes of two closely related species of East Asian cyprinid directly using magnetic separation

Background: Short and long interspersed elements (SINEs and LINEs, respectively), two types of retroposons, are active in shaping the architecture of genomes and powerful tools for studies of phylogeny and population biology. Here we developed special protocol to apply biotin-streptavidin bead system into isolation of interspersed repeated sequences rapidly and efficiently, in which SINEs and LINEs were captured directly from digested genomic DNA by hybridization to bead-probe complex in solution instead of traditional strategy including genomic library construction and screening. Results: A new couple of SINEs and LINEs that shared an almost identical 3'tail was isolated and characterized in silver carp and bighead carp of two closely related species. These SINEs (34 members), designated HAmo SINE family, were little divergent in sequence and flanked by obvious TSD indicated that HAmo SINE was very young family. The copy numbers of this family was estimated to 2 x 10(5) and 1.7 x 10(5) per haploid genome by Real-Time qPCR, respectively. The LINEs, identified as the homologs of LINE2 in other fishes, had a conserved primary sequence and secondary structures of the 3'tail region that was almost identical to that of HAmo SINE. These evidences suggest that HAmo SINEs are active and amplified recently utilizing the enzymatic machinery for retroposition of HAmoL2 through the recognition of higher-order structures of the conserved 42-tail region. We analyzed the possible structures of HAmo SINE that lead to successful amplification in genome and then deduced that HAmo SINE, SmaI SINE and FokI SINE that were similar in sequence each other, were probably generated independently and created by LINE family within the same lineage of a LINE phylogeny in the genomes of different hosts. Conclusion: The presented results show the advantage of the novel method for retroposons isolation and a pair of young SINE family and its partner LINE family in two carp fishes, which strengthened the hypotheses containing the slippage model for initiation of reverse transcription, retropositional parasitism of SINEs on LINEs, the formation of the stem loop structure in 3'tail region of some SINEs and LINEs and the mechanism of template switching in generating new SINE family.

Rana grylio virus thymidine kinase gene: an early gene of iridovirus encoding for a cytoplasmic protein

The presence of thymidine kinase (TK) is a feature of many large DNA viruses. Here, a TK gene homologue was cloned and characterized from Rana grylio virus (RGV), a member of family Iridoviridae. RGV TK encodes a protein of 195 aa with a predicted molecular mass of 22.1 kDa. Homologues of the protein were present in all the currently sequenced iridoviruses, and phylogenetic analysis showed that it was much close to cellular TK type 2 (TK2), deoxycytidine kinase (dCK) and deoxyguanosine kinase (dGK). Subsequently, Western blotting revealed TK expression increased with time from 6 h post-infection in RGV-infected cells. Using drug inhibition analysis by protein synthesis inhibitor (cycloheximide) and DNA replication inhibitor (cytosine arabinofuranoside), RGV TK was classified as the early expression gene during in vitro infection. Subcellular localization by TK-GFP fusion protein expression and immunofluorescence staining showed RGV TK was an exclusively cytoplasmic protein in fish cells. Collectively, current data indicate that RGV TK was an early gene of iridovirus which encoded a cytoplasmic protein in fish cells.

Comparative study and expression analysis of the interferon gamma gene locus cytokines in Xenopus tropicalis

Using bioinformatics approach, the genome locus containing interleukin (IL)-22, IL-26, and interferon gamma (IFN-gamma) genes has been identified in the amphibian, Xenopus tropicalis. Like that in other vertebrates such as fish, birds, and mammals, the Xenopus IL-22, IL-26, and IFN-gamma are clustered in the same chromosome and the adjacent genes are conserved. The genomic structures of the Xenopus IL-22, IL-26, and IFN-gamma gene were identical to that of their mammalian counterparts. The Xenopus IL-22 and IL-26 genes contained five exons and four introns while the Xenopus IFN-gamma gene consisted of four exons and three introns. The Xenopus IL-22, IL-26, and IFN-gamma share 14.1-41.6%, 14.6-31.2%, and 23.7-36.5% identity to their counterparts in other species, respectively. Reverse-transcription polymerase chain reaction (PCR) and real-time quantitative PCR analyses revealed that the expression of IL-22, IL-26, and IFN-gamma genes was significantly upregulated after simulation with bacterial polyliposaccharide and/or synthetic double-stranded poly(I:C), suggesting these cytokines like those in other vertebrates play an important role in regulating immune response in Xenopus.

Phylogeography of the freshwater catfish Hemibagrus guttatus (Siluriformes, Bagridae): Implications for South China biogeography and influence of sea-level changes

National Natural Science Foundation of China (NSFC) ; [2007CB411600]; [30530120]

Cloning, characterization, and gene expression analysis of a novel cadmium metallothionein gene in Tetrahymena pigmentosa

A novel cadmium-inducible metallothionein (MT) gene (Tpig-MT1) was cloned and sequenced from the ciliate Tetrahymena pigmentosa. The number of deduced amino acids is 118. The polypeptide possesses CCC and CC clusters characteristic of typical Tetrahymena Cd-inducible MTs. The structure of Tpig-MT1 is different from the reported Cd-MT in T. pyriformis, T. thermophila and T. pigmentosa. Tpig-MT1 contains two intragenic tandem repeats with 72.9% identity described as Tpig-MT1 (repeat A1) and Tpig-MT1 (repeat A2). The transcriptional response of Tpig-MT1 gene to different heavy metals (Cd, Cu, Zn, Hg, Pb) and oxidative stress (H2O2) was measured using real-time quantitative PCR. The results showed that the gene was quickly induced (1 h) by the five heavy metals and the order of expression level was Hg>Pb>Cd>Cu>Zn. The induction effect of H2O2 was 5-fold after about 15 min, but soon decreased to a non-significant level (30 min). The genetic diversity of Tetrahymena MT genes is discussed in relation to the unique structure of the Tpig-MT1 gene and other reported Cd-MT isoforms. (C) 2008 Elsevier B.V. All rights reserved.

Molecular characterization and expression pattern of fetuin-B in gibel carp (Carassius auratus gibelio)

Fetuin-B has recently been cloned and identified from rats, mice, and humans; their expression patterns, however, have not been elucidated yet. In this study, Cagfetuin-B has been cloned in gibel carp. RT-PCR and Western blot detection revealed that Cagfetuin-B is first transcribed from the blastula stage and at a relatively stable level afterward during embryogenesis and the larval stage. Cagfetuin-B transcripts are predominantly distributed over the yolk syncytial layer in the early embryos and later restricted to the cells of liver and brain in newly hatched larvae. Moreover, a dynamic distribution of Cagfetuin-B protein was observed in brain, kidney, liver, and skin during morphogenesis. In adult fish, Cagfetuin-B transcripts are restricted in liver and ovary. Our work, for the first time, revealed the extra-hepatic transcription and a dynamic distribution of fetuin-B during embryogenesis and in adults, which indicates the potential roles of fetuin-B in fish organogenesis.

Population differentiation of Procypris rabaudi (Tchang) in the upper reaches of the Yangtze River revealed by mitochondrial DNA and RAPD markers

Procypris rabaudi (Tchang) is a cyprinid fish endermic to middle and upper reaches of the Yangtze River. Besides in main stream and large tributaries, there exists an early matured, small-sized ecological type in a small tributary, Tang River. In this study, mitochondrial DNA cytochrome h (cyt b) gene sequence analysis and randomly amplified polymorphic DNA (RAPD) analysis were performed to investigate the differentiation of the Tang River population from the Mudong reach population of the Yangtze River, with the purpose of conservation and exploitation of this fish. In the 1140 bps of cyt b gene sequence surveyed, 20 sites were found polymorphic, which defined 23 haplotypes. Among them, four haplotypes accounted for 54.4% of all individuals, while population-specific haplotypes occurred in low frequencies. Analysis of molecular variation on cyt b data revealed no significant partition existing between Tang River population and Mudong reach population. Analyses of 132 RAPD loci suggested that genetic variation between populations was significant, though values of different F-ST were not very high. The results revealed low genetic diversity and the beginning of population differentiation, suggesting that Tang River population should be designated as a separate Management Unit.

Development and characterization of polymorphic microsatellite loci in the endangered Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientalis)

The finless porpoise (Neophocaena phocaenoides) is one of the smallest cetacean species widely distributed in the shallow coastal waters of the Indo-Pacific Oceans. The population size of the Yangtze subspecies (N. p. asiaeorientalis) has sharply decreased in the last two decades and access to objective data on its population structure and genetic diversity would be of great assistance for their proper management. Here we report on the isolation of nine polymorphic microsatellite using the "Fast Isolation by AFLP of Sequences Containing repeats" (FIASCO) protocol. Polymorphism was assessed using 30 porpoise individuals randomly sampled in the Yangtze River. The number of alleles per locus varies from 2 to 9, with an average value of 5.56, whereas the ranges of observed and expected heterozygosi ties were 0.300-0.633 (mean 0.496) and 0.473-0.804 (mean 0.659), respectively.

Variation patterns of the mitochondrial 16S rRNA gene with secondary structure constraints and their application to phylogeny of cyprinine fishes (Teleostei : Cypriniformes)

The mitochondrial 16S ribosomal RNA (rRNA) gene sequences from 93 cyprinid fishes were examined to reconstruct the phylogenetic relationships within the diverse and economically important subfamily Cyprininae. Within the subfamily a biased nucleotide composition (A > T, C > G) was observed in the loop regions of the gene, and in stem regions apparent selective pressures of base pairing showed a bias in favor of G over C and T over A. The bias may be associated with transition-transversion bias. Rates of nucleotide substitution were lower in stems than in loops. Analysis of compensatory substitutions across these taxa demonstrates 68% covariation in the gene and a logical weighting factor to account for dependence in mutations for phylogenetic inference should be 0.66. Comparisons of varied stem-loop weighting schemes indicate that the down-weightings for stem regions could improve the phylogenetic analysis and the degree of non-independence of stem substitutions was not as important as expected. Bayesian inference under four models of nucleotide substitution indicated that likelihood-based phylogenetic analyses were more effective in improving the phylogenetic performance than was weighted parsimony analysis. In Bayesian analyses, the resolution of phylogenies under the 16-state models for paired regions, incorporating GTR + G + I models for unpaired regions was better than those under other models. The subfamily Cyprininae was resolved as a monophyletic group, as well as tribe Labein and several genera. However, the monophyly of the currently recognized tribes, such as Schizothoracin, Barbin, Cyprinion + Onychostoma lineages, and some genera was rejected. Furthermore, comparisons of the parsimony and Bayesian analyses and results of variable length bootstrap analysis indicates that the mitochondrial 16S rRNA gene should contain important character variation to recover well-supported phylogeny of cyprinid taxa whose divergences occurred within the recent 8 MY, but could not provide resolution power for deep phylogenies spanning 10-19 MYA. (c) 2008 Published by Elsevier Inc.