Phylogenetic position of Psalteriomonas lanterna deduced from the SSU rDNA sequence

The small subunit ribosomal DNA sequence (SSU rDNA) of the microaerophilic free-living amoeboflagellate Psalteriomonas lanterna has been sequenced and analyzed. The gene is 1,945 bp long and has a G + C content of 33.4%. Based upon ultrastructural studies, P. lanterna has been placed in the class Lyromonadea within the phylum Percolozoa Cavalier-Smith, 1991. However, based upon cytological characteristics, this microaerophilic free-living amoeboflagellate appears to be very primitive. It shares certain characteristics in common with some archezoans, i.e. it lacks mitochondria and dictyosomes but contains hydrogenosomes. Despite sharing these characteristics with the amitochondriate taxa, P. lanterna is not related to any of these taxa but instead to the Vahlkampfiidae. Therefore, we used primary sequence data and the secondary structure of the SSU rDNA gene to determine the placement P. lanterna in the phylogenetic tree. Our analyses showed that P. lanterna groups as a sister taxon to the Vahlkampfiidae but probably diverged from them quite early.

The AIDS Pandemic is New, but is HIV Not New?

The determinations made by Mindell, Shultz and Ewald regarding the ancestral host for immunodeficiency retroviruses, and their conclusion that monkeys acquired their infections as a result of a host-switch from humans, do not withstand rigorous scrutiny. Their hypothesis requires the complete uniformativeness of third position transitions and of gapped regions in the alignment. When all of the data are permitted to corroborate or refute relationships, optimizing hosts on the viral phylogeny renders either equivocal statements or an unequivocal simian ancestry. However, merely optimizing hosts as characters on the viral phylogeny is illogical. Not only does this treat hosts as dependent on the viruses (instead of the reverse) but it ignores 15 years of methodological developments specifically designed to answer questions regarding cospeciation or host-switching.

Deciphering posttranslational processing events in the pituitary of a neopterygian fish: cloning of a gar proopiomelanocortin cDNA

A cDNA that codes for the polypeptide hormone precursor proopiomelanocortin (POMC) was cloned and sequenced from a gar (Lepisosteus osseus) pituitary cDNA library. The gar POMC cDNA is 1237 bp and contains a 780-bp open reading frame. The deduced amino acid sequence for gar POMC is 259 amino acids in length. The general organization of gar POMC is very similar to that of other gnathostome POMC sequences. The beta-endorphin sequence had 91% sequence identity with sockeye A beta-endorphin and 71% sequence identity with Xenopus laevis beta-endorphin. Three melanocyte-stimulating hormone (MSH) core sequences [HFR(W)] were detected. The gar alpha-MSH sequence was identical to the alpha-MSH sequence in rat POMC. The gar beta-MSH sequence had 77% sequence identity with salmonid forms of beta-MSH and 53% sequence identity with tetrapod forms of beta-MSH. The gamma-MSH region of gar POMC only had 26% primary sequence identity with tetrapod gamma-MSH sequences. Gar gamma-MSH had an incomplete MSH core sequence (HRF), an apparent internal deletion of five amino acids, and lacked flanking paired basic amino acids essential for proteolytic cleavage. The apparent degenerate nature of gar gamma-MSH is discussed in light of the absence of this sequence in salmonid fish.

Primary structure of insulin from the african lungfish, Protopterus annectens

Among the extant Sarcopterygii, the interrelationship between the Dipnoi (lungfishes), Actinistia (coelacanths), and Tetrapoda (tetrapods) is controversial. Insulin has been purified from an extract of the pancreas of the African lungfish Protopterus annectens and its primary structure established as A-chain, Gly-Ile-Val-Glu-Gln-Cys-Cys-His-Lys-Pro10-Cys-Ser-Leu- Tyr -Glu-Leu-Glu-Asn-Tyr-Cys20-Asn-Val-Pro; and B-chain, Ala-Val-Leu-Asn-Gln-His-Leu-Cys-Gly-Ser10-His-Leu-Val- Glu- Ala-Leu-Tyr-Leu-Val-Cys20-Ala-Asp-Asn-Gly-Phe- Phe-Tyr-Lys-Pro-Ser30-Gly. Lungfish insulin contains unusual structural features, such as the dipeptide extension to the C-terminus of the A-chain and the substitution Arg --> Asn at position B-23 in the putative receptor binding region of insulin, which may be expected to influence appreciably its biological potency relative to mammalian insulins. Lungfish insulin also contains amino acid substitutions such as Gly --> Ala at position B-21, Glu --> Asp at position B-22, and a Lys --> Ser residue at position B-30, previously found in insulins from amphibia. This observation is consistent with paleontological data suggesting that lungfish and amphibia share a close phylogenetic relationship.

East Gondwana ancestry of the sunflower alliance of families

The sunflower alliance of families comprises nearly 10% of all flowering plant species and includes the largest of all plant families, the sunflower family Asteraceae, which has 23,000 species, and the bellflower family Campanulaceae. Both are worldwide in distribution, but the majority of their species occur in the northern hemisphere. Recently it has been shown that a number of small, woody families from the Australian-Southwest Pacific area also belong in this relationship. Here we add yet another such family and present phylogenetic, biogeographic, and chronological analyses elucidating the origin of this large group of plants. We show that the ancestral lineages are confined to Malesia, Australia, New Guinea, and New Zealand and that the sunflower and bellflower families represent phylogenetically derived lineages within a larger group with a Cretaceous and southern-hemisphere, presumably East Gondwana, ancestry. Their highly derived position in the flowering plant phylogeny makes this significant for understanding the evolution of flowering plants in general.

The phylogeny of the Caelifera (Insecta, Orthoptera) as deduced from mtrRNA gene sequences

Fragments of both mitochondrial ribosomal RNA genes of 32 caeliferan taxa (representing six of the seven superfamilies) and six outgroup Orthopteroids were sequenced. The combined alignment length was 630 bp after removal of all ambiguously aligned positions. Separation for the basal taxa was problematic and analysis using the LogDet transformation indicated that shared base composition biases were a confounding factor. The suborder Caelifera and all traditional caeliferan superfamilies except the Pamphagoidea are retrieved as monophyletic groups, though the Eumastacoidea lack significant bootstrap support. Of the traditional pamphagoid taxa, the Pamphagidae is embedded between classically acridid subfamilies, whereas Pyrgomorphidae is placed close to the Pneumoroidea. The morphological similarities of the Pyrgomorphidae and Pamphagidae may thus be homoplasic. A consensus tree based on five different methods of analysis indicated the following order: (Tridactyloidea, Tetrigoidea (Eumastacidae, Proscopiidae (Pneumoridae, Pyrgomorphidae (Acrididae + Pamphagidae)))).

Molecular characterization and phylogenetic analyses of a new, highly divergent simian T-cell lymphotropic virus type 1 (STLV-1marc1) in Macaca arctoides

A serological survey of a captive colony of Asian monkeys indicated that six Macaca arctoides had antibodies to human T-cell leukemia/lymphotropic viruses (HTLV). Over a 4-year interval, sera from these animals continued to exhibit a peculiar Western blot (WB) pattern resembling an HTLV-2 pattern (p24gag reactivity of equal or greater intensity than that of p19gag and a strong reactivity to recombinant gp21) but also exhibiting, in five of six cases, a reactivity against MTA-1, an HTLV-1 gp46 peptide. PCR experiments on DNA extracted from peripheral blood mononuclear cells using HTLV-1- or HTLV-2-specific long terminal repeat, gag, pol, env, and tax primers yielded negative results. However, highly conserved primers successfully amplified three different gene segments of env, tax, and env-tax. The results of comparative sequence analysis demonstrated that STLV-1marc1 was not closely related to any known STLV-1 strain, was the most divergent strain of the HTLV-1-STLV-1 group, and lacked the ATG initiation codons corresponding to the p12 and p13 proteins of HTLV-1. Phylogenetic analyses incorporating representative strains of all known HTLV-STLV clades consistently depicted STLV-1marc1 within the HTLV-1-STLV-1 type 1 lineage, but it probably diverged early, since its position is clearly different from all known viral strains of this group and it had a bootstrap resampling value of 100%. Genetic distance estimates between STLV-1marc1 and all other type 1 viruses were of the same order of magnitude as those between STLV-2PanP and all other type 2 viruses. In light of the recent demonstration of interspecies transmission of some STLV-1 strains, our results suggest the existence in Asia of HTLV-1 strains related to this new divergent STLV-1marc1 strain, which may be derived from a common ancestor early in the evolution of the type 1 viruses and could be therefore considered a prototype of a new HTLV-STLV clade.

Phylogeny of the side-blotched lizards (Phrynosomatidae:Uta) based on mtDNA sequences: support for midpeninsular seaway in Baja California

We investigated the phylogenetic relationships of western North American Side-blotched lizards, genus Uta, using mitochondrial DNA partial gene sequences from cytochorme b and ATPase 6. Uta stejnegeri appeared basal in our tree followed by U. palmeri. Uta stansburiana from the islands of Angel de la Guarda, Mejia, and Raza formed the next clade, followed by U. antiqua and other populations of U. stansburiana. These relationships suggest that the populations on the Angel de la Guarda island block should be recognized as a distinct species. Remaining populations of U. stansburiana formed two clades, corresponding to the northern and southern Baja California peninsula. Uta stellata and U. squamata occurred within the northern and southern clades, respectively. The discontinuity requires a long-lasting isolation, the only reasonable explanation being a former midpeninsular seaway. Correlations between our cladogram and magnetic anomalies in the Gulf of California date formation of the seaway at 1 million years ago.

Evolution of the entire arthropod Hox gene set predated the origin and radiation of the onychophoran/arthropod clade

BACKGROUND

Dramatic changes in body size and pattern occurred during the radiation of many taxa in the Cambrian, and these changes are best documented for the arthropods. The sudden appearance of such diverse body plans raises the fundamental question of when the genes and the developmental control systems that regulate these designs evolved. As Hox genes regulate arthropod body patterns, the evolution of these genes may have played a role in the origin and diversification of the arthropod body plan from a homonomous ancestor. To trace the origin of arthropod Hox genes, we examined their distribution in a myriapod and in the Onychophora, a sister group to the arthropods.

RESULTS

Despite the limited segmental diversity within myriapods and Onychophora, all insect Hox genes are present in both taxa, including the trunk Hox genes Ultrabithorax and abdominal-A as well as an ortholog of the fushi tarazu gene. Comparative analysis of Hox gene deployment revealed that the anterior boundary of expression of trunk Hox genes has shifted dramatically along the anteroposterior axis between Onychophora and different arthropod classes. Furthermore, we found that repression of expression of the Hox target gene Distal-less is unique to the insect lineage.

CONCLUSIONS

A complete arthropod Hox gene family existed in the ancestor of the onychophoran/arthropod clade. No new Hox genes were therefore required to catalyze the arthropod radiation; instead, arthropod body-plan diversity arose through changes in the regulation of Hox genes and their downstream targets.

The biological relevance of testing for perfect symmetry

No abstract

Hadamard conjugations and modeling sequence evolution with unequal rates across sites

This paper considers the many different distributions that may approximate the distribution of site rates in DNA sequences and shows how the Hadamard conjugation may be modified to take these into account. This is done for both 2-state and 4-state data. Distributions which give simple closed forms include the gamma (gamma) distribution, the inverse Gaussian distribution (which is similar to the lognormal), and a mixture of either of these with a proportion of sites which cannot change (invariant sites). It is seen that the tail of a distribution can have major effects upon the coefficient of variation of site rates. Because the Hadamard conjugation can be used to either correct data or predict the data given the model (i.e., the likelihood of site patterns), light is shed on properties of maximum likelihood tree selection with unequal site rates. Analysis of rRNA shows how unequal rates across sites can change the optimal tree. Maximum likelihood analysis also shows that distinct distributions fit each data set, with the gamma often not being the best. Analyzing both these data and a long stretch of primate mtDNA reveals evidence of many "hidden" multiple substitutions, while signals not corresponding to the preferred biological tree generally decrease an unequal rates are allowed for. Last, we discuss the expected behavior of sequences evolving by models where stabilizing selection alone explains unequal site rates. Such models do not explain "synapomorphies" or informative changes in ancient molecules, because while stabilizing selection can vastly decrease change at a site, it will also vastly accelerate back-substitution (leaving only a covarion model to explain old synapomorphies). When and why models allowing a continuous distribution of site rates (e.g., gamma) will approximate covarion evolution requires further study.

Genes, peoples, and languages

The genetic history of a group of populations is usually analyzed by reconstructing a tree of their origins. Reliability of the reconstruction depends on the validity of the hypothesis that genetic differentiation of the populations is mostly due to population fissions followed by independent evolution. If necessary, adjustment for major population admixtures can be made. Dating the fissions requires comparisons with paleoanthropological and paleontological dates, which are few and uncertain. A method of absolute genetic dating recently introduced uses mutation rates as molecular clocks; it was applied to human evolution using microsatellites, which have a sufficiently high mutation rate. Results are comparable with those of other methods and agree with a recent expansion of modern humans from Africa. An alternative method of analysis, useful when there is adequate geographic coverage of regions, is the geographic study of frequencies of alleles or haplotypes. As in the case of trees, it is necessary to summarize data from many loci for conclusions to be acceptable. Results must be independent from the loci used. Multivariate analyses like principal components or multidimensional scaling reveal a number of hidden patterns and evaluate their relative importance. Most patterns found in the analysis of human living populations are likely to be consequences of demographic expansions, determined by technological developments affecting food availability, transportation, or military power. During such expansions, both genes and languages are spread to potentially vast areas. In principle, this tends to create a correlation between the respective evolutionary trees. The correlation is usually positive and often remarkably high. It can be decreased or hidden by phenomena of language replacement and also of gene replacement, usually partial, due to gene flow.

Genetics and the origin of bird species

External (environmental) factors affecting the speciation of birds are better known than the internal (genetic) factors. The opposite is true for several groups of invertebrates, Drosophila being the outstanding example. Ideas about the genetics of speciation in general trace back to Dobzhansky who worked with Drosophila. These ideas are an insufficient guide for reconstructing speciation in birds for two main reasons. First, speciation in birds proceeds with the evolution of behavioral barriers to interbreeding; postmating isolation usually evolves much later, perhaps after gene exchange has all but ceased. As a consequence of the slow evolution of postmating isolating factors the scope for reinforcement of premating isolation is small, whereas the opportunity for introgressive hybridization to influence the evolution of diverging species is large. Second, premating isolation may arise from nongenetic, cultural causes; isolation may be affected partly by song, a trait that is culturally inherited through an imprinting-like process in many, but not all, groups of birds. Thus the genetic basis to the origin of bird species is to be sought in the inheritance of adult traits that are subject to natural and sexual selection. Some of the factors involved in premating isolation (plumage, morphology, and behavior) are under single-gene control, most are under polygenic control. The genetic basis of the origin of postmating isolating factors affecting the early development of embryos (viability) and reproductive physiology (sterility) is almost completely unknown. Bird speciation is facilitated by small population size, involves few genetic changes, and occurs relatively rapidly.

Molecular evidence that echiurans and pogonophorans are derived annelids

The Annelida, which includes the polychaetes and the clitellates, has long held the taxonomic rank of phylum. The unsegmented, mud-dwelling echiuran spoon worms and the gutless, deep-sea pogonophoran tube worms (including vestimentiferans) share several embryological and morphological features with annelids, but each group also has been considered as a separate metazoan phylum based on the unique characters each group displays. Phylogenetic analyses of DNA sequences from the nuclear gene elongation factor-1alpha place echiurans and pogonophorans within the Annelida. This result, indicating the derived loss of segmentation in echiurans, has profound implications for our understanding of the evolution of metazoan body plans and challenges the traditional view of the phylum-level diversity and evolutionary relationships of protostome worms.

Phylogenetics and the origin of species

A recent criticism that the biological species concept (BSC) unduly neglects phylogeny is examined under a novel modification of coalescent theory that considers multiple, sex-defined genealogical pathways through sexual organismal pedigrees. A competing phylogenetic species concept (PSC) also is evaluated from this vantage. Two analytical approaches are employed to capture the composite phylogenetic information contained within the braided assemblages of hereditary pathways of a pedigree: (i) consensus phylogenetic trees across allelic transmission routes and (ii) composite phenograms from quantitative values of organismal coancestry. Outcomes from both approaches demonstrate that the supposed sharp distinction between biological and phylogenetic species concepts is illusory. Historical descent and reproductive ties are related aspects of phylogeny and jointly illuminate biotic discontinuity.

Incipient species formation in salamanders of the Ensatina complex

The Ensatina eschscholtzii complex of plethodontid salamanders, a well-known "ring species," is thought to illustrate stages in the speciation process. Early research, based on morphology and coloration, has been extended by the incorporation of studies of protein variation and mitochondrial DNA sequences. The new data show that the complex includes a number of geographically and genetically distinct components that are at or near the species level. The complex is old and apparently has undergone instances of range contraction, isolation, differentiation, and then expansion and secondary contact. While the hypothesis that speciation is retarded by gene flow around the ring is not supported by molecular data, the general biogeographical hypothesis is supported. There is evidence of a north to south range expansion along two axes, with secondary contact and completion of the ring in southern California. Current research targets regions once thought to show primary intergradation, but which molecular markers reveal to be zones of secondary contact. Here emphasis is on the subspecies E. e. xanthoptica, which is involved in four distinct secondary contacts in central California. There is evidence of renewed genetic interactions upon recontact, with greater genetic differentiation within xanthoptica than between it and some of the interacting populations. The complex presents a full array of intermediate conditions between well-marked species and geographically variable populations. Geographically differentiated segments represent a diversity of depths of time of isolation and admixture, reflecting the complicated geomorphological history of California. Ensatina illustrates the continuing difficulty in making taxonomic assignments in complexes studied during species formation.

Cleavage patterns and the topology of the metazoan tree of life

Several major alliances of metazoan phyla have been identified by small subunit rRNA sequence comparisons. It is possible to arrange the phyla to produce a parsimonious distribution of cleavage types, requiring only one change from a radial ancestral condition to spiral cleavage and one other to "idiosyncratic" cleavage; this arrangement is consistent with most of the recent molecular phylogenies. The cleavage shifts are correlated with changes in many of the features that once were used to distinguish Protostomia and Deuterostomia. It is hypothesized that changes in cleavage direction are causally associated with changes in blastomere fates and thus that cleavage type correlates with such features as the identity of mesoderm founder cells, which in turn can constrain the mode of origination of the eucelom. Cleavage changes may also affect the timing of cell fate specification. In a tree that emphasizes cleavage parsimony, radial cleavage, regulative development, and enterocely are ancestral within the Bilateria, and spiral or idiosyncratic cleavages, mosaic development, and schizocely are associated with a change in cleavage direction. Deuterostomy is presumably ancestral and is correlated with radial cleavage for this reason, rather than mechanistically.

Path costs in evolutionary tree reconstruction

This paper describes a dynamic programming algorithm to solve a family of problems in the reconstruction of evolutionary trees from protein sequence data, that of constructing "minimal" colorings. This dynamic programming formulation can be modified to efficiently enumerate the number of minimal colorings and thereby be used to calculate the average cost of any given edge, where the average is taken over the entire set of minimal colorings. An extension of our dynamic programming formulation allows for the calculation of average path costs amongst all minimal colorings. our results resolve questions raised in (Hendy and penny, 1987); in particular, we develop polynomial time procedures to find the minimum, maximum, and average (expected) cost of an edge, and more generally of a path, for a minimal coloring. Our algorithm is distinguished in its flexibility to address further distribution and statistical questions relating to the minimal colorings. Furthermore, the more general concept of calculating statistics describing the set of optimal solutions may be of interest in other combinatorial problems.

On a Mirkin-Muchnik-Smith conjecture for comparing molecular phylogenies

A conjecture of Mirkin, Muchnik, and Smith is answered affirmatively which connects the inconsistency function, a biologically meaningful similarity/dissimilarity measure for a gene tree and a species tree, to the mutation cost function, a combinatorial measure based on the mapping of trees. A linear-time algorithm for computing the mutation cost function is also derived from the conjecture.

Evolution of the mitochondrial DNA control region in the mbuna (Cichlidae) species flock of Lake Malawi, East Africa

Considerable controversy has surrounded the application of mitochondrial DNA data to reconstruction of evolutionary relationships among the endemic cichlids of Lake Malawi. Central to this debate has been the issue of whether lineage sorting is complete, and thus whether these data actually reflect species phylogeny, or simply gene genealogy. Review of all mtDNA control region sequences available for members of one monophyletic subset of this species flock, the Malawi rockfishes, or mbuna, strongly indicates that lineage sorting is incomplete: Character-based analyses of these sequences reconstruct gene, not species, interrelationships. Analysis of the pattern of nucleotide substitutions differentiating these mtDNA alleles suggests that pyrimidine residues undergo transition substitutions more often than do purines. Estimation of the magnitude of derived sequence differentiation in light of the reconstructed gene genealogy suggests that the mbuna may be of considerably more recent vintage than previous molecular characterizations have indicated.

Preliminary results on the genetic variability of mitochondrial DNA in the signal crayfish, Pacifastacus leniusculus Dana

This study presents a revised method for the extraction of mitochondrial genome (mtDNA) from crayfish species of aquacultural interest, Pacifastacus leniusculus. The mean size of the mitochondrial genome is approximately 18,000 bp. Restriction fragment length polymorphism (RFLP) analysis performed on the mtDNA has revealed extensive genetic variability within this species. The estimated percentage of nucleotide sequence divergence among the four haplotypes found for P. leniusculus ranges from 1.16 to 3.99%. These data suggest that RFLP analysis of mtDNA may provide greater resolution than protein electrophoresis for population identification among signal crayfish populations. This marker offers new avenues for aquaculture and for understanding the genetic structure and evolutionary history of crayfish populations.

DNA phylogeny of the marsupial wolf resolved

The phylogenetic position of the recently extinct marsupial 'wolf', or thylacine (Thylacinus cynocephalus), has been a source of contention in mammalian systematics for nearly a century. Thylacines were endemic to Australasia, but possessed striking anatomical similarities to Oligo-Miocene borhyaenid marsupials of South America. At issue has been whether these features are indicative of common ancestry or convergent adaptation to carnivory. Recent morphological studies have supported both conclusions. Although current marsupial classifications group thylacines with Australian dasyuromorphians, this putative clade is characterized by mostly primitive morphological features. Attempts to determine thylacine affinities with ancient protein and DNA analyses have supported, but not resolved, a dasyuromorphian placement. We report 1546 bp of mitochondrial DNA sequence (from cytochrome b and 12S rRNA genes) and 841 bp of nuclear protamine gene sequence from the thylacine and representatives of all or most other marsupial orders. Phylogenetic analysis of these sequences shows unambiguously that thylacines are members of Dasyuromorphia, and suggests a late Oligocene or very early Miocene divergence of familial lineages.

Molecular systematics of Middle American cichlid fishes and the evolution of trophic-types in 'Cichlasoma (Amphilophus)' and 'C. (Thorichthys)'

The majority of Middle American cichlids are placed in the informal assemblage 'Cichlasoma.' The group is divided into eight sections which appear to be based primarily on trophic morphology. Although several members of 'Cichlasoma' have been used in ecomorphological, behavioral, and biogeographic studies, no phylogenetic hypotheses for the group exist. In an attempt to develop a better understanding of the phylogenetic relationships of 'cichlasomine' cichlids, we examined the evolution of the trophic specialization, substratum-sifting, in two sections, 'Cichlasoma (Thorichthys)' and 'C. (Amphilophus),' to determine whether the trait reflects common ancestry. We sequenced the complete mitochondrial cytochrome b gene for 19 cichlids representing six sections of 'Cichlasoma,' and representatives of other Neotropical Cichlidae. Additional cichlid, and noncichlid outgroup sequences were included for a total of 22 taxa. The molecular phylogeny supports the recognition of the section 'C. (Thoricthys)' as a natural group, and we place those cichlids in the genus Thorichthys. The phylogeny also depicts 'C. (Amphilophus)' as paraphyletic, with substratum-sifters and generalized predators forming separate nonsister clades. We recommend that the substratum-sifting clade of the section 'C. (Amphilophus)' be placed in the resurrected genus Astatheros. The generalized predator clade of 'C. (Amphilophus)' contains only two species, 'C. (A.) citrinellum' and 'C. (A.) labiatum,' which we place in the genus Amphilophus. The phylogenetic hypotheses generated indicate that the substratum-sifting genera Thorichthys and Amphilophus do not share a common ancestor. Reconstruction of the evolution of substratum-sifting is equivocal, requiring either the independent evolution of the trait on two separate occasions or its presence in a more inclusive clade and subsequent loss in nonsubstratum sifting species.

Phylogenetic relationships of cottontails (Sylvilagus, Lagomorpha): congruence of 12S rDNA and cytogenetic data

The genus Sylvilagus, which comprises the New World cottontail rabbits, contains several commercially important as well as endangered (or threatened) species. Understanding the evolution of this group is pertinent to their management and conservation. The purpose of this study was to examine the evolutionary history of the cottontails using sequence data from the mitochondrial 12S rRNA gene. The 12S data provide a robust phylogeny which was supported under a variety of phylogenetic approaches and transition/transversion (Ti/Tv) weighting schemes. Stem and loop regions of the gene were analyzed separately and two different methods of estimating Ti/Tv ratios were employed. The phylogeny obtained was consistent with available cytogenetic information. The 12S data indicate that separate generic status for the pygmy rabbit, Brachylagus idahoensis, is warranted based on its phylogenetic position and sequence divergence values. Additionally, the taxa which are geographically adjacent are also phylogenetically closely related; for example, the marsh rabbit, S. palustris, and the swamp rabbit, S. aquaticus, are sister taxa, as are the mountain cottontail, S. nuttallii, and desert cottontail, S. audubonii. This finding suggests that recent vicariance events might explain the diversification of several cottontail lineages.

Correlation of functional domains and rates of nucleotide substitution in cytochrome b

Distinguishing noise from signal presents a problem when DNA sequences are used for phylogeny reconstruction. Multiple substitutions at sites are a primary cause of noise and this is compounded by variation in substitution rates among sites. For protein-coding genes, one method used to determine if data are noisy is to assess levels of saturation of substitutions by codon position. However, this procedure may not be a fine enough filter for assessing noise. Variation in substitution rates may also be caused by constraints on change imposed by the function of the protein product. Using a structural model of the cytochrome b protein as a template, I divided cyt b sequence data for species within the avian family Falconidae (falcons and caracaras) into three functional domains. Saturation of substitutions of sequences within these regions was assessed graphically. This qualitative determination of saturation was then used to differentially weight phylogenetic analysis, resulting in an hypothesis congruent with existing cladistic analyses and traditional morphology. These results demonstrate that saturation of substitutions is correlated with functional regions of cytochrome b and that using this information improves phylogenetic inference.

Consistency, characters, and the likelihood of correct phylogenetic inference

Computer simulations of character-state evolution in 8, 16, 32, and 64 ingroup taxa with a known set of relationships demonstrate that the maximum probability of correct phylogenetic inference increases with the number of variable (or informative) characters and their consistency index and decreases with the number of taxa, when the consistency index has been standardized to eliminate its dependence on the number of taxa. Equations for the probability of correct phylogenetic inference and for the standardized consistency indices (including or excluding autapomorphies) are derived. Given that actual studies based on DNA restriction sites and sequences generate more characters with a higher level of consistency than comparable studies based on morphology, calculations suggest that such molecular studies may often provide a more precise guide to phylogenetic relationships.

A survey of homeobox genes in Chaetopterus variopedatus and analysis of polychaete homeodomains

A survey of genomic DNA from the polychaete Chaetopterus variopedatus was conducted using the polymerase chain reaction. Twelve unique homeobox-containing gene fragments were recovered. Phylogenetic analysis indicates that seven of the fragments are from genes belonging to Hox homeobox classes. Other fragments show orthology with Xlox, caudal, and Prh homeobox classes, with two fragments not definitely assignable to a homeobox class by our analysis. Orthology with gene sequences reported for the polychaete Ctenodrilus serratus, by Dick and Buss (1994), was calculated and indicated that at least eight of the C. variopedatus fragments are homologous to these previously reported sequences. Tabulation of the Hox gene relationships suggest that polychaetes have representative genes of each of the Hox cognate groups except Abd-B. This conclusion further suggests that the Hox cluster in the basal protostome ancestor had already undergone the gene duplications leading to the complete complement of homeotic genes known in Drosophila, with the possible loss of Abd-B in the polychaete lineage.

A molecular phylogeny of the Aphidiinae (Hymenoptera: Braconidae)

Phylogenetic relationships within the Aphidiinae, and between this and other subfamilies of Braconidae (Hymenoptera), were investigated using sequence data from three genes: elongation factor-1alpha, cytochrome b, and the second expansion segment of the 28S ribosomal subunit. Variation in both protein-coding genes was characterized by a high level of homoplasy, but analysis of the expansion segment--robust over a range of alignment methods and parameters-resolved some of the older divergences. Parsimony analysis of the combined data suggests the following tribal relationships: (Ephedrini + (Praini + (Aphidiini + Trioxini))). In addition, the cyclostome subfamilies were found to form a clade separate from the Aphidiinae, but relationships between the Aphidiinae and the noncyclostome braconids could not be resolved. The inferred phylogeny also supported a secondary loss of internal pupation within the Praini and a polyphyletic origin of endoparasitism within the Braconidae.

Phylogenetic relationships of the liverworts (Hepaticae), a basal embryophyte lineage, inferred from nucleotide sequence data of the chloroplast gene rbcL

Sequence data from the chloroplast-encoded gene rbcL were obtained for 24 liverworts, a basal group of embryophytes. Maximum likelihood and parsimony analyses of these data, along with data from other major green plant lineages, confirm hypotheses based on morphological data, such as the paraphyly of bryophytes, and the basal position of liverworts. Molecular data corroborate the deep separation between the complex thalloid and leafy/simple thalloid liverworts implied by morphological data, but the monophyly of liverworts could not be rejected. The effects of accounting for site-to-site rate heterogeneity in these data were examined using maximum likelihood methods. Comparison of trees obtained with and without rate heterogeneity showed that simply allowing for heterogeneity had a greater improvement on likelihood score than optimization of transition/transversion bias. Incorporation of site-to-site rate heterogeneity in the larger analysis, however, did not necessarily change which topology was favored. Properties of rbcL sequences from the two liverwort groups were compared. Significantly different substitution rates were found between leafy/simple thalloid and complex thalloid liverwort taxa, with rates of rbcL sequence evolution in leafy/simple thalloid taxa being higher and more indicative of those of vascular plants, and with those of complex thalloid taxa (such as Marchantia) being slower. Codon usage in rbcL in complex thalloid liverworts was biased toward NNU and NNA, compared to the leafy/simple thalloid liverworts. Although base composition and relative substitution rates differed between the two groups, no significant differences were detected within each of the two groups of liverworts. The signal present in first and second codon sites versus third codon sites was compared. While the third codon positions in rbcL across this taxon sampling are highly variable (with only 15 constant sites of 439), the trees obtained were in general agreement with trees from the entire data set and with trees obtained from independent sources of data. The presence of signal in third codon positions across greater than 400 MY of plant evolution means that definitions of saturation based on pair-wise comparisons of sequences inadequately assess phylogenetic signal.

A cladistic analysis of mitochondrial ribosomal DNA from the Bovidae

There is a huge data base of genetic information for the domestic artiodactyl species Bos taurus (cow), Ovis aries (sheep), and Capra hircus (goat). However, the phylogenetic relationships of these economically critical taxa and their close relatives, family Bovidae, remain for the most part unresolved. In this report, we aligned new mitochondrial (mt) 12S and 16S ribosomal (r) DNA sequences from 26 bovid taxa with published sequences. Phylogenetic analyses of the more than 64 kilobases of mt rDNA from 57 taxa support a basal division in the Bovidae that separates Bos and its close relatives from Capra, Ovis, and their kin. As suggested by previous molecular and morphological studies, "antelopes" are a paraphyletic assemblage. Caprinae (sheep, goats, goat antelopes, and musk oxen) groups consistently with hippotragine and alcelaphine antelopes, while Bovini (cattle and buffaloes) clusters with tragelaphine and boselaphine antelopes. The traditional tribal subdivisions of Bovidae are supported in most cases, but there are exceptions within Caprinae and Antilopinae (gazelles and close relatives). The rDNA data consistently place the enigmatic genera Pelea, Pantholops, and Saiga, but the origin of Aepyceros, the impala, remains obscure. Combined phylogenetic analyses of the rDNA data with the skeletal characters of Gentry (1992) were used to assess the stability of the molecular results.

Evidence for a clade of nematodes, arthropods and other moulting animals

The arthropods constitute the most diverse animal group, but, despite their rich fossil record and a century of study, their phylogenetic relationships remain unclear. Taxa previously proposed to be sister groups to the arthropods include Annelida, Onychophora, Tardigrada and others, but hypotheses of phylogenetic relationships have been conflicting. For example, onychophorans, like arthropods, moult periodically, have an arthropod arrangement of haemocoel, and have been related to arthropods in morphological and mitochondrial DNA sequence analyses. Like annelids, they possess segmental nephridia and muscles that are a combination of smooth and obliquely striated fibres. Our phylogenetic analysis of 18S ribosomal DNA sequences indicates a close relationship between arthropods, nematodes and all other moulting phyla. The results suggest that ecdysis (moulting) arose once and support the idea of a new clade, Ecdysozoa, containing moulting animals: arthropods, tardigrades, onychophorans, nematodes, nematomorphs, kinorhynchs and priapulids. No support is found for a clade of segmented animals, the Articulata, uniting annelids with arthropods. The hypothesis that nematodes are related to arthropods has important implications for developmental genetic studies using as model systems the nematode Caenorhabditis elegans and the arthropod Drosophila melanogaster, which are generally held to be phylogenetically distant from each other.

Old World fruitbat phylogeny: evidence for convergent evolution and an endemic African clade

Knud Andersen (1912, Catalogue of the Chiroptera in the Collections of the British Museum: I. Megachiroptera, British Museum of Natural History, London) divided Old World fruitbats (family Pteropodidae) into the rousettine, cynopterine, epomophorine, eonycterine, and notopterine sections. The latter two sections comprise the subfamily Macroglossinae; members of this subfamily exhibit specializations for nectarivory (e.g., elongated, protrusible, brushy tongues) and cluster together in cladistic analyses based on anatomical characters. Other evidence, including single-copy DNA hybridization, suggests that macroglossines are either paraphyletic or polyphyletic; this implies that adaptations for pollen and nectar feeding evolved independently in different macroglossine lineages or were lost in nonmacroglossines after evolving in a more basal common ancestor. Hybridization data also contradict Andersen's phylogeny in providing support for an endemic African clade that includes representatives of three of Andersen's sections. Here, we present complete mitochondrial 12S rRNA and valine tRNA gene sequences for 20 pteropodids, including representatives of all of Andersen's sections, and examine the aforementioned controversies. Maximum likelihood, minimum evolution, and maximum parsimony analyses all contradict macroglossine monophyly and provide support for an African clade that associates Megaloglossus and Lissonycteris and those two with Epomophorus. In conjunction with the DNA hybridization results, there are now independent lines of molecular evidence suggesting: (i) convergent evolution of specializations for nectarivory, at least in Megaloglossus versus other macroglossines, and (ii) a previously unrecognized clade of endemic Africa taxa. Estimates of divergence time based on 12S rRNA and DNA hybridization data are also in good agreement and suggest that extant fruitbats trace back to a common ancestor 25 million to 36 million years ago.

Reassessing embryogenesis in the Ctenophora: the inductive role of e1 micromeres in organizing ctene row formation in the ‘mosaic’ embryo, Mnemiopsis leidyi

Ctenophores are a phylum of diploblastic marine animals displaying biradial symmetry organized along an oral-aboral axis. One of the apomorphic sets of adult structures in ctenophores are the eight external comb rows, which run along the oral-aboral axis. Comb rows consist of serial arrays of individual comb plates of cilia, which beat in a coordinated fashion for locomotory behavior. Classical cell lineage experiments using chalk particles indicated that comb rows are derived exclusively from the four e1 micromeres at the 16-cell stage. This conclusion was also supported by the fact that no ctene rows (or their underlying endodermal canals) form when all four e1 micromeres were deleted. We have used intracellular diI cell lineage tracing to determine that, in addition to e1 micromeres, the four m1 micromeres also make significant contributions to the ctene rows. Thus, e1 micromere derivatives not only generate comb plates but are required for ctene row formation by m1 derivatives. These results demonstrate that inductive interactions are an important component of early development in ctenophores and indicate that e1 micromeres influence the development of adjacent cell lineages (both m1 and endodermal lineages) during ctenophore embryogenesis. In addition, intracellular labeling has revealed that there are subtle variations in the composition of clones derived from identified embryonic blastomeres. Together these findings reveal a picture of ctenophore embryogenesis, which is in marked contrast to the former rigid ‘mosaic’ reputation of ctenophore development, and invite speculation as to the role of the cleavage program in embryonic patterning in the lower Metazoa.

The origin and evolution of animal appendages

Animals have evolved diverse appendages adapted for locomotion, feeding and other functions. The genetics underlying appendage formation are best understood in insects and vertebrates. The expression of the Distal-less (Dll) homeoprotein during arthropod limb outgrowth and of Dll orthologs (Dlx) in fish fin and tetrapod limb buds led us to examine whether expression of this regulatory gene may be a general feature of appendage formation in protostomes and deuterostomes. We find that Dll is expressed along the proximodistal axis of developing polychaete annelid parapodia, onychophoran lobopodia, ascidian ampullae, and even echinoderm tube feet. Dll/Dlx expression in such diverse appendages in these six coelomate phyla could be convergent, but this would have required the independent co-option of Dll/Dlx several times in evolution. It appears more likely that ectodermal Dll/Dlx expression along proximodistal axes originated once in a common ancestor and has been used subsequently to pattern body wall outgrowths in a variety of organisms. We suggest that this pre-Cambrian ancestor of most protostomes and the deuterostomes possessed elements of the genetic machinery for and may have even borne appendages.