Mitochondrial phylogeny of Anura (Amphibia): a case study of congruent phylogenetic reconstruction using amino acid and nucleotide characters

We explore whether phylogenetic analyses of the same sequence data set at the amino acid and nucleotide level are able to recover congruent topologies, as well as the advantages and limitations of both alternative approaches. As a case study, mitochondrial protein-coding genes were used to discern among competing hypotheses on the phylogenetic relationships of major anuran amphibian lineages. To properly address this phylogenetic question, the complete nucleotide sequences of the mitochondrial genomes of two archaeobatrachian species, Ascaphus truei and Pelobates cultripes, were determined anew. Bayesian and maximum likelihood phylogenetic inferences of the same sequence data set were performed based on both amino acid and nucleotide characters, with the latter analysed either as codons or as a reduced data set of first+second (P12) codon positions. In addition, likelihood-based ratio tests were performed to evaluate the support of alternative topologies. The different data sets arrived at congruent and highly supported topologies, suggesting a similar phylogenetic resolving power of the two character types provided that correctly selected sites and appropriate evolutionary models are used. The reconstructed anuran mitochondrial phylogeny supports the paraphyly of Archaeobatrachia, with Ascaphus as sister group to all the remaining anurans, and Pelobates as sister group of Neobatrachia. However, the employed tree reconstruction methods and likelihood-based ratio tests seemed to be negatively affected by the fast evolving sequences of neobatrachians, suggesting that the phylogeny of Anura here presented is not definitive, and needs further investigation using an extended taxon sampling.     

Endochitinase gene-based phylogenetic analysis of Trichoderma

DNA sequences of the single-copy gene coding for the 42 kDa endochitinase enzyme (EC 3.2.1.14) were used for phylogenetic analysis in Trichoderma. A set of 12 primers was developed and the entire gene was sequenced for 16 strains, and nucleotide and deduced amino acid sequences were compared to data from GenBank for additional Trichoderma strains. Analysis of the sequences revealed parsimony informative variation from 2.4 to 43.6% depending on the part of the gene (exons/introns) and the taxonomic level considered. Results are discussed in comparison to previous data from ITS-1 and ITS-2 rDNA sequencing and suggest the 42 kDa endochitinase gene as a potential molecular marker for reconstructing phylogenetic relationships in the genus Trichoderma at species level.     

The complete mitochondrial DNA sequence of the Mekong giant catfish (Pangasianodon gigas), and the phylogenetic relationships among Siluriformes

The Mekong giant catfish (Pangasianodon gigas) is the largest scale-less freshwater fish of the world, and a critically endangered species. We determined the complete nucleotide sequence (16,533 bp) of the mitochondrial genome of the Mekong giant catfish, and conducted phylogenetic analyses based on mitochondrial protein (the combined amino acid sequences of all 13 mitochondrial protein coding genes) and rRNA (the combined nucleotide sequences of mitochondrial 12S and 16S rRNA genes) data sets in order to further clarify the relative phylogenetic position of P. gigas, and to recover phylogenetic relationships among 15 out of the 33 families of Siluriformes. Phylogenetic analyses (maximum parsimony, minimum evolution, maximum likelihood, Bayesian inference) of the protein data set were congruent with a basal split of the order into Loricarioidei and Siluroidei, and supported a closer relationship of the Mekong giant catfish (family Pangasiidae) to Siluridae than to Bagridae. The rRNA-based Bayesian phylogeny recovered Callichthyidae as the sister group of all other analyzed non-diplomystid catfish families, rendering Loricarioidei paraphyletic. In addition, Loricariidae were recovered as paraphyletic due to the inclusion of Astroblepidae. However, none of the two relationships received bootstrap support in the maximum parsimony, minimum evolution, and maximum likelihood analyses, and should be interpreted with caution. The derived position of Cetopsidae within Siluroidei, the sister group relationship of Pseudopimelodidae and Pimelodidae, and a close relationship of Doradidae and Auchenipteridae to the exclusion of Mochokidae were strongly supported. Pangasiidae was placed as a single lineage without clear affinities.     

Phylogenetic relationships among families of the Scaphopoda (Mollusca)

Phylogenetic relationships among families in the molluscan class Scaphopoda were analysed using morphological characters and cladistic parsimony methods. A maximum parsimony analysis of 34 discrete characters, treated as unordered and equally weighted, from nine ingroup terminal taxa produced a single most parsimonious tree; supplementary analyses of tree length frequency distribution and Bremer support indices indicate a strong phylogenetic signal from the data and moderate to minimally supported clades. The traditional major division of the class, the orders Dentaliida and Gadilida, is supported as both taxa are confirmed as monophyletic clades. Within the Dentaliida, two clades are recognized, the first comprised of the families Dentaliidae and Fustiariidae, the second of the Rhabdidae and Calliodentaliidae; together, these groups comprise a third clade, which has the Gadilinidae as sister. Within the Gadilida, a nested series of relationships is found among [Entalinidae, [Pulsellidae, [Wemersoniellidae, Gadilidae]]]. These results lend cladistic support to earlier hypotheses of shared common ancestry for some families, but are at variance with other previous hypotheses of evolution in the Scaphopoda. Furthermore, analysis of constituent Gadilinidae representatives provide evidence for paraphyly of this family. The relationships supported here provide a working hypothesis that the development of new characters and greater breadth of taxonomic sampling can test, with a suggested primary goal of establishing monophyly at the family level.     

Interspecies sequence differences in the Mip protein from the genus Legionella: implications for function and evolutionary relatedness

The nucleotide sequence of the mip genes and their inferred amino acid sequences were determined from 35 Legionella species and compared with the published sequences for L . pneumophila , L . micdadei and L . longbeachae . The sequences were 69–97% conserved at the nucleotide level and 82–99% at the amino acid level, with total conservation of amino acids determined to be associated with sites known to be involved in peptidyl prolyl cis–trans isomerase activity. No apparent difference could be determined in the arrangement of amino acids that would predict a functional difference in Mip from species associated with disease and Mip from species isolated only from the environment. Additionally, a phylogenetic comparison of the sequences with published 16S RNA sequences, using both genetic distance and maximum parsimony methods, was performed. Few relationships were apparent that were well supported by both data sets, the most robust being a clade comprising {[( cincinnatiensis , longbeachae , sainthelensi , santicrucis ) gratiana ] ( moravica , quateirensis , shakespearei , worsleiensis ) anisa , bozemanii , cherrii , dumoffii , gormanii , jordanis , parisiensis , pneumophila , steigerwaltii , tucsonensis , and wadsworthii  }.     

ON COMBINING PROTEIN SEQUENCES AND NUCLEIC ACID SEQUENCES IN PHYLOGENETIC ANALYSIS: THE HOMEOBOX PROTEIN CASE

Abstract— Amino acid encoding genes contain character state information that may be useful for phylogenetic analysis on at least two levels. The nucleotide sequence and the translated amino acid sequences have both been employed separately as character states for cladistic studies of various taxa, including studies of the genealogy of genes in multigene families. In essence, amino acid sequences and nucleic acid sequences are two different ways of character coding the information in a gene. Silent positions in the nucleotide sequence (first or third positions in codons that can accrue change without changing the identity of the amino acid that the triplet codes for) may accrue change relatively rapidly and become saturated, losing the pattern of historical divergence. On the other hand, non-silent nucleotide alterations and their accompanying amino acid changes may evolve too slowly to reveal relationships among closely related taxa. In general, the dynamics of sequence change in silent and non-silent positions in protein coding genes result in homoplasy and lack of resolution, respectively. We suggest that the combination of nucleic acid and the translated amino acid coded character states into the same data matrix for phylogenetic analysis addresses some of the problems caused by the rapid change of silent nucleotide positions and overall slow rate of change of non-silent nucleotide positions and slowly changing amino acid positions. One major theoretical problem with this approach is the apparent non-independence of the two sources of characters. However, there are at least three possible outcomes when comparing protein coding nucleic acid sequences with their translated amino acids in a phylogenetic context on a codon by codon basis. First, the two character sets for a codon may be entirely congruent with respect to the information they convey about the relationships of a certain set of taxa. Second, one character set may display no information concerning a phylogenetic hypothesis while the other character set may impart information to a hypothesis. These two possibilities are cases of non-independence, however, we argue that congruence in such cases can be thought of as increasing the weight of the particular phylogenetic hypothesis that is supported by those characters. In the third case, the two sources of character information for a particular codon may be entirely incongruent with respect to phylogenetic hypotheses concerning the taxa examined. In this last case the two character sets are independent in that information from neither can predict the character states of the other. Examples of these possibilities are discussed and the general applicability of combining these two sources of information for protein coding genes is presented using sequences from the homeobox region of 46 homeobox genes fromDrosophila melanogasterto develop a hypothesis of genealogical relationship of these genes in this large multigene family.     

A large‐scale survey of heterochrony in anuran cranial ossification patterns

Most anurans have a biphasic life cycle, which includes metamorphosis from a tadpole stage to an adult frog. This process involves extensive transformations of the cranial skeleton, which have been of long‐standing interest with respect to anuran skeletal evolution and taxonomy. In this study, large‐scale patterns of anuran skeletal ossification are assessed by collecting the most comprehensive data set on anuran cranial ossification to date from the literature, including data for 45 anuran and one caudate outgroup species. Ossification sequences were translated into event‐pair matrices for explorative phylogenetic analysis and phylogenetically informed parsimony search for heterochrony using the Parsimov algorithm. Rank variability of single bones across species was also analysed. Little phylogenetic signal was retrieved from a parsimony‐based phylogenetic analysis of event‐pairs, and only a few species that are generally agreed to be closely related are placed close to each other (e.g. some Pipidae and Costata). Parsimov analysis revealed some clade‐specific heterochrony in anuran clades of varying inclusiveness. Our results show that relating heterochronic changes in anuran cranial ontogeny to parameters such as direct development or miniaturization is problematic because of the high evolvability of cranial ossification sequences. Rank variation analysis suggests that anuran cranial bones are highly variable in their sequence positioning, possibly because tadpole and adult cranial morphology do not co‐evolve. Elements which are lost in some species ossify at the end of the sequence, providing evidence for the notion that failure of anuran cranial elements to ossify is due to processes of paedomorphosis.     

Molecular differentiation and phylogeny of entomopathogenic nematodes (rhabditida: heterorhabditidae) based on ND4 gene sequences of Mitochondrial DNA.

We determined partial ND4 gene sequences of mitochondrial DNA from 15 heterorhabditid nematode isolates, representing 5 species collected from different regions of the world, by using polymerase chain reaction (PCR) and direct-sequencing of PCR products. Aligned nucleotide as well as amino acid sequences were used to differentiate nematode species by comparing sequence divergence and to infer phylogeny of the nematodes by using maximum parsimony and likelihood methods. Robustness of our phylogenetic trees was checked by bootstrap tests. The 15 nematode isolates can be divided into 7 haplotypes based on DNA sequences. On a larger scale, the sequence divergence revealed 4 distinct groups corresponding to 4 described species. No sequence divergence was detected from 5 isolates of Heterorhabditis bacteriophora or between Heterorhabditis marelatus to Heterorhabditis hepialius. Our sequence data yielded phylogenetic trees with identical topologies when different tree-building methods were used. Most relationships were also confirmed by using amino acid sequences in maximum parsimony analysis. Our molecular phylogeny of Heterorhabditis species support an existing taxonomy that is based largely on morphology and the sequence divergence of the ND4 gene permits species identification.     

Phylogenetic analysis of the vertebrate glycoprotein hormone family including new sequences of sturgeon (Acipenser baeri) beta subunits of the two gonadotropins and the thyroid-stimulating hormone

The beta subunits of the two gonadotropins (GTH1 and GTH2) and of the thyroid-stimulating hormone (TSH) of a chondrostean fish, Acipenser baeri, were cloned. These new sequences and selected representative members of beta subunits of vertebrate glycoprotein hormones, including tetrapod follicle-stimulating hormones (FSH) and luteinizing hormones (LH), allowed us to infer the phylogenetic relationships within this family. Both distance matrix and maximum parsimony methods were used on both nucleotide and amino acid sequences, with bootstrapping evaluation over 1000 replicates. The four trees obtained had highly similar topologies. In each case, three monophylogenetic lineages, TSH, GTH1-FSH, and GTH2-LH were clearly identified. The three monophylogenetic lineages were supported by 21-23 specific characters at the amino acid level, out of a total of 121 characters. The resolved topologies within each monophyletic hormone cluster were congruent with the known phylogenetic relationships between the related species. The inferred parental relationships within gonadotropins are in agreement with data concerning their biological functions. The present study demonstrates that GTH1 and GTH2 are the actinopterygian homologues of tetrapod FSH and LH, respectively.     

Granule-bound starch synthase: structure, function, and phylogenetic utility

Interest in the use of low-copy nuclear genes for phylogenetic analyses of plants has grown rapidly, because highly repetitive genes such as those commonly used are limited in number. Furthermore, because low- copy genes are subject to different evolutionary processes than are plastid genes or highly repetitive nuclear markers, they provide a valuable source of independent phylogenetic evidence. The gene for granule-bound starch synthase (GBSSI or waxy) exists in a single copy in nearly all plants examined so far. Our study of GBSSI had three parts: (1) Amino acid sequences were compared across a broad taxonomic range, including grasses, four dicotyledons, and the microbial homologs of GBSSI. Inferred structural information was used to aid in the alignment of these very divergent sequences. The informed alignments highlight amino acids that are conserved across all sequences, and demonstrate that structural motifs can be highly conserved in spite of marked divergence in amino acid sequence. (2) Maximum-likelihood (ML) analyses were used to examine exon sequence evolution throughout grasses. Differences in probabilities among substitution types and marked among-site rate variation contributed to the observed pattern of variation. Of the parameters examined in our set of likelihood models, the inclusion of among-site rate variation following a gamma distribution caused the greatest improvement in likelihood score. (3) We performed cladistic parsimony analyses of GBSSI sequences throughout grasses, within tribes, and within genera to examine the phylogenetic utility of the gene. Introns provide useful information among very closely related species, but quickly become difficult to align among more divergent taxa. Exons are variable enough to provide extensive resolution within the family, but with low bootstrap support. The combined results of amino acid sequence comparisons, maximum-likelihood analyses, and phylogenetic studies underscore factors that might affect phylogenetic reconstruction. In this case, accommodation of the variable rate of evolution among sites might be the first step in maximizing the phylogenetic utility of GBSSI.      

A cladistic analysis of chloroplast DNA restriction site variation and morphology for the genera of the Juglandaceae

The phylogenetic relationships of the genera of the Juglandaceae are examined with cladistic analyses of chloroplast DNA (cpDNA) restriction site variation and morphology. Rates of evolution of the chloroplast genome are slower than in many other groups of plants, enabling the entire genome to be utilized at the intergeneric level. The trees resulting from the two independent analyses were completely congruent. The combined analysis of the two data sets produced a tree completely congruent with the cladogram from the two data sets analyzed independently. The cladogram is compared with previous classifications, cladistic analyses, and fossil history for the family. Although the topology resulting from the cladistic analyses of this study was strongly congruent with previous estimates of relationships within the family, the fossil record indicates that the basal-most lineages in the cladistic trees arose later than the more terminal lineages. This reversed order of origin indicates that perhaps the rooting of the trees is erroneous.     

Phylogeny of West AfricanCaryedon(Coleoptera: Bruchidae): Congruence between Molecular and Morphological Data

Seed beetles belonging to the Old World genusCaryedonfeed in the seeds of various Caesalpinioideae, Mimosoideae, and Combretaceae. In an attempt to resolve broad phylogenetic relationships within the genus, we obtained 332 base pair sequences of mitochondrial 12S ribosomal DNA and morphological data for the 16 West AfricanCaryedonspecies. Morphological characters were analyzed under maximum parsimony and sequences were compared under maximum parsimony, maximum likelihood, and neighbor joining. Using a partition homogeneity test, we determined that morphological and molecular data sets were combinable. Combined data were analyzed under maximum parsimony. Morphological and molecular trees were congruent at the species group level and total evidence analyses yielded the same topologies as molecular data with each of the three outgroups used. Four main terminal clades are recognized, each corresponding with a group of species generally feeding on the same host plant family, subfamily, genus, or species. The monophyly of legume feedingCaryedonis supported by both data sets, and Combretaceae feeders split in two monophyletic assemblages.     

Molecular phylogeny of ateline new world monkeys (Platyrrhini, atelinae) based on gamma-globin gene sequences: evidence that brachyteles is the sister group of lagothrix

Nucleotide sequences, each spanning approximately 7 kb of the contiguous gamma1 and gamma2 globin genomic loci, were determined for seven species representing all extant genera (Ateles, Lagothrix, Brachyteles, and Alouatta) of the New World monkey subfamily Atelinae. After aligning these seven ateline sequences with outgroup sequences from several other primate (non-ateline) genera, they were analyzed by maximum parsimony, maximum likelihood, and neighbor-joining algorithms. All three analyzes estimated the same phylogenetic relationships: [Alouatta [Ateles (Brachyteles, Lagothrix)]]. Brachyteles and Lagothrix are sister-groups supported by 100% of bootstrap replications in the parsimony analyses. Ateles joins this clade, followed by the basal genus Alouatta; these joinings were strongly supported, again with 100% bootstrap values. This cladistic pattern for the four ateline genera is congruent with that obtained in previous studies utilizing epsilon-globin, IRBP, and G6PD nuclear genomic sequences as well as mitochondrial COII sequences. Because the number of aligned nucleotide positions is much larger in the present datasetoff than in any of these other datasets, much stronger support was obtained for the cladistic classification that divides subfamily Atelinae into tribes Alouattini (Alouatta) and Atelini, while the latter divides into subtribes Atelina (Ateles) and Brachytelina (Brachyteles and Lagothrix).     

Molecular phylogeny of the New World monkeys (Platyrrhini,primates) based on two unlinked nuclear genes: IRBP intron 1 and ϵ-globin sequences

Nuclear sequences of the 1.8 kilobase (kb) long intron 1 of the interstitial retinol-binding protein gene (IRBP), previously determined for 11 of the 16 extant genera of New World monkeys (superfamily Ceboidea, infraorder Platyrrhini), have now been determined for the remaining 5 genera. The maximum parsimony trees found, first with IRBP sequences alone and then with tandemly combined IRBP and ϵ-globin gene sequences from the same species, supported a provisional cladistic classification with the following clusters. Subtribes Callitrichina (Callithrix, Cebuella), Callimiconina (Callimico), Leontopithecina (Leontopithecus) and Saguina (Saguinus) constitute subfamily Callitrichinae, and subfamilies Callitrichinae, Aotinae (Aotus), and Cebinae (Cebus, Saimiri) constitute family Cebidae. Subtribes Chiropotina (Chiropotes, Cacajao) and Pitheciina (Pithecia) constitute tribe Pitheciini; and tribes Pitheciini and Callicebini (Callicebus) constitute subfamily Pitheciinae. Subtribes Brachytelina (Brachyteles, Lagothrix) and Atelina (Ateles) constitute tribe Atelini, and tribes Atelini and Alouattini (Alouatta) constitute subfamily Atelinae. The parsimony results were equivocal as to whether Pitheciinae should be grouped with Atelinae in family Atelidae or have its own family Pitheciidae. The cladistic groupings of extant ceboids were also examined by different stochastic evolutionary models that employed the same stochastic process of nucleotide substitutions but alternative putative phylogenetic trees on which the nucleotide substitutions occurred. Each model, i.e., each different tree, predicted a different multinomial distribution of nucleotide character patterns for the contemporary sequences. The predicted distributions that were closest to the actual observed distributions identified the best fitting trees. The cladistic relationships depicted in these best fitting trees agreed in almost all cases with those depicted in the maximum parsimony trees. © 1996 Wiley-Liss, Inc.     

Plethodontid salamander mitochondrial genomics: A parsimony evaluation of character conflict and implications for historical biogeography

A new parsimony analysis of 27 complete mitochondrial genomic sequences is conducted to investigate the phylogenetic relationships of plethodontid salamanders. This analysis focuses on the amount of character conflict between phylogenetic trees recovered from newly conducted parsimony searches and the Bayesian and maximum likelihood topology reported by Mueller et al. (2004 ; PNAS, 101, 13820–13825). Strong support for Hemidactylium as the sister taxon to all other plethodontids is recovered from parsimony analyses. Plotting area relationships on the most parsimonious phylogenetic tree suggests that eastern North America is the origin of the family Plethodontidae supporting the “Out of Appalachia” hypothesis. A new taxonomy that recognizes clades recovered from phylogenetic analyses is proposed. © The Willi Hennig Society 2005.     

Phylogeny of the Asian spiny frog tribe Paini (Family Dicroglossidae) sensu Dubois

The anuran tribe Paini, family Dicroglossidae, is known in this group only from Asia. The phylogenetic relationships and often the taxonomic recognition of species are controversial. In order to stabilize the classification, we used approximately 2100 bp of nuclear (rhodopsin, tyrosinase) and mitochondrial (12S, 16S rRNA) DNA sequence data to infer the phylogenetic relationships of these frogs. Phylogenetic trees reconstructed using Bayesian inference and maximum parsimony methods supported a monophyletic tribe Paini. Two distinct groups (I,II) were recovered with the mtDNA alone and the total concatenated data (mtDNA+nuDNA). The recognition of two genera, Quasipaa and Nanorana, was supported. Group I, Quasipaa, is widespread east of the Hengduan Mountain Ranges and consists of taxa from relatively low elevations in southern China, Vietnam and Laos. Group II, Nanorana, contains a mix of species occurring from high to low elevation predominantly in the Qinghai-Tibetan Plateau and Hengduan Mountain Ranges. The occurrence of frogs at high elevations appears to be a derived ecological condition. The composition of some major species groups based on morphological characteristics strongly conflicts with the molecular analysis. Some possible cryptic species are indicated by the molecular analyses. The incorporation of genetic data from type localities helped to resolve some of the taxonomic problems, although further combined analyses of morphological data from type specimens are required. The two nuDNA gene segments proved to be very informative for resolving higher phylogenetic relationships and more nuclear data should be explored to be more confident in the relationships.     

Phylogenetic relationships of discoglossid frogs (Amphibia:Anura:Discoglossidae) based on complete mitochondrial genomes and nuclear genes

The complete nucleotide sequence of the mitochondrial (mt) genome was determined for three species of discoglossid frogs (Amphibia:Anura:Discoglossidae), representing three of the four recognized genera: Alytes obstetricans, Bombina orientalis, and Discoglossus galganoi. The organization and size of these newly determined mt genomes are similar to those previously reported for other vertebrates. Phylogenetic analyses (maximum likelihood, Bayesian inference, minimum evolution, and maximum parsimony) of mt protein-coding genes at the amino acid level were performed in combination with already published mt genome sequence data of three species of Neobatrachia, one of Pipoidea, and four of Caudata. Phylogenetic analyses based on the deduced amino acid sequences of all mt protein-coding genes arrived at the same topology. The monophyly of Discoglossidae is strongly supported. Within the Discoglossidae, Alytes is consistently recovered as sister group of Discoglossus, to the exclusion of Bombina. The three species representing Neobatrachia exhibited extremely long branches irrespective of the phylogenetic inference method used, and hence their relative position with respect to Discoglossidae and Xenopus may be artefactual due to a severe long branch attraction effect. To further investigate the phylogenetic intrarelationships of discoglossids, nucleotide sequences of four nuclear protein-coding genes (CXCR4, RAG1, RAG2, and Rhodopsin) with sequences available for the three discoglossid genera and Xenopus were retrieved from GenBank, and together with a concatenated nucleotide sequence data set containing all mt protein-coding genes except ND6 were subjected to separate and combined phylogenetic analyses. In all cases, a sister group relationship between Alytes and Discoglossus was recovered with high statistical support.     

Protein coding gene trees in ciliates: comparison with rRNA-based phylogenies

We have reexamined the phylogeny of the ciliates using alpha-tubulin and phosphoglycerate kinase gene sequences. For alpha-tubulin, we have compared the amino acid and nucleotide sequences of 20 species representing seven of the nine classes of the phylum (Karyorelictea, Heterotrichea, Hypotrichea, Oligohymenophorea, Colpodea, Nassophorea, and Litostomatea). The phylogenetic tree resembles a bush from which three monophyletic lineages can be distinguished which correspond to the three classes Hypotrichea, Oligohymenophorea, and Litostomatea. For phosphoglycerate kinase, we have compared the amino acid sequences from 7 species representing three classes (Heterotrichea, Hypotrichea, and Oligohymenophorea). The branching pattern is resolved in three deeply separated branches with an early emergence of the heterotrich. Our comparative analysis shows that if alpha-tubulin phylogeny is not informative at the interclass level, the preliminary data from the phosphoglycerate kinase molecule appear more promising. Nevertheless, at low taxonomic level and at the class level, the resolved phylogenetic relationships inferred from both protein and rRNA sequence data are congruent.     

Relationships among didelphid marsupials based on sequence variation in the mitochondrial cytochrome B gene

Variation in the mitochondrial cytochrome b gene (nucleotide and amino acid sequences) is evaluated for 9 genera and 15 species of American opossums in the family Didelphidae, using the American caenolestid rat opossumLestoros and the New Guinean peroryctid bandicootEchimypera as outgroups. Phylogenetic analyses (parsimony and distance) strongly support the monophyly of the Didelphidae and delineate two major clades; (1)Didelphis andPhilander are strongly aligned sister taxa, withMetachirus weakly but consistently associated with them, and (2)Marmosa plusMicoureus, withMonodelphis falling outside that pair. The generaMarmosops, Caluromys, andGlironia exhibit varied relationships, depending upon the method of analysis and data (DNA or amino acid sequences) used, but generally are placed individually or in combinations near or at the base of the didelphid radiation. Some aspects of these relationships are consistent with current taxonomic views, but others are in marked contrast. Specifically, a clade comprised of the mouse opossumsMarmosa, Micoureus, andMarmosops is strongly rejected by log-likelihood analysis, contrary to expectations from some current classifications. Also, the woolly opossumsCaluromys andGlironia also do not form a sister-taxon relationship, as suggested by their placement in a subfamily separate from the remaining didelphids examined. However, such a relationship cannot be rejected from log-likelihood analyses. The relationships suggested fromcyt-b sequences are strongly concordant with those based on DNA-DNA hybridization analyses. In addition to systematic and phylogenetic properties, molecular evolution of the didelphid cytochrome b gene sequence is characterized according to nucleotide bias and rate differentials at each codon position and across the entire sequence.To whom correspondence should be addressed.