Set Of Orthologous Genes Research Articles

Orthopoxvirus Genome Evolution: The Role of Gene Loss

Poxviruses are highly successful pathogens, known to infect a variety of hosts. The family Poxviridae includes Variola virus, the causative agent of smallpox, which has been eradicated as a public health threat but could potentially reemerge as a bioterrorist threat. The risk scenario includes other animal poxviruses and genetically engineered manipulations of poxviruses. Studies of orthologous gene sets have established the evolutionary relationships of members within the Poxviridae family. It is not clear, however, how variations between family members arose in the past, an important issue in understanding how these viruses may vary and possibly produce future threats. Using a newly developed poxvirus-specific tool, we predicted accurate gene sets for viruses with completely sequenced genomes in the genus Orthopoxvirus. Employing sensitive sequence comparison techniques together with comparison of syntenic gene maps, we established the relationships between all viral gene sets. These techniques allowed us to unambiguously identify the gene loss/gain events that have occurred over the course of orthopoxvirus evolution. It is clear that for all existing Orthopoxvirus species, no individual species has acquired protein-coding genes unique to that species. All existing species contain genes that are all present in members of the species Cowpox virus and that cowpox virus strains contain every gene present in any other orthopoxvirus strain. These results support a theory of reductive evolution in which the reduction in size of the core gene set of a putative ancestral virus played a critical role in speciation and confining any newly emerging virus species to a particular environmental (host or tissue) niche.

Comparisons of early transcriptome responses to low-oxygen environments in three dicotyledonous plant species

Waterlogging is a serious impediment to crop productivity worldwide which acts to reduce oxygen levels in the rhizosphere due to the low diffusion rate of molecular oxygen in water. Plants respond to low oxygen through rapid and specific changes at both the transcriptional and translational levels. Transcriptional changes to low-oxygen (hypoxia) stress have been studied in a number of plant species using whole genome microarrays. Using transcriptome data from root tissue from early time points (4-5 h) from cotton (Gossypium hirsutum), Arabidopsis and gray poplar (Populus x canescens), we have identified a core set of orthologous genes that responded to hypoxia in similar ways between species, and others that showed species specific responses . Responses to hypoxia were most similar between Arabidopsis and cotton, while the waterlogging tolerant poplar species exhibited some significant differences.

Genome-wide expression analysis of Saccharomyces pastorianus orthologous genes using oligonucleotide microarrays

The lager brewing yeast, Saccharomyces pastorianus, an allopolyploid species hybrid, contains 2 diverged sub-genomes; one derived from Saccharomyces cerevisiae (Sc-type) and the other from Saccharomyces bayanus (Sb-type). We analyzed the functional roles of these orthologous genes in determining the phenotypic features of S. pastorianus. We used a custom-made oligonucleotide microarray containing probes designed for both Sc-type and Sb-type ORFs for a comprehensive expression analysis of S. pastorianus in a pilot-scale fermentation. We showed a high degree of correlation between the expression levels and the expression changes for a majority of orthologous gene sets during the fermentation process. We screened the functional categories and metabolic pathways where Sc- or Sb-type genes have higher expression levels than the corresponding orthologous genes. Our data showed that, for example, pathways for sulfur metabolism, cellular import, and production of branched amino acids are dominated by Sb-type genes. This comprehensive expression analysis of orthologous genes can provide valuable insights on understanding the phenotype of S. pastorianus.

Evolution of reduced prokaryotic genomes and the minimal cell concept: Variations on a theme

Prokaryotic genomes of endosymbionts and parasites are examples of naturally evolved minimal cells, the study of which can shed light on life in its minimum form. Their diverse biology, their lack of a large set of orthologous genes and the existence of essential linage (and environmentally) specific genes all illustrate the diversity of genes building up naturally evolved minimal cells. This conclusion is reinforced by the fact that sometimes the same essential function is performed by genes from different evolutionary origins. Nevertheless, all cells perform a set of life-essential functions however different their cell machinery and environment in which they thrive. An upcoming challenge for biologists will be to discern, by studying differences and similarities in current biodiversity, how cells with reduced genomes have adapted while retaining all basic life-supporting functions.

Systematic discovery of nonobvious human disease models through orthologous phenotypes

Biologists have long used model organisms to study human diseases, particularly when the model bears a close resemblance to the disease. We present a method that quantitatively and systematically identifies nonobvious equivalences between mutant phenotypes in different species, based on overlapping sets of orthologous genes from human, mouse, yeast, worm, and plant (212,542 gene-phenotype associations). These orthologous phenotypes, or phenologs, predict unique genes associated with diseases. Our method suggests a yeast model for angiogenesis defects, a worm model for breast cancer, mouse models of autism, and a plant model for the neural crest defects associated with Waardenburg syndrome, among others. Using these models, we show that SOX13 regulates angiogenesis, and that SEC23IP is a likely Waardenburg gene. Phenologs reveal functionally coherent, evolutionarily conserved gene networks-many predating the plant-animal divergence-capable of identifying candidate disease genes.

The Evolution of tRNA Genes in Drosophila

The structure and function of transfer RNA (tRNA) genes have been extensively studied for several decades, yet the general mechanisms controlling tRNA gene family evolution remain unclear, primarily because previous phylogenetics-based methods fail to distinguish between paralogs and orthologs that are highly similar in sequence. We have developed a system for identifying orthologs of tRNAs using flanking sequences to identify regions of conserved synteny and used it to annotate sets of orthologous tRNA genes across the 12 sequenced species of Drosophila. These data have allowed us to place the gains and losses of individual tRNA genes on each branch of the Drosophila tree and estimate rates of tRNA gene turnover. Our results show extensive rearrangement of the Drosophila tRNA gene complement over the last 60 My. We estimate a combined average rate of 2.18 ± 0.10 tRNA gene gains and losses per million years across the Drosophila lineage. We have identified 192 tRNAs that are ancestral to the genus, of which 157 are “core” tRNAs conserved in at least 11 of 12 extant species. We provide evidence that the core set of tRNA genes encode a nearly complete set of anticodons and have different properties from other “peripheral” tRNA genes, such as preferential location outside large tRNA clusters and higher sequence conservation. We also demonstrate that tRNA isoacceptor and alloacceptor changes by anticodon shifts have occurred several times in Drosophila, annotating 16 such events in functional tRNAs during the evolution of the genus.

Positively Selected Disease Response Orthologous Gene Sets in the Cereals Identified Using Sorghum bicolor L. Moench Expression Profiles and Comparative Genomics

Disease response genes (DRGs) diverge under recurrent positive selection as a result of a molecular arms race between hosts and pathogens. Most of these studies were conducted in animals, and few defense genes have been shown to evolve adaptively in plants. To test for adaptation in the molecules mediating disease resistance in the cereals, we first combined information from the expression pattern of Sorghum bicolor genes and from divergence to the full genome of rice to identify candidate DRGs. We then used evolutionary analyses of orthologous gene sets from several grass species, to determine whether the DRGs show signals of positive selection and the residues targeted. We found 140 divergent genes upregulated under biotic stress in S. bicolor by evaluating the relative abundance of expressed sequence tags in different libraries and comparing them with rice genes. For 10 of these genes, we found sets of orthologs including sequences from rice and three other cereals; six genes showed a pattern of substitution that was consistent with positive selection. Three of these genes, a thaumatin, a peroxidase, and a barley mlo homolog, are known antifungal proteins. The other three genes with evidence of positive selection were a MCM-1 agamous deficiens SRF- (MADS) box transcription factor, an eIF5 translation initiation factor, and a gene of unknown function but with evidence of expression during stress. Permutation analyses, using different ortholog and paralog sequences, consistently identified five positively selected codons in the peroxidase, a member of a cluster of genes and a large gene family. We mapped the positively selected residues onto the structure of the peroxidase and thaumatin and found that all sites are on the surface of these proteins and several are close to biochemically determined active sites. Identifying new positively selected plant disease resistance genes and the critical amino acid sites provides a basis for functional studies that may increase our understanding of their underlying molecular mechanisms of action. Additionally, it may lead to the identification of individuals having variation at functionally important sites, as well as eventually using this information in the rational design and engineering of proteins involved in plant disease resistance.

Databases of homologous gene families for comparative genomics

BackgroundComparative genomics is a central step in many sequence analysis studies, from gene annotation and the identification of new functional regions in genomes, to the study of evolutionary processes at the molecular level (speciation, single gene or whole genome duplications, etc.) and phylogenetics. In that context, databases providing users high quality homologous families and sequence alignments as well as phylogenetic trees based on state of the art algorithms are becoming indispensable.MethodsWe developed an automated procedure allowing massive all-against-all similarity searches, gene clustering, multiple alignments computation, and phylogenetic trees construction and reconciliation. The application of this procedure to a very large set of sequences is possible through parallel computing on a large computer cluster.ResultsThree databases were developed using this procedure: HOVERGEN, HOGENOM and HOMOLENS. These databases share the same architecture but differ in their content. HOVERGEN contains sequences from vertebrates, HOGENOM is mainly devoted to completely sequenced microbial organisms, and HOMOLENS is devoted to metazoan genomes from Ensembl. Access to the databases is provided through Web query forms, a general retrieval system and a client-server graphical interface. The later can be used to perform tree-pattern based searches allowing, among other uses, to retrieve sets of orthologous genes. The three databases, as well as the software required to build and query them, can be used or downloaded from the PBIL (Pôle Bioinformatique Lyonnais) site at .

Are Protein Domains Modules of Lateral Genetic Transfer?

BackgroundIn prokaryotes and some eukaryotes, genetic material can be transferred laterally among unrelated lineages and recombined into new host genomes, providing metabolic and physiological novelty. Although the process is usually framed in terms of gene sharing (e.g. lateral gene transfer, LGT), there is little reason to imagine that the units of transfer and recombination correspond to entire, intact genes. Proteins often consist of one or more spatially compact structural regions (domains) which may fold autonomously and which, singly or in combination, confer the protein's specific functions. As LGT is frequent in strongly selective environments and natural selection is based on function, we hypothesized that domains might also serve as modules of genetic transfer, i.e. that regions of DNA that are transferred and recombined between lineages might encode intact structural domains of proteins.Methodology/Principal FindingsWe selected 1,462 orthologous gene sets representing 144 prokaryotic genomes, and applied a rigorous two-stage approach to identify recombination breakpoints within these sequences. Recombination breakpoints are very significantly over-represented in gene sets within which protein domain-encoding regions have been annotated. Within these gene sets, breakpoints significantly avoid the domain-encoding regions (domons), except where these regions constitute most of the sequence length. Recombination breakpoints that fall within longer domons are distributed uniformly at random, but those that fall within shorter domons may show a slight tendency to avoid the domon midpoint. As we find no evidence for differential selection against nucleotide substitutions following the recombination event, any bias against disruption of domains must be a consequence of the recombination event per se.Conclusions/SignificanceThis is the first systematic study relating the units of LGT to structural features at the protein level. Many genes have been interrupted by recombination following inter-lineage genetic transfer, during which the regions within these genes that encode protein domains have not been preferentially preserved intact. Protein domains are units of function, but domons are not modules of transfer and recombination. Our results demonstrate that LGT can remodel even the most functionally conservative modules within genomes.

Gene Order Phylogeny of the Genus Prochlorococcus

BackgroundUsing gene order as a phylogenetic character has the potential to resolve previously unresolved species relationships. This character was used to resolve the evolutionary history within the genus Prochlorococcus, a group of marine cyanobacteria.Methodology/Principal FindingsOrthologous gene sets and their genomic positions were identified from 12 species of Prochlorococcus and 1 outgroup species of Synechococcus. From this data, inversion and breakpoint distance-based phylogenetic trees were computed by GRAPPA and FastME. Statistical support of the resulting topology was obtained by application of a 50% jackknife resampling technique. The result was consistent and congruent with nucleotide sequence-based and gene-content based trees. Also, a previously unresolved clade was resolved, that of MIT9211 and SS120.Conclusions/SignificanceThis is the first study to use gene order data to resolve a bacterial phylogeny at the genus level. It suggests that the technique is useful in resolving the Tree of Life.

Trees of Trees: An Approach to Comparing Multiple Alternative Phylogenies

Phylogenetic analysis very commonly produces several alternative trees for a given fixed set of taxa. For example, different sets of orthologous genes may be analyzed, or the analysis may sample from a distribution of probable trees. This article describes an approach to comparing and visualizing multiple alternative phylogenies via the idea of a "tree of trees" or "meta-tree." A meta-tree clusters phylogenies with similar topologies together in the same way that a phylogeny clusters species with similar DNA sequences. Leaf nodes on a meta-tree correspond to the original set of phylogenies given by some analysis, whereas interior nodes correspond to certain consensus topologies. The construction of meta-trees is motivated by analogy with construction of a most parsimonious tree for DNA data, but instead of using DNA letters, in a meta-tree the characters are partitions or splits of the set of taxa. An efficient algorithm for meta-tree construction is described that makes use of a known relationship between the majority consensus and parsimony in terms of gain and loss of splits. To illustrate these ideas meta-trees are constructed for two datasets: a set of gene trees for species of yeast and trees from a bootstrap analysis of a set of gene trees in ray-finned fish. A software tool for constructing meta-trees and comparing alternative phylogenies is available online, and the source code can be obtained from the author.

The genome sequence of the model ascomycete fungus Podospora anserina

A 10X draft sequence of Podospora anserina genome shows highly dynamic evolution since its divergence from Neurospora crassa.

Complex evolution of the DAL5 transporter family

BackgroundGenes continuously duplicate and the duplicated copies remain in the genome or get deleted. The DAL5 subfamily of transmembrane transporter genes has eight known members in S. cerevisiae. All are putative anion:cation symporters of vitamins (such as allantoate, nicotinate, panthotenate and biotin). The DAL5 subfamily is an old and important group since it is represented in both Basidiomycetes ("mushrooms") and Ascomycetes ("yeast"). We studied the complex evolution of this group in species from the kingdom of fungi particularly the Ascomycetes.ResultsWe identified numerous gene duplications creating sets of orthologous and paralogous genes. In different lineages the DAL5 subfamily members expanded or contracted and in some lineages a specific member could not be found at all. We also observed a close relationship between the gene YIL166C and its homologs in the Saccharomyces sensu stricto species and two "wine spoiler" yeasts, Dekkera bruxellensis and Candida guilliermondi, which could possibly be the result of horizontal gene transfer between these distantly related species. In the analyses we detect several well defined groups without S. cerevisiae representation suggesting new gene members in this subfamily with perhaps altered specialization or function.ConclusionThe transmembrane DAL5 subfamily was found to have a very complex evolution in yeast with intra- and interspecific duplications and unusual relationships indicating specialization, specific deletions and maybe even horizontal gene transfer. We believe that this group will be important in future investigations of evolution in fungi and especially the evolution of transmembrane proteins and their specialization.

Gramene: a growing plant comparative genomics resource

Gramene (www.gramene.org) is a curated resource for genetic, genomic and comparative genomics data for the major crop species, including rice, maize, wheat and many other plant (mainly grass) species. Gramene is an open-source project. All data and software are freely downloadable through the ftp site (ftp.gramene.org/pub/gramene) and available for use without restriction. Gramene's core data types include genome assembly and annotations, other DNA/mRNA sequences, genetic and physical maps/markers, genes, quantitative trait loci (QTLs), proteins, ontologies, literature and comparative mappings. Since our last NAR publication 2 years ago, we have updated these data types to include new datasets and new connections among them. Completely new features include rice pathways for functional annotation of rice genes; genetic diversity data from rice, maize and wheat to show genetic variations among different germplasms; large-scale genome comparisons among Oryza sativa and its wild relatives for evolutionary studies; and the creation of orthologous gene sets and phylogenetic trees among rice, Arabidopsis thaliana, maize, poplar and several animal species (for reference purpose). We have significantly improved the web interface in order to provide a more user-friendly browsing experience, including a dropdown navigation menu system, unified web page for markers, genes, QTLs and proteins, and enhanced quick search functions.

A Tree of Life Based on Protein Domain Organizations

It is desirable to estimate a tree of life, a species tree including all available species in the 3 superkingdoms, Archaea, Bacteria, and Eukaryota, using not a limited number of genes but full-scale genome information. Here, we report a new method for constructing a tree of life based on protein domain organizations, that is, sequential order of domains in a protein, of all proteins detected in a genome of an organism. The new method is free from the identification of orthologous gene sets and therefore does not require the burdensome and error-prone computation. By pairwise comparisons of the repertoires of protein domain organizations of 17 archaeal, 136 bacterial, and 14 eukaryotic organisms, we computed evolutionary distances among them and constructed a tree of life. Our tree shows monophyly in Archaea, Bacteria, and Eukaryota and then monophyly in each of eukaryotic kingdoms and in most bacterial phyla. In addition, the branching pattern of the bacterial phyla in our tree is consistent with the widely accepted bacterial taxonomy and is very close to other genome-based trees. A couple of inconsistent aspects between the traditional trees and the genome-based trees including ours, however, would perhaps urge to revise the conventional view, particularly on the phylogenetic positions of hyperthermophiles.

Phase-dependent nucleotide substitution in protein-coding sequences

It is well known that due to the degeneracy of genetic code, most of the silent substitutions appear in the third codon position, so the mutation frequency of the third codon position is much higher than that of the first two positions. However, it remains unknown whether the directionality of point mutation in three codon positions is similar or not. In this paper, through analyzing 15 sets of orthologous genes, it is revealed that most of the substitution types are significantly different between any two codon positions, especially between the 2nd and the 3rd phases. Furthermore, the average frequencies of each type of substitution calculated from the fifteen sets of orthologous genes are similar to those identified in single nucleotide polymorphisms (SNPs) of human and mouse genome. The present analyses suggest that the nucleotide substitution in protein-coding sequences is not only context-dependent (so called neighboring-nucleotide effects), but also phase-dependent, which is of significance to improving the prevalent nucleotide-evolution models.

Detecting non-coding selective pressure in coding regions

BackgroundComparative genomics approaches, where orthologous DNA regions are compared and inter-species conserved regions are identified, have proven extremely powerful for identifying non-coding regulatory regions located in intergenic or intronic regions. However, non-coding functional elements can also be located within coding region, as is common for exonic splicing enhancers, some transcription factor binding sites, and RNA secondary structure elements affecting mRNA stability, localization, or translation. Since these functional elements are located in regions that are themselves highly conserved because they are coding for a protein, they generally escaped detection by comparative genomics approaches.ResultsWe introduce a comparative genomics approach for detecting non-coding functional elements located within coding regions. Codon evolution is modeled as a mixture of codon substitution models, where each component of the mixture describes the evolution of codons under a specific type of coding selective pressure. We show how to compute the posterior distribution of the entropy and parsimony scores under this null model of codon evolution. The method is applied to a set of growth hormone 1 orthologous mRNA sequences and a known exonic splicing elements is detected. The analysis of a set of CORTBP2 orthologous genes reveals a region of several hundred base pairs under strong non-coding selective pressure whose function remains unknown.ConclusionNon-coding functional elements, in particular those involved in post-transcriptional regulation, are likely to be much more prevalent than is currently known. With the numerous genome sequencing projects underway, comparative genomics approaches like that proposed here are likely to become increasingly powerful at detecting such elements.

OrthoMaM: A database of orthologous genomic markers for placental mammal phylogenetics

BackgroundMolecular sequence data have become the standard in modern day phylogenetics. In particular, several long-standing questions of mammalian evolutionary history have been recently resolved thanks to the use of molecular characters. Yet, most studies have focused on only a handful of standard markers. The availability of an ever increasing number of whole genome sequences is a golden mine for modern systematics. Genomic data now provide the opportunity to select new markers that are potentially relevant for further resolving branches of the mammalian phylogenetic tree at various taxonomic levels.DescriptionThe EnsEMBL database was used to determine a set of orthologous genes from 12 available complete mammalian genomes. As targets for possible amplification and sequencing in additional taxa, more than 3,000 exons of length > 400 bp have been selected, among which 118, 368, 608, and 674 are respectively retrieved for 12, 11, 10, and 9 species. A bioinformatic pipeline has been developed to provide evolutionary descriptors for these candidate markers in order to assess their potential phylogenetic utility. The resulting OrthoMaM (Orthologous Mammalian Markers) database can be queried and alignments can be downloaded through a dedicated web interface .ConclusionThe importance of marker choice in phylogenetic studies has long been stressed. Our database centered on complete genome information now makes possible to select promising markers to a given phylogenetic question or a systematic framework by querying a number of evolutionary descriptors. The usefulness of the database is illustrated with two biological examples. First, two potentially useful markers were identified for rodent systematics based on relevant evolutionary parameters and sequenced in additional species. Second, a complete, gapless 94 kb supermatrix of 118 orthologous exons was assembled for 12 mammals. Phylogenetic analyses using probabilistic methods unambiguously supported the new placental phylogeny by retrieving the monophyly of Glires, Euarchontoglires, Laurasiatheria, and Boreoeutheria. Muroid rodents thus do not represent a basal placental lineage as it was mistakenly reasserted in some recent phylogenomic analyses based on fewer taxa. We expect the OrthoMaM database to be useful for further resolving the phylogenetic tree of placental mammals and for better understanding the evolutionary dynamics of their genomes, i.e., the forces that shaped coding sequences in terms of selective constraints.

Conservation of alternative polyadenylation patterns in mammalian genes

BackgroundAlternative polyadenylation is a widespread mechanism contributing to transcript diversity in eukaryotes. Over half of mammalian genes are alternatively polyadenylated. Our understanding of poly(A) site evolution is limited by the lack of a reliable identification of conserved, equivalent poly(A) sites among species. We introduce here a working definition of conserved poly(A) sites as sites that are both (i) properly aligned in human and mouse orthologous 3' untranslated regions (UTRs) and (ii) supported by EST or cDNA data in both species.ResultsWe identified about 4800 such conserved poly(A) sites covering one third of the orthologous gene set studied. Characteristics of conserved poly(A) sites such as processing efficiency and tissue-specificity were analyzed. Conserved sites show a higher processing efficiency but no difference in tissular distribution when compared to non-conserved sites. In general, alternative poly(A) sites are species-specific and involve minor, non-conserved sites that are unlikely to play essential roles. However, there are about 500 genes with conserved tandem poly(A) sites. A significant fraction of these conserved tandems display a conserved arrangement of major/minor sites in their 3' UTR, suggesting that these alternative 3' ends may be under selection.ConclusionThis analysis allows us to identify potential functional alternative poly(A) sites and provides clues on the selective mechanisms at play in the appearance of multiple poly(A) sites and their maintenance in the 3' UTRs of genes.

On conditioned reconstruction, gene content data, and the recovery of fusion genomes

Conditioned reconstruction (CR)1 represents a new phylogenetic method that has been presented as a means of utilizing vast amounts of gene absence/presence data to reconstruct phylogenetic relationships and to directly study the inXuence of genome fusion on evolution (Lake and Rivera, 2004; Rivera and Lake, 2004; Simonson et al., 2005). In the Wrst direct application of CR, the results were stated to unambiguously support the role of an archaeal and eubacterial genome fusion event in the ancestry of the eukaryotic genome (McInerney and Wilkinson, 2005; Rivera and Lake, 2004; Simonson et al., 2005). In the present manuscript, we outline the basic components of CR, discuss the logic behind their application, and subsequently identify concerns with aspects of data interpretation and the current use of conditioning genomes and aligned networks in the study of genome fusion. Conditioned reconstruction begins with the selection of a conditioning genome (CG) that will represent the full set of orthologous genes coded during matrix development (Lake and Rivera, 2004). Strong emphasis has been placed on the strict use of relatively small genomes as conditioning genomes (Lake and Rivera, 2004). This recommendation is based on the observation that the use of relatively large conditioning genomes can induce an artifact referred to as “big genome attraction,” which can mislead phylogenetic inference. All other genomes that will be used in the analysis are scored for the absence or presence of a putative ortholog to each gene in the CG. Orthologs present in other genomes included in the