Primary Structure and Phylogenetic Relationships of a Malate Dehydrogenase Gene from Giardia lamblia

The lactate and malate dehydrogenases comprise a complex protein superfamily with multiple enzyme homologues found in eubacteria, archaebacteria, and eukaryotes. In this study we describe the sequence and phylogenetic relationships of a malate dehydrogenase (MDH) gene from the amitochondriate diplomonad protist, Giardia lamblia. Parsimony, distance, and maximum-likelihood analyses of the MDH protein family solidly position G. lamblia MDH within a eukaryote cytosolic MDH clade, to the exclusion of chloroplast, mitochondrial, and peroxisomal homologues. Furthermore, G. lamblia MDH is specifically related to a homologue from Trichomonas vaginalis. This MDH topology, together with published phylogenetic analyses of β-tubulin, chaperonin 60, valyl-tRNA synthetase, and EF-1α, suggests a sister-group relationship between diplomonads and parabasalids. Since these amitochondriate lineages contain genes encoding proteins which are characteristic of mitochondria and α-proteobacteria, their shared ancestry suggests that mitochondrial properties were lost in the common ancestor of both groups. Received: 14 September 1998 / Accepted: 29 December 1998     

Actin Phylogeny Identifies Mesostigma viride as a Flagellate Ancestor of the Land Plants

Green algae and land plants trace their evolutionary history to a unique common ancestor. This ``green lineage' is phylogenetically subdivided into two distinct assemblages, the Chlorophyta and the Streptophyta. The Chlorophyta includes the Chlorophyceae, Trebouxiophyceae, Ulvophyceae, and Prasinopohyceae, whereas the Streptophyta includes the Charophyceae plus the bryophytes, ferns, and all other multicellular land plants (Embryophyta). The Prasinophyceae is believed to contain the earliest divergences within the green lineage. Phylogenetic analyses using rDNA sequences identify the prasinophytes as a paraphyletic taxon that diverges at the base of the Chlorophyta. rDNA analyses, however, provide ambiguous results regarding the identity of the flagellate ancestor of the Streptophyta. We have sequenced the actin-encoding cDNAs from Scherffelia dubia (Prasinophyceae), Coleochaete scutata, Spirogyra sp. (Charophyceae), and the single-copy actin gene from Mesostigma viride (Prasinophyceae). Phylogenetic analyses show Mesostigma to be the earliest divergence within the Streptophyta and provide direct evidence for a scaly, biflagellate, unicellular ancestor for this lineage. This result is supported by the existence of two conserved actin-coding region introns (positions 20-3, 152-1), and one intron in the 5′-untranslated region of the actin gene shared by Mesostigma and the embryophytes. Received: 10 July 1997 / Accepted: 9 April 1998     

Apolipophorin II/I, Apolipoprotein B, Vitellogenin, and Microsomal Triglyceride Transfer Protein Genes Are Derived from a Common Ancestor

Large lipid transfer proteins (LLTP) are nonexchangeable apolipoproteins and intracellular lipid-exchange proteins involved in the assembly, secretion, and metabolism of lipoproteins. We have identified contiguous conserved sequence motifs in alignments of insect apolipophorin II/I precursor (apoLp-II/I), human apolipoprotein B (apoB), invertebrate and vertebrate vitellogenins (VTG), and the large subunit of mammalian microsomal triglyceride transfer protein (MTP). Conserved motifs present in the N-terminal part of nonexchangeable apolipoproteins encompass almost completely the large subunit of MTP, suggesting a derivation from a common ancestral functional unit, termed large lipid transfer (LLT) module. Divergence of LLTP from a common ancestor is supported by (1) the statistical significance of the combined match scores obtained after motif-based database searches, (2) the presence of several identical amino acid residues in all LLTP sequences currently available, (3) the conservation of hydrophobic clusters in an α-helical domain, (4) the phylogenetic analysis of the conserved sequences related to the von Willebrand factor D (VWD) module identified in nonexchangeable apolipoproteins, and (5) the presence of four and one ancestral exon boundaries in the LLT and VWD modules, respectively. Our data indicate that the genes coding for apoLp-II/I, apoB, VTG, and the MTP large subunit are members of the same multigene superfamily. LLTP have emerged from an ancestral molecule designed to ensure a pivotal event in the intracellular and extracellular transfer of lipids and liposoluble substances. Received: 8 June 1998 / Accepted: 15 February 1999     

The Evolutionary History of Carbamoyltransferases: A Complex Set of Paralogous Genes Was Already Present in the Last Universal Common Ancestor

Forty-four sequences of ornithine carbamoyltransferases (OTCases) and 33 sequences of aspartate carbamoyltransferases (ATCases) representing the three domains of life were multiply aligned and a phylogenetic tree was inferred from this multiple alignment. The global topology of the composite rooted tree (each enzyme family being used as an outgroup to root the other one) suggests that present-day genes are derived from paralogous ancestral genes which were already of the same size and argues against a mechanism of fusion of independent modules. A closer observation of the detailed topology shows that this tree could not be used to assess the actual order of organismal descent. Indeed, this tree displays a complex topology for many prokaryotic sequences, with polyphyly for Bacteria in both enzyme trees and for the Archaea in the OTCase tree. Moreover, representatives of the two prokaryotic Domains are found to be interspersed in various combinations in both enzyme trees. This complexity may be explained by assuming the occurrence of two subfamilies in the OTCase tree (OTC α and OTC β) and two other ones in the ATCase tree (ATC I and ATC II). These subfamilies could have arisen from duplication and selective losses of some differentiated copies during the successive speciations. We suggest that Archaea and Eukaryotes share a common ancestor in which the ancestral copies giving the present-day ATC II/OTC β combinations were present, whereas Bacteria comprise two classes: one containing the ATC II/OTC α combination and the other harboring the ATC I/OTC β combination. Moreover, multiple horizontal gene transfers could have occurred rather recently amongst prokaryotes. Whichever the actual history of carbamoyltransferases, our data suggest that the last common ancestor to all extant life possessed differentiated copies of genes coding for both carbamoyltransferases, indicating it as a rather sophisticated organism.     

Molecular Data from the Chloroplast rpoC1 Gene Suggest a Deep and Distinct Dichotomy of Contemporary Spermatophytes into Two Monophyla: Gymnosperms (Including Gnetales) and Angiosperms

Partial sequences of the rpoC1 gene from two species of angiosperms and three species of gymnosperms (8330 base pairs) were determined and compared. The data obtained support the hypothesis that angiosperms and gymnosperms are monophyletic and none of the recent groups of the latter is sister to angiosperms. Received: 20 November 1998 / Accepted: 26 April 1999     

Dynamic Diversification from a Putative Common Ancestor of Scorpion Toxins Affecting Sodium, Potassium, and Chloride Channels

Scorpions have survived successfully over millions of years without detectable changes in their morphology. Instead, they have developed an efficient alomonal machinery and a stinging device supporting their needs for prey and defense. They produce a large variety of polypeptidic toxins that bind and modulate ion channel conductance in excitable tissues. The binding site, mode of action, and chemical properties of many toxins have been studied extensively, but little is known about their genomic organization and diversity. Genes representing each of the major classes of Buthidae scorpion toxins, namely, ``long' toxins, affecting sodium channels (alpha, depressant, and excitatory), and ``short' toxins, affecting potassium and chloride channels, were isolated from a single scorpion segment and analyzed. Each toxin type was found to be encoded by a gene family. Regardless of toxin length, 3-D structure, and site of action, all genes contain A+T-rich introns that split, at a conserved location, an amino acid codon of the signal sequence. The introns vary in length and sequence but display identical boundaries, agree with the GT/AG splice junctions, and contain T-runs downstream of a putative branch point, 5′-TAAT-3′. Despite little sequence similarity among all toxin classes, the conserved gene organization, intron features, and common cysteine-stabilized α-helical (CSH) core connecting an α-helix to a three-stranded β-sheet suggest, that they all evolved from an ancestral common progenitor. Furthermore, the vast diversity found among genomic copies, cDNAs, and their protein products for each toxin suggests an extensive evolutionary process of the scorpion ``pharmaceutical factory,' whose success is due, most likely, to the inherent permissiveness of the toxin exterior to structural alterations. Received: 16 March 1998 / Accepted: 30 July 1998     

Limitations of Metazoan 18S rRNA Sequence Data: Implications for Reconstructing a Phylogeny of the Animal Kingdom and Inferring the Reality of the Cambrian Explosion

We document the phylogenetic behavior of the 18S rRNA molecule in 67 taxa from 28 metazoan phyla and assess the effects of among-site rate variation on reconstructing phylogenies of the animal kingdom. This empirical assessment was undertaken to clarify further the limits of resolution of the 18S rRNA gene as a phylogenetic marker and to address the question of whether 18S rRNA phylogenies can be used as a source of evidence to infer the reality of a Cambrian explosion. A notable degree of among-site rate variation exists between different regions of the 18S rRNA molecule, as well as within all classes of secondary structure. There is a significant negative correlation between inferred number of nucleotide substitutions and phylogenetic information, as well as with the degree of substitutional saturation within the molecule. Base compositional differences both within and between taxa exist and, in certain lineages, may be associated with long branches and phylogenetic position. Importantly, excluding sites with different degrees of nucleotide substitution significantly influences the topology and degree of resolution of maximum-parsimony phylogenies as well as neighbor-joining phylogenies (corrected and uncorrected for among-site rate variation) reconstructed at the metazoan scale. Together, these data indicate that the 18S rRNA molecule is an unsuitable candidate for reconstructing the evolutionary history of all metazoan phyla, and that the polytomies, i.e., unresolved nodes within 18S rRNA phylogenies, cannot be used as a single or reliable source of evidence to support the hypothesis of a Cambrian explosion. Received: 9 December 1997 / Accepted: 23 March 1998     

Phylogenetic analysis of viroid and viroid-like satellite RNAs from plants: a reassessment

The proposed monophyletic origin of a group of subviral plant pathogens (viroids and viroid-like satellite RNAs), as well as the phylogenetic relationships and the resulting taxonomy of these entities, has been recently questioned. The criticism comes from the (apparent) lack of sequence similarity among these RNAs necessary to reliably infer a phylogeny. Here we show that, despite their low overall sequence similarity, a sequence alignment manually adjusted to take into account all the local similarities and the insertions/deletions and duplications/rearrangements described in the literature for viroids and viroid-like satellite RNA, along with the use of an appropriate estimator of genetic distances, constitutes a data set suitable for a phylogenetic reconstruction. When the likelihood-mapping method was applied to this data set, the tree-likeness obtained was higher than that corresponding to a sequence alignment that does not take into consideration the local similarities. In addition, bootstrap analysis also supports the major groups previously proposed and the reconstruction is consistent with the biological properties of this RNAs. Received: 17 January 2001 / Accepted: 16 March 2001     

Unusual and Strongly Structured Sequence Variation in a Complex Satellite DNA Family from the Nematode Meloidogyne chitwoodi

An AluI satellite DNA family has been isolated in the genome of the root-knot nematode Meloidogyne chitwoodi. This repeated sequence was shown to be present at approximately 11,400 copies per haploid genome, and represents about 3.5% of the total genomic DNA. Nineteen monomers were cloned and sequenced. Their length ranged from 142 to 180 bp, and their A + T content was high (from 65.7 to 79.1%), with frequent runs of As and Ts. An unexpected heterogeneity in primary structure was observed between monomers, and multiple alignment analysis showed that the 19 repeats could be unambiguously clustered in six subfamilies. A consensus sequence has been deduced for each subfamily, within which the number of positions conserved is very high, ranging from 86.7% to 98.6%. Even though blocks of conserved regions could be observed, multiple alignment of the six consensus sequences did not enable the establishment of a general unambiguous consensus sequence. Screening of the six consensus sequences for evidence of internal repeated subunits revealed a 6-bp motif (AAATTT), present in both direct and inverted orientation. This motif was found up to nine times in the consensus sequences, also with the occurrence of degenerated subrepeats. Along with the meiotic parthenogenetic mode of reproduction of this nematode, such structural features may argue for the evolution of this satellite DNA family either (1) from a common ancestral sequence by amplification followed by mechanisms of sequence divergence, or (2) through independent mutations of the ancestral sequence in isolated amphimictic nematode populations and subsequent hybridization events. Overall, our results suggest the ancient origin of this satellite DNA family, and may reflect for M. chitwoodi a phylogenetic position close to the ancestral amphimictic forms of root-knot nematodes. Received: 23 April 1997 / Accepted: 9 July 1997     

Anopheles stephensi Dox-A2 Shares Common Ancestry with Genes from Distant Groups of Eukaryotes Encoding a 26S Proteasome Subunit and Is in a Conserved Gene Cluster

The sequence of a cloned Anopheles stephensi gene showed 72% inferred amino acid identity with Drosophila melanogaster Dox-A2 and 93% with its putative ortholog in Anopheles gambiae. Dox-A2 is the reported but herein disputed structural locus for diphenol oxidase A2. Database searches identified Dox-A2 related gene sequences from 15 non-insect species from diverse groups. Phylogenetic trees based on alignments of inferred protein sequences, DNA, and protein motif searches and protein secondary structure predictions produced results consistent with expectations for genes that are orthologous. The only inconsistency was that the C-terminus appears to be more primitive in the yeasts than in plants. In mammals, plants, and yeast these genes have been shown to code for a non-ATPase subunit of the PA700 (19S) regulatory complex of 26S proteasome. The analyses indicated that the insect genes contain no divergent structural features, which taken within an appraisal of all available data, makes the reported alternative function highly improbable. A plausible additional role, in which the 26S proteasome is implicated in regulation of phenol oxidase, would also apply to at least the mammalian genes. No function has yet been reported for the other included sequences. These were from genome projects and included Caenorhabiditus elegans, Arabidopsis thaliana, Fugu rubripes, and Toxoplasma gondii. A consensus of the results predicts a protein containing exceptionally long stretches of helix with a hydrophilic C-terminus. Phosphorylation site motifs were identified at two conserved positions. Possible SRY and GATA-1 binding motifs were found at conserved positions upstream of the mosquito genes. The location of A. stephensi Dox-A2 was determined by in situ hybridization at 34D on chromosome arm 3R. It is in a conserved gene cluster with respect to the other insects. However, the A. stephensi cluster contains a gene showing significant sequence identity to human and pigeon carnitine acetyltransferase genes, therefore showing divergence with the distal end of the D. melanogaster cluster. Received: 3 July 1998 / Accepted: 22 December 1999     

Partial Sequence of a Sponge Mitochondrial Genome Reveals Sequence Similarity to Cnidaria in Cytochrome Oxidase Subunit II and the Large Ribosomal RNA Subunit

A 2550-bp portion of the mitochondrial genome of a Demosponge, genus Tetilla, was amplified from whole genomic DNA extract and sequenced. The sequence was found to code for the 3′ end of the 16S rRNA gene, cytochrome c oxidase subunit II, a lysine tRNA, ATPase subunit 8, and a 5′ portion of ATPase subunit 6. The Porifera cluster distinctly within the eumetazoan radiation, as a sister group to the Cnidaria. Also, the mitochondrial genetic code of this sponge is likely identical to that found in the Cnidaria. Both the full COII DNA and protein sequences and a portion of the 16S rRNA gene were found to possess a striking similarity to published Cnidarian mtDNA sequences, allying the Porifera more closely to the Cnidaria than to any other metazoan phylum. The gene arrangement, COII—tRNA^Lys—ATP8—ATP6, is observed in many Eumetazoan phyla and is apparently ancestral in the metazoa. Received: 24 November 1997 / Accepted: 14 September 1998     

A unique fungal lysine biosynthesis enzyme shares a common ancestor with tricarboxylic acid cycle and leucine biosynthetic enzymes found in diverse organisms

Fungi have evolved a unique α-aminoadipate pathway for lysine biosynthesis. The fungal-specific enzyme homoaconitate hydratase from this pathway is moderately similar to the aconitase-family proteins from a diverse array of taxonomic groups, which have varying modes of obtaining lysine. We have used the similarity of homoaconitate hydratase to isopropylmalate isomerase (serving in leucine biosynthesis), aconitase (from the tricarboxylic acid cycle), and iron-responsive element binding proteins (cytosolic aconitase) from fungi and other eukaryotes, eubacteria, and archaea to evaluate possible evolutionary scenarios for the origin of this pathway. Refined sequence alignments show that aconitase active site residues are highly conserved in each of the enzymes, and intervening sequence sites are quite dissimilar. This pattern suggests strong purifying selection has acted to preserve the aconitase active site residues for a common catalytic mechanism; numerous other substitutions occur due to adaptive evolution or simply lack of functional constraint. We hypothesize that the similarities are the remnants of an ancestral gene duplication, which may not have occurred within the fungal lineage. Maximum likelihood, neighbor joining, and maximum parsimony phylogenetic comparisons show that the α-aminoadipate pathway enzyme is an outgroup to all aconitase family proteins for which sequence is currently available. Received: 7 October 1997