Orthologous Pairs Research Articles

Discrimination of thermophilic and mesophilic proteins

Background There is a considerable literature on the source of the thermostability of proteins from thermophilic organisms. Understanding the mechanisms for this thermostability would provide insights into proteins generally and permit the design of synthetic hyperstable biocatalysts.Results We have systematically tested a large number of sequence and structure derived quantities for their ability to discriminate thermostable proteins from their non-thermostable orthologs using sets of mesophile-thermophile ortholog pairs. Most of the quantities tested correspond to properties previously reported to be associated with thermostability. Many of the structure related properties were derived from the Delaunay tessellation of protein structures.Conclusions Carefully selected sequence based indices discriminate better than purely structure based indices. Combined sequence and structure based indices improve performance somewhat further. Based on our analysis, the strongest contributors to thermostability are an increase in ion pairs on the protein surface and a more strongly hydrophobic interior.

Functional Shifts in Insect microRNA Evolution

MicroRNAs (miRNAs) are short endogenous RNA molecules that regulate gene expression at the posttranscriptional level and have been shown to play critical roles during animal development. The identification and comparison of miRNAs in metazoan species are therefore paramount for our understanding of the evolution of body plans. We have characterized 203 miRNAs from the red flour beetle Tribolium castaneum by deep sequencing of small RNA libraries. We can conclude, from a single study, that the Tribolium miRNA set is at least 15% larger than that in the model insect Drosophila melanogaster (despite tens of high-throughput sequencing experiments in the latter). The rate of birth and death of miRNAs is high in insects. Only one-third of the Tribolium miRNA sequences are conserved in D. melanogaster, and at least 18 Tribolium miRNAs are conserved in vertebrates but lost in Drosophila. More than one-fifth of miRNAs that are conserved between Tribolium and Drosophila exhibit changes in the transcription, genomic organization, and processing patterns that lead to predicted functional shifts. For example, 13% of conserved miRNAs exhibit seed shifting, and we describe arm-switching events in 11% of orthologous pairs. These shifts fundamentally change the predicted targets and therefore function of orthologous miRNAs. In general, Tribolium miRNAs are more representative of the insect ancestor than Drosophila miRNAs and are more conserved in vertebrates.

The CAZyome of Phytophthora spp.: A comprehensive analysis of the gene complement coding for carbohydrate-active enzymes in species of the genus Phytophthora

BackgroundEnzymes involved in carbohydrate metabolism include Carbohydrate esterases (CE), Glycoside hydrolases (GH), Glycosyl transferases (GT), and Polysaccharide lyases (PL), commonly referred to as carbohydrate-active enzymes (CAZymes). The CE, GH, and PL superfamilies are also known as cell wall degrading enzymes (CWDE) due to their role in the disintegration of the plant cell wall by bacterial and fungal pathogens. In Phytophthora infestans, penetration of the plant cells occurs through a specialized hyphal structure called appressorium; however, it is likely that members of the genus Phytophthora also use CWDE for invasive growth because hyphal forces are below the level of tensile strength exhibited by the plant cell wall. Because information regarding the frequency and distribution of CAZyme coding genes in Phytophthora is currently unknown, we have scanned the genomes of P. infestans, P. sojae, and P. ramorum for the presence of CAZyme-coding genes using a homology-based approach and compared the gene collinearity in the three genomes. In addition, we have tested the expression of several genes coding for CE in cultures grown in vitro.ResultsWe have found that P. infestans, P. sojae and P. ramorum contain a total of 435, 379, and 310 CAZy homologs; in each genome, most homologs belong to the GH superfamily. Most GH and PL homologs code for enzymes that hydrolyze substances present in the pectin layer forming the middle lamella of the plant cells. In addition, a significant number of CE homologs catalyzing the deacetylation of compounds characteristic of the plant cell cuticle were found. In general, a high degree of gene location conservation was observed, as indicated by the presence of sequential orthologous pairs in the three genomes. Such collinearity was frequently observed among members of the GH superfamily. On the other hand, the CE and PL superfamilies showed less collinearity for some of their putative members. Quantitative PCR experiments revealed that all genes are expressed in P. infestans when this pathogen grown in vitro. However, the levels of expression vary considerably and are lower than the expression levels observed for the constitutive control.ConclusionsIn conclusion, we have identified a highly complex set of CAZy homologs in the genomes of P. infestans, P. sojae, and P. ramorum, a significant number of which could play roles critical for pathogenicity, by participating in the degradation of the plant cell wall.

Identification of mammalian orthologs using local synteny

BackgroundAccurate determination of orthology is central to comparative genomics. For vertebrates in particular, very large gene families, high rates of gene duplication and loss, multiple mechanisms of gene duplication, and high rates of retrotransposition all combine to make inference of orthology between genes difficult. Many methods have been developed to identify orthologous genes, mostly based upon analysis of the inferred protein sequence of the genes. More recently, methods have been proposed that use genomic context in addition to protein sequence to improve orthology assignment in vertebrates. Such methods have been most successfully implemented in fungal genomes and have long been used in prokaryotic genomes, where gene order is far less variable than in vertebrates. However, to our knowledge, no explicit comparison of synteny and sequence based definitions of orthology has been reported in vertebrates, or, more specifically, in mammals.ResultsWe test a simple method for the measurement and utilization of gene order (local synteny) in the identification of mammalian orthologs by investigating the agreement between coding sequence based orthology (Inparanoid) and local synteny based orthology. In the 5 mammalian genomes studied, 93% of the sampled inter-species pairs were found to be concordant between the two orthology methods, illustrating that local synteny is a robust substitute to coding sequence for identifying orthologs. However, 7% of pairs were found to be discordant between local synteny and Inparanoid. These cases of discordance result from evolutionary events including retrotransposition and genome rearrangements.ConclusionsBy analyzing cases of discordance between local synteny and Inparanoid we show that local synteny can distinguish between true orthologs and recent retrogenes, can resolve ambiguous many-to-many orthology relationships into one-to-one ortholog pairs, and might be used to identify cases of non-orthologous gene displacement by retroduplicated paralogs.

Time warping of evolutionary distant temporal gene expression data based on noise suppression

BackgroundComparative analysis of genome wide temporal gene expression data has a broad potential area of application, including evolutionary biology, developmental biology, and medicine. However, at large evolutionary distances, the construction of global alignments and the consequent comparison of the time-series data are difficult. The main reason is the accumulation of variability in expression profiles of orthologous genes, in the course of evolution.ResultsWe applied Pearson distance matrices, in combination with other noise-suppression techniques and data filtering to improve alignments. This novel framework enhanced the capacity to capture the similarities between the temporal gene expression datasets separated by large evolutionary distances. We aligned and compared the temporal gene expression data in budding (Saccharomyces cerevisiae) and fission (Schizosaccharomyces pombe) yeast, which are separated by more then ~400 myr of evolution. We found that the global alignment (time warping) properly matched the duration of cell cycle phases in these distant organisms, which was measured in prior studies. At the same time, when applied to individual ortholog pairs, this alignment procedure revealed groups of genes with distinct alignments, different from the global alignment.ConclusionOur alignment-based predictions of differences in the cell cycle phases between the two yeast species were in a good agreement with the existing data, thus supporting the computational strategy adopted in this study. We propose that the existence of the alternative alignments, specific to distinct groups of genes, suggests presence of different synchronization modes between the two organisms and possible functional decoupling of particular physiological gene networks in the course of evolution.

Conserved repertoire of orthologous vomeronasal type 1 receptor genes in ruminant species

BackgroundIn mammals, pheromones play an important role in social and innate reproductive behavior within species. In rodents, vomeronasal receptor type 1 (V1R), which is specifically expressed in the vomeronasal organ, is thought to detect pheromones. The V1R gene repertoire differs dramatically between mammalian species, and the presence of species-specific V1R subfamilies in mouse and rat suggests that V1R plays a profound role in species-specific recognition of pheromones. In ruminants, however, the molecular mechanism(s) for pheromone perception is not well understood. Interestingly, goat male pheromone, which can induce out-of-season ovulation in anestrous females, causes the same pheromone response in sheep, and vice versa, suggesting that there may be mechanisms for detecting "inter-species" pheromones among ruminant species.ResultsWe isolated 23 goat and 21 sheep intact V1R genes based on sequence similarity with 32 cow V1R genes in the cow genome database. We found that all of the goat and sheep V1R genes have orthologs in their cross-species counterparts among these three ruminant species and that the sequence identity of V1R orthologous pairs among these ruminants is much higher than that of mouse-rat V1R orthologous pairs. Furthermore, all goat V1Rs examined thus far are expressed not only in the vomeronasal organ but also in the main olfactory epithelium.ConclusionOur results suggest that, compared with rodents, the repertoire of orthologous V1R genes is remarkably conserved among the ruminants cow, sheep and goat. We predict that these orthologous V1Rs can detect the same or closely related chemical compound(s) within each orthologous set/pair. Furthermore, all identified goat V1Rs are expressed in the vomeronasal organ and the main olfactory epithelium, suggesting that V1R-mediated ligand information can be detected and processed by both the main and accessory olfactory systems. The fact that ruminant and rodent V1Rs have distinct features suggests that ruminant and rodent V1Rs have evolved distinct functions.

Anopheles Fibrinogen-related Proteins Provide Expanded Pattern Recognition Capacity against Bacteria and Malaria Parasites

The fibrinogen-related protein family (FREP, also known as FBN) is an evolutionarily conserved immune gene family found in mammals and invertebrates. It is the largest pattern recognition receptor gene family in Anopheles gambiae, with as many as 59 putative members, while the Drosophila melanogaster genome has only 14 known FREP members. Our sequence and phylogenetic analysis suggest that this remarkable gene expansion in the mosquito is the result of tandem duplication of the fibrinogen domain. We found that the majority of the FREP genes displayed immune-responsive transcription after challenge with bacteria, fungi, or Plasmodium, and these expression patterns correlated strongly with gene phylogeny and chromosomal location. Using RNAi-mediated gene-silencing assays, we further demonstrated that some FREP members are essential factors of the mosquito innate immune system that are required for maintaining immune homeostasis, and members of this family have complementary and synergistic functions. One of the most potent anti-Plasmodium FREP proteins, FBN9, was found to interact with both Gram-negative and Gram-positive bacteria and strongly co-localized with both rodent and human malaria parasites in the mosquito midgut epithelium, suggesting that its defensive activity involves direct interaction with the pathogen. Interestingly, FBN9 formed dimers that bound to the bacterial surfaces with different affinities. Our findings indicate that the A. gambiae FREP gene family plays a central role in the mosquito innate immune system and provides an expanded pattern recognition and anti-microbial defense repertoire.

Dynamic evolution of V1R putative pheromone receptors between Mus musculus and Mus spretus

BackgroundThe mammalian vomeronasal organ (VNO) expresses two G-protein coupled receptor gene families that mediate pheromone responses, the V1R and V2R receptor genes. In rodents, there are ~150 V1R genes comprising 12 subfamilies organized in gene clusters at multiple chromosomal locations. Previously, we showed that several of these subfamilies had been extensively modulated by gene duplications, deletions, and gene conversions around the time of the evolutionary split of the mouse and rat lineages, consistent with the hypothesis that V1R repertoires might be involved in reinforcing speciation events. Here, we generated genome sequence for one large cluster containing two V1R subfamilies in Mus spretus, a closely related and sympatric species to Mus musculus, and investigated evolutionary change in these repertoires along the two mouse lineages.ResultsWe describe a comparison of spretus and musculus with respect to genome organization and synteny, as well as V1R gene content and phylogeny, with reference to previous observations made between mouse and rat. Unlike the mouse-rat comparisons, synteny seems to be largely conserved between the two mouse species. Disruption of local synteny is generally associated with differences in repeat content, although these differences appear to arise more from deletion than new integrations. Even though unambiguous V1R orthology is evident, we observe dynamic modulation of the functional repertoires, with two of seven V1Rb and one of eleven V1Ra genes lost in spretus, two V1Ra genes becoming pseudogenes in musculus, two additional orthologous pairs apparently subject to strong adaptive selection, and another divergent orthologous pair that apparently was subjected to gene conversion.ConclusionTherefore, eight of the 18 (~44%) presumptive V1Ra/V1Rb genes in the musculus-spretus ancestor appear to have undergone functional modulation since these two species diverged. As compared to the rat-mouse split, where modulation is evident by independent expansions of these two V1R subfamilies, divergence between musculus and spretus has arisen more by mutations within coding sequences. These results support the hypothesis that adaptive changes in functional V1R repertoires contribute to the delineation of very closely related species.

Phylogenetic and functional assessment of orthologs inference projects and methods.

Accurate genome-wide identification of orthologs is a central problem in comparative genomics, a fact reflected by the numerous orthology identification projects developed in recent years. However, only a few reports have compared their accuracy, and indeed, several recent efforts have not yet been systematically evaluated. Furthermore, orthology is typically only assessed in terms of function conservation, despite the phylogeny-based original definition of Fitch. We collected and mapped the results of nine leading orthology projects and methods (COG, KOG, Inparanoid, OrthoMCL, Ensembl Compara, Homologene, RoundUp, EggNOG, and OMA) and two standard methods (bidirectional best-hit and reciprocal smallest distance). We systematically compared their predictions with respect to both phylogeny and function, using six different tests. This required the mapping of millions of sequences, the handling of hundreds of millions of predicted pairs of orthologs, and the computation of tens of thousands of trees. In phylogenetic analysis or in functional analysis where high specificity is required, we find that OMA and Homologene perform best. At lower functional specificity but higher coverage level, OrthoMCL outperforms Ensembl Compara, and to a lesser extent Inparanoid. Lastly, the large coverage of the recent EggNOG can be of interest to build broad functional grouping, but the method is not specific enough for phylogenetic or detailed function analyses. In terms of general methodology, we observe that the more sophisticated tree reconstruction/reconciliation approach of Ensembl Compara was at times outperformed by pairwise comparison approaches, even in phylogenetic tests. Furthermore, we show that standard bidirectional best-hit often outperforms projects with more complex algorithms. First, the present study provides guidance for the broad community of orthology data users as to which database best suits their needs. Second, it introduces new methodology to verify orthology. And third, it sets performance standards for current and future approaches.

The G protein-coupled receptor subset of the dog genome is more similar to that in humans than rodents

BackgroundThe dog is an important model organism and it is considered to be closer to humans than rodents regarding metabolism and responses to drugs. The close relationship between humans and dogs over many centuries has lead to the diversity of the canine species, important genetic discoveries and an appreciation of the effects of old age in another species. The superfamily of G protein-coupled receptors (GPCRs) is one of the largest gene families in most mammals and the most exploited in terms of drug discovery. An accurate comparison of the GPCR repertoires in dog and human is valuable for the prediction of functional similarities and differences between the species.ResultsWe searched the dog genome for non-olfactory GPCRs and obtained 353 full-length GPCR gene sequences, 18 incomplete sequences and 13 pseudogenes. We established relationships between human, dog, rat and mouse GPCRs resolving orthologous pairs and species-specific duplicates. We found that 12 dog GPCR genes are missing in humans while 24 human GPCR genes are not part of the dog GPCR repertoire. There is a higher number of orthologous pairs between dog and human that are conserved as compared with either mouse or rat. In almost all cases the differences observed between the dog and human genomes coincide with other variations in the rodent species. Several GPCR gene expansions characteristic for rodents are not found in dog.ConclusionThe repertoire of dog non-olfactory GPCRs is more similar to the repertoire in humans as compared with the one in rodents. The comparison of the dog, human and rodent repertoires revealed several examples of species-specific gene duplications and deletions. This information is useful in the selection of model organisms for pharmacological experiments.

Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx.

BackgroundNew methods are needed for genomic-scale analysis of emerging model organisms that exemplify important biological questions but lack fully sequenced genomes. For example, there is an urgent need to understand the potential for corals to adapt to climate change, but few molecular resources are available for studying these processes in reef-building corals. To facilitate genomics studies in corals and other non-model systems, we describe methods for transcriptome sequencing using 454, as well as strategies for assembling a useful catalog of genes from the output. We have applied these methods to sequence the transcriptome of planulae larvae from the coral Acropora millepora.ResultsMore than 600,000 reads produced in a single 454 sequencing run were assembled into ~40,000 contigs with five-fold average sequencing coverage. Based on sequence similarity with known proteins, these analyses identified ~11,000 different genes expressed in a range of conditions including thermal stress and settlement induction. Assembled sequences were annotated with gene names, conserved domains, and Gene Ontology terms. Targeted searches using these annotations identified the majority of genes associated with essential metabolic pathways and conserved signaling pathways, as well as novel candidate genes for stress-related processes. Comparisons with the genome of the anemone Nematostella vectensis revealed ~8,500 pairs of orthologs and ~100 candidate coral-specific genes. More than 30,000 SNPs were detected in the coral sequences, and a subset of these validated by re-sequencing.ConclusionThe methods described here for deep sequencing of the transcriptome should be widely applicable to generate catalogs of genes and genetic markers in emerging model organisms. Our data provide the most comprehensive sequence resource currently available for reef-building corals, and include an extensive collection of potential genetic markers for association and population connectivity studies. The characterization of the larval transcriptome for this widely-studied coral will enable research into the biological processes underlying stress responses in corals and evolutionary adaptation to global climate change.

Ortholog-based protein-protein interaction prediction and its application to inter-species interactions

BackgroundThe rapid growth of protein-protein interaction (PPI) data has led to the emergence of PPI network analysis. Despite advances in high-throughput techniques, the interactomes of several model organisms are still far from complete. Therefore, it is desirable to expand these interactomes with ortholog-based and other methods.ResultsOrthologous pairs of 18 eukaryotic species were expanded and merged with experimental PPI datasets. The contributions of interologs from each species were evaluated. The expanded orthologous pairs enable the inference of interologs for various species. For example, more than 32,000 human interactions can be predicted. The same dataset has also been applied to the prediction of host-pathogen interactions. PPIs between P. falciparum calmodulin and several H. sapiens proteins are predicted, and these interactions may contribute to the maintenance of host cell Ca2+ concentration. Using comparisons with Bayesian and structure-based approaches, interactions between putative HSP40 homologs of P. falciparum and the H. sapiens TNF receptor associated factor family are revealed, suggesting a role for these interactions in the interference of the human immune response to P. falciparum.ConclusionThe PPI datasets are available from POINT and POINeT . Further development of methods to predict host-pathogen interactions should incorporate multiple approaches in order to improve sensitivity, and should facilitate the identification of targets for drug discovery and design.

Eukaryotic-like protein kinases in the prokaryotes and the myxobacterial kinome

Ser/Thr/Tyr kinases, which together comprise a major class of regulatory proteins in eukaryotes, were not believed to play an important role in prokaryotes until recently. However, our analysis of 626 prokaryotic genomes reveals that eukaryotic-like protein kinases (ELKs) are found in nearly two-thirds of the sequenced strains. We have identified 2697 ELKs, most of which are encoded by multicellular strains of the phyla Proteobacteria (Myxococcales), Actinobacteria, Cyanobacteria, and Chloroflexi, and 2 Acidobacteria and 1 Planctomycetes. Astonishingly, 7 myxobacterial strains together encode 892 ELKs, with 4 of the strains exhibiting a genomic ELK density similar to that observed in eukaryotes. Most myxobacterial ELKs show a modular organization in which the kinase domain is located at the N terminus. The C-terminal portion of the ELKs is highly diverse and often contains sequences with similarity to characterized domains, most of them involved in signaling mechanisms or in protein-protein interactions. However, many of these architectures are unique to the myxobacteria, an observation that suggests that this group exploits sophisticated and novel signal transduction systems. Phylogenetic reconstruction using the kinase domains revealed many orthologous sequence pairs and a huge number of gene duplications that probably occurred after speciation. Furthermore, studies of the microsynteny in the ELK-encoding regions reveal only low levels of synteny among Myxococcus xanthus, Plesiocystis pacifica, and Sorangium cellulosum. However, extensive similarities between M. xanthus, Stigmatella aurantiaca, and 3 Anaeromyxobacter strains were observed, indicating that they share regulatory pathways involving various ELKs.

Genomics of industrial Aspergilli and comparison with toxigenic relatives

Aspergillus oryzae has been used in Japanese fermentation industries for more than a thousand years. The species produces large amounts of various hydrolytic enzymes and has been successfully applied to modern biotechnology. The size of the A. oryzae genome (37.5 Mb) is very close to that of A. flavus and A. niger, and 20–30% larger than that of either A. nidulans or A. fumigatus. A. oryzae and A. flavus have exactly the same number of aspartic proteinase genes, of which each orthologous pair shares highly conserved amino acid sequences. Synteny analysis with A. fumigatus and A. nidulans showed that the A. oryzae genome has a mosaic structure consisting of syntenic and non-syntenic blocks. In the microorganisms to be compared, the density of the genes having homologs was obviously higher on the syntenic than on the non-syntenic blocks. Expression analysis by the DNA microarray supported the significantly lower expression of genes on the non-syntenic than on the syntenic blocks.

Conservation of peptide acceptor preferences between Drosophila and mammalian polypeptide-GalNAc transferase ortholog pairs

UDP-GalNAc:polypeptide alpha-N-acetylgalactosaminyltrans- ferases (ppGalNAc Ts) comprise a large family of glycosyltransferases that initiate mucin-type protein O-glycosylation, transferring alpha-GalNAc to Thr and Ser residues of polypeptide acceptors. Families of ppGalNAc Ts are found across diverse eukaryotes with orthologs identifiable from mammals to single-cell organisms. The peptide substrate specificity and specific protein targets of the individual ppGalNAc T family members remain poorly understood. Previously, we reported a series of oriented random peptide substrate libraries for quantitatively determining the peptide substrate specificities of the mammalian ppGalNAc T1 and T2 (Gerken TA, Raman J, Fritz TA, Jamison O. 2006. Identification of common and unique peptide substrate preferences for the UDP-GalNAc:polypeptide alpha-N-acetylgalactosaminyltransferases T1 & T2 (ppGalNAc T1 & T2) derived from oriented random peptide substrates. J Biol Chem. 281:32403-32416). With these substrates, previously unknown features of the transferases were revealed. Utilizing these and a new lengthened set of random peptides, studies have now been performed on PGANT5 and PGANT2, the Drosophila orthologs of T1 and T2. The results from these studies suggest that the major peptide substrate determinants for these transferases are contained within 2 to 3 residues flanking the site of glycosylation. It is further found that the mammalian and fly T1 orthologs display very similar peptide substrate preferences, while the T2 orthologs are nearly indistinguishable, suggesting similar peptide preferences amongst orthologous pairs have been maintained across evolution. This conclusion is further supported by sequence homology comparisons of each of the transferase orthologs, showing that the peptide substrate and UDP binding site residues are more highly conserved between species relative to their remaining catalytic and lectin domain residues.

Preferential subfunctionalization of slow-evolving genes after allopolyploidization in Xenopus laevis

As paleopolyploid genomes evolve, the expression profiles of retained gene pairs are expected to diverge. To examine this divergence process on a large scale in a vertebrate system, we compare Xenopus laevis, which has retained approximately 40% of loci in duplicate after a recent whole-genome duplication (WGD), with its unduplicated relative Silurana (Xenopus) tropicalis. This comparison of ingroup pairs to an outgroup allows the direction of change in expression profiles to be inferred for a set of 1,300 X. laevis pairs, relative to their single orthologs in S. tropicalis, across 11 tissues. We identify 68 pairs in which X. laevis is inferred to have undergone a significant reduction of expression in at least two tissues since WGD. Of these pairs, one-third show evidence of subfunctionalization, with decreases in the expression levels of different gene copies in two different tissues. Surprisingly, we find that genes with slow rates of evolution are particularly prone to subfunctionalization, even when the tendency for highly expressed genes to evolve slowly is controlled for. We interpret this result to be an effect of allopolyploidization. We then compare the outcomes of this WGD with an independent one that happened in the teleost fish lineage. We find that if a gene pair was retained in duplicate in X. laevis, the orthologous pair is more likely to have been retained in duplicate in zebrafish, suggesting that similar factors, among them subfunctionalization, determined which gene pairs survived in duplicate after the two WGDs.

CpG dinucleotides and the mutation rate of non-CpG DNA

The neutral mutation rate is equal to the base substitution rate when the latter is not affected by natural selection. Differences between these rates may reveal that factors such as natural selection, linkage, or a mutator locus are affecting a given sequence. We examined the neutral base substitution rate by measuring the sequence divergence of approximately 30,000 pairs of inactive orthologous L1 retrotransposon sequences interspersed throughout the human and chimpanzee genomes. In contrast to other studies, we related ortholog divergence to the time (age) that the L1 sequences resided in the genome prior to the chimpanzee and human speciation. As expected, the younger orthologs contained more hypermutable CpGs than the older ones because of their conversion to TpGs (and CpAs). Consequently, the younger orthologs accumulated more CpG mutations than the older ones during the approximately 5 million years since the human and chimpanzee lineages separated. But during this same time, the younger orthologs also accumulated more non-CpG mutations than the older ones. In fact, non-CpG and CpG mutations showed an almost perfect (R2 = 0.98) correlation for approximately 97% of the ortholog pairs. The correlation is independent of G + C content, recombination rate, and chromosomal location. Therefore, it likely reflects an intrinsic effect of CpGs, or mutations thereof, on non-CpG DNA rather than the joint manifestation of the chromosomal environment. The CpG effect is not uniform for all regions of non-CpG DNA. Therefore, the mutation rate of non-CpG DNA is contingent to varying extents on local CpG content. Aside from their implications for mutational mechanisms, these results indicate that a precise determination of a uniform genome-wide neutral mutation rate may not be attainable.

Measuring Global Credibility with Application to Local Sequence Alignment

Computational biology is replete with high-dimensional (high-D) discrete prediction and inference problems, including sequence alignment, RNA structure prediction, phylogenetic inference, motif finding, prediction of pathways, and model selection problems in statistical genetics. Even though prediction and inference in these settings are uncertain, little attention has been focused on the development of global measures of uncertainty. Regardless of the procedure employed to produce a prediction, when a procedure delivers a single answer, that answer is a point estimate selected from the solution ensemble, the set of all possible solutions. For high-D discrete space, these ensembles are immense, and thus there is considerable uncertainty. We recommend the use of Bayesian credibility limits to describe this uncertainty, where a (1−α)%, 0≤α≤1, credibility limit is the minimum Hamming distance radius of a hyper-sphere containing (1−α)% of the posterior distribution. Because sequence alignment is arguably the most extensively used procedure in computational biology, we employ it here to make these general concepts more concrete. The maximum similarity estimator (i.e., the alignment that maximizes the likelihood) and the centroid estimator (i.e., the alignment that minimizes the mean Hamming distance from the posterior weighted ensemble of alignments) are used to demonstrate the application of Bayesian credibility limits to alignment estimators. Application of Bayesian credibility limits to the alignment of 20 human/rodent orthologous sequence pairs and 125 orthologous sequence pairs from six Shewanella species shows that credibility limits of the alignments of promoter sequences of these species vary widely, and that centroid alignments dependably have tighter credibility limits than traditional maximum similarity alignments.

Quantitative evidence for conserved longevity pathways between divergent eukaryotic species

Studies in invertebrate model organisms have been a driving force in aging research, leading to the identification of many genes that influence life span. Few of these genes have been examined in the context of mammalian aging, however, and it remains an open question as to whether and to what extent the pathways that modulate longevity are conserved across different eukaryotic species. Using a comparative functional genomics approach, we have performed the first quantitative analysis of the degree to which longevity genes are conserved between two highly divergent eukaryotic species, the yeast Saccharomyces cerevisiae and the nematode Caenorhabditis elegans. Here, we report the replicative life span phenotypes for single-gene deletions of the yeast orthologs of worm aging genes. We find that 15% of these yeast deletions are long-lived. In contrast, only 3.4% of a random set of deletion mutants are long-lived-a statistically significant difference. These data suggest that genes that modulate aging have been conserved not only in sequence, but also in function, over a billion years of evolution. Among the longevity determining ortholog pairs, we note a substantial enrichment for genes involved in an evolutionarily conserved pathway linking nutrient sensing and protein translation. In addition, we have identified several conserved aging genes that may represent novel longevity pathways. Together, these findings indicate that the genetic component of life span determination is significantly conserved between divergent eukaryotic species, and suggest pathways that are likely to play a similar role in mammalian aging.

Investigation of Loss and Gain of Introns in the Compact Genomes of Pufferfishes (Fugu and Tetraodon)

We have investigated intron evolution in the compact genomes of 2 closely related species of pufferfishes, Fugu rubripes and Tetraodon nigroviridis, that diverged about 32 million years ago (MYA). Analysis of 148,028 aligned intron positions in 13,547 gene pairs using human as an outgroup identified 57 and 24 intron losses in Tetraodon and fugu lineages, respectively, and no gain in either lineage. For comparison, we analyzed 144,545 intron positions in 12,866 orthologous pairs of genes in human and mouse that diverged about 61 MYA using dog as an outgroup and identified 51 intron losses in mouse and 3 losses in human and no gain. The rate of intron loss in Tetraodon is higher than that in fugu, mouse, and human but lower than the previous estimates for other eukaryotes. The introns lost in pufferfishes and mammals are significantly shorter than the mean size of introns in the genome. One intron deleted in fugu and another in Tetraodon have left behind 6 and 3 nucleotides, respectively, suggesting that they were lost due to genomic deletions. Such losses of introns are likely to be the result of a higher rate of DNA deletions experienced by the genomes of pufferfishes compared with mammals. The shorter generation time of Tetraodon compared with fugu, and the rich diversity and higher activity of transposable elements in pufferfishes compared with mammals, may be responsible for the higher rate of intron loss in Tetraodon. Our findings indicate that overall very little intron turnover has occurred in pufferfishes and mammals during recent evolution and that intron gain is an extremely rare event in vertebrate evolution.