Prediction and analysis of paralogous proteins in Trichomonas vaginalis genome

Trichomonas vaginalis causes trichomoniasis, second most sexually transmitted disease. The genome sequence draft of T. vaginalis was published by The Institute of Genomic Research reveals an abnormally large genome size of 160 Mb. It was speculated that a significant portion of the proteome contains paralogous proteins. The present study was aimed at identification and analysis of the paralogous proteins. The all against all search approach is used to identify the paralogous proteins. The dataset of proteins was retrieved from TIGR and TrichDB FTP server. The BLAST-P program performed all against all database searches against the protein database of Trichomonas vaginalis available at NCBI genome database. In the present study about 50,000 proteins were searched where 2,700 proteins were found to be paralogous under the rigid selection criteria. The Pfam database search has identified significant number of paralogous proteins which were further categorized among different 1496 paralogous protein in pfam families, 1027 paralogous protein contains domain, 60 proteins were having different repeats and 1092 paralogous protein sequences of clans. Such identification and functional annotation of paralogous proteins will also help in removing paralogous proteins from possible drug targets in future. Presence of huge number of paralogous proteins across wide range of gene families and domains may be one of the possible mechanisms involved in the T. vaginalis genome expansion and evolution.     

KEGG: Kyoto Encyclopedia of Genes and Genomes.

Kyoto Encyclopedia of Genes and Genomes (KEGG) is a knowledge base for systematic analysis of gene functions in terms of the networks of genes and molecules. The major component of KEGG is the PATHWAY database that consists of graphical diagrams of biochemical pathways including most of the known metabolic pathways and some of the known regulatory pathways. The pathway information is also represented by the ortholog group tables summarizing orthologous and paralogous gene groups among different organisms. KEGG maintains the GENES database for the gene catalogs of all organisms with complete genomes and selected organisms with partial genomes, which are continuously re-annotated, as well as the LIGAND database for chemical compounds and enzymes. Each gene catalog is associated with the graphical genome map for chromosomal locations that is represented by Java applet. In addition to the data collection efforts, KEGG develops and provides various computational tools, such as for reconstructing biochemical pathways from the complete genome sequence and for predicting gene regulatory networks from the gene expression profiles. The KEGG databases are daily updated and made freely available (http://www.genome.ad.jp/kegg/).     

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome‐scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence‐level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.     

4种模式植物LRR VIII-2亚家族基因的鉴定和进化历史分析

全基因组重复与串联重复是发生基因重复的重要机制, 也是基因组和遗传系统多样化的重要动力。LRR-RLK编码富含亮氨酸重复的类受体蛋白激酶, 是被子植物进化史上发生大规模扩张而形成的多基因家族。拟南芥(Arabidopsis thaliana) AtLRR-RLK包含15个亚家族, AtLRR VIII-2是其中发生串联重复比例最高的亚家族。通过分析拟南芥、杨树(Populus trichocarpa)、葡萄(Vitis vinifera)和番木瓜(Carica papaya) 4种模式植物中LRR VIII-2亚家族基因的扩张及差异保留情况, 结果显示, LRR VIII-2在杨树中的扩张程度最高, 在拟南芥和葡萄中的扩张程度居中, 但在番木瓜中发生丢失。拟南芥、杨树和葡萄LRR VIII-2亚家族具有旁系同源基因对, 但在番木瓜中未发现旁系同源基因。除杨树中的1对旁系同源基因外, 4种模式植物中LRR VIII-2亚家族的旁系和直系同源基因都受到较强的纯化选择作用。对LRR VIII-2亚家族进化历史的深入分析有助于理解基因重复在植物进化中的作用和意义, 可为预测同源基因功能及解析其它基因家族进化历史提供参考。     

New genes originated via multiple recombinational pathways in the beta-globin gene family of rodents

Species differences in the size or membership composition of multigene families can be attributed to lineage-specific additions of new genes via duplication, losses of genes via deletion or inactivation, and the creation of chimeric genes via domain shuffling or gene fusion. In principle, it should be possible to infer the recombinational pathways responsible for each of these different types of genomic change by conducting detailed comparative analyses of genomic sequence data. Here, we report an attempt to unravel the complex evolutionary history of the beta-globin gene family in a taxonomically diverse set of rodent species. The main objectives were: 1) to characterize the genomic structure of the beta-globin gene cluster of rodents; 2) to assign orthologous and paralogous relationships among duplicate copies of beta-like globin genes; and 3) to infer the specific recombinational pathways responsible for gene duplications, gene deletions, and the creation of chimeric fusion genes. Results of our comparative genomic analyses revealed that variation in gene family size among rodent species is mainly attributable to the differential gain and loss of later expressed beta-globin genes via unequal crossing-over. However, two distinct recombinational mechanisms were implicated in the creation of chimeric fusion genes. In muroid rodents, a chimeric gamma/epsilon fusion gene was created by unequal crossing-over between the embryonic epsilon- and gamma-globin genes. Interestingly, this gamma/epsilon fusion gene was generated in the same fashion as the "anti-Lepore" 5'-delta-(beta/delta)-beta-3' duplication mutant in humans (the reciprocal exchange product of the pathological hemoglobin Lepore deletion mutant). By contrast, in the house mouse, Mus musculus, a chimeric beta/delta fusion pseudogene was created by a beta-globin --> delta-globin gene conversion event. Although the gamma/epsilon and beta/delta fusion genes share a similar chimeric gene structure, they originated via completely different recombinational pathways.     

Rapid rates of lineage-specific gene duplication and deletion in the alpha-globin gene family

Phylogeny reconstructions of the globin gene families have revealed that paralogous genes within species are often more similar to one another than they are to their orthologous counterparts in closely related species. This pattern has been previously attributed to mechanisms of concerted evolution such as interparalog gene conversion that homogenize sequence variation between tandemly duplicated genes and therefore create the appearance of recent common ancestry. Here we report a comparative genomic analysis of the alpha-globin gene family in mammals that reveal a surprisingly high rate of lineage-specific gene duplication and deletion via unequal crossing-over. Results of our analysis reveal that patterns of sequence similarity between paralogous alpha-like globin genes from the same species are only partly explained by concerted evolution between preexisting gene duplicates. In a number of cases, sequence similarity between paralogous sequences from the same species is attributable to recent ancestry between the products of de novo gene duplications. As a result of this surprisingly rapid rate of gene gain and loss, many mammals possess alpha-like globin genes that have no orthologous counterparts in closely related species. The resultant variation in gene copy number among species may represent an important source of regulatory variation that affects physiologically important aspects of blood oxygen transport and aerobic energy metabolism.     

Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes

Taste reception is fundamental to diet selection in many animals. The genetic basis underlying the evolution and diversity of taste reception, however, is not well understood. Recent discoveries of T1R sweet/umami receptor genes and T2R bitter receptor genes in humans and mice provided an opportunity to address this question. Here, we report the identification of 20 putatively functional T1R genes and 167 T2R genes from the genome sequences of nine vertebrates, including three fishes, one amphibian, one bird, and four mammals. Our comparative genomic analysis shows that orthologous T1R sequences are relatively conserved in evolution and that the T1R gene repertoire remains virtually constant in size across most vertebrates, except for the loss of the T1R2 sweet receptor gene in the sweet-insensitive chicken and the absence of all T1R genes in the tongueless western clawed frog. In contrast, orthologous T2R sequences are more variable, and the T2R repertoire diverges tremendously among species, from only three functional genes in the chicken to 49 in the frog. These evolutionary patterns suggest the relative constancy in the number and type of sweet and umami tastants encountered by various vertebrates or low binding specificities of T1Rs but a large variation in the number and type of bitter compounds detected by different species. Although the rate of gene duplication is much lower in T1Rs than in T2Rs, signals of positive selection are detected during the functional divergences of paralogous T1Rs, as was previously found among paralogous T2Rs. Thus, functional divergence and specialization of taste receptors generally occurred via adaptive evolution.     

Campylobacter fetus subspecies: Comparative genomics and prediction of potential virulence targets

The genus Campylobacter contains pathogens causing a wide range of diseases, targeting both humans and animals. Among them, the Campylobacter fetus subspecies fetus and venerealis deserve special attention, as they are the etiological agents of human bacterial gastroenteritis and bovine genital campylobacteriosis, respectively. We compare the whole genomes of both subspecies to get insights into genomic architecture, phylogenetic relationships, genome conservation and core virulence factors. Pan-genomic approach was applied to identify the core- and pan-genome for both C. fetus subspecies and members of the genus. The C. fetus subspecies conserved (76%) proteome were then analyzed for their subcellular localization and protein functions in biological processes. Furthermore, with pathogenomic strategies, unique candidate regions in the genomes and several potential core-virulence factors were identified. The potential candidate factors identified for attenuation and/or subunit vaccine development against C. fetus subspecies contain: nucleoside diphosphate kinase (Ndk), type IV secretion systems (T4SS), outer membrane proteins (OMP), substrate binding proteins CjaA and CjaC, surface array proteins, sap gene, and cytolethal distending toxin (CDT). Significantly, many of those genes were found in genomic regions with signals of horizontal gene transfer and, therefore, predicted as putative pathogenicity islands. We found CRISPR loci and dam genes in an island specific for C. fetus subsp. fetus, and T4SS and sap genes in an island specific for C. fetus subsp. venerealis. The genomic variations and potential core and unique virulence factors characterized in this study would lead to better insight into the species virulence and to more efficient use of the candidates for antibiotic, drug and vaccine development.     

Ancestral paralogs and pseudoparalogs and their role in the emergence of the eukaryotic cell

Gene duplication is a crucial mechanism of evolutionary innovation. A substantial fraction of eukaryotic genomes consists of paralogous gene families. We assess the extent of ancestral paralogy, which dates back to the last common ancestor of all eukaryotes, and examine the origins of the ancestral paralogs and their potential roles in the emergence of the eukaryotic cell complexity. A parsimonious reconstruction of ancestral gene repertoires shows that 4137 orthologous gene sets in the last eukaryotic common ancestor (LECA) map back to 2150 orthologous sets in the hypothetical first eukaryotic common ancestor (FECA) [paralogy quotient (PQ) of 1.92]. Analogous reconstructions show significantly lower levels of paralogy in prokaryotes, 1.19 for archaea and 1.25 for bacteria. The only functional class of eukaryotic proteins with a significant excess of paralogous clusters over the mean includes molecular chaperones and proteins with related functions. Almost all genes in this category underwent multiple duplications during early eukaryotic evolution. In structural terms, the most prominent sets of paralogs are superstructure-forming proteins with repetitive domains, such as WD-40 and TPR. In addition to the true ancestral paralogs which evolved via duplication at the onset of eukaryotic evolution, numerous pseudoparalogs were detected, i.e. homologous genes that apparently were acquired by early eukaryotes via different routes, including horizontal gene transfer (HGT) from diverse bacteria. The results of this study demonstrate a major increase in the level of gene paralogy as a hallmark of the early evolution of eukaryotes.     

Prediction of missing enzyme genes in a bacterial metabolic network. Reconstruction of the lysine-degradation pathway of Pseudomonas aeruginosa

The metabolic network is an important biological network which consists of enzymes and chemical compounds. However, a large number of metabolic pathways remains unknown, and most organism-specific metabolic pathways contain many missing enzymes. We present a novel method to identify the genes coding for missing enzymes using available genomic and chemical information from bacterial genomes. The proposed method consists of two steps: (a) estimation of the functional association between the genes with respect to chromosomal proximity and evolutionary association, using supervised network inference; and (b) selection of gene candidates for missing enzymes based on the original candidate score and the chemical reaction information encoded in the EC number. We applied the proposed methods to infer the metabolic network for the bacteria Pseudomonas aeruginosa from two genomic datasets: gene position and phylogenetic profiles. Next, we predicted several missing enzyme genes to reconstruct the lysine-degradation pathway in P. aeruginosa using EC number information. As a result, we identified PA0266 as a putative 5-aminovalerate aminotransferase (EC 2.6.1.48) and PA0265 as a putative glutarate semialdehyde dehydrogenase (EC 1.2.1.20). To verify our prediction, we conducted biochemical assays and examined the activity of the products of the predicted genes, PA0265 and PA0266, in a coupled reaction. We observed that the predicted gene products catalyzed the expected reactions; no activity was seen when both gene products were omitted from the reaction.     

Going nuclear: gene family evolution and vertebrate phylogeny reconciled

Gene duplications have been common throughout vertebrate evolution, introducing paralogy and so complicating phylogenetic inference from nuclear genes. Reconciled trees are one method capable of dealing with paralogy, using the relationship between a gene phylogeny and the phylogeny of the organisms containing those genes to identify gene duplication events. This allows us to infer phylogenies from gene families containing both orthologous and paralogous copies. Vertebrate phylogeny is well understood from morphological and palaeontological data, but studies using mitochondrial sequence data have failed to reproduce this classical view. Reconciled tree analysis of a database of 118 vertebrate gene families supports a largely classical vertebrate phylogeny.     

The effect of paralogous lineages on the application of reconciliation analysis by cophylogeny mapping

Paralogy defines similarity caused by duplication rather than common descent and is well known in the case of paralogous gene copies within a single genome. The term is here extended to paralogous lineages of associates within a single host. The phylogenies of four genera within the Herpesviridae were reconciled with host phylogenies using cophylogenetic mapping. The observed correspondence for each pair of phylogenies was evaluated through randomization of the viral phylogeny and demonstrated to be greater than expected by chance. A simulation study was then carried out to assess the influence of paralogous lineages on the efficacy of reconciliation analysis. Combining viral taxa from different genera that infected common hosts introduced incongruence into the cophylogenies and reduced both the minimum and maximum observed number of codivergence events relative to the initial analysis of orthologous clades. However, at an average sample size this did not alter the fundamental significance of observed correspondence. With smaller sample sizes, the number of orthologous taxa selected at random from the pool of taxa was reduced. False-negative results then increased in proportion from 0.02 to 0.33. These results demonstrated that reconciliation analysis is robust under conditions of paralogy at "normal" sample sizes but is adversely affected by a combination of paralogy and low sample size. Consideration of phylogenies for Papillomavirus, Atadenovirus, and Mastadenovirus suggest that paralogous lineages may be a widespread phenomenon among DNA viruses and that duplication irrespective of host speciation is an important cause of viral diversification.     

Box H/ACA snoRNAs are preferred substrates for the trimethylguanosine synthase in the divergent unicellular eukaryote Trichomonas vaginalis

The 2,2,7-trimethylguanosine caps of eukaryal snRNAs and snoRNA are formed by the enzyme Tgs1, which catalyzes sequential guanine-N2 methylations of m(7)G caps. Atypically, in the divergent unicellular eukaryote Trichomonas vaginalis, spliceosomal snRNAs lack a guanosine cap and the recombinant T. vaginalis trimethylguanosine synthase (TvTgs) produces only m(2,7)G in vitro. Here, we show by direct metabolic labeling that endogenous T. vaginalis RNAs contain m(7)G, m(2,7)G, and m(2,2,7)G caps. Immunodepletion of TvTgs from cell extracts and TvTgs add-back experiments demonstrate that TvTgs produces m(2,7)G and m(2,2,7)G caps. Expression of TvTgs in yeast tgs1Δ cells leads to the formation of m(2,7)G and m(2,2,7)G caps and complementation of the lethality of a tgs1Δ mud2Δ strain. Whereas TvTgs is present in the nucleus and cytosol of T. vaginalis cells, TMG-containing RNAs are localized primarily in the nucleolus. Molecular cloning of anti-TMG affinity-purified T. vaginalis RNAs identified 16 box H/ACA snoRNAs, which are implicated in guiding RNA pseudouridylation. The ensemble of new T. vaginalis H/ACA snoRNAs allowed us to predict and partially validate an extensive map of pseudouridines in T. vaginalis rRNA.     

Identification of antigenic proteins in Trichomonas vaginalis

Trichomoniasis is a sexually transmitted disease due to infection with Trichomonas vaginalis, and it can cause serious consequences for women's health. To study the virulence factors of this pathogen, T. vaginalis surface proteins were investigated using polyclonal antibodies specific to the membrane fractions of T. vaginalis. The T. vaginalis expression library was constructed by cloning the cDNA derived from mRNA of T. vaginalis into a phage λ Uni-ZAP XR vector, and then used for immunoscreening with the anti-membrane proteins of T. vaginalis antibodies. The immunoreactive proteins identified included adhesion protein AP65-1, α-actinin, kinesin-associated protein, teneurin, and 2 independent hypothetical proteins. Immunofluorescence assays showed that AP65-1, one of the identified immunogenic clones, is prevalent in the whole body of T. vaginalis. This study led us to identify T. vaginalis proteins which may stimulate immune responses by human cells.     

Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Reconstruction of phylogenetic relationships among recently diverged species is complicated by three general problems: segregation of polymorphisms that pre-date species divergence, gene flow during and after speciation, and intra-locus recombination. In light of these difficulties, the Y chromosome offers several important advantages over other genomic regions as a source of phylogenetic information. These advantages include the absence of recombination, rapid coalescence, and reduced opportunity for interspecific introgression due to hybrid male sterility. In this report, we test the phylogenetic utility of Y-chromosomal sequences in two groups of closely related and partially inter-fertile Drosophila species. In the D. bipectinata species complex, Y-chromosomal loci unambiguously recover the phylogeny most consistent with previous multi-locus analysis and with reproductive relationships, and show no evidence of either post-speciation gene flow or persisting ancestral polymorphisms. In the D. simulans species complex, the situation is complicated by the duplication of at least one Y-linked gene region, followed by intrachromosomal recombination between the duplicate genes that scrambles their genealogy. We suggest that Y-chromosomal sequences are a useful tool for resolving phylogenetic relationships among recently diverged species, especially in male-heterogametic organisms that conform to Haldane's rule. However, duplication of Y-linked genes may not be uncommon, and special care should be taken to distinguish between orthologous and paralogous sequences.     

Comparative genomics of two closely related unicellular thermo-acidophilic red algae, Galdieria sulphuraria and Cyanidioschyzon merolae, reveals the molecular basis of the metabolic flexibility of Galdieria sulphuraria and significant differences in carbohydrate metabolism of both algae

Unicellular algae serve as models for the study and discovery of metabolic pathways, for the functional dissection of cell biological processes such as organellar division and cell motility, and for the identification of novel genes and gene functions. The recent completion of several algal genome sequences and expressed sequence tag collections and the establishment of nuclear and organellar transformation methods has opened the way for functional genomics approaches using algal model systems. The thermo-acidophilic unicellular red alga Galdieria sulphuraria represents a particularly interesting species for a genomics approach owing to its extraordinary metabolic versatility such as heterotrophic and mixotrophic growth on more than 50 different carbon sources and its adaptation to hot acidic environments. However, the ab initio prediction of genes required for unknown metabolic pathways from genome sequences is not trivial. A compelling strategy for gene identification is the comparison of similarly sized genomes of related organisms with different physiologies. Using this approach, candidate genes were identified that are critical to the metabolic versatility of Galdieria. Expressed sequence tags and high-throughput genomic sequence reads covering >70% of the G. sulphuraria genome were compared to the genome of the unicellular, obligate photoautotrophic red alga Cyanidioschyzon merolae. More than 30% of the Galdieria sequences did not relate to any of the Cyanidioschyzon genes. A closer inspection of these sequences revealed a large number of membrane transporters and enzymes of carbohydrate metabolism that are unique to Galdieria. Based on these data, it is proposed that genes involved in the uptake of reduced carbon compounds and enzymes involved in their metabolism are crucial to the metabolic flexibility of G. sulphuraria.     

Hydrogenosomal activity of Trichomonas vaginalis cultivated under different iron conditions

To evaluate whether iron concentration in TYM medium influence on hydrogenosomal enzyme gene expression and hydrogenosomal membrane potential of Trichomonas vaginalis, trophozoites were cultivated in irondepleted, normal and iron-supplemented TYM media. The mRNA of hydrogenosomal enzymes, such as pyruvate ferredoxin oxidoreductase (PFOR), hydrogenase, ferredoxin and malic enzyme, was increased with iron concentrations in T. vaginalis culture media, measured by RT-PCR. Hydrogenosomal membrane potentials measured with DiOC6 also showed similar tendency, e.g. T. vaginalis cultivated in iron-depleted and iron-supplemented media for 3 days showed a significantly reduced and enhanced hydrogenosomal membrane potential compared with that of normal TYM media, respectively. Therefore, it is suggested that iron may regulate hydrogenosomal activity through hydrogenosomal enzyme expression and hydrogenosomal membrane potential.     

On the incidence of intron loss and gain in paralogous gene families

Understanding gene duplication and gene structure evolution are fundamental goals of molecular evolutionary biology. A previous study by Babenko et al. (2004. Prevalence of intron gain over intron loss in the evolution of paralogous gene families. Nucleic Acids Res. 32:3724-3733) employed Dollo parsimony to infer spliceosomal intron losses and gains in paralogous gene families and concluded that there was a general excess of gains over losses. This result contrasts with patterns in orthologous genes, in which most lineages show an excess of intron losses over gains, suggesting the possibility of fundamentally different modes of intron evolution between orthologous and paralogous genes. We further studied the data and found a low level of intron position conservation with outgroups, and this led to problems with using Dollo parsimony to analyze the data. Statistical reanalysis of the data suggests, instead, that intron losses have outnumbered intron gains in paralogous gene families.