Against expectation: a short sequence with high signal elucidates cone snail phylogeny

A short (259 nucleotide) conserved intronic sequence (CIS) is surprisingly informative for delineating deep phylogenetic relationships in cone snails. Conus species previously have been assigned to clades based on the evidence from mitochondrial 12S and 16S rRNA gene sequences (1129 bp). Despite their length, these genes lack the phylogenetic information necessary to resolve the relationships among the clades. Here we show that the relationships can be inferred from just 46 sites in the very short CIS sequence (a portion of "intron 9" of the γ-glutamyl carboxylase gene). This is counterintuitive because in short sequences sampling error (noise) often drowns out phylogenetic signal. The intron 9 CIS is rich in synapomorphies that define the divergence patterns among eight clades of worm- and fish-hunting Conus, and it contains almost no homoplasy. Parsimony, maximum likelihood and Bayesian analyses of the combined sequences (mt rRNA+CIS) confirm most of the relationships among 23 Conus sequences. This phylogeny implies that fish-hunting behavior evolved at least twice during the history of Conus-once among New World species and independently in the Indo-Pacific clades.     

Evaluation of atpB nucleotide sequences for phylogenetic studies of ferns and other pteridophytes

Inferring basal relationships among vascular plants poses a major challenge to plant systematists. The divergence events that describe these relationships occurred long ago and considerable homoplasy has since accrued for both molecular and morphological characters. A potential solution is to examine phylogenetic analyses from multiple data sets. Here I present a new source of phylogenetic data for ferns and other pteridophytes. I sequenced the chloroplast gene atpB from 23 pteridophyte taxa and used maximum parsimony to infer relationships. A 588-bp region of the gene appeared to contain a statistically significant amount of phylogenetic signal and the resulting trees were largely congruent with similar analyses of nucleotide sequences from rbcL. However, a combined analysis of atpB plus rbcL produced a better resolved tree than did either data set alone. In the shortest trees, leptosporangiate ferns formed a monophyletic group. Also, I detected a well-supported clade of Psilotaceae (Psilotum and Tmesipteris) plus Ophioglossaceae (Ophioglossum and Botrychium). The demonstrated utility of atpB suggests that sequences from this gene should play a role in phylogenetic analyses that incorporate data from chloroplast genes, nuclear genes, morphology, and fossil data.     

Phylogeography of Aglais urticae (Lepidoptera) based on DNA sequences of the mitochondrial COI gene and control region

Mitochondrial DNA (mtDNA) sequences of the COI gene and the control region were used to examine the genetic population structure of Aglais urticae L. (Lepidoptera) over its entire geographic range, i.e., the Palaearctic. The phylogenetic relationships within and between A. urticae subspecies were determined and patterns of mtDNA divergence and ecological differentiation were compared. High gene flow together with a recent and sudden population expansion characterise the genetic population structure of this species. No geographically induced differentiation was observed, nor were subspecies identified as separate evolutionary units. The discrepancy between the genetic and ecological variation is most likely due to the slower rate of mtDNA evolution compared to ecological differentiation. The control region proved to be a less useful molecular marker for the population genetics and the phylogenetic reconstruction of closely related taxa in A. urticae than it has for other species. The extreme bias in adenine and thymine content (A+T=90.91%) probably renders this region highly susceptible to homoplasy, resulting in a less informative molecular marker.     

Molecular phylogeny of the nasuta subgroup of Drosophila based on 12S rRNA, 16S rRNA and CoI mitochondrial genes, RAPD and ISSR polymorphisms

The nasuta subgroup is a cluster of morphologically almost similar forms with a wide range of geographic distribution. During the last three decades nature of inter-relationship among the members has been investigated at different levels of organization. The phylogenetic relationships of the members of the nasuta subgroup of the immigrans species group of Drosophila was made by employing Random Amplified Polymorphic DNA (RAPD), Inter Simple Sequence Repeats-PCR (ISSR-PCR) polymorphisms, mitochondrial 12S rRNA, 16S rRNA and Cytochrome C Oxidase subunit I (CoI) gene sequences. The phylogenetic tree generated by RAPD analysis is in nearly complete congruence with the classification based on morphophenotypic characters. The 12S and 16S rRNA genes were highly conserved across the nasuta subgroup and revealed only 3 and 4 variable sites respectively, of which only one site was informative. The CoI gene, on the other hand, revealed 57 variable sites of which 25 sites were informative. All the three species of orbital sheen complex were included in a major cluster in the phylogenetic trees derived from mitochondrial gene sequence data consistent with the morphophenotypic classification. The CoI analysis placed two species of frontal sheen complex, D. n. nasuta and D. n. albomicans in two different clades and this is inconsistent with morphological classification. The molecular clock suggested that divergence between the kohkoa complex and the albomicans complex occurred approximately 2.2 MYA, indicating recent evolution of the nasuta subgroup. The higher transition bias in the mitochondrial genes reported in the present study also suggested recent evolution of the nasuta subgroup.     

5S rRNA sequences of 12 species of flatworms: implications for the phylogeny of the Platyhelminthes

5S rRNAs from 12 species of free living and parasitic platyhelminthes were sequenced. In the phylogenetic analysis, attention was focused on the statistical estimates of the trees corresponding to existing phylogenetic hypotheses. The available 5S rRNA data agree well with widely accepted views on the relationships between the Acoela, Polycladida, Tricladida, and Neorhabdocoela; our analysis of the published 18S rRNA sequences also demonstrated good correspondence between these views and molecular data. With available 5S rRNA data the hypothesis that the dalyellioid turbellarians is the sister group of the Neodermata is less convincing than the hypotheses proposing the Neodermata as the sister group of the Neorhabdocoela, or of the Seriata, or of the branch uniting them. A relatively low rate of base replacement in parasitic flatworms, probably, accounts for the uncertain position of the Neodermata, while a relatively high rate in planarians may explain a relatively too early divergence of the Tricladida in several published phylogenetic trees constructed from various rRNA data.     

COASTAL BOTTLENOSE DOLPHINS FROM SOUTHEASTERN AUSTRALIA ARE TURSIOPS ADUNCUS ACCORDING TO SEQUENCES OF THE MITOCHONDRIAL DNA CONTROL REGION

Sequence analysis of the mitochondrial DNA control region was used to clarify the taxonomic status of two coastal bottlenose dolphin populations from southeastern Australia currently classified as Tursiops truncatus . A 368-bp segment of the control region of 57 biopsy-sampled, photo-identified dolphins of Jervis Bay and Port Stephens was compared to published sequences of T. truncatus and T. aduncus from different oceanic regions. Sequence divergence between haplotypes from southeastern Australia and T. aduncus was much lower than that from T. truncatus . Analyses using two different methods of phylogenetic reconstruction unambiguously placed all haplotypes from southeastern Australia in a group composed exclusively of T. aduncus . The results strongly indicated that these two bottlenose dolphin populations belong to T. aduncus , extending the range of the species to subtropical waters of the Western South Pacific Ocean.     

Amplification and sequence analysis of the full length toxR gene in Vibrio harveyi

This study was focused on obtaining the complete gene sequence of the toxR gene in V. harveyi by using toxR-targeted PCR to amplify 5' and 3' regions flanking the 576-bp Vibrio harveyi (NBRC 15634) toxR gene fragment previously amplified using degenerate PCR. To obtain the 5' flanking sequences, a forward PCR primer (VhtoxRpv) was designed based on known sequences upstream of toxR in V. parahaemolyticus and V. vulnificus. The reverse primer (VctoxR2R) was based on the sequence of the 576-bp Vibrio harveyi toxR fragment. The resulting 750-bp amplicon was sequenced, providing the 5' sequences of the V. harveyi (NBRC 15634) toxR gene. The 3' flanking region was amplified using a primer pair toxRS1 and toxRS2 based on V. parahaemolyticus and V. vulnificus toxR and toxS, resulting in a 900-bp amplicon that contained the remaining 3' sequences of the V. harveyi NBRC 15634 toxR. This paper reports, for the first time, a complete 882-bp nucleotide sequence for toxR in Vibrio harveyi. Sequence analysis and alignment revealed that the complete toxR gene in V. harveyi shares 87% sequence similarity with toxR of V. parahaemolyticus, 84% similarity with V. fluvialis, 83% with V. vulnificus and partial sequence of V. campbellii. The phylogenetic trees revealed wider divergence in toxR compared to 16S rRNA genes, so that V. harveyi could easily be distinguished from V. campbellii and V. parahaemolyticus.     

Molecular phylogenetic relationships of pond frogs distributed in the Palearctic region inferred from DNA sequences of mitochondrial 12S ribosomal RNA and cytochrome b genes

The evolutionary relationships of pond frogs distributed in the Far East and Europe were investigated by analyses of nucleotide sequences of mitochondrial 12S ribosomal RNA (12S rRNA) and cytochrome b (cyt b) genes. The nucleotide sequences of a 412-bp segment of the 12S rRNA gene and a 534-bp segment of the cyt b gene were determined by the PCR-direct sequencing method using 19 frogs belonging to six species and one subspecies distributed in the Palearctic region. Phylogenetic trees were constructed by the neighbor-joining and maximum-likelihood methods using Rana catesbeiana or Xenopus laevis as an outgroup. The 412-bp segment of the 12S rRNA gene contained 65 variable sites including gap sites, and the 534-bp segment of the cyt b gene contained 160 variable sites. The nucleotide sequence divergences of the 12S rRNA gene were 0.25-4.83% within the Far Eastern frogs, 0.25-6.22% within the European frogs, and 8.74-11.24% between the Far Eastern and the European frogs, whereas those of the cyt b gene were 3.64-14.73% within the Far Eastern frogs, 0.38-14.42% within the European frogs, and 16.53-23.58% between the Far Eastern and the European frogs. Although most nucleotide substitutions were at the third codon position of the cyt b gene and were silent mutations, 4 amino acid replacements occurred within the Far Eastern frogs, 4 within the European frogs, and 11 between the Far Eastern and the European frogs. The phylogenetic trees constructed from the nucleotide sequence divergences showed slightly different topologies for the 12S rRNA and cyt b genes. R. esculenta from Ukraine was closely related to R. lessonae from Luxembourg in both the 12S rRNA and the cyt b gene sequences.     

Conserved flanking microsatellite sequences (ReFS) differentiate between Lepidoptera species,and provide insight into microsatellite evolution

High rates of mutation and homoplasy mean that microsatellites generally are not considered to be useful molecular markers for inferring systematic relationships between species. However, an earlier pilot study suggested that conserved flanking microsatellite sequences, also known as repetitive flanking sequences (ReFS), may form a basis for a dominant marker that can differentiate between species of Lepidoptera. We present data that demonstrate that ReFS are quick and easy to use, and generate highly repeatable banding patterns from a range of Lepidoptera species. Sequence data from a subset of ReFS‐amplified bands revealed microsatellite families with flanking sequences that are more conserved within than among species: this is probably attributable to recombination‐mediated events, transposition of mobile elements or a combination of the two. Our data support the use of ReFS as dominant interspecific molecular markers, and add to the growing literature on the evolution of microsatellites in Lepidoptera.     

Phylogenetic position of someChlorella species within the chlorococcales based upon complete small-subunit ribosomal RNA sequences

Summary Complete small-subunit rRNA (16S-like rRNA) coding region sequences were determined for eight species of the Chlorococcales (Chlorophyceae). The genera investigated includePrototheca, Ankistrodesmus, Scenedesmus, and fiveChlorella species. Distance matrix methods were used to infer a phylogenetic tree that describes evolutionary relationships between several plant and green algal groups. The tree exhibits a bifurcation within the Chlorococcales consistent with the division into Oocystaceae and Scenedesmaceae, but three of the fiveChlorella species are more similar to other algae than toChlorella vulgaris. All of the sequences contain primary and secondary structural features that are characteristic of 16S-like rRNAs of chlorophytes and higher plants.Anikstrodesmus stipitatus, however, contains a 394-bp group I intervening sequence in its 16S-like rRNA coding region.     

Phylogenetic analysis of the mitochondrial cytochrome c oxidase subunit 1 gene from 13 sipunculan genera: intra- and interphylum relationships

Abstract. Sipunculans are a phylum of non-segmented, marine worms. Although they are well characterized morphologically, relationships within the phylum and the relationship of Sipuncula to other spiralian phyla have been strongly debated. I analyzed representatives of 13 of 17 described genera using a 654-bp fragment of the mitochondrial gene, cytochrome c oxidase subunit I, to construct the first intraphylum phylogenetic hypothesis for sipunculans based on molecular sequence data. Within the phylum, tree topologies are loosely congruent with a previously published morphological analysis, except that the monotypic genus Phascolopsis occurred within the Golfingiaformes as a sister group to, or nested within, the Themistidae. Phylogenetic analyses, including 30 sequences from additional invertebrate taxa, suggest that sipunculans are most closely related to the Annelida (including Echiura). A previously proposed sipunculan-molluscan relationship is not supported. While not universally accepted, this hypothesis is consistent with other recent and past data on phylum-level relationships.     

Discordant phylogenies within the rrn loci of Rhizobia

It is evident from complete genome sequencing results that lateral gene transfer and recombination are essential components in the evolutionary process of bacterial genomes. Since this has important implications for bacterial systematics, the primary objective of this study was to compare estimated evolutionary relationships among a representative set of alpha-Proteobacteria by sequencing analysis of three loci within their rrn operons. Tree topologies generated with 16S rRNA gene sequences were significantly different from corresponding trees assembled with 23S rRNA gene and internally transcribed space region sequences. Besides the incongruence in tree topologies, evidence that distinct segments along the 16S rRNA gene sequences of bacteria currently classified within the genera Bradyrhizobium, Mesorhizobium and Sinorhizobium have a reticulate evolutionary history was also obtained. Our data have important implications for bacterial taxonomy, because currently most taxonomic decisions are based on comparative 16S rRNA gene sequence analysis. Since phylogenetic placement based on 16S rRNA gene sequence divergence perhaps is questionable, we suggest that the proposals of bacterial nomenclature or changes in their taxonomy that have been made may not necessarily be warranted. Accordingly, a more conservative approach should be taken in the future, in which taxonomic decisions are based on the analysis of a wider variety of loci and comparative analytical methods are used to estimate phylogenetic relationships among the genomes under consideration.     

An early origin of plastids within the cyanobacterial divergence is suggested by evolutionary trees based on complete 16S rRNA sequences

It is generally accepted that the plastids arose from a cyanobacterial ancestor, but the exact phylogenetic relationships between cyanobacteria and plastids are still controversial. Most studies based on partial 16S rRNA sequences suggested a relatively late origin of plastids within the cyanobacterial divergence. In order to clarify the exact relationship and divergence order of cyanobacteria and plastids, we studied their phylogeny on the basis of nearly complete 16S rRNA gene sequences. The data set comprised 15 strains of cyanobacteria from different morphological groups, 1 prochlorophyte, and plastids belonging to 8 species of plants and 12 species of diverse algae. This set included three cyanobacterial sequences determined in this study. This is the most comprehensive set of complete cyanobacterial and plastidial 16S rRNA sequences used so far. Phylogenetic trees were constructed using neighbor joining and maximum parsimony, and the reliability of the tree topologies was tested by different methods. Our results suggest an early origin of plastids within the cyanobacterial divergence, preceded only by the divergence of two cyanobacterial genera, Gloeobacter and Pseudanabaena.      

Evolution of Microsatellite Alleles in Four Species of Mice (Genus Apodemus)

Microsatellite length variation was investigated at a highly variable microsatellite locus in four species of Apodemus. Information obtained from microsatellite allele sequences was contrasted with allele sizes, which included 18 electromorphs. Additional analysis of a 400-bp unique sequence in the flanking region identified 26 different haplotype sequences or ``true' alleles in the sample. Three molecular mechanisms, namely, (1) addition/deletion of repeats, (2) substitutions and indels in the flanking region, and (3) mutations interrupting the repeat, contributed to the generation of allelic variation. Size homoplasy can be inferred for alleles within populations, from different populations of the same species, and from different species. We propose that microsatellite flanking sequences may be informative markers for investigating mutation processes in microsatellite repeats as well as phylogenetic relationships among alleles, populations, and species. Received: 3 November 1999 / Accepted: 2 May 2000     

RNA SEQUENCE DATA IN PHYLOGENETIC RECONSTRUCTION: TESTING THE LIMITS OF ITS RESOLUTION

Abstract— 18S ribosomal RNA sequences from 11 echinoderms are analysed using parsimony to investigate phylogenetic relationships. Their estimated divergence limes range from less than 20 Ma to more than 550 Ma before present. Phylogenies based on 18S rRNA sequence data are compared with well-established morphological phylogenies to discover at what evolutionary distance the two approaches start to produce incongruent results. Three regions of the 18S rRNA molecule are analysed separately and together, and paired and unpaired sites are also treated separately and combined.
Results show that a parsimony analysis of sequence data produces reliable results only when taxa have diverged more recently than about 100 Ma. At greater evolutionary distances (up to 250 Ma), paired nucleotides produce more reliable results than unpaired, while paired and unpaired data combined produce intermediate results. All trees within about 1% of the most parsimonious solution ought to be accepted. Transversions give results almost as reliable as paired regions though there were relatively few informative sites. The relationships of echinoderm classes, which diverged 450–550 Ma ago, are unresolved by 18S rRNA data.     

Multiple nuclear pseudogenes of mitochondrial DNA exist in the canine genome

Many copies of nuclear counterparts of mitochondrial DNA (mtDNA) were found in nuclear DNA from sperm heads of the domestic dog, Canis familiaris, by DNA-DNA hybridization and DNA sequencing. Nuclear counterparts homologous to the mtDNA D-loop region were cloned into lambda phage vectors (EMBL4 and lambda gt11), and nucleotide sequences of seven different mtDNA pseudogenes were then determined. The seven pseudogenes were E3 (474 bp; 82% homology with canine mtDNA), E13 (1867 bp; 67%), 8B (2375 bp; 78%), 12A (2650 bp; 79%), 33 (4131 bp; 86%), 47 (4251 bp; 86%), and E17 (5721 bp; 71%). These seven mtDNA pseudogenes corresponded to portions of cytoplasmic mtDNA containing the genes ile, ND1, leu, 16S rRNA, val, 12S rRNA, phe, D-loop, pro, thr, cytb, and glu. A neighbor-joining phylogenetic tree constructed from 12S rRNA sequences in mtDNA pseudogenes 8B, 33, 47, and E17 and in 10 mtDNA fragments from other species showed that these four pseudogenes form a monophyletic clade with canine mtDNA. A neighbor-joining phylogenetic tree based on the 318-bp cytb region showed that the canine pseudogenes existed before the divergence of 17 related canids, and their divergence dates were calculated at around 4.4 to 8.6 million years ago.     

Extreme intraspecific mitochondrial DNA sequence divergence in Galaxias maculatus (Osteichthys: Galaxiidae), one of the world's most widespread freshwater fish

Biogeographic controversies surrounding the widespread freshwater fish, Galaxias maculatus, were addressed with DNA sequence data. Mitochondrial cytochrome b and 16S rRNA sequences were obtained from representatives of six populations of this species. Substantial levels of cytochrome b (maximum 14.6%) and 16S rRNA sequence divergence (maximum 6.0%) were detected between western Pacific (Tasmania-New Zealand) and South American (Chile-Falkland Islands) haplotypes. A considerable level of divergence was also detected between Tasmanian and New Zealand haplotypes (maximum 5.1%) and within and among Chilean and Falkland Island G. maculatus (maximum 3. 8%). The phylogenetic structure of haplotypes conflicts with the accepted pattern of continental fragmentation. Molecular clock calibrations suggest that haplotype divergences postdate the fragmentation of Gondwana. These findings point to marine dispersal rather than ancient vicariance as an explanation for the wide distribution. The phylogenetic structure of South American haplotypes was not consistent with their geographic distribution. We consider factors such as population divergence, population size, dispersal, secondary contact, and philopatry as potential causes of the high level of mtDNA nucleotide diversity in this species.     

Mitochondrial COI sequences of brachiopods: genetic code shared with protostomes and limits of utility for phylogenetic reconstruction

A 1230-bp region of the cytochrome c oxidase subunit I (COI) gene of mitochondrial DNA of each of 16 brachiopod species, representing all five living orders, was amplified by polymerase chain reaction and sequenced. Pairwise comparisons of sequence differences plotted against divergence times estimated from the brachiopod fossil record revealed that, although there are considerable variations in the expected substitution rate among different lineages, amino acid substitutions of the COI sequences may largely become saturated in 100 Ma, due mostly to multiple substitutions at the same site. Coinciding with this result, phylogenetic analysis indicated low bootstrap values for nodes corresponding to divergence events that occurred before 100 Ma, suggesting that COI sequences are suitable only for inference of phylogenetic events subsequent to the Mesozoic. Examination of brachiopod codons corresponding to invariant amino acids in the COI of various other animals suggest the nonuniversal codon relationships UGA = Trp, AUA = Met, AAA/G = Lys, and AGA/G = Ser. These are identical to those in mollusks, annelids, and arthropods, consistent with the conclusion that brachiopods are protostomes, as indicated by previous molecular analyses.     

Diversity and relationships of Bradyrhizobia from Amphicarpaea bracteata based on partial nod and ribosomal sequences.

Partial sequences of three nod genes (nodC, nodD1, and nodA 5' flanking region) and of 16S and 23S rDNA were obtained from isolates of Bradyrhizobium sp. associated with the native North American legume Amphicarpaea bracteata. Isolates from Amphicarpaea had identical sequences in the three nod gene regions, but differed from all other Bradyrhizobium taxa at > 10% of nucleotide sites. Parsimony analysis of all nod gene segments indicated a phylogenetic relationship of these bacteria to B. elkanii, with B. japonicum diverging prior to the diversification of these taxa. All Bradyrhizobium isolates from Amphicarpaea were also identical to B. elkanii in the size of the intervening sequence (IVS) in the 5' region of the 23S rRNA gene, while B. japonicum had an IVS length variant with 29 additional nucleotides. Parsimony analysis of both 16S and 23S partial rDNA sequences grouped Bradyrhizobium sp. isolates from Amphicarpaea into a clade together with B. elkanii, consistent with the relationships inferred from nod sequences.