Phylogenetic analysis of tmRNA genes within a bacterial subgroup reveals a specific structural signature

Bacterial tmRNA mediates a trans-translation reaction, which permits the recycling of stalled ribosomes and probably also contributes to the regulated expression of a subset of genes. Its action results in the addition of a small number of C-terminal amino acids to protein whose synthesis had stalled and these constitute a proteolytic recognition tag for the degradation of these incompletely synthesized proteins. Previous work has identified pseudoknots and stem–loops that are widely conserved in divergent bacteria. In the present work an alignment of tmRNA gene sequences within 13 β-proteobacteria reveals an additional sub-structure specific for this bacterial group. This sub-structure is in pseudoknot Pk2, and consists of one to two additional stem–loop(s) capped by stable GNRA tetraloop(s). Three-dimensional models of tmRNA pseudoknot 2 (Pk2) containing various topological versions of the additional sub-structure suggest that the sub-structures likely point away from the core of the RNA, containing both the tRNA and the mRNA domains. A putative tertiary interaction has also been identified.

Multiple maternal origins and weak phylogeographic structure in domestic goats

Domestic animals have played a key role in human history. Despite their importance, however, the origins of most domestic species remain poorly understood. We assessed the phylogenetic history and population structure of domestic goats by sequencing a hypervariable segment (481 bp) of the mtDNA control region from 406 goats representing 88 breeds distributed across the Old World. Phylogeographic analysis revealed three highly divergent goat lineages (estimated divergence >200,000 years ago), with one lineage occurring only in eastern and southern Asia. A remarkably similar pattern exists in cattle, sheep, and pigs. These results, combined with recent archaeological findings, suggest that goats and other farm animals have multiple maternal origins with a possible center of origin in Asia, as well as in the Fertile Crescent. The pattern of goat mtDNA diversity suggests that all three lineages have undergone population expansions, but that the expansion was relatively recent for two of the lineages (including the Asian lineage). Goat populations are surprisingly less genetically structured than cattle populations. In goats only ≈10% of the mtDNA variation is partitioned among continents. In cattle the amount is ≥50%. This weak structuring suggests extensive intercontinental transportation of goats and has intriguing implications about the importance of goats in historical human migrations and commerce.

Effects of passage history and sampling bias on phylogenetic reconstruction of human influenza A evolution

In this paper we determine the extent to which host-mediated mutations and a known sampling bias affect evolutionary studies of human influenza A. Previous phylogenetic reconstruction of influenza A (H3N2) evolution using the hemagglutinin gene revealed an excess of nonsilent substitutions assigned to the terminal branches of the tree. We investigate two hypotheses to explain this observation. The first hypothesis is that the excess reflects mutations that were either not present or were at low frequency in the viral sample isolated from its human host, and that these mutations increased in frequency during passage of the virus in embryonated eggs. A set of 22 codons known to undergo such “host-mediated” mutations showed a significant excess of mutations assigned to branches attaching sequences from egg-cultured (as opposed to cell-cultured) isolates to the tree. Our second hypothesis is that the remaining excess results from sampling bias. Influenza surveillance is purposefully biased toward sequencing antigenically dissimilar strains in an effort to identify new variants that may signal the need to update the vaccine. This bias produces an excess of mutations assigned to terminal branches simply because an isolate with no close relatives is by definition attached to the tree by a relatively long branch. Simulations show that the magnitude of excess mutations we observed in the hemagglutinin tree is consistent with expectations based on our sampling protocol. Sampling bias does not affect inferences about evolution drawn from phylogenetic analyses. However, if possible, the excess caused by host-mediated mutations should be removed from studies of the evolution of influenza viruses as they replicate in their human hosts.

From symmetry to asymmetry: Phylogenetic patterns of asymmetry variation in animals and their evolutionary significance

Phylogenetic analyses of asymmetry variation offer a powerful tool for exploring the interplay between ontogeny and evolution because (i) conspicuous asymmetries exist in many higher metazoans with widely varying modes of development, (ii) patterns of bilateral variation within species may identify genetically and environmentally triggered asymmetries, and (iii) asymmetries arising at different times during development may be more sensitive to internal cytoplasmic inhomogeneities compared to external environmental stimuli. Using four broadly comparable asymmetry states (symmetry, antisymmetry, dextral, and sinistral), and two stages at which asymmetry appears developmentally (larval and postlarval), I evaluated relations between ontogenetic and phylogenetic patterns of asymmetry variation. Among 140 inferred phylogenetic transitions between asymmetry states, recorded from 11 classes in five phyla, directional asymmetry (dextral or sinistral) evolved directly from symmetrical ancestors proportionally more frequently among larval asymmetries. In contrast, antisymmetry, either as an end state or as a transitional stage preceding directional asymmetry, was confined primarily to postlarval asymmetries. The ontogenetic origin of asymmetry thus significantly influences its subsequent evolution. Furthermore, because antisymmetry typically signals an environmentally triggered asymmetry, the phylogenetic transition from antisymmetry to directional asymmetry suggests that many cases of laterally fixed asymmetries evolved via genetic assimilation.

Retroposon analysis of major cetacean lineages: The monophyly of toothed whales and the paraphyly of river dolphins

SINE (short interspersed element) insertion analysis elucidates contentious aspects in the phylogeny of toothed whales and dolphins (Odontoceti), especially river dolphins. Here, we characterize 25 informative SINEs inserted into unique genomic loci during evolution of odontocetes to construct a cladogram, and determine a total of 2.8 kb per taxon of the flanking sequences of these SINE loci to estimate divergence times among lineages. We demonstrate that: (i) Odontocetes are monophyletic; (ii) Ganges River dolphins, beaked whales, and ocean dolphins diverged (in this order) after sperm whales; (iii) three other river dolphin taxa, namely the Amazon, La Plata, and Yangtze river dolphins, form a monophyletic group with Yangtze River dolphins being the most basal; and (iv) the rapid radiation of extant cetacean lineages occurred some 28–33 million years B.P., in strong accord with the fossil record. The combination of SINE and flanking sequence analysis suggests a topology and set of divergence times for odontocete relationships, offering alternative explanations for several long-standing problems in cetacean evolution.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.

Extreme morphological and ecological homoplasy in tropical salamanders

Fossorial salamanders typically have elongate and attenuated heads and bodies, diminutive limbs, hands and feet, and extremely elongate tails. Batrachoseps from California, Lineatriton from eastern México, and Oedipina from southern México to Ecuador, all members of the family Plethodontidae, tribe Bolitoglossini, resemble one another in external morphology, which has evolved independently. Whereas Oedipina and Batrachoseps are elongate because there are more trunk vertebrae, a widespread homoplasy (parallelism) in salamanders, the genus Lineatriton is unique in having evolved convergently by an alternate “giraffe-neck” developmental program. Lineatriton has the same number of trunk vertebrae as related, nonelongated taxa, but individual trunk vertebrae are elongated. A robust phylogenetic hypothesis, based on sequences of three mtDNA genes, finds Lineatriton to be deeply nested within a clade characterized by generalized ecology and morphology. Lineatriton lineolus, the only currently recognized taxon in the genus, shows unanticipated genetic diversity. Surprisingly, geographically separated populations of L. lineolus are not monophyletic, but are sister taxa of different species of the morphologically generalized genus Pseudoeurycea. Lineatriton, long thought to be a unique monospecific lineage, is polyphyletic. Accordingly, the specialized morphology of Lineatriton displays homoplasy at two hierarchical levels: (i) with respect to other elongate lineages in the family (convergence), and (ii) within what is currently recognized as a single taxon (parallelism). These evolutionary events are of adaptive significance because to invade the lowland tropics salamanders must be either arboreal or fossorial; the repeated evolution of elongation and attenuation has led to multiple lowland invasions.

Comparison of circumsporozoite proteins from avian and mammalian malarias: biological and phylogenetic implications.

The circumsporozoite (CS) protein of malaria parasites (Plasmodium) covers the surface of sporozoites that invade hepatocytes in mammalian hosts and macrophages in avian hosts. CS genes have been characterized from many Plasmodium that infect mammals; two domains of the corresponding proteins, identified initially by their conservation (region I and region II), have been implicated in binding to hepatocytes. The CS gene from the avian parasite Plasmodium gallinaceum was characterized to compare these functional domains to those of mammalian Plasmodium and for the study of Plasmodium evolution. The P. gallinaceum protein has the characteristics of CS proteins, including a secretory signal sequence, central repeat region, regions of charged amino acids, and an anchor sequence. Comparison with CS signal sequences reveals four distinct groupings, with P. gallinaceum most closely related to the human malaria Plasmodium falciparum. The 5-amino acid sequence designated region I, which is identical in all mammalian CS and implicated in hepatocyte invasion, is different in the avian protein. The P. gallinaceum repeat region consists of 9-amino acid repeats with the consensus sequence QP(A/V)GGNGG(A/V). The conserved motif designated region II-plus, which is associated with targeting the invasion of liver cells, is also conserved in the avian protein. Phylogenetic analysis of the aligned Plasmodium CS sequences yields a tree with a topology similar to the one obtained using sequence data from the small subunit rRNA gene. The phylogeny using the CS gene supports the proposal that the human malaria P. falciparum is significantly more related to avian parasites than to other parasites infecting mammals, although the biology of sporozoite invasion is different between the avian and mammalian species.

Paternal and maternal DNA lineages reveal a bottleneck in the founding of the Finnish population.

An analysis of Y-chromosomal haplotypes in several European populations reveals an almost monomorphic pattern in the Finns, whereas Y-chromosomal diversity is significantly higher in other populations. Furthermore, analyses of nucleotide positions in the mitochondrial control region that evolve slowly show a decrease in genetic diversity in Finns. Thus, relatively few men and women have contributed the genetic lineages that today survive in the Finnish population. This is likely to have caused the so-called "Finnish disease heritage"-i.e., the occurrence of several genetic diseases in the Finnish population that are rare elsewhere. A preliminary analysis of the mitochondrial mutations that have accumulated subsequent to the bottleneck suggests that it occurred about 4000 years ago, presumably when populations using agriculture and animal husbandry arrived in Finland.

Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.

Phylogenetic analyses are increasingly used in attempts to clarify transmission patterns of human immunodeficiency virus type 1 (HIV-1), but there is a continuing discussion about their validity because convergent evolution and transmission of minor HIV variants may obscure epidemiological patterns. Here we have studied a unique HIV-1 transmission cluster consisting of nine infected individuals, for whom the time and direction of each virus transmission was exactly known. Most of the transmissions occurred between 1981 and 1983, and a total of 13 blood samples were obtained approximately 2-12 years later. The p17 gag and env V3 regions of the HIV-1 genome were directly sequenced from uncultured lymphocytes. A true phylogenetic tree was constructed based on the knowledge about when the transmissions had occurred and when the samples were obtained. This complex, known HIV-1 transmission history was compared with reconstructed molecular trees, which were calculated from the DNA sequences by several commonly used phylogenetic inference methods [Fitch-Margoliash, neighbor-joining, minimum-evolution, maximum-likelihood, maximum-parsimony, unweighted pair group method using arithmetic averages (UPGMA), and a Fitch-Margoliash method assuming a molecular clock (KITSCH)]. A majority of the reconstructed trees were good estimates of the true phylogeny; 12 of 13 taxa were correctly positioned in the most accurate trees. The choice of gene fragment was found to be more important than the choice of phylogenetic method and substitution model. However, methods that are sensitive to unequal rates of change performed more poorly (such as UPGMA and KITSCH, which assume a constant molecular clock). The rapidly evolving V3 fragment gave better reconstructions than p17, but a combined data set of both p17 and V3 performed best. The accuracy of the phylogenetic methods justifies their use in HIV-1 research and argues against convergent evolution and selective transmission of certain virus variants.

The root of the universal tree and the origin of eukaryotes based on elongation factor phylogeny.

The genes for the protein synthesis elongation factors Tu (EF-Tu) and G (EF-G) are the products of an ancient gene duplication, which appears to predate the divergence of all extant organismal lineages. Thus, it should be possible to root a universal phylogeny based on either protein using the second protein as an outgroup. This approach was originally taken independently with two separate gene duplication pairs, (i) the regulatory and catalytic subunits of the proton ATPases and (ii) the protein synthesis elongation factors EF-Tu and EF-G. Questions about the orthology of the ATPase genes have obscured the former results, and the elongation factor data have been criticized for inadequate taxonomic representation and alignment errors. We have expanded the latter analysis using a broad representation of taxa from all three domains of life. All phylogenetic methods used strongly place the root of the universal tree between two highly distinct groups, the archaeons/eukaryotes and the eubacteria. We also find that a combined data set of EF-Tu and EF-G sequences favors placement of the eukaryotes within the Archaea, as the sister group to the Crenarchaeota. This relationship is supported by bootstrap values of 60-89% with various distance and maximum likelihood methods, while unweighted parsimony gives 58% support for archaeal monophyly.