Evolutionary Sequence Analysis Research Articles

WebPRANK: a phylogeny-aware multiple sequence aligner with interactive alignment browser

BackgroundPhylogeny-aware progressive alignment has been found to perform well in phylogenetic alignment benchmarks and to produce superior alignments for the inference of selection on codon sequences. Its implementation in the PRANK alignment program package also allows modelling of complex evolutionary processes and inference of posterior probabilities for sequence sites evolving under each distinct scenario, either simultaneously with the alignment of sequences or as a post-processing step for an existing alignment. This has led to software with many advanced features, and users may find it difficult to generate optimal alignments, visualise the full information in their alignment results, or post-process these results, e.g. by objectively selecting subsets of alignment sites.ResultsWe have created a web server called webPRANK that provides an easy-to-use interface to the PRANK phylogeny-aware alignment algorithm. The webPRANK server supports the alignment of DNA, protein and codon sequences as well as protein-translated alignment of cDNAs, and includes built-in structure models for the alignment of genomic sequences. The resulting alignments can be exported in various formats widely used in evolutionary sequence analyses. The webPRANK server also includes a powerful web-based alignment browser for the visualisation and post-processing of the results in the context of a cladogram relating the sequences, allowing (e.g.) removal of alignment columns with low posterior reliability. In addition to de novo alignments, webPRANK can be used for the inference of ancestral sequences with phylogenetically realistic gap patterns, and for the annotation and post-processing of existing alignments. The webPRANK server is freely available on the web at http://tinyurl.com/webprank .ConclusionsThe webPRANK server incorporates phylogeny-aware multiple sequence alignment, visualisation and post-processing in an easy-to-use web interface. It widens the user base of phylogeny-aware multiple sequence alignment and allows the performance of all alignment-related activity for small sequence analysis projects using only a standard web browser.

Association of hepatitis B virus mutations in basal core promoter and precore regions with severity of liver disease: an investigation of 793 Chinese patients with mild and severe chronic hepatitis B and acute-on-chronic liver failure

ObjectiveTo investigate the features of hepatitis B virus (HBV) basal core promoter/precore (BCP/PC) mutations and genotypes in a large number of mild/severe chronic hepatitis B (CHB-M/CHB-S), and acute-on-chronic liver failure (ACLF) patients and analyze the clinical implications of the virologic features.Patients and methodsSera of 793 (325 CHB-M, 170 CHB-S, and 298 ACLF) patients admitted to or who had visited Beijing 302 Hospital from January 2005 to December 2008 were collected and successfully amplified for the HBV BCP/PC and a 1225-bp-long S/Pol (nt 54–1278) gene regions. Biochemical and serological parameters and HBV DNA level were routinely performed. Viral DNA was extracted and subjected to a nested PCR. Genotypes/subgenotypes were determined based on complete genomic sequence or on analysis of the 1225-bp-long S/Pol-gene sequence. HBV genotyping was performed by direct PCR sequencing followed by molecular evolutionary analysis of the viral sequences. A P value of <0.05 (two-sided) was considered to be statistically significant.ConclusionsOur findings suggest that CHB patients infected with BCP/PC mutant viruses are more susceptible to severe hepatitis and ACLF than those with the BCP/PC wild-type virus and that ACLF patients with PC mutant viruses have an increased risk of death. As such, the HBV PC mutation is a potential predictive indicator of ACLF outcome.

Evolutionary analysis of regulatory sequences (EARS) in plants

Identification of regulatory sequences within non-coding regions of DNA is an essential step towards elucidation of gene networks. This approach constitutes a major challenge, however, as only a very small fraction of non-coding DNA is thought to contribute to gene regulation. The mapping of regulatory regions traditionally involves the laborious construction of promoter deletion series which are then fused to reporter genes and assayed in transgenic organisms. Bioinformatic methods can be used to scan sequences for matches for known regulatory motifs, however these methods are currently hampered by the relatively small amount of such motifs and by a high false-discovery rate. Here, we demonstrate a robust and highly sensitive, in silico method to identify evolutionarily conserved regions within non-coding DNA. Sequence conservation within these regions is taken as evidence for evolutionary pressure against mutations, which is suggestive of functional importance. We test this method on a small set of well characterised promoters, and show that it successfully identifies known regulatory regions. We further show that these evolutionarily conserved sequences contain clusters of transcription binding sites, often described as regulatory modules. A version of the tool optimised for the analysis of plant promoters is available online at http://wsbc.warwick.ac.uk/ears/main.php.

Analysis of hepatitis C virus strains circulating in Republic of the Congo

The aim of this study was to assess the seroprevalence, viremia, genotype distribution, and demographic history of hepatitis C virus (HCV) in the Republic of the Congo. Testing was carried out on sera samples collected in 2005 from 807 Bantus belonging to the Kongo, Teke, and Ngala subgroups and 80 Pygmies. Positive HCV serology was found in 50 (5.6%) individuals including 31 (60%) who were viremic. Seroprevalence increased with age with a cutoff at 50 years: 2.8% <50 versus 12% >50. Twenty-one strains belonged to four described subtypes, that is, 4c in eight cases, 4h in two, 4k in three, and 4r in eight. Ten strains could not be assigned to any known subtype and may represent six new variants, that is, subtype 4 in five cases and subtype 2 in one. Evolutionary analysis of subtype 4c and 4r sequences indicated a period of enhanced transmission in the mid-twentieth century probably due to iatrogenic causes. This study underlines the high genetic diversity of strains in the Republic of the Congo with nine subtypes 4 and one subtype 2.

Evolutionary analysis of sequences expressed in tumors

Earlier we have identified a new class of human genomic sequences expressed preferentially in tumors. Here we use a comparative genomics approach and conservation analysis to study evolutionary specificity of nine human tumor-specific sequences, described previously. Three sequences had originated in the primate lineage. The other three had mammalian orthologs, but conservation analysis has shown that these sequences evolved neutrally. Three sequences were conservative. These data confirm previously formulated hypothesis that evolutionarily new genes are expressed in tumors.

Statistical Comparison of Nucleotide, Amino Acid, and Codon Substitution Models for Evolutionary Analysis of Protein-Coding Sequences

Statistical models for the evolution of molecular sequences play an important role in the study of evolutionary processes. For the evolutionary analysis of protein-coding sequences, 3 types of evolutionary models are available: 1) nucleotide, 2) amino acid, and 3) codon substitution models. Selecting appropriate models can greatly improve the estimation of phylogenies and divergence times and the detection of positive selection. Although much attention has been paid to the comparisons among the same types of models, relatively little attention has been paid to the comparisons among the different types of models. Additionally, because such models have different data structures, comparison of those models using conventional model selection criteria such as Akaike information criterion (AIC) or Bayesian information criterion (BIC) is not straightforward. Here, we suggest new procedures to convert models of the above-mentioned 3 types to 64-dimensional models with nucleotide triplet substitution. These conversion procedures render it possible to statistically compare the models of these 3 types by using AIC or BIC. By analyzing divergent and conserved interspecific mammalian sequences and intraspecific human population data, we show the superiority of the codon substitution models and discuss the advantages and disadvantages of the models of the 3 types.

Evolution of multiple phosphodiesterase isoforms in stickleback involved in cAMP signal transduction pathway

BackgroundDuplicate genes are considered to have evolved through the partitioning of ancestral functions among duplicates (subfunctionalization) and/or the acquisition of novel functions from a beneficial mutation (neofunctionalization). Additionally, an increase in gene dosage resulting from duplication may also confer an advantageous effect, as has been suggested for histone, tRNA, and rRNA genes. Currently, there is little understanding of the effect of increased gene dosage on subcellular networks like signal transduction pathways. Addressing this issue may provide further insights into the evolution by gene duplication.ResultsWe analyzed the evolution of multiple stickleback phosphodiesterase (PDE, EC: 3.1.4.17) 1C genes involved in the cyclic nucleotide signaling pathway. Stickleback has 8–9 copies of this gene, whereas only one or two loci exist in other model vertebrates. Our phylogenetic and synteny analyses suggested that the multiple PDE1C genes in stickleback were generated by repeated duplications of >100-kbp chromosome segments. Sequence evolution analysis did not provide strong evidence for neofunctionalization in the coding sequences of stickleback PDE1C isoforms. On the other hand, gene expression analysis suggested that the derived isoforms acquired expression in new organs, implying their neofunctionalization in terms of expression patterns. In addition, at least seven isoforms of the stickleback PDE1C were co-expressed with olfactory-type G-proteins in the nose, suggesting that PDE1C dosage is increased in the stickleback olfactory transduction (OT) pathway. In silico simulations of OT implied that the increased PDE1C dosage extends the longevity of the depolarization signals of the olfactory receptor neuron.ConclusionThe predicted effect of the increase in PDE1C products on the OT pathway may play an important role in stickleback behavior and ecology. However, this possibility should be empirically examined. Our analyses imply that an increase in gene product sometimes has a significant, yet unexpected, effect on the functions of subcellular networks.

Phylogenetic Analysis and Molecular Evolution of Guanine Deaminases: From Guanine to Dendrites

Guanine deaminase (GDA; guanase) is a ubiquitous enzyme that catalyzes the first step of purine metabolism by hydrolytic deamination of guanine, resulting in the production of xanthine. This hydrolase subfamily member plays an essential role in maintaining homeostasis of cellular triphosphate nucleotides for energy, signal transduction pathways, and nitrogen sources. In mammals, GDA protein levels can play a role in neuronal development by regulating dendritic arborization. We previously demonstrated that the most abundant alternative splice form of GDA in mammals, termed cypin (cytosolic PSD-95 interactor), interacts with postsynaptic density proteins, regulates microtubule polymerization, and increases dendrite number. Since purine metabolism and dendrite development were previously thought to be independent cellular processes, this multifunctional protein serves as a new target for the treatment of cognitive disorders characterized by aberrant neuronal morphology and purine metabolism. Although the enzymatic activity of GDA has been conserved during evolution from prokaryotes to higher eukaryotes, a detailed evolutionary assessment of the principal domains in GDA proteins has not yet been put forward. In this study, we perform a complete evolutionary analysis of the full-length sequences and the principal domains in guanine deaminases. Furthermore, we reconstruct the molecular phylogeny of guanine deaminases with neighbor-joining, maximum-likelihood, and UPGMA methods of phylogenetic inference. This study can act as a model whereby a universal housekeeping enzyme may be adapted to act also as a key regulator of a developmental process.

Molecular characterization and expression of a novel homolog of uncoupling protein 5 (UCP5) from the eastern oyster Crassostrea virginica (Bivalvia: Ostreidae)

Uncoupling proteins (UCPs) belong to the mitochondrial anion carrier gene family which has been implicated in diverse physiological functions ranging from thermoregulation to antioxidant defense. In mammals, the UCP family is well characterized and contains five members (UCP1-5). In contrast, invertebrate homologues of uncoupling proteins are much less studied both from the viewpoints of structure and function. In this study we report nucleotide and predicted protein structure of an important member of UCP family, UCP5 from eastern oysters Crassostrea virginica. UCP5 from oysters appears to be a close homolog of the mammalian brain mitochondrial carrier protein (BMCP1, or UCP5) and is the first full-length UCP described from a Lophotrochozoan invertebrate. Evolutionary analysis of UCP sequences indicates at least three monophyletic UCP branches (UCP1-3, UCP4 and UCP5) that have diverged early in the evolution, prior to the divergence of vertebrates and invertebrates. In oysters, two forms of UCP5 transcript are found (UCP5S and UCP5L) that differ by 152 bp in length due to the presence of an intron in UCP5L. UCP5 was expressed in all studied oyster tissues, unlike mammals, where UCP5 is predominantly expressed in brains and male gonads. Hypoxia-reoxygenation stress, sublethal Cd exposure (50 μg L − 1 Cd for 56 days) and acclimation to different temperatures (12 and 20 °C) had no significant effect on UCP5 mRNA expression in oysters indicative of its relative unimportance in antioxidant defense and temperature adaptation of oyster mitochondria. These data suggest that despite the relatively high degree of evolutionary conservation of the UCP5 amino acid sequence, its functional significance in mitochondria changed in the course of evolution of mollusks and vertebrates.

Time-dependent ARMA modeling of genomic sequences

BackgroundOver the past decade, many investigators have used sophisticated time series tools for the analysis of genomic sequences. Specifically, the correlation of the nucleotide chain has been studied by examining the properties of the power spectrum. The main limitation of the power spectrum is that it is restricted to stationary time series. However, it has been observed over the past decade that genomic sequences exhibit non-stationary statistical behavior. Standard statistical tests have been used to verify that the genomic sequences are indeed not stationary. More recent analysis of genomic data has relied on time-varying power spectral methods to capture the statistical characteristics of genomic sequences. Techniques such as the evolutionary spectrum and evolutionary periodogram have been successful in extracting the time-varying correlation structure. The main difficulty in using time-varying spectral methods is that they are extremely unstable. Large deviations in the correlation structure results from very minor perturbations in the genomic data and experimental procedure. A fundamental new approach is needed in order to provide a stable platform for the non-stationary statistical analysis of genomic sequences.ResultsIn this paper, we propose to model non-stationary genomic sequences by a time-dependent autoregressive moving average (TD-ARMA) process. The model is based on a classical ARMA process whose coefficients are allowed to vary with time. A series expansion of the time-varying coefficients is used to form a generalized Yule-Walker-type system of equations. A recursive least-squares algorithm is subsequently used to estimate the time-dependent coefficients of the model. The non-stationary parameters estimated are used as a basis for statistical inference and biophysical interpretation of genomic data. In particular, we rely on the TD-ARMA model of genomic sequences to investigate the statistical properties and differentiate between coding and non-coding regions in the nucleotide chain. Specifically, we define a quantitative measure of randomness to assess how far a process deviates from white noise. Our simulation results on various gene sequences show that both the coding and non-coding regions are non-random. However, coding sequences are "whiter" than non-coding sequences as attested by a higher index of randomness.ConclusionWe demonstrate that the proposed TD-ARMA model can be used to provide a stable time series tool for the analysis of non-stationary genomic sequences. The estimated time-varying coefficients are used to define an index of randomness, in order to assess the statistical correlations in coding and non-coding DNA sequences. It turns out that the statistical differences between coding and non-coding sequences are more subtle than previously thought using stationary analysis tools: Both coding and non-coding sequences exhibit statistical correlations, with the coding regions being "whiter" than the non-coding regions. These results corroborate the evolutionary periodogram analysis of genomic sequences and revoke the stationary analysis' conclusion that coding DNA behaves like random sequences.

High prevalence of hepatitis C virus infection and predominance of genotype 4 in rural gabon

Hepatitis C (HCV) molecular epidemiology is documented poorly in central African countries. In response to this, a population-based study of 319 consenting adults resident in a remote village of Gabon was undertaken (mean age: 38 years; age range: 13-85+; sex ratio: 0.74). Screening for anti-HCV antibodies was performed using ELISA and recombinant immunoblot assay. Seropositive samples were assessed further with viral load and genotyping techniques. Sixty-six (20.7%) individuals were HCV seropositive. Viral loads ranged from 600 to 24.9 million IU/ml (median: 372,500). Seroprevalence and viral loads increased significantly with age (P < 10(-5) and P < 0.003, respectively). HCV sequences of the 5'UTR genome region were obtained from 60 (90.9%) samples and NS5B region sequences were obtained from 22 (36.6%) samples. All strains belonged to subtypes of genotype 4: 4e (72.7%), 4c (13.6%), 4p (4.5%), 4r (4.5%) and one unclassified genotype 4 strain. Evolutionary analysis of the subtype 4e sequences indicates a period of raised transmission during the early twentieth century.

Hadamard Conjugation for the Kimura 3ST Model: Combinatorial Proof Using Path Sets

Under a stochastic model of molecular sequence evolution the probability of each possible pattern of a characters is well defined. The Kimura's three-substitution-types (K3ST) model of evolution, allows analytical expression for these probabilities of by means of the Hadamard conjugation as a function of the phylogeny T and the substitution probabilities on each edge of TM . In this paper we produce a direct combinatorial proof of these results, using pathset distances which generalise pairwise distances between sequences. This interpretation provides us with tools that were proved useful in related problems in the mathematical analysis of sequence evolution.

Mitochondrial p32 Is a Critical Mediator of ARF-Induced Apoptosis

The shared exon 2 of the p14ARF-p16INK4a locus is frequently mutated in human cancers. However, in contrast to the exon 1beta-encoded N-terminal half of ARF, the function of the exon 2-encoded C-terminal half of ARF has been elusive. Here, we report that the mitochondrial protein p32/C1QBP binds the ARF C terminus. We show that p32 is required for ARF to localize to mitochondria and induce apoptosis, and that ARF mutations specifically disrupting p32 binding can impair both of these functions. Wild-type ARF, but not a p32-binding-deficient ARF mutant, localizes to mitochondria, reduces mitochondrial membrane potential, and sensitizes cells to p53-induced apoptosis. These findings provide a potential explanation for the frequent human cancer mutations targeting the ARF C terminus.

Synonymous Substitutions Substantially Improve Evolutionary Inference from Highly Diverged Proteins

Codon-and amino acid-substitution models are widely used for the evolutionary analysis of protein-coding DNA sequences. Using codon models, the amounts of both nonsynonymous and synonymous DNA substitutions can be estimated. The ratio of these amounts represents the strength of selective pressure. Using amino acid models, the amount of nonsynonymous substitutions is estimated, but that of synonymous substitutions is ignored. Although amino acid models lose any information regarding synonymous substitutions, they explicitly incorporate the information for amino acid replacement, which is empirically derived from databases. It is often presumed that when the protein-coding sequences are highly divergent, synonymous substitutions might be saturated and the evolutionary analysis may be hampered by synonymous noise. However, there exists no quantitative procedure to verify whether synonymous substitutions can be ignored; therefore, amino acid models have been arbitrarily selected. In this study, we investigate the issue of a statistical comparison between codon-and amino acid-substitution models. For this purpose, we propose a new procedure to transform a 20-dimensional amino acid model to a 61-dimensional codon model. This transformation reveals that amino acid models belong to a subset of the codon models and enables us to test whether synonymous substitutions can be ignored by using the likelihood ratio. Our theoretical results and analyses of real data indicate that synonymous substitutions are very informative and substantially improve evolutionary inference, even when the sequences are highly divergent. Therefore, we note that amino acid models should be adopted only after carefully investigating and discarding the possibility that synonymous substitutions can reveal important evolutionary information.

Evolution of the redox function in mammalian apurinic/apyrimidinic endonuclease

Human apurinic/apyrimidinic endonuclease (hApe1) encodes two important functional activities: an essential base excision repair (BER) activity and a redox activity that regulates expression of a number of genes through reduction of their transcription factors, AP-1, NFκB, HIF-1α, CREB, p53 and others. The BER function is highly conserved from prokaryotes ( E. coli exonuclease III) to humans (hApe1). Here, we provide evidence supporting a redox function unique to mammalian Apes. An evolutionary analysis of Ape sequences reveals that, of the 7 Cys residues, Cys 93, 99, 208, 296, and 310 are conserved in both mammalian and non-mammalian vertebrate Apes, while Cys 65 is unique to mammalian Apes. In the zebrafish Ape (zApe), selected as the vertebrate sequence most distant from human, the residue equivalent to Cys 65 is Thr 58. The wild-type zApe enzyme was tested for redox activity in both in vitro EMSA and transactivation assays and found to be inactive, similar to C65A hApe1. Substitution of Thr 58 with Cys in zApe, however, resulted in a redox active enzyme, suggesting that a Cys residue in this position is indeed critical for redox function. In order to further probe differences between redox active and inactive enzymes, we have determined the crystal structures of vertebrate redox inactive enzymes, the C65A human Ape1 enzyme and the zApe enzyme at 1.9 and 2.3 Å, respectively. Our results provide new insights on the redox function and highlight a dramatic gain-of-function activity for Ape1 in mammals not found in non-mammalian vertebrates or lower organisms.

MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences

The Molecular Evolutionary Genetics Analysis (MEGA) software is a desktop application designed for comparative analysis of homologous gene sequences either from multigene families or from different species with a special emphasis on inferring evolutionary relationships and patterns of DNA and protein evolution. In addition to the tools for statistical analysis of data, MEGA provides many convenient facilities for the assembly of sequence data sets from files or web-based repositories, and it includes tools for visual presentation of the results obtained in the form of interactive phylogenetic trees and evolutionary distance matrices. Here we discuss the motivation, design principles and priorities that have shaped the development of MEGA. We also discuss how MEGA might evolve in the future to assist researchers in their growing need to analyze large data set using new computational methods.

The Average Mutual Information Profile as a Genomic Signature

BackgroundOccult organizational structures in DNA sequences may hold the key to understanding functional and evolutionary aspects of the DNA molecule. Such structures can also provide the means for identifying and discriminating organisms using genomic data. Species specific genomic signatures are useful in a variety of contexts such as evolutionary analysis, assembly and classification of genomic sequences from large uncultivated microbial communities and a rapid identification system in health hazard situations.ResultsWe have analyzed genomic sequences of eukaryotic and prokaryotic chromosomes as well as various subtypes of viruses using an information theoretic framework. We confirm the existence of a species specific average mutual information (AMI) profile. We use these profiles to define a very simple, computationally efficient, alignment free, distance measure that reflects the evolutionary relationships between genomic sequences. We use this distance measure to classify chromosomes according to species of origin, to separate and cluster subtypes of the HIV-1 virus, and classify DNA fragments to species of origin.ConclusionAMI profiles of DNA sequences prove to be species specific and easy to compute. The structure of AMI profiles are conserved, even in short subsequences of a species' genome, rendering a pervasive signature. This signature can be used to classify relatively short DNA fragments to species of origin.

Inferring horizontal transfers in the presence of rearrangements by the minimum evolution criterion

The evolution of viruses is very rapid and in addition to local point mutations (insertion, deletion, substitution) it also includes frequent recombinations, genome rearrangements and horizontal transfer of genetic materials (HGTS). Evolutionary analysis of viral sequences is therefore a complicated matter for two main reasons: First, due to HGTs and recombinations, the right model of evolution is a network and not a tree. Second, due to genome rearrangements, an alignment of the input sequences is not guaranteed. These facts encourage developing methods for inferring phylogenetic networks that do not require aligned sequences as input. In this work, we present the first computational approach which deals with both genome rearrangements and horizontal gene transfers and does not require a multiple alignment as input. We formalize a new set of computational problems which involve analyzing such complex models of evolution. We investigate their computational complexity, and devise algorithms for solving them. Moreover, we demonstrate the viability of our methods on several synthetic datasets as well as four biological datasets. The code is available from the authors upon request.

The molecular evolution of sperm zonadhesin

Based on pioneering work of Hardy and Garbers, zonadhesin has become one of the best studied sperm ligands in boreoeutherian mammals, both from a biochemical and evolutionary perspective. Zonadhesin is a mosaic-type protein that localizes to the apical head of spermatozoa. In pig, cattle, rabbit and primates, zonadhesin precursor essentially consists of two or three MAM (meprin/A5 antigen/mu receptor tyrosine phosphatase) domains, one mucin-like domain, one incomplete and four complete D domains (homologous to vWFD). Mouse zonadhesin is distinguished from this general pattern by 20 extra partial D3 domains. While concerted evolution drives the divergence of the mucin-like domain in the ortholog comparison, MAM and D domains mainly diverge under the influence of drift and positive selection, both in the paralog and ortholog comparison. As can be seen particularly well within a putative binding region in the most C-terminal MAM domain, positive selection not only causes amino acid exchanges, but also promotes changes in the pattern of predicted posttranslational modification. Moving window and correlation analyses of sequence evolution and sexual body dimorphism further suggest that sexual selection, especially sperm competition, drives zonadhesin divergence. However, considering its zona pellucida avidity, female cryptic choice might as well contribute to zonadhesin evolution. Despite the general tendency for divergence of zonadhesin, conservation by negative selection dominates the evolution of most codon sites. In accordance, the distribution of EGF (epidermal growth factor)-like motifs, DP-doublets, single cysteines and CGLC motifs suggests a wide conservation of processing, folding and oligomerization of zonadhesin in pig, rabbit and primates.

OCPAT: an online codon-preserved alignment tool for evolutionary genomic analysis of protein coding sequences.

BackgroundRapidly accumulating genome sequence data from multiple species offer powerful opportunities for the detection of DNA sequence evolution. Phylogenetic tree construction and codon-based tests for natural selection are the prevailing tools used to detect functionally important evolutionary change in protein coding sequences. These analyses often require multiple DNA sequence alignments that maintain the correct reading frame for each collection of putative orthologous sequences. Since this feature is not available in most alignment tools, codon reading frames often must be checked manually before evolutionary analyses can commence.ResultsHere we report an online codon-preserved alignment tool (OCPAT) that generates multiple sequence alignments automatically from the coding sequences of any list of human gene IDs and their putative orthologs from genomes of other vertebrate tetrapods. OCPAT is programmed to extract putative orthologous genes from genomes and to align the orthologs with the reading frame maintained in all species. OCPAT also optimizes the alignment by trimming the most variable alignment regions at the 5' and 3' ends of each gene. The resulting output of alignments is returned in several formats, which facilitates further molecular evolutionary analyses by appropriate available software. Alignments are generally robust and reliable, retaining the correct reading frame. The tool can serve as the first step for comparative genomic analyses of protein-coding gene sequences including phylogenetic tree reconstruction and detection of natural selection. We aligned 20,658 human RefSeq mRNAs using OCPAT. Most alignments are missing sequence(s) from at least one species; however, functional annotation clustering of the ~1700 transcripts that were alignable to all species shows that genes involved in multi-subunit protein complexes are highly conserved.ConclusionThe OCPAT program facilitates large-scale evolutionary and phylogenetic analyses of entire biological processes, pathways, and diseases.