Phylogenetic positions of several amitochondriate protozoa

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group―Archezoa. The main evidence for this is their ‘lacking mitochondria’ and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of ‘long-branch attraction’ appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplo-monads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be poly-phyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

Phylogenetic positions of several amitochondriate protozoa-Evidence from phylogenetic analysis of DNA topoisomerase II

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group―Archezoa. The main evidence for this is their 'lacking mitochondria' and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of 'long-branch attraction' appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplomonads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be polyphyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

Spliceosomal introns as tools for genomic and evolutionary analysis

Over the past 5 years, the availability of dozens of whole genomic sequences from a wide variety of eukaryotic lineages has revealed a very large amount of information about the dynamics of intron loss and gain through eukaryotic history, as well as the evolution of intron sequences. Implicit in these advances is a great deal of information about the structure and evolution of surrounding sequences. Here, we review the wealth of ways in which structures of spliceosomal introns as well as their conservation and change through evolution may be harnessed for evolutionary and genomic analysis. First, we discuss uses of intron length distributions and positions in sequence assembly and annotation, and for improving alignment of homologous regions. Second, we review uses of introns in evolutionary studies, including the utility of introns as indicators of rates of sequence evolution, for inferences about molecular evolution, as signatures of orthology and paralogy, and for estimating rates of nucleotide substitution. We conclude with a discussion of phylogenetic methods utilizing intron sequences and positions.     

Phylogenetic positions of several amitochondriate protozoa

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, En-tamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group—Archezoa. The main evidence for this is their ‘lacking mitochondria’ and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of ‘long-branch attraction’ appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplomonads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be polyphyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

Phylogenetic analysis of 70 kD heat shock protein sequences suggests a chimeric origin for the eukaryotic cell nucleus

Background: The evolutionary relationships between archaebacteria, eubacteria and eukaryotic cells are of central importance in biology. The current view is that each of these three groups of organisms constitutes a monophyletic domain, and that eukaryotic cells have evolved from an archaebacterial ancestor. Recent studies on a number of highly conserved protein sequences do not, however, support this view and raise important questions concerning the evolutionary relationships between all extant organisms, particularly regarding the origin of eukaryotic cells.Results We have used sequences of 70 kD heat shock protein (hsp70) — the most conserved protein found to date in all species — to examine the evolutionary relationship between various species. We have obtained two new archaebacterial hsp70 sequences from the species, Thermoplasma acidophilum and Halobacterium cutirubrum. A global comparison of hsp70 sequences, including our two new sequences, shows that all known archaebacterial homologs share a number of sequence signatures with the Gram-positive group of bacteria that are not found in any other prokaryotic or eukaryotic species. In contrast, the eukaryotic homologs are shown to share a number of unique sequence features with the Gram-negative bacteria that are not present in any archaebacteria. Detailed phylogenetic analyses of hsp70 sequences strongly support a specific evolutionary relationship between archaebacteria and Gram-positive bacteria on the one hand, and Gram-negative bacteria and eukaryotes on the other. The phylogenetic analyses also indicate a polyphyletic branching of archaebacteria within the Gram-positive species. The possibility that the observed relationships are due to horizontal gene transfers can be excluded on the basis of sequence characteristics of different groups of homologs.Conclusion Our results do not support the view that archaebacteria constitute a monophyletic domain, but instead suggest a close evolutionary linkage between archaebacteria and Gram-positive bacteria. Furthermore, in contrast to the presently accepted view, eukaryotic hsp70s show a close and specific relationship to those from Gram-negative species. To explain the phylogenies based on different gene sequences, a chimeric model for the origin of the eukaryotic cell nucleus involving fusion between an archaebacterium and a Gram-negative eubacterium is proposed. Several predictions from the chimeric model are discussed.     

Comparative analysis of hemagglutinin of 2009 H1N1 influenza A pandemic indicates its evolution to 1918 H1N1 pandemic

To gain insight into the possible origin of the hemagglutinin of 2009 outbreak, we performed its comparative analysis with hemagglutinin of influenza viral strains from 2005 to 2008 and the past pandemics of 1977, 1968, 1957 and 1918. This insilico analysis showed a maximum sequence similarity between 2009 and 1918 pandemics. Primary structure analysis, antigenic and glycosylation site analyses revealed that this protein has evolved from 1918 pandemic. Phylogenetic analysis of HA amino acid sequence of 2009 influenza A(H1N1) viruses indicated that this virus possesses a distinctive evolutionary trait with 1918 influenza A virus. Although the disordered sequences are different among all the isolates, the disordered positions and sequences between 2009 and 1918 isolates show a greater similarity. Thus these analyses contribute to the evidence of the evolution of 2009 pandemic from 1918 influenza pandemic. This is the first computational evolutionary analysis of HA protein of 2009 H1N1 pandemic.     

Evolution and Horizontal Transfer of dUTPase-Encoding Genes in Viruses and Their Hosts

dUTPase is a ubiquitous and essential enzyme responsible for regulating cellular levels of dUTP. The dut gene exists as single, tandemly duplicated, and tandemly triplicated copies. Crystallized single-copy dUTPases have been shown to assemble as homotrimers. dUTPase is encoded as an auxiliary gene in a number of virus genomes. The origin of viral dut genes has remained unresolved since their initial discovery. A comprehensive analysis of dUTPase amino acid sequence relationships was performed to explore the evolutionary dynamics of dut in viruses and their hosts. Our data set, comprised of 24 host and 51 viral sequences, includes representative sequences from available eukaryotes, archaea, eubacteria cells, and viruses, including herpesviruses. These amino acid sequences were aligned by using a hidden Markov model approach developed to align divergent data. Known secondary structures from single-copy crystals were mapped onto the aligned duplicate and triplicate sequences. We show how duplicated dUTPases might fold into a monomer, and we hypothesize that triplicated dUTPases also assemble as monomers. Phylogenetic analysis revealed at least five viral dUTPase sequence lineages in well-supported monophyletic clusters with eukaryotic, eubacterial, and archaeal hosts. We have identified all five as strong examples of horizontal transfer as well as additional potential transfer of dut genes among eubacteria, between eubacteria and viruses, and between retroviruses. The evidence for horizontal transfers is particularly interesting since eukaryotic dut genes have introns, while DNA virus dut genes do not. This implies that an intermediary retroid agent facilitated the horizontal transfer process between host mRNA and DNA viruses.     

Ab initio modeling led annotation suggests nucleic acid binding function for many DUFs

Expansions of sequence databases driven by new sequencing technology continue apace. These result in a continuous supply of protein sequences and domains that cannot be straightforwardly annotated by simple homology methods. For these, structure-based function prediction may contribute to an improved annotation. Here, short Domains of Unknown Function (DUFs) are ab initio modeled with ROSETTA and screened for likely nucleic acid binding function. Thirty-two DUFs are thereby predicted to have a nucleic acid binding function. In most cases, additional evidence supporting that function could be obtained from structure comparison, domain architectures, distant evolutionary relationships, genome context or protein-protein interaction data. These predictions contribute to the function annotation of thousands of proteins.     

Human prolactin. cDNA structural analysis and evolutionary comparisons

Prolactin (Prl), growth hormone, and chorionic sommatomammotropin form a set (the "Prl set") of hormones which is thought to have evolved from a common ancestral gene. This assumption is based on several lines of evidence: overlap in their biological and immunological properties, similarities in their amino acid sequences, and homologies in the nucleic acid sequences of their structural genes. In the current study we report the cloning, amplification in bacteria, and sequence analysis of DNA complementary to Prl mRNA isolated from human pituitary Prl-secreting adenomas. The cloned DNA contains 914 bases, which includes the entire coding sequence of human prePrl as well as portions of the 5- and 3'-untranslated regions of the mRNA. The amino acid sequence predicted by our data differs from a previously reported amino acid sequence in 8 positions. With the results of this study we can now compare in one species the nucleotide sequences of the structural gene coding for each of the hormones of the Prl set. The sequence divergence at replacement sites is used to establish an evolutionary clock for the Prl set of genes. Using this clock, we postulate that the chromosomal segregation of human Prl and human growth hormone occurred about 392 million years ago and that growth hormone and chorionic sommatomammotropin underwent an intrachromosomal recombination within the last 10 million years.     

Prokaryotic and eukaryotic pyridoxal-dependent decarboxylases are homologous

Summary A database search has revealed significant and extensive sequence similarities among prokaryotic and eukaryotic pyridoxal phosphate (PLP)-dependent decarboxylases, includingDrosophila glutamic acid decarboxylase (GAD) and bacterial histidine decarboxylase (HDC). Based on these findings, the sequences of seven PLP-dependent decarboxylases from five different organisms have been aligned to derive a consensus sequence for this family of enzymes. In addition, quantitative methods have been employed to calculate the relative evolutionary distances between pairs of the decarboxylases comprising this family. The multiple sequence analysis together with the quantitative results strongly suggest an ancient and common origin for all PLP-dependent decarboxylases. This analysis also indicates that prokaryotic and eukaryotic HDC activities evolved independently. Finally, a sensitive search algorithm (PROFILE) was unable to detect additional members of this decarboxylase family in protein sequence databases.     

Molecular evolution of the carbonic anhydrase genes: Calculation of divergence time for mouse carbonic anhydrase I and II

Summary A cDNA clone in pBR322 that cross-hybridizes with a mouse carbonic anhydrase form II (CAII) probe has been sequenced and identified as mouse carbonic anhydrase form I (CAI). The 1224-base-pair clone encodes the entire 260-amino-acid protein and appears to contain an Alu-like element in the 3 untranslated region. The deduced amino acid sequence exhibits 77% homology to human CAI and contains 17 of the 20 residues that are considered unique to and invariant for all mammalian CAI isozymes. The results of a detailed comparison of the nucleic acid sequences spanning the coding regions of mouse CAI and rabbit CAI have been used to calibrate an evolutionary clock for the carbonic anhydrases (CAs). These data have been applied to a comparison of the mouse CAI and CAII nucleic acid sequences to calculate the divergence time between the two genes. The divergence-time calculation provides the first estimation of the evolutionary relationship between CAs based entirely on nucleotide sequence comparison.     

Locked nucleic acid flow cytometry-fluorescence in situ hybridization (LNA flow-FISH): a method for bacterial small RNA detection

Fluorescence in situ hybridization (FISH) is a powerful technique that is used to detect and localize specific nucleic acid sequences in the cellular environment. In order to increase throughput, FISH can be combined with flow cytometry (flow-FISH) to enable the detection of targeted nucleic acid sequences in thousands of individual cells. As a result, flow-FISH offers a distinct advantage over lysate/ensemble-based nucleic acid detection methods because each cell is treated as an independent observation, thereby permitting stronger statistical and variance analyses. These attributes have prompted the use of FISH and flow-FISH methods in a number of different applications and the utility of these methods has been successfully demonstrated in telomere length determination, cellular identification and gene expression, monitoring viral multiplication in infected cells, and bacterial community analysis and enumeration. Traditionally, the specificity of FISH and flow-FISH methods has been imparted by DNA oligonucleotide probes. Recently however, the replacement of DNA oligonucleotide probes with nucleic acid analogs as FISH and flow-FISH probes has increased both the sensitivity and specificity of each technique due to the higher melting temperatures (T(m)) of these analogs for natural nucleic acids. Locked nucleic acid (LNA) probes are a type of nucleic acid analog that contain LNA nucleotides spiked throughout a DNA or RNA sequence. When coupled with flow-FISH, LNA probes have previously been shown to outperform conventional DNA probes and have been successfully used to detect eukaryotic mRNA and viral RNA in mammalian cells. Here we expand this capability and describe a LNA flow-FISH method which permits the specific detection of RNA in bacterial cells (Figure 1). Specifically, we are interested in the detection of small non-coding regulatory RNA (sRNA) which have garnered considerable interest in the past few years as they have been found to serve as key regulatory elements in many critical cellular processes. However, there are limited tools to study sRNAs and the challenges of detecting sRNA in bacterial cells is due in part to the relatively small size (typically 50-300 nucleotides in length) and low abundance of sRNA molecules as well as the general difficulty in working with smaller biological cells with varying cellular membranes. In this method, we describe fixation and permeabilzation conditions that preserve the structure of bacterial cells and permit the penetration of LNA probes as well as signal amplification steps which enable the specific detection of low abundance sRNA (Figure 2).     

Directed Molecular Evolution

We propose the existence of a relationship of stereochemical complementarity between gene sequences that code for interacting components: nucleic acid-nucleic acid, protein-protein and protein-nucleic acid. Such a relationship would impose evolutionary constraints on the DNA sequences themselves, thus retaining these sequences and governing the direction of the evolutionary process. Therefore, we propose that prebiotic, template-directed autocatalytic synthesis of mutally cognate peptides and polynucleotides resulted in their amplification and evolutionary conservation in contemporary prokaryotic and eukaryotic organisms as a genetic regulatory apparatus. If this proposal is correct, then the relationships between the sequences in DNA coding for these interactions constitute a life code of which the genetic code is only one aspect of the many related interactions encoded in DNA.     

Looking through the Lens of the Ribosome Biogenesis Evolutionary History: Possible Implications for Archaeal Phylogeny and Eukaryogenesis

Our understanding of microbial diversity and its evolutionary relationships has increased substantially over the last decade. Such an understanding has been greatly fueled by culture-independent metagenomics analyses. However, the outcome of some of these studies and their biological and evolutionary implications, such as the origin of the eukaryotic lineage from the recently discovered archaeal Asgard superphylum, is debated. The sequences of the ribosomal constituents are amongst the most used phylogenetic markers. However, the functional consequences underlying the analysed sequence diversity and their putative evolutionary implications are essentially not taken into consideration. Here, we propose to exploit additional functional hallmarks of ribosome biogenesis to help disentangle competing evolutionary hypotheses. Using selected examples, such as the multiple origins of halophily in archaea or the evolutionary relationship between the Asgard archaea and Eukaryotes, we illustrate and discuss how function-aware phylogenetic framework can contribute to refining our understanding of archaeal phylogeny and the origin of eukaryotic cells.     

Evolution of HSP70 gene and its implications regarding relationships between archaebacteria,eubacteria, and eukaryotes

The 70-kDa heat-shock protein (HSP70) constitutes the most conserved protein present in all organisms that is known to date. Based on global alignment of HSP70 sequences from organisms representing all three domains, numerous sequence signatures that are specific for prokaryotic and eukaryotic homologs have been identified. HSP70s from the two archaebacterial species examined (viz., Halobacterium marismortui and Methanosarcina mazei) have been found to contain all eubacterial but no eukaryotic signature sequences. Based on several novel features of the HSP70 family of proteins (viz., presence of tandem repeats of a 9-amino-acid [a.a.] polypeptide sequence and structural similarity between the first and second quadrants of HSP70, homology of the N-terminal half of HSP70 to the bacterial MreB protein, presence of a conserved insert of 23–27 a.a. in all HSP70s except those from archaebacteria and gram-positive eubacteria) a model for the evolution of HSP70 gene from an early stage is proposed. The HSP70 homologs from archaebacteria and gram-positive bacteria lacking the insert in the N-terminal quadrants are indicated to be the ancestral form of the protein. Detailed phylogenetic analyses of HSP70 sequence data (viz., by bootstrap analyses, maximum parsimony, and maximum likelihood methods) provide evidence that archaebacteria are not monophyletic and show a close evolutionary linkage with the gram-positive eubacteria. These results do not support the traditional archaebacterial tree, where a close relationship between archaebacterial and eukaryotic homologs is observed. To explain the phylogenies based on HSP70 and other gene sequences, a model for the origin of eukaryotic cells involving fusion between archaebacteria and gram-negative eubacteria is proposed. Correspondence to: R. S. Gupta     

Methods for the analysis of developmental sequence data

SUMMARY Heterochrony, evolutionary changes in developmental timing, can be studied either by examining changes in growth or changes in the sequence of developmental events. Developmental sequence data has the potential to address many important questions in the field of developmental evolution, but methodological challenges remain due to the biological and logical dependence of events in a ranked sequence. In the past 10 years, the study of sequence heterochrony has undergone a rebirth, with the creation of several new methods for the analysis of this type of data. These methods can be divided into two broad categories: phenetic comparisons between terminal taxa that strive to uncover integrations within the developmental sequences and putative shared sequence heterochronies, and phylogeny-based methods that derive ancestor-descendent sequence heterochronies and establish statements of sequence evolution. In this review, we will discuss the strengths and weaknesses of the methodologies that have been proposed to quantitatively examine developmental sequence data, and studies that have attempted to implement them in an evolutionary context.     

Conformation and structural transitions in the EF-hands of calmodulin

Calcium plays a key role in cellular signal transduction. Calmodulin, a protein binding four calcium ions, is found in all eukaryotic cells and is believed to activate such processes. The calcium binding loop found in this protein, the canonical EF-hand, is also found in a large number of other proteins such as troponins, parvalbumins, calbindins etc. Earlier analysis of the amino acid sequences of these proteins with a view of understanding evolution of protein families and signaling mechanisms have provided extensive evidence for a characteristic double gene duplication event in this family of proteins. These analyses have been extended here to the three dimensional structures and the biophysical properties of the sequence segments of calmodulin EF-hands. The clear evolutionary history that shows up in sequences is not reflected as clearly in the conformation of individual EF-hands, which may be a consequence of the much higher conservation pressure on the structure. Some evidence for the proposed gene duplication is implicit in the apo-holo structural transitions of the EF-hands. The profile of amino acid properties that might be significant for calcium binding, however, clearly reflects the gene duplication. These profiles might also provide insightful information on the calcium affinity of the EF-hand motifs and the nature of amino acid residues that constitute them.     

Nitric oxide synthase sequences in the marine fish Stenotomus chrysops and the sea urchin Arbacia punctulata, and phylogenetic analysis of nitric oxide synthase calmodulin-binding domains

The phylogenetic distribution and structural diversity of the nitric oxide synthases (NOS) remain important and issues that are little understood. We present sequence information, as well as phylogenetic analysis, for three NOS cDNAs identified in two non-mammalian species: the vertebrate marine teleost fish Stenotomus chrysops (scup) and the invertebrate echinoderm Arbacia punctulata (sea urchin). Partial gene sequences containing the well-conserved calmodulin (CaM)-binding domain were amplified by RT-PCR. Identical 375-bp cDNAs were amplified from scup brain, heart, liver and spleen; this sequence shares 82% nucleic acid and 91% predicted amino acid identity with the corresponding region of human neuronal NOS. A 387-bp cDNA was amplified from sea urchin ovary and testes; this sequence shares 72% nucleic acid identity and 65% deduced amino acid identity with human neuronal NOS. A second cDNA of 381 bp was amplified from sea urchin ovary and it shares 66% nucleic acid and 57% deduced amino acid identity with the first sea urchin sequence. Together with earlier reports of neuronal and inducible NOS sequences in fish, these data indicate that multiple NOS isoforms exist in non-mammalian species. Phylogenetic analysis of these sequences confirms the conserved nature of NOS, particularly of the calmodulin-binding domains.