Nuclear counterparts of the cytoplasmic mitochondrial 12S rRNA gene: A problem of ancient DNA and molecular phylogenies

Monkey mummy bones and teeth originating from the North Saqqara Baboon Galleries (Egypt), soft tissue from a mummified baboon in a museum collection, and nineteenth/twentieth-century skin fragments from mangabeys were used for DNA extraction and PCR amplification of part of the mitochondrial 12S rRNA gene. Sequences aligning with the 12S rRNA gene were recovered but were only distantly related to contemporary monkey mitochondrial 12S rRNA sequences. However, many of these sequences were identical or closely related to human nuclear DNA sequences resembling mitochondrial 12S rRNA (isolated from a cell line depleted in mitochondria) and therefore have to be considered contamination. Subsequently in a separate study we were able to recover genuine mitochondrial 12S rRNA sequences from many extant species of nonhuman Old World primates and sequences closely resembling the human nuclear integrations. Analysis of all sequences by the neighbor-joining (NJ) method indicated that mitochondrial DNA sequences and their nuclear counterparts can be divided into two distinct clusters. One cluster contained all temporary cytoplasmic mitochondrial DNA sequences and approximately half of the monkey nuclear mitochondriallike sequences. A second cluster contained most human nuclear sequences and the other half of monkey nuclear sequences with a separate branch leading to human and gorilla mitochondrial and nuclear sequences. Sequences recovered from ancient materials were equally divided between the two clusters. These results constitute a warning for when working with ancient DNA or performing phylogenetic analysis using mitochondrial DNA as a target sequence: Nuclear counterparts of mitochondrial genes may lead to faulty interpretation of results.Correspondence to: A.C. van der Kuyl     

Variation of 3'-terminal fragment of 16S rRNA gene in closely related species of Drosophila virilis group

Primary sequence of 3-terminal fragment of the mitochondrial 16S rRNA gene has been determined in 12 Drosophila species of the virilis group. The functionally important elements in secondary structure of the RNA product were defined. The region corresponding to the peptidyltransferase center has been localized. Variation of the 3'-terminal region of 16S rRNA gene has been described in 12 species of the virilis group. Phylogeny of the Drosophila virilis species group is discussed.     

Discordant 16S and 23S rRNA gene phylogenies for the genus Helicobacter: implications for phylogenetic inference and systematics

Analysis of 16S rRNA gene sequences has become the primary method for determining prokaryotic phylogeny. Phylogeny is currently the basis for prokaryotic systematics. Therefore, the validity of 16S rRNA gene-based phylogenetic analyses is of fundamental importance for prokaryotic systematics. Discrepancies between 16S rRNA gene analyses and DNA-DNA hybridization and phenotypic analyses have been noted in the genus Helicobacter. To clarify these discrepancies, we sequenced the 23S rRNA genes for 55 helicobacter strains representing 41 taxa (>2,700 bases per sequence). Phylogenetic-tree construction using neighbor-joining, parsimony, and maximum likelihood methods for 23S rRNA gene sequence data yielded stable trees which were consistent with other phenotypic and genotypic methods. The 16S rRNA gene sequence-derived trees were discordant with the 23S rRNA gene trees and other data. Discrepant 16S rRNA gene sequence data for the helicobacters are consistent with the horizontal transfer of 16S rRNA gene fragments and the creation of mosaic molecules with loss of phylogenetic information. These results suggest that taxonomic decisions must be supported by other phylogenetically informative macromolecules, such as the 23S rRNA gene, when 16S rRNA gene-derived phylogeny is discordant with other credible phenotypic and genotypic methods. This study found Wolinella succinogenes to branch with the unsheathed-flagellum cluster of helicobacters by 23S rRNA gene analyses and whole-genome comparisons. This study also found intervening sequences (IVSs) in the 23S rRNA genes of strains of 12 Helicobacter species. IVSs were found in helices 10, 25, and 45, as well as between helices 31' and 27'. Simultaneous insertion of IVSs at three sites was found in H. mesocricetorum.     

Huxley's line demarcates extensive genetic divergence between eastern and western forms of the giant freshwater prawn, Macrobrachium rosenbergii

Phylogenetic analysis of representatives from 18 wild populations of the giant freshwater prawn, Macrobrachium rosenbergii, utilising a fragment of the 16S rRNA mitochondrial gene, identified two major reciprocally monophyletic clades either side of a well-known biogeographic barrier, Huxley's line. The level of divergence between the two clades (maximum 6.2%) far exceeds divergence levels within either clade (maximum 0.9%), and does not concord with geographical distance among sites. 'Eastern' and 'western' M. rosenbergii clades have probably been separated since Miocene times. Within-clade diversity appears to have been shaped by dispersal events influenced by eustatic change.     

Isolation of cDNA clones coding for mitochondrial 16S ribosomal RNA from the crustacean Artemia

cDNA clones coding for Artemia mitochondrial 16S ribosomal RNA (rRNA) have been isolated. The clones cover from nucleotide 650 of the RNA molecule to its 3' end. The comparison of Artemia sequence with both vertebrate and invertebrate mitochondrial 16S rRNA sequences has shown the existence of regions of high similarity between them. A model for the secondary structure of the 3' half of Artemia mitochondrial 16S rRNA is proposed. The size of the rRNA molecule has been estimated at 1.35 kb. Despite the similarity of the Artemia gene to insect rRNA in size, sequence and secondary structure, the G + C content of the Artemia gene (42%) is closer to that of mammals than to the insect genes. The number of mitochondria in Artemia has been estimated at 1500 per diploid genome in the cyst and 4000 in the nauplius. In contrast, the amount of mt 16S rRNA is constant at all stages of Artemia development.     

Biogeography and phylogeny of the New Zealand cicada genera (Hemiptera: Cicadidae) based on nuclear and mitochondrial DNA data

Aim Determine the geographical and temporal origins of New Zealand cicadas. Location New Zealand, eastern Australia and New Caledonia. Methods DNA sequences from 14 species of cicadas from New Zealand, Australia, and New Caledonia were examined. A total of 4628 bp were analysed from whole genome extraction of four mitochondrial genes (cytochrome oxidase subunits I and II, and ribosomal 12S and 16S subunits) and one nuclear gene (elongation factor‐1 alpha). These DNA sequences were aligned and analysed using standard phylogenetic methods based primarily on the maximum likelihood optimality criterion. Dates of divergences between clades were determined using several molecular clock methods. Results New Zealand cicadas form two well‐defined clades. One clade groups with Australian taxa, the other with New Caledonian taxa. The molecular clock analyses indicate that New Zealand genera diverged from the Australian and New Caledonian genera within the last 11.6 Myr. Main conclusions New Zealand was likely colonized by two or more invasions. One NZ lineage has its closest relatives in Australia and the other in New Caledonia. These invasions occurred well after New Zealand became isolated from other land masses, therefore cicadas must have crossed large bodies of water to reach New Zealand.     

Evolution on a volcanic conveyor belt: using phylogeographic reconstructions and K–Ar-based ages of the Hawaiian Islands to estimate molecular evolutionary rates

The Hawaiian Islands form as the Pacific Plate moves over a 'hot spot' in the earth's mantle where magma extrudes through the crust to build huge shield volcanos. The islands subside and erode as the plate carries them to the north-west, eventually to become coral atolls and seamounts. Thus islands are ordered linearly by age, with the oldest islands in the north-west (e.g. Kauai at 5.1 Ma) and the youngest in the south-east (e.g. Hawaii at 0.43 Ma). K–Ar estimates of the date of an island's formation provide a maximum age for the taxa inhabiting the island. These ages can be used to calibrate rates of molecular change under the following assumptions: (i) K–Ar dates are accurate; (ii) tree topologies show that derivation of taxa parallels the timing of island formation; (iii) populations do not colonize long after island emergence; (iv) the coalescent point for sister taxa does not greatly predate the formation of the colonized younger island; (v) saturation effects and (vi) among-lineage rate variation are minimal or correctable; and (vii) unbiased standard errors of distances and regressions can be estimated from multiple pairwise comparisons. We use the approach to obtain overall corrected rate calibrations for: (i) part of the mitochondrial cytochrome b gene in Hawaiian drepanidines (0.016 sequence divergence/Myr); (ii) the Yp1 gene in Hawaiian Drosophila (0.019/Myr Kambysellis et al. 1995 ); and (iii) parts of the mitochondrial 12S and 16S rRNA and tRNA^val in Laupala crickets (0.024–0.102/Myr, Shaw 1996 ). We discuss the reliability of the estimates given the assumptions (i–vii) above and contrast the results with previous calibrations of Adh in Hawaiian Drosophila and chloroplast DNA in lobeliods.     

Nucleotide Sequence and Gene Organization of the Starfish Asterina Pectinifera Mitochondrial Genome

The 16,260-bp mitochondrial DNA (mtDNA) from the starfish Asterina pectinifera has been sequenced. The genes for 13 proteins, two rRNAs and 22 tRNAs are organized in an extremely economical fashion, similar to those of other animal mtDNAs, with some of the genes overlapping each other. The gene organization is the same as that for another echinoderm, sea urchin, except for the inversion of a 4.6-kb segment that contains genes for two proteins, 13 tRNAs and the 16S rRNA. Judging from the organization of the protein coding genes, mammalian mtDNAs resemble the sea urchin mtDNA more than that of the starfish. The region around the 3' end of the 12S rRNA gene of the starfish shows a high similarity with those for vertebrates. This region encodes a possible stem and loop structure; similar potential structures occur in this region of vertebrate mtDNAs and also in nonmitochondrial small subunit rRNA. A similar stem and loop structure is also found at the 3' end of the 16S rRNA genes in A. pectinifera, in another starfish Pisaster ochraceus, in vertebrates and in Drosophila, but not in sea urchins. The full sequence data confirm the presumption that AGA/AGG, AUA and AAA codons, respectively, code for serine, isoleucine, and asparagine in the starfish mitochondria, and that AGA/AGG codons are read by tRNA(GCU)(Ser), which possesses a truncated dihydrouridine arm, that was previously suggested from a partial mtDNA sequence. The structural characteristics of tRNAs and possible mechanisms for the change in the mitochondrial genetic code are also discussed.     

Mitochondrial phylogeny of the Chrysisignita (Hymenoptera: Chrysididae) species group based on simultaneous Bayesian alignment and phylogeny reconstruction

The ignita species group within the genus Chrysis includes over 100 cuckoo wasp species, which all lead a parasitic lifestyle and exhibit very similar morphology. The lack of robust, diagnostic morphological characters has hindered phylogenetic reconstructions and contributed to frequent misidentification and inconsistent interpretations of species in this group. Therefore, molecular phylogenetic analysis is the most suitable approach for resolving the phylogeny and taxonomy of this group. We present a well-resolved phylogeny of the Chrysis ignita species group based on mitochondrial sequence data from 41 ingroup and six outgroup taxa. Although our emphasis was on European taxa, we included samples from most of the distribution range of the C. ignita species group to test for monophyly. We used a continuous mitochondrial DNA sequence consisting of 16S rRNA, tRNA(Val), 12S rRNA and ND4. The location of the ND4 gene at the 3' end of this continuous sequence, following 12S rRNA, represents a novel mitochondrial gene arrangement for insects. Due to difficulties in aligning rRNA genes, two different Bayesian approaches were employed to reconstruct phylogeny: (1) using a reduced data matrix including only those positions that could be aligned with confidence; or (2) using the full sequence dataset while estimating alignment and phylogeny simultaneously. In addition maximum-parsimony and maximum-likelihood analyses were performed to test the robustness of the Bayesian approaches. Although all approaches yielded trees with similar topology, considerably more nodes were resolved with analyses using the full data matrix. Phylogenetic analysis supported the monophyly of the C. ignita species group and divided its species into well-supported clades. The resultant phylogeny was only partly in accordance with published subgroupings based on morphology. Our results suggest that several taxa currently treated as subspecies or names treated as synonyms may in fact constitute separate species. Our study provides a solid basis for further systematic investigations of this enigmatic insect group.     

Differential susceptibility of tree species to oviposition by periodical cicadas

Abstract 1. Periodical cicadas occur widely in deciduous forests of eastern North America and represent the largest insect emergences in the world. Their effects on forest community dynamics, via oviposition damage and root feeding, are potentially great but not well understood. If tree species vary in their susceptibility to cicada oviposition and root colonization, then periodical cicadas may affect ecological interactions among tree species and forest composition. 2. We measured oviposition damage on 63 tree species at 52 sites in southern Indiana, U.S.A. for two broods of periodical cicadas. We compared oviposition damage among different sites, habitats, and tree species. We also measured cicada emergence rates in mature and successional forest habitats and from beneath different tree species. 3. Certain tree species received significantly more oviposition damage than others, and the most and least damaged species were consistent between the two broods. For some species, oviposition damage varied between mature and successional sites. Despite the differences in oviposition rates among species, there was no difference in emergence rates from beneath different tree species, suggesting that initial differences in cicada density among tree species are not maintained through the nymphal life cycle. 4. Cicadas oviposited more at successional sites, but emerged more in mature forest stands, suggesting that successive generations of cicadas represent a shifting mosaic between forests of different ages. As tree species consistently vary in their susceptibility to oviposition damage, periodical cicadas could alter local competitive interactions and therefore affect forest dynamics.     

Mapping of the RNA recognition site of Escherichia coli ribosomal protein S7

Bacterial ribosomal protein S7 initiates the folding of the 3' major domain of 16S ribosomal RNA by binding to its lower half. The X-ray structure of protein S7 from thermophilic bacteria was recently solved and found to be a modular structure, consisting of an alpha-helical domain with a beta-ribbon extension. To gain further insights into its interaction with rRNA, we cloned the S7 gene from Escherichia coli K12 into a pET expression vector and introduced 4 deletions and 12 amino acid substitutions in the protein sequence. The binding of each mutant to the lower half of the 3' major domain of 16S rRNA was assessed by filtration on nitrocellulose membranes. Deletion of the N-terminal 17 residues or deletion of the B hairpins (residues 72-89) severely decreased S7 affinity for the rRNA. Truncation of the C-terminal portion (residues 138-178), which includes part of the terminal alpha-helix, significantly affected S7 binding, whereas a shorter truncation (residues 148-178) only marginally influenced its binding. Severe effects were also observed with several strategic point mutations located throughout the protein, including Q8A and F17G in the N-terminal region, and K35Q, G54S, K113Q, and M115G in loops connecting the alpha-helices. Our results are consistent with the occurrence of several sites of contact between S7 and the 16S rRNA, in line with its role in the folding of the 3' major domain.     

Drosophila Melanogaster Mitochondrial DNA: Gene Organization and Evolutionary Considerations

The sequence of a 8351-nucleotide mitochondrial DNA (mtDNA) fragment has been obtained extending the knowledge of the Drosophila melanogaster mitochondrial genome to 90% of its coding region. The sequence encodes seven polypeptides, 12 tRNAs and the 3' end of the 16S rRNA and CO III genes. The gene organization is strictly conserved with respect to the Drosophila yakuba mitochondrial genome, and different from that found in mammals and Xenopus. The high A + T content of D. melanogaster mitochondrial DNA is reflected in a reiterative codon usage, with more than 90% of the codons ending in T or A, G + C rich codons being practically absent. The average level of homology between the D. melanogaster and D. yakuba sequences is very high (roughly 94%), although insertion and deletions have been detected in protein, tRNA and large ribosomal genes. The analysis of nucleotide changes reveals a similar frequency for transitions and transversions, and reflects a strong bias against G + C on both strands. The predominant type of transition is strand specific.     

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.     

Molecular structures of mitochondrial-DNA-like sequences in human nuclear DNA.

Two lambda phage clones carrying mitochondrial-DNA-like (mtDNA-like) sequences isolated from a human gene library were named Lm E-1 and Lm C-2, and their DNA structures were characterized. Lm E-1 contains about 0.4 kb DNA homologous to the 5' portion of the mitochondrial 16S ribosomal RNA (rRNA) gene and Lm C-2, a 1.6 kb DNA homologous to the 3' portion of the 12S rRNA gene and to almost all of the 16S rRNA gene. Comparisons of their nucleotide sequences with those of the corresponding regions of the human mtDNA revealed no detectable DNA rearrangement and their homologies to the human mtDNA are 84% and 80%, respectively. There are neither terminal repeats in the nuclear mtDNA-like sequences nor duplications of the nuclear DNAs flanking the mtDNA-like sequences. Evolutionary relationship between these two human nuclear mtDNA-like sequences and the human and bovine mtDNAs is discussed.     

Unraveling evolutionary lineages in the limbless fossorial skink genus Acontias (Sauria: Scincidae): are subspecies equivalent systematic units?

Subspecies in the limbless, endemic African fossorial skink genus Acontias constitute ill-defined operational taxonomic units, consequently considerable systematic debate has lingered on the systematic diversity within Acontias. In the present study, the systematic affinities among acontine taxa are explored with the utility of partial sequence data from two mitochondrial gene loci (16S rRNA and cytochrome oxidase subunit 1 (COI)) for all taxa, while two additional loci (12S rRNA, cytochrome b) were used to investigate relationships within the Acontias meleagris complex. Phylogenetic results, derived from the combined analysis, revealed two monophyletic clades. Clade 1 is comprised of small-bodied skinks while clade 2 comprised the medium bodied skinks. Within clade 2 none of the traditionally recognized subspecies formed reciprocally monophyletic groups. Furthermore, constraining the topology and enforcing sister taxa relationships between the assumed subspecies, consistently recovered a topology that was statistically significant worse, indicating that the traditionally designated subspecies groupings probably represent invalid taxonomic units, thus clearly reflecting considerable discord with current taxonomy. The burrowing life style of these lizards has probably led to marked convergent evolution and constrained the development of diagnostic morphological characters among these species. Morphological similarities in color as well as scale architecture within Acontias are labile and highly homoplaseous and do not reflect the evolutionary history of the group. Taxonomic implications of these results are discussed.