PubMed Health⌕ Search

Biomedical subjects

Andrew J Roger

Publications and source records attributed to Andrew J Roger.

At least 37 records · Page 2Linked to original sources
PubMed2004

On inconsistency of the neighbor-joining, least squares, and minimum evolution estimation when substitution processes are incorrectly modeled.

Using analytical methods, we show that under a variety of model misspecifications, Neighbor-Joining, minimum evolution, and least squares estimation procedures are statistically inconsistent. Failure to correctly account for differing rates-across-sites processes, failure to correctly model rate matrix parameters, and failure to adjust for parallel rates-across-sites changes (a rates-across-subtrees process) are all shown to lead to a "long branch attraction" form of inconsistency. In addition, failure to account for rates-across-sites processes is also shown to result in underestimation of evolutionary distances for a wide variety of substitution models, generalizing an earlier analytical result for the Jukes-Cantor model reported in Golding and a similar bias result for the GTR or REV model in Kelly and Rice (1996). Although standard rates-across-sites models can be employed in many of these cases to restore consistency, current models cannot account for other kinds of misspecification. We examine an idealized but biologically relevant case, where parallel changes in rates at sites across subtrees is shown to give rise to inconsistency. This changing rates-across-subtrees type model misspecification cannot be adjusted for with conventional methods or without carefully considering the rate variation in the larger tree. The results are presented for four-taxon trees, but the expectation is that they have implications for larger trees as well. To illustrate this, a simulated 42-taxon example is given in which the microsporidia, an enigmatic group of eukaryotes, are incorrectly placed at the archaebacteria-eukaryotes split because of incorrectly specified pairwise distances. The analytical nature of the results lend insight into the reasons that long branch attraction tends to be a common form of inconsistency and reasons that other forms of inconsistency like "long branches repel" can arise in some settings. In many of the cases of inconsistency presented, a particular incorrect topology is estimated with probability converging to one, the implication being that measures of uncertainty like bootstrap support will be unable to detect that there is a problem with the estimation. The focus is on distance methods, but previous simulation results suggest that the zones of inconsistency for distance methods contain the zones of inconsistency for maximum likelihood methods as well.

Archaea↗
PubMed2004

Estimating and comparing the rates of gene discovery and expressed sequence tag (EST) frequencies in EST surveys.

MOTIVATION: Expressed sequence tag (EST) surveys are an efficient way to characterize large numbers of genes from an organism. The rate of gene discovery in an EST survey depends on the degree of redundancy of the cDNA libraries from which sequences are obtained. However, few statistical methods have been developed to assess and compare redundancies of various libraries from preliminary EST surveys. RESULTS: We consider statistics for the comparison of EST libraries based upon the frequencies with which genes occur in subsamples of reads. These measures are useful in determining which one of several libraries is more likely to yield new genes in future reads and what proportion of additional reads one might want to take from the libraries in order to be likely to obtain new genes. One approach is to compare single sample measures that have been successfully used in species estimation problems, such as coverage of a library, defined as the proportion of the library that is represented in the given sample of reads. Another single library measure is an estimate of the expected number of additional genes that will be found in a new sample of reads. We also propose statistics that jointly use data from all the libraries. Analogous formulas for coverage and the expected numbers of new genes are presented. These measures consider coverage in a single library based upon reads from all libraries and similarly, the expected numbers of new genes that will be discovered by taking reads from all libraries with fixed proportions. Together, the statistics presented provide useful comparative measures for the libraries that can be used to guide sampling from each of the libraries to maximize the rate of gene discovery. Finally, we present tests for whether genes are equally represented or expressed in a set of libraries. Binomial and chi2 tests are presented for gene-by-gene comparisons of expression. Overall tests of the equality of proportional representation are presented and multiple comparisons issues are addressed. These methods can be used to evaluate changes in gene expression reflected in the composition of EST libraries prepared from different tissue types or cells exposed to different environmental conditions. AVAILABILITY: Software will be made available at http://www.mathstat.dal.ca/~tsusko

Algorithms↗
PubMed2004

Covarion shifts cause a long-branch attraction artifact that unites microsporidia and archaebacteria in EF-1alpha phylogenies.

Microsporidia branch at the base of eukaryotic phylogenies inferred from translation elongation factor 1alpha (EF-1alpha) sequences. Because these parasitic eukaryotes are fungi (or close relatives of fungi), it is widely accepted that fast-evolving microsporidian sequences are artifactually "attracted" to the long branch leading to the archaebacterial (outgroup) sequences ("long-branch attraction," or "LBA"). However, no previous studies have explicitly determined the reason(s) why the artifactual allegiance of microsporidia and archaebacteria ("M + A") is recovered by all phylogenetic methods, including maximum likelihood, a method that is supposed to be resistant to classical LBA. Here we show that the M + A affinity can be attributed to those alignment sites associated with large differences in evolutionary site rates between the eukaryotic and archaebacterial subtrees. Therefore, failure to model the significant evolutionary rate distribution differences (covarion shifts) between the ingroup and outgroup sequences is apparently responsible for the artifactual basal position of microsporidia in phylogenetic analyses of EF-1alpha sequences. Currently, no evolutionary model that accounts for discrete changes in the site rate distribution on particular branches is available for either protein or nucleotide level phylogenetic analysis, so the same artifacts may affect many other "deep" phylogenies. Furthermore, given the relative similarity of the site rate patterns of microsporidian and archaebacterial EF-1alpha proteins ("parallel site rate variation"), we suggest that the microsporidian orthologs may have lost some eukaryotic EF-1alpha-specific nontranslational functions, exemplifying the extreme degree of reduction in this parasitic lineage.

Animals↗
PubMed2004

Protein phylogenies robustly resolve the deep-level relationships within Euglenozoa.

The deepest-level relationships amongst Euglenozoa remain poorly resolved, despite a rich history of morphological examination and numerous molecular phylogenetic studies of small subunit ribosomal RNA (SSU rRNA) data. We address this question using two nuclear-encoded proteins, the cytosolic isoforms of heat shock protein 90 (hsp90) and heat shock protein 70 (hsp70). For both proteins we examined sequences from the three primary groups within Euglenozoa (euglenids, diplonemids, and kinetoplastids), and from their close relatives, Heterolobosea. Maximum likelihood (ML) and ML distance analyses of these proteins support a close relationship between diplonemids and kinetoplastids to the exclusion of the euglenid Euglena gracilis. In hsp90 and combined protein analyses bootstrap support is very strong and alternative topologies are generally rejected by 'approximately unbiased' (AU) tests. This result is consistent with recent molecular biological and morphological data, but contradicts early structural accounts and many SSU rRNA analyses that favour a closer relationship between diplonemids and euglenids. However, a re-examination of an important SSU rRNA data set highlights the instability of the inferences from this marker. The protein analyses also suggest that bodonids are paraphyletic, with trypanosomatids grouping with 'clade 2' and 'clade 3' bodonids to the exclusion of 'clade 1' bodonids.

Animals↗
PubMed2004

Early evolution within kinetoplastids (euglenozoa), and the late emergence of trypanosomatids.

Many important relationships amongst kinetoplastids, including the position of trypanosomatids, remain uncertain, with limited taxon sampling of markers other than small subunit ribosomal RNA (SSUrRNA). We report gene sequences for cytosolic heat shock proteins 90 and/or 70 (HSP90, HSP70) from the potentially early-diverging kinetoplastids Ichthyobodo necator and Rhynchobodo sp., and from bodonid clades '2' (Parabodonidae) and '3' (Eubodonidae). Some of the new cytosolic HSP70 sequences represent a distinct paralog family (HSP70-B), which is related to yet another paralog known from trypanosomatids (HSP70-C). The (HSP70-B, HSP70-C) clade seemingly diverged before the separation between kinetoplastids and diplonemids. Protein phylogenies support the basal placement of Ichthyobodo within kinetoplastids. Unexpectedly, Rhynchobodo usually forms the next most basal group, separated from the clade '1' bodonids with which it has been allied. Bootstrap support is often weak, but the possibility that Rhynchobodo represents a separate early-diverging lineage within core kinetoplastids deserves further testing. Trypanosomatids always fall remote from the root of kinetoplastids, forming a specific relationship with bodonid clades 2 (and 3), generally with strong bootstrap support. These protein trees with improved taxon sampling provide the best evidence to date for a 'late' emergence of trypanosomatids, contradicting recent SSUrRNA-based proposals for a relatively early divergence of this group.

Amino Acid Sequence↗
PubMed2003

Evidence for Golgi bodies in proposed 'Golgi-lacking' lineages.

Golgi bodies are nearly ubiquitous in eukaryotic cells. The apparent lack of such structures in certain eukaryotic lineages might be taken to mean that these protists evolved prior to the acquisition of the Golgi, and it raises questions of how these organisms function in the absence of this crucial organelle. Here, we report gene sequences from five proposed 'Golgi-lacking' organisms (Giardia intestinalis, Spironucleus barkhanus, Entamoeba histolytica, Naegleria gruberi and Mastigamoeba balamuthi). BLAST and phylogenetic analyses show these genes to be homologous to those encoding components of the retromer, coatomer and adaptin complexes, all of which have Golgi-related functions in mammals and yeast. This is, to our knowledge, the first molecular evidence for Golgi bodies in two major eukaryotic lineages (the pelobionts and heteroloboseids). This substantiates the suggestion that there are no extant primitively 'Golgi-lacking' lineages, and that this apparatus was present in the last common eukaryotic ancestor, but has been altered beyond recognition several times.

Amino Acid Sequence↗
PubMed2003

Assessing functional divergence in EF-1alpha and its paralogs in eukaryotes and archaebacteria.

A number of methods have recently been published that use phylogenetic information extracted from large multiple sequence alignments to detect sites that have changed properties in related protein families. In this study we use such methods to assess functional divergence between eukaryotic EF-1alpha (eEF-1alpha), archaebacterial EF-1alpha (aEF-1alpha) and two eukaryote-specific EF-1alpha paralogs-eukaryotic release factor 3 (eRF3) and Hsp70 subfamily B suppressor 1 (HBS1). Overall, the evolutionary modes of aEF-1alpha, HBS1 and eRF3 appear to significantly differ from that of eEF-1alpha. However, functionally divergent (FD) sites detected between aEF-1alpha and eEF-1alpha only weakly overlap with sites implicated as putative EF-1beta or aminoacyl-tRNA (aa-tRNA) binding residues in EF-1alpha, as expected based on the shared ancestral primary translational functions of these two orthologs. In contrast, FD sites detected between eEF-1alpha and its paralogs significantly overlap with the putative EF-1beta and/or aa-tRNA binding sites in EF-1alpha. In eRF3 and HBS1, these sites appear to be released from functional constraints, indicating that they bind neither eEF-1beta nor aa-tRNA. These results are consistent with experimental observations that eRF3 does not bind to aa-tRNA, but do not support the 'EF-1alpha-like' function recently proposed for HBS1. We re-assess the available genetic data for HBS1 in light of our analyses, and propose that this protein may function in stop codon-independent peptide release.

Amino Acid Sequence↗
PubMed2003

Evolution of glutamate dehydrogenase genes: evidence for lateral gene transfer within and between prokaryotes and eukaryotes.

BACKGROUND: Lateral gene transfer can introduce genes with novel functions into genomes or replace genes with functionally similar orthologs or paralogs. Here we present a study of the occurrence of the latter gene replacement phenomenon in the four gene families encoding different classes of glutamate dehydrogenase (GDH), to evaluate and compare the patterns and rates of lateral gene transfer (LGT) in prokaryotes and eukaryotes. RESULTS: We extend the taxon sampling of gdh genes with nine new eukaryotic sequences and examine the phylogenetic distribution pattern of the various GDH classes in combination with maximum likelihood phylogenetic analyses. The distribution pattern analyses indicate that LGT has played a significant role in the evolution of the four gdh gene families. Indeed, a number of gene transfer events are identified by phylogenetic analyses, including numerous prokaryotic intra-domain transfers, some prokaryotic inter-domain transfers and several inter-domain transfers between prokaryotes and microbial eukaryotes (protists). CONCLUSION: LGT has apparently affected eukaryotes and prokaryotes to a similar extent within the gdh gene families. In the absence of indications that the evolution of the gdh gene families is radically different from other families, these results suggest that gene transfer might be an important evolutionary mechanism in microbial eukaryote genome evolution.

Algal Proteins↗
PubMed2003

Phylogenetic analyses of diplomonad genes reveal frequent lateral gene transfers affecting eukaryotes.

BACKGROUND: Lateral gene transfer (LGT) is an important evolutionary mechanism among prokaryotes. The situation in eukaryotes is less clear; the human genome sequence failed to give strong support for any recent transfers from prokaryotes to vertebrates, yet a number of LGTs from prokaryotes to protists (unicellular eukaryotes) have been documented. Here, we perform a systematic analysis to investigate the impact of LGT on the evolution of diplomonads, a group of anaerobic protists. RESULTS: Phylogenetic analyses of 15 genes present in the genome of the Atlantic Salmon parasite Spironucleus barkhanus and/or the intestinal parasite Giardia lamblia show that most of these genes originated via LGT. Half of the genes are putatively involved in processes related to an anaerobic lifestyle, and this finding suggests that a common ancestor, which most probably was aerobic, of Spironucleus and Giardia adapted to an anaerobic environment in part by acquiring genes via LGT from prokaryotes. The sources of the transferred diplomonad genes are found among all three domains of life, including other eukaryotes. Many of the phylogenetic reconstructions show eukaryotes emerging in several distinct regions of the tree, strongly suggesting that LGT not only involved diplomonads, but also involved other eukaryotic groups. CONCLUSIONS: Our study shows that LGT is a significant evolutionary mechanism among diplomonads in particular and protists in general. These findings provide insights into the evolution of biochemical pathways in early eukaryote evolution and have important implications for studies of eukaryotic genome evolution and organismal relationships. Furthermore, "fusion" hypotheses for the origin of eukaryotes need to be rigorously reexamined in the light of these results.

Adaptation, Physiological↗
PubMed2003

Inferring functional constraints and divergence in protein families using 3D mapping of phylogenetic information.

Comparative sequence analysis has been used to study specific questions about the structure and function of proteins for many years. Here we propose a knowledge-based framework in which the maximum likelihood rate of evolution is used to quantify the level of constraint on the identity of a site. We demonstrate that site-rate mapping on 3D structures using datasets of rhodopsin-like G-protein receptors and alpha- and beta-tubulins provides an excellent tool for pinpointing the functional features shared between orthologous and paralogous proteins. In addition, functional divergence within protein families can be inferred by examining the differences in the site rates, the differences in the chemical properties of the side chains or amino acid usage between aligned sites. Two novel analytical methods are introduced to characterize rate- independent functional divergence. These are tested using a dataset of two classes of HMG-CoA reductases for which only one class can perform both the forward and reverse reaction. We show that functionally divergent sites occur in a cluster of sites interacting with the catalytic residues and that this information should facilitate the design of experimental strategies to directly test functional properties of residues.

Animals↗
PubMed2003

Estimation of rates-across-sites distributions in phylogenetic substitution models.

Previous work has shown that it is often essential to account for the variation in rates at different sites in phylogenetic models in order to avoid phylogenetic artifacts such as long branch attraction. In most current models, the gamma distribution is used for the rates-across-sites distributions and is implemented as an equal-probability discrete gamma. In this article, we introduce discrete distribution estimates with large numbers of equally spaced rate categories allowing us to investigate the appropriateness of the gamma model. With large numbers of rate categories, these discrete estimates are flexible enough to approximate the shape of almost any distribution. Likelihood ratio statistical tests and a nonparametric bootstrap confidence-bound estimation procedure based on the discrete estimates are presented that can be used to test the fit of a parametric family. We applied the methodology to several different protein data sets, and found that although the gamma model often provides a good parametric model for this type of data, rate estimates from an equal-probability discrete gamma model with a small number of categories will tend to underestimate the largest rates. In cases when the gamma model assumption is in doubt, rate estimates coming from the discrete rate distribution estimate with a large number of rate categories provide a robust alternative to gamma estimates. An alternative implementation of the gamma distribution is proposed that, for equal numbers of rate categories, is computationally more efficient during optimization than the standard gamma implementation and can provide more accurate estimates of site rates.

Evolution, Molecular↗
PubMed2002

Eukaryotic evolution: getting to the root of the problem.

Comparative analyses of multiple genes suggest most known eukaryotes can be classified into half a dozen 'super-groups'. A new investigation of the distribution of a fused gene pair amongst these 'super-groups' has greatly narrowed the possible positions of the root of the eukaryote tree, clarifying the broad outlines of early eukaryote evolution.

Animals↗
PubMed2002

Eukaryotic evolution: early origin of canonical introns.

Spliceosomal introns, one of the hallmarks of eukaryotic genomes, were thought to have originated late in evolution and were assumed not to exist in eukaryotes that diverged early -- until the discovery of a single intron with an aberrant splice boundary in the primitive 'protozoan' Giardia. Here we describe introns from a close relative of Giardia, Carpediemonas membranifera, that have boundary sequences of the normal eukaryotic type, indicating that canonical introns are likely to have arisen very early in eukaryotic evolution.

Animals↗
PubMed2002

Lateral transfer of an EF-1alpha gene: origin and evolution of the large subunit of ATP sulfurylase in eubacteria.

It is generally accepted that new genes arise via duplication and functional divergence of existing genes, in accordance with Ohno's model, now called "Mutation During Redundancy," or MDR. In this model, one of the two gene copies is free to acquire novel (although likely related) activities through mutation, since only one copy is required for its original function. However, duplication within a genome is not the only process that might give rise to this situation: acquisition of a functionally redundant gene by lateral gene transfer (LGT) could also initiate the MDR process. Here we describe a probable instance, involving LGT of an archaeal or eukaryotic elongation factor 1alpha (EF-1alpha) gene. The large subunit of ATP sulfurylase (CysN or the N-terminal portion of NodQ), found mainly in proteobacteria, is clearly related to translation elongation factors. However, our analyses show that cysN arose from an EF-1alpha gene initially acquired by LGT, not from a within-genome duplication of the resident EF-Tu gene. To our knowledge, this is the first unequivocal case of LGT followed by functional modification to be described; this mechanism could be a potentially important force in establishing genes with novel functions in genomes.

Amino Acid Sequence↗
PubMed2002

Gene duplication and gene conversion shape the evolution of archaeal chaperonins.

Chaperonins are multi-subunit double-ring complexes that mediate the folding of nascent or denatured proteins. Gene duplication has been a potent force in the evolution of chaperonins in Archaea. Here we show that gene conversion has also been an important factor. We utilized a novel maximum likehood-based phylogenetic method for scanning DNA sequence alignments for regions of anomalous phylogenetic signal, such as those affected by gene conversion. Our results suggest that in crenarchaeotes, where an ancient gene duplication producing alpha and beta subunits took place in the common ancestor of the Pyrodictium, Aeropyrum, Pyrobaculum and Sulfolobus lineages, multiple independent gene conversions have occurred between the alpha and beta genes independently in each of these groups. Significantly, the conversions have repeatedly homogenized the region of the gene encoding the substrate-binding domain. This suggests that while the alpha and beta subunits in crenarchaeotes share only 50-60% overall amino acid sequence identity, they do not possess distinct roles in the binding of substrate. Cryptic gene conversion between distantly related paralogs may be more common than is currently appreciated, and could be a significant factor in slowing the functional differentiation of proteins encoded by duplicate genes long after their duplication.

Amino Acid Sequence↗
PubMed2002

A cyanobacterial gene in nonphotosynthetic protists--an early chloroplast acquisition in eukaryotes?

Since the incorporation of mitochondria and chloroplasts (plastids) into the eukaryotic cell by endosymbiosis, genes have been transferred from the organellar genomes to the nucleus of the host, via an ongoing process known as endosymbiotic gene transfer. Accordingly, in photosynthetic eukaryotes, nuclear genes with cyanobacterial affinity are believed to have originated from endosymbiotic gene transfer from chloroplasts. Analysis of the Arabidopsis thaliana genome has shown that a significant fraction (2%-9%) of the nuclear genes have such an endosymbiotic origin. Recently, it was argued that 6-phosphogluconate dehydrogenase (gnd)-the second enzyme in the oxidative pentose phosphate pathway-was one such example. Here we show that gnd genes with cyanobacterial affinity also are present in several nonphotosynthetic protistan lineages, such as Heterolobosea, Apicomplexa, and parasitic Heterokonta. Current data cannot definitively resolve whether these groups acquired the gnd gene by primary and/or secondary endosymbiosis or via an independent lateral gene transfer event. Nevertheless, our data suggest that chloroplasts were introduced into eukaryotes much earlier than previously thought and that several major groups of heterotrophic eukaryotes have secondarily lost photosynthetic plastids.

Arabidopsis↗