Replication slippage may cause parallel evolution in the secondary structures of mitochondrial transfer RNAs

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

The ultrastructural features of the infundibulum of the green iguana, Iguana iguana

The purpose of this study is to describe, in detail, the ultrastructure of the infundibulum of the sexually mature and active female green iguana, Iguana iguana. The infundibulum of five iguanas was remarkably distinct from the uterus, and was also clearly demarcated into cranial (expanded v-shaped) and caudal (tubular) divisions. Tissue samples obtained from five portions (three from the cranial division and two from the caudal division) of the infundibulum were processed conventionally for light and electron microscopy. The epithelial lining of the most anterior, middle, and posterior, parts of the cranial division displayed nonciliated cells predominantly, and occasionally ciliated cells. The numerous secretory granules in nonciliated type 1 cell found in the fimbrial aspect of the infundibulum were homogenous and deeply electron-dense, but those in the other two regions were variants of this cell type because they contained variably electron-dense secretory granules. Two main types of nonciliated cells (type 2 and its variant, type 3, as well as type 4) occurred in the epithelial lining of the caudal division of the infundibulum, but they, clearly, showed no dense secretory granules. Whereas the nonciliated type 2 cell and its variant (type 3 cell) contained large glycogen deposits, the type 4 cell lacked these deposits but its apical part contained large lipid-like droplets and, remarkably, blebbed into the duct lumen. The nonciliated cells lining the mucosal tubular glands contained highly electron-dense secretory granules, which were similar to those found in the nonciliated type 1 cell in the epithelial lining of the fimbrial part of the cranial division of the infundibulum.

Rolling Circles as a Means of Encoding Genes in the RNA World

The rolling circle mechanism found in viroids and some RNA viruses is a likely way that replication could have begun in the RNA World. Here, we consider simulations of populations of protocells, each containing multiple copies of rolling circle RNAs that can replicate non-enzymatically. The mechanism requires the presence of short self-cleaving ribozymes such as hammerheads, which can cleave and re-circularize RNA strands. A rolling circle must encode a hammerhead and the complement of a hammerhead, so that both plus and minus strands can cleave. Thus, the minimal functional length is twice the length of the hammerhead sequence. Selection for speed of replication will tend to reduce circles to this minimum length. However, if sequence errors occur when copying the hammerhead sequence, this prevents cleavage at one point, but still allows cleavage on the next passage around the rolling circle. Thus, there is a natural doubling mechanism that creates strands that are multiple times the length of the minimal sequence. This can provide space for the origin of new genes with beneficial functions. We show that if a beneficial gene appears in this new space, the longer sequence with the beneficial function can be selected, even though it replicates more slowly. This provides a route for the evolution of longer circles encoding multiple genes.

Describing novel mitochondrial genomes of Antarctic amphipods

To date, only one mitogenome from an Antarctic amphipod has been published. Here, novel complete mitochondrial genomes (mitogenomes) of two morphospecies are assembled, namely, Charcotia amundseni and Eusirus giganteus. For the latter species, we have assembled two mitogenomes from different genetic clades of this species. The lengths of Eusirus and Charcotia mitogenomes range from 15,534 to 15,619 base pairs and their mitogenomes are composed of 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and 1 putative control region CR. Some tRNAs display aberrant structures suggesting that minimalization is also ongoing in amphipod mitogenomes. The novel mitogenomes of the two Antarctic species have features distinguishing them from other amphipod mitogenomes such as a lower AT-richness in the whole mitogenomes and a negative GC- skew in both strands of protein coding genes. The genetically most variable mitochondrial regions of amphipods are nad6 and atp8, while cox1 shows low nucleotide diversity among closely and more distantly related species. In comparison to the pancrustacean mitochondrial ground pattern, E. giganteus shows a translocation of the nad1 gene, while cytb and nad6 genes are translocated in C. amundseni. Phylogenetic analysis based on mitogenomes illustrates that Eusirus and Charcotia cluster together with other species belonging to the same amphipod superfamilies. In the absence of reference nuclear genomes, mitogenomes can be useful to develop markers for studying population genetics or evolutionary relationships at higher taxonomic levels.

Evidence of two deeply divergent co-existing mitochondrial genomes in the Tuatara reveals an extremely complex genomic organization

Animal mitochondrial genomic polymorphism occurs as low-level mitochondrial heteroplasmy and deeply divergent co-existing molecules. The latter is rare, known only in bivalvian mollusks. Here we show two deeply divergent co-existing mt-genomes in a vertebrate through genomic sequencing of the Tuatara (Sphenodon punctatus), the sole-representative of an ancient reptilian Order. The two molecules, revealed using a combination of short-read and long-read sequencing technologies, differ by 10.4% nucleotide divergence. A single long-read covers an entire mt-molecule for both strands. Phylogenetic analyses suggest a 7-8 million-year divergence between genomes. Contrary to earlier reports, all 37 genes typical of animal mitochondria, with drastic gene rearrangements, are confirmed for both mt-genomes. Also unique to vertebrates, concerted evolution drives three near-identical putative Control Region non-coding blocks. Evidence of positive selection at sites linked to metabolically important transmembrane regions of encoded proteins suggests these two mt-genomes may confer an adaptive advantage for an unusually cold-tolerant reptile.

Acanthosaura aurantiacrista (Squamata: Agamidae), a new long horn lizard from northern Thailand

BACKGROUND

In Thailand, five species of Acanthosaura have been recorded so far, including Acanthosaura armata from the southern region, A. cardamomensis from the eastern region, A. crucigera from the western region, A. lepidogaster from the northern region and A. phuketensis from the Phuket Island and south-western region. However, comprehensive studies of diversity patterns and distribution of Acanthosaura are still lacking in some areas and need further information for designating areas of special conservation importance and nature protection planning in Thailand.

NEW INFORMATION

Acanthosaura aurantiacrista is a new species of long-horned lizard of the genus Acanthosaura from northern Thailand. It is distinguished from all other species of Acanthosaura by a dagger-like nuchal spine with yellowish-orange colouration in females, bright yellow colouration in males and a combination of other morphological characters: a greater tail length to snout-vent length ratio; a larger postorbital spine, nuchal spine, dorsal spine and occipital spine compared to its head length; a smaller diastema to snout-vent length ratio; a greater number of subdigital lamellae on the fourth finger and fourth toe; and a larger gular pouch than other Acanthosaura species. Analysis of mitochondrial ND2 gene sequences revealed a sister clade between the A. aurantiacrista lineage and the A. crucigera lineage with a 100% probability of divergence, according to Bayesian analysis and strong support value for Maximum Likelihood analysis. The pairwise distance ranged from 13.8-15.0% between A. aurantiacrista and A. cardamomensis, 10.9-14.5% between A. aurantiacrista and A. crucigera and 0-1.2% amongst A. aurantiacrista populations. The discovery of this lizard increases the known endemic herpetological diversity and underscores the importance of conservation in the mountain rainforest region of northern Thailand.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora)

BACKGROUND

Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages.

RESULTS

Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses.

CONCLUSIONS

In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

An Ecologist’s Guide to Mitochondrial DNA Mutations and Senescence

Longevity plays a key role in the fitness of organisms, so understanding the processes that underlie variance in senescence has long been a focus of ecologists and evolutionary biologists. For decades, the performance and ultimate decline of mitochondria have been implicated in the demise of somatic tissue, but exactly why mitochondrial function declines as individual’s age has remained elusive. A possible source of decline that has been of intense debate is mutations to the mitochondrial DNA. There are two primary sources of such mutations: oxidative damage, which is widely discussed by ecologists interested in aging, and mitochondrial replication error, which is less familiar to most ecologists. The goal of this review is to introduce ecologists and evolutionary biologists to the concept of mitochondrial replication error and to review the current status of research on the relative importance of replication error in senescence. We conclude by detailing some of the gaps in our knowledge that currently make it difficult to deduce the relative importance of replication error in wild populations and encourage organismal biologists to consider this variable both when interpreting their results and as viable measure to include in their studies.

Evolution of mitochondrial genomes in Baikalian amphipods

BACKGROUND

Amphipods (Crustacea) of Lake Baikal are a very numerous and diverse group of invertebrates generally believed to have originated by adaptive radiation. The evolutionary history and phylogenetic relationships in Baikalian amphipods still remain poorly understood. Sequencing of mitochondrial genomes is a relatively feasible way for obtaining a set of gene sequences suitable for robust phylogenetic inferences. The architecture of mitochondrial genomes also may provide additional information on the mechanisms of evolution of amphipods in Lake Baikal.

RESULTS

Three complete and four nearly complete mitochondrial genomes of Baikalian amphipods were obtained by high-throughput sequencing using the Illumina platform. A phylogenetic inference based on the nucleotide sequences of all mitochondrial protein coding genes revealed the Baikalian species to be a monophyletic group relative to the nearest non-Baikalian species with a completely sequenced mitochondrial genome - Gammarus duebeni. The phylogeny of Baikalian amphipods also suggests that the shallow-water species Eulimnogammarus has likely evolved from a deep-water ancestor, however many other species have to be added to the analysis to test this hypothesis. The gene order in all mitochondrial genomes of studied Baikalian amphipods differs from the pancrustacean ground pattern. Mitochondrial genomes of four species possess 23 tRNA genes, and in three genomes the extra tRNA gene copies have likely undergone remolding. Widely varying lengths of putative control regions and other intergenic spacers are typical for the mitochondrial genomes of Baikalian amphipods.

CONCLUSIONS

The mitochondrial genomes of Baikalian amphipods display varying organization suggesting an intense rearrangement process during their evolution. Comparison of complete mitochondrial genomes is a potent approach for studying the amphipod evolution in Lake Baikal.

The complete mitochondrial genome of the geophilomorph centipede Strigamia maritima

Strigamia maritima (Myriapoda; Chilopoda) is a species from the soil-living order of geophilomorph centipedes. The Geophilomorpha is the most speciose order of centipedes with over a 1000 species described. They are notable for their large number of appendage bearing segments and are being used as a laboratory model to study the embryological process of segmentation within the myriapods. Using a scaffold derived from the recently published genome of Strigamia maritima that contained multiple mitochondrial protein-coding genes, here we report the complete mitochondrial genome of Strigamia, the first from any geophilomorph centipede. The mitochondrial genome of S. maritima is a circular molecule of 14,938 base pairs, within which we could identify the typical mitochondrial genome complement of 13 protein-coding genes and 2 ribosomal RNA genes. Sequences resembling 16 of the 22 transfer RNA genes typical of metazoan mitochondrial genomes could be identified, many of which have clear deviations from the standard ‘cloverleaf’ secondary structures of tRNA. Phylogenetic trees derived from the concatenated alignment of protein-coding genes of S. maritima and >50 other metazoans were unable to resolve the Myriapoda as monophyletic, but did support a monophyletic group of chilopods: Strigamia was resolved as the sister group of the scolopendromorph Scolopocryptos sp. and these two (Geophilomorpha and Scolopendromorpha), along with the Lithobiomorpha, formed a monophyletic group the Pleurostigmomorpha. Gene order within the S. maritima mitochondrial genome is unique compared to any other arthropod or metazoan mitochondrial genome to which it has been compared. The highly unusual organisation of the mitochondrial genome of Strigamia maritima is in striking contrast with the conservatively evolving nuclear genome: sampling of more members of this order of centipedes will be required to see whether this unusual organization is typical of the Geophilomorpha or results from a more recent reorganisation in the lineage leading to Strigamia.

Mendelian-mutationism: the forgotten evolutionary synthesis

According to a classical narrative, early geneticists, failing to see how Mendelism provides the missing pieces of Darwin's theory, rejected gradual changes and advocated an implausible yet briefly popular view of evolution-by-mutation; after decades of delay (in which synthesis was prevented by personal conflicts, disciplinary rivalries, and anti-Darwinian animus), Darwinism emerged on a new Mendelian basis. Based on the works of four influential early geneticists - Bateson, de Vries, Morgan and Punnett -, and drawing on recent scholarship, we offer an alternative that turns the classical view on its head. For early geneticists, embracing discrete inheritance and the mutation theory (for the origin of hereditary variation) did not entail rejection of selection, but rejection of Darwin's non-Mendelian views of heredity and variation, his doctrine of naturanon facitsaltum, and his conception of "natural selection" as a creative force that shapes features out of masses of infinitesimal differences. We find no evidence of a delay in synthesizing mutation, rules of discrete inheritance, and selection in a Mendelian-Mutationist Synthesis. Instead, before 1918, early geneticists had conceptualized allelic selection, the Hardy-Weinberg equilibrium, the evolution of a quantitative trait under selection, the probability of fixation of a new mutation, and other key innovations. Contemporary evolutionary thinking seems closer to their more ecumenical view than to the restrictive mid-twentieth-century consensus known as the Modern Synthesis.

A comparative study of nemertean complete mitochondrial genomes, including two new ones for Nectonemertes cf. mirabilis and Zygeupolia rubens, may elucidate the fundamental pattern for the phylum Nemertea

BACKGROUND

The mitochondrial genome is important for studying genome evolution as well as reconstructing the phylogeny of organisms. Complete mitochondrial genome sequences have been reported for more than 2200 metazoans, mainly vertebrates and arthropods. To date, from a total of about 1275 described nemertean species, only three complete and two partial mitochondrial DNA sequences from nemerteans have been published. Here, we report the entire mitochondrial genomes for two more nemertean species: Nectonemertes cf. mirabilis and Zygeupolia rubens.

RESULTS

The sizes of the entire mitochondrial genomes are 15365 bp for N. cf. mirabilis and 15513 bp for Z. rubens. Each circular genome contains 37 genes and an AT-rich non-coding region, and overall nucleotide composition is AT-rich. In both species, there is significant strand asymmetry in the distribution of nucleotides, with the coding strand being richer in T than A and in G than C. The AT-rich non-coding regions of the two genomes have some repeat sequences and stem-loop structures, both of which may be associated with the initiation of replication or transcription. The 22 tRNAs show variable substitution patterns in nemerteans, with higher sequence conservation in genes located on the H strand. Gene arrangement of N. cf. mirabilis is identical to that of Paranemertes cf. peregrina, both of which are Hoplonemertea, while that of Z. rubens is the same as in Lineus viridis, both of which are Heteronemertea. Comparison of the gene arrangements and phylogenomic analysis based on concatenated nucleotide sequences of the 12 mitochondrial protein-coding genes revealed that species with closer relationships share more identical gene blocks.

CONCLUSION

The two new mitochondrial genomes share many features, including gene contents, with other known nemertean mitochondrial genomes. The tRNA families display a composite substitution pathway. Gene order comparison to the proposed ground pattern of Bilateria and some lophotrochozoans suggests that the nemertean ancestral mitochondrial gene order most closely resembles the heteronemertean type. Phylogenetic analysis proposes a sister-group relationship between Hetero- and Hoplonemertea, which supports one of two recent alternative hypotheses of nemertean phylogeny.

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.

Evolutionary origins and diversification of dragon lizards in Australia's tropical savannas

Australia's monsoonal tropics are dominated by the largest and least modified savanna woodlands in the world, and they are globally significant for their high biodiversity and regional endemism. Despite this, there have been very few molecular studies of the evolutionary origins and diversification of vertebrates in this region. The semi-arboreal dragon lizards of Lophognathus and Amphibolurus are widely distributed in the savanna and dry sclerophyll woodlands of Australasia, including the monsoon tropics. We sequenced a ~1400 bp region of mitochondrial DNA and a ~1400 bp nuclear gene (RAG1) to investigate the phylogenetic relationships and phylogeographic structuring of all seven species of Lophognathus and Amphibolurus. Our analyses show that there is a higher level of species and generic diversity in the monsoon tropics than previously thought, and a full morphological review and taxonomic revision of these genera is required. Relaxed molecular clock analyses indicate that species across both genera originated in the late Miocene and early Pliocene, with significant phylogeographic structure within species. We did not find any evidence that the monsoon tropics species were a monophyletic group that had diversified within the region; instead Amphibolurus and Lophognathus represent at least three independent evolutionary colonizations of the monsoon tropics. It is probable that the origins and phylogeographic patterns of the northern Lophognathus species have evolved under the climatic influence of the Australian monsoon, rather than being either an ancient Gondwanan lineage that pre-dates the monsoon or the result of a more recent dispersal event across Wallace's Line.

The mitochondrial genomes of two nemerteans, Cephalothrix sp. (Nemertea: Palaeonemertea) and Paranemertes cf. peregrina (Nemertea: Hoplonemertea)

The mitochondrial genome sequences were determined for two species of nemerteans, Cephalothrix sp. (15,800 bp sequenced, near-complete) and Paranemertes cf. peregrina (14,558 bp, complete). As seen in most metazoans, the genomes encode 13 protein, 2 ribosomal RNA and 22 transfer RNA genes. The nucleotide composition is strongly biased toward A and T, as is typical for metazoan mtDNAs. There is also a significant strand skew in the distribution of these nucleotides, with the coding strand being richer in T than A and in G than C. All genes are transcribed in the same direction except for trnP and trnT, which is consistent with that reported for Cephalothrix hongkongiensis and Lineus viridis. Gene arrangement of Cephalothrix sp. is identical to that of C. hongkongiensis, while in P. cf. peregrina it is similar to L. viridis, but differs significantly from the three Cephalothrix species in the position of four protein-coding genes and seven tRNAs. Some protein-coding genes have 3' end stem-loop structures, which may allow mRNA processing without flanking tRNAs. The major non-coding regions observed in the two genomes with potential to form stem-loop structures may be involved in the initiation of replication or transcription. The average Ka/Ks values, varying from 0.12 to 0.89, are markedly different among the 13 mitochondrial protein-coding genes, suggesting that there may exist different selective pressure among mitochondrial genes of nemerteans.

ND6 gene "lost" and found: evolution of mitochondrial gene rearrangement in Antarctic notothenioids

Evolution of Antarctic notothenioids in the frigid and oxygen-rich Southern Ocean had led to remarkable genomic changes, most notably the gain of novel antifreeze glycoproteins and the loss of oxygen-binding hemoproteins in the icefish family. Recently, the mitochondrial (mt) NADH dehydrogenase subunit 6 (ND6) gene and the adjacent transfer RNA(Glu) (tRNA(Glu)) were also reportedly lost. ND6 protein is crucial for the assembly and function of Complex I of the mt electron transport chain that produces adenosine triphosphate (ATP) essential for life; thus, ND6 absence would be irreconcilable with Antarctic notothenioids being thriving species. Here we report our discovery that the ND6 gene and tRNA(Glu) were not lost but had been translocated to the control region (CR) from their canonical location between ND5 and cytochrome b genes. We characterized the CR and adjacent sequences of 22 notothenioid species representing all eight families of Notothenioidei to elucidate the mechanism and evolutionary history of this mtDNA rearrangement. Species of the three basal non-Antarctic families have the canonical vertebrate mt gene order, whereas species of all five Antarctic families have a rearranged CR bearing the embedded ND6 (ND6(CR)) and tRNA(Glu), with additional copies of tRNA(Thr), tRNA(Pro), and noncoding region in various lineages. We hypothesized that an initial duplication of the canonical mt region from ND6 through CR occurred in the common ancestor to the Antarctic clade, and we deduced the succession of loss or modification of the duplicated region leading to the extant patterns of mt DNA reorganization that is consistent with notothenioid evolutionary history. We verified that the ND6(CR) gene in Antarctic notothenioids is transcribed and therefore functional. However, ND6(CR)-encoded protein sequences differ substantially from basal non-Antarctic notothenioid ND6, and we detected lineage-specific positive selection on the branch leading to the Antarctic clade of ND6(CR) under the branch-site model. Collectively, the novel mt ND6(CR) genotype of the Antarctic radiation represents another major molecular change in Antarctic notothenioid evolution and may reflect an adaptive change conducive to the functioning of the protein (Complex I) machinery of mt respiration in the polar environment, driven by the advent of freezing, oxygen-rich conditions in the Southern Ocean.

Phylogeography of the Pygmy Rain Frog (Pristimantis ridens) across the lowland wet forests of isthmian Central America

We used a phylogeographic approach to elucidate the evolutionary history of a lineage of frogs, known as Pristimantis (formerly Eleutherodactylus) ridens (Anura: Brachycephalidae), restricted to the wet forests occurring along the Caribbean versant of isthmian Central America as well as the disjunct wet forest on the Pacific slope of Costa Rica. We placed our phylogeographic study of P. ridens within a larger molecular phylogenetic analysis of Central American Pristimantis. All phylogenetic inferences were based on a 1455 base pair fragment of mitochondrial DNA, containing the complete ND2 gene and five flanking tRNA genes. Our reconstruction of the intraspecific phylogeny of P. ridens yielded a basal trichotomy dating to an estimated 12+ million years ago (Ma), consisting of central Panama, western Panama, and Costa Rica plus Honduras. Thus, the presence of P. ridens appears to predate the completion of the Isthmus 3.1Ma. Using a parametric bootstrap (SOWH) test, we evaluated four a priori zoogeographic hypotheses for the origin and spread of P. ridens. This analysis suggested that the P. ridens populations on the Caribbean versant of Costa Rica were established by Pacific versant ancestors only recently, in contrast to the very old lineages found in Panama. Our results support a model of Miocene colonization, long-term geographic stasis, followed by rapid dispersal across the Caribbean lowlands during the Pliocene or Pleistocene.

Phylogenetic relationships of the heath dragons (Rankinia adelaidensis and R. parviceps) from the south-western Australian biodiversity hotspot

Although the south-western Australian region is recognised as a global biodiversity hotspot, there are still significant gaps in our understanding of the biodiversity of this region. We present a phylogenetic study of the heath dragons (Rankinia adelaidensis and R. parviceps) from this region, incorporating a 1612-bp section of mtDNA and two nuclear introns, Gapdh (~244 bp) and Enol (~330 bp). In addition, we present a generic-level analysis of three gene regions (mtDNA, Gapdh, BDNF), which provides clear evidence that Rankinia adelaidensis and R. parviceps are not closely related to Rankinia diemensis from eastern Australia. Instead, the heath dragons are strongly supported as forming a clade with the genus Ctenophorus. In addition, we find that there are significant levels of haplotype divergence between currently recognised subspecies of the heath dragons (R. a. adelaidensis, R. a. chapmani, R. p. parviceps, R. p. butleri). We suggest that the genetic divergences between subspecies result from geographic isolation in allopatry owing to habitat preferences, followed by drift and/or selection. On the basis of these deep divergences and consistent morphological differences between subspecies, we recommend elevating all taxa to full species, and provide a taxonomic revision of the genera Rankinia and Ctenophorus.

A manually curated database of tetrapod mitochondrially encoded tRNA sequences and secondary structures

BACKGROUND

Mitochondrial tRNAs have been the subject of study for structural biologists interested in their secondary structure characteristics, evolutionary biologists have researched patterns of compensatory and structural evolution and medical studies have been directed towards understanding the basis of human disease. However, an up to date, manually curated database of mitochondrially encoded tRNAs from higher animals is currently not available.

DESCRIPTION

We obtained the complete mitochondrial sequence for 277 tetrapod species from GenBank and re-annotated all of the tRNAs based on a multiple alignment of each tRNA gene and secondary structure prediction made independently for each tRNA. The mitochondrial (mt) tRNA sequences and the secondary structure based multiple alignments are freely available as Supplemental Information online.

CONCLUSION

We compiled a manually curated database of mitochondrially encoded tRNAs from tetrapods with completely sequenced genomes. In the course of our work, we reannotated more than 10% of all tetrapod mt-tRNAs and subsequently predicted the secondary structures of 6060 mitochondrial tRNAs. This carefully constructed database can be utilized to enhance our knowledge in several different fields including the evolution of mt-tRNA secondary structure and prediction of pathogenic mt-tRNA mutations. In addition, researchers reporting novel mitochondrial genome sequences should check their tRNA gene annotations against our database to ensure a higher level of fidelity of their annotation.

The mitochondrial genome of the screamer louse Bothriometopus (phthiraptera: ischnocera): effects of extensive gene rearrangements on the evolution of the genome

Mitochondrial (mt) genome rearrangement has generally been studied with respect to the phenomenon itself, focusing on their phylogenetic distribution and causal mechanisms. Rearrangements have additional significance through effects on substitution, transcription, and mRNA processing. Lice are an ideal group in which to study the interactions between rearrangements and these factors due to the heightened rearrangement rate within this group. The entire mt genome of the screamer louse Bothriometopus was sequenced and compared to previously sequenced louse genomes. The mt genome is 15,564 bp, circular, and all genes are encoded on the same strand. The gene arrangement differs radically from both other louse species and the ancestral insect. Nucleotide composition is A+T biased, but there is no skew which may be due to reversal of replication direction or a transcriptional effect. Bothriometopus has both tRNA duplication and concerted evolution which has not been observed previously. Eleven of the 13 protein-coding genes have 3' end stem-loop structures which may allow mRNA processing without flanking tRNAs and so facilitate gene rearrangements. There are five candidate control regions capable of forming stem-loop structures. Two are structurally more similar to the control regions of other insect species than those of other lice. Analyses of Bothriometopus demonstrate that louse mt genomes, in addition to being extensively rearranged, differ significantly from most insect species in nucleotide composition biases, tRNA evolution, protein-coding gene structures and putative signaling sites such as the control region. These may be either a cause or a consequence of gene rearrangements.

A phylogenetic and biogeographic perspective on the evolution of poeciliid fishes

Phylogenetic relationships of members of the subfamily Poeciliinae (Cyprinodontiformes) are investigated to test alternate hypotheses of diversification resulting from the assembly of the Central America and the Caribbean from the Cretaceous period onwards. We use 4333 aligned base pairs of mitochondrial DNA and 1549 aligned base pairs of nuclear DNA from 55 samples representing 48 ingroup and seven outgroup species to test this hypothesis. Mitochondrial genes analyzed include those encoding the 12S and 16S ribosomal RNAs; transfer RNAs coding for valine, leucine, isoleucine, glutamine, methionine, tryptophan, alanine, asparagine, cysteine and tyrosine; and complete cytochrome b and NADH dehydrogenase subunit I and II; nuclear gene analyzed included the third exon of the recombination activation gene 1 (RAG1). Analyses of combined mtDNA and nuclear DNA data sets result in a well-supported phylogenetic hypothesis. This hypothesis is in conflict with the classical taxonomic assignment of genera into tribes and phylogenetic hypotheses based on the taxonomy; however, the molecular hypothesis defines nine clades that are geographically restricted and consistent with the geological evolution of Central America and the Caribbean. Our analyses support multiple colonization events of Middle America followed by a mix of vicariance and dispersal events.

Mutation-biased adaptation in a protein NK model

Evolutionary trends responsible for systematic differences in genome and proteome composition have been attributed to GC:AT mutation bias in the context of neutral evolution or to selection acting on genome composition. A possibility that has been ignored, presumably because it is part of neither the Modern Synthesis nor the Neutral Theory, is that mutation may impose a directional bias on adaptation. This possibility is explored here with simulations of the effect of a GC:AT bias on amino acid composition during adaptive walks on an abstract protein fitness landscape called an "NK" model. The results indicate that adaptation does not preclude mutation-biased evolution. In the complete absence of neutral evolution, a modest GC:AT bias of realistic magnitude can displace the trajectory of adaptation in a mutationally favored direction, to such a degree that amino acid composition is biased substantially and persistently. Thus, mutational explanations for evolved patterns need not presuppose neutral evolution.

Resolving deep phylogenetic relationships in salamanders: analyses of mitochondrial and nuclear genomic data

Phylogenetic relationships among salamander families illustrate analytical challenges inherent to inferring phylogenies in which terminal branches are temporally very long relative to internal branches. We present new mitochondrial DNA sequences, approximately 2,100 base pairs from the genes encoding ND1, ND2, COI, and the intervening tRNA genes for 34 species representing all 10 salamander families, to examine these relationships. Parsimony analysis of these mtDNA sequences supports monophyly of all families except Proteidae, but yields a tree largely unresolved with respect to interfamilial relationships and the phylogenetic positions of the proteid genera Necturus and Proteus. In contrast, Bayesian and maximum-likelihood analyses of the mtDNA data produce a topology concordant with phylogenetic results from nuclear-encoded rRNA sequences, and they statistically reject monophyly of the internally fertilizing salamanders, suborder Salamandroidea. Phylogenetic simulations based on our mitochondrial DNA sequences reveal that Bayesian analyses outperform parsimony in reconstructing short branches located deep in the phylogenetic history of a taxon. However, phylogenetic conflicts between our results and a recent analysis of nuclear RAG-1 gene sequences suggest that statistical rejection of a monophyletic Salamandroidea by Bayesian analyses of our mitochondrial genomic data is probably erroneous. Bayesian and likelihood-based analyses may overestimate phylogenetic precision when estimating short branches located deep in a phylogeny from data showing substitutional saturation; an analysis of nucleotide substitutions indicates that these methods may be overly sensitive to a relatively small number of sites that show substitutions judged uncommon by the favored evolutionary model.

Molecular Mechanisms of Extensive Mitochondrial Gene Rearrangement in Plethodontid Salamanders

Extensive gene rearrangement is reported in the mitochondrial genomes of lungless salamanders (Plethodontidae). In each genome with a novel gene order, there is evidence that the rearrangement was mediated by duplication of part of the mitochondrial genome, including the presence of both pseudogenes and additional, presumably functional, copies of duplicated genes. All rearrangement-mediating duplications include either the origin of light-strand replication and the nearby tRNA genes or the regions flanking the origin of heavy-strand replication. The latter regions comprise nad6, trnE, cob, trnT, an intergenic spacer between trnT and trnP and, in some genomes, trnP, the control region, trnF, rrnS, trnV, rrnL, trnL1, and nad1. In some cases, two copies of duplicated genes, presumptive regulatory regions, and/or sequences with no assignable function have been retained in the genome following the initial duplication; in other genomes, only one of the duplicated copies has been retained. Both tandem and nontandem duplications are present in these genomes, suggesting different duplication mechanisms. In some of these mitochondrial DNAs, up to 25% of the total length is composed of tandem duplications of noncoding sequence that includes putative regulatory regions and/or pseudogenes of tRNAs and protein-coding genes along with the otherwise unassignable sequences. These data indicate that imprecise initiation and termination of replication, slipped-strand mispairing, and intramolecular recombination may all have played a role in generating repeats during the evolutionary history of plethodontid mitochondrial genomes.

The mitochondrial sequences of Heptathela hangzhouensis and Ornithoctonus huwena reveal unique gene arrangements and atypical tRNAs

We have sequenced the complete mitochondrial genomes of the spiders Heptathela hangzhouensis and Ornithoctonus huwena. Both genomes encode 13 protein-coding genes, 22 tRNA genes, and 2 ribosomal RNA genes. H. hangzhouensis, a species of the suborder Mesothelae and a representative of the most basal clade of Araneae, possesses a gene order identical to that of Limulus polyphemus of Xiphosura. On the other hand, O. huwena, a representative of suborder Opisthothelae, infraorder Mygalomorphae, was found to have seven tRNA genes positioned differently from those of Limulus. The rrnL-trnL1-nad1 arrangement shared by the araneomorph families Salticidae, Nesticidae, and Linyphiidae and the mygalomorph family Theraphosidae is a putative synapomorphy joining the mygalomorph with the araneomorph. Between the two species examined, base compositions also differ significantly. The lengths of most protein-coding genes in H. hangzhouensis and O. huwena mtDNA are either identical to or slightly shorter than their Limulus counterparts. Usage of initiation and termination codons in these protein-coding genes seems to follow patterns conserved among most arthropod and some other metazoan mitochondrial genomes. The sequences of the 3' ends of rrnS and rrnL in the two species are similar to those reported for Limulus, and the entire genes are shortened by about 100-250 nucleotides with respect to Limulus. The lengths of most tRNA genes from the two species are distinctly shorter than those of Limulus and the sequences reveal unusual inferred tRNA secondary structures. Our finding provides new molecular evidence supporting that the suborder Mesothelae is basal to opisthothelids.

Cenozoic biogeography and evolution in direct-developing frogs of Central America (Leptodactylidae: Eleutherodactylus) as inferred from a phylogenetic analysis of nuclear and mitochondrial genes

We report the first phylogenetic analysis of DNA sequence data for the Central American component of the genus Eleutherodactylus (Anura: Leptodactylidae: Eleutherodactylinae), one of the most ubiquitous, diverse, and abundant components of the Neotropical amphibian fauna. We obtained DNA sequence data from 55 specimens representing 45 species. Sampling was focused on Central America, but also included Bolivia, Brazil, Jamaica, and the USA. We sequenced 1460 contiguous base pairs (bp) of the mitochondrial genome containing ND2 and five neighboring tRNA genes, plus 1300 bp of the c-myc nuclear gene. The resulting phylogenetic inferences were broadly concordant between data sets and among analytical methods. The subgenus Craugastor is monophyletic and its initial radiation was potentially rapid and adaptive. Within Craugastor, the earliest splits separate three northern Central American species groups, milesi, augusti, and alfredi, from a clade comprising the rest of Craugastor. Within the latter clade, the rhodopis group as formerly recognized comprises three deeply divergent clades that do not form a monophyletic group; we therefore restrict the content of the rhodopis group to one of two northern clades, and use new names for the other northern (mexicanus group) and one southern clade (bransfordii group). The new rhodopis and bransfordii groups together form the sister taxon to a clade comprising the biporcatus, fitzingeri, mexicanus, and rugulosus groups. We used a Bayesian MCMC approach together with geological and biogeographic assumptions to estimate divergence times from the combined DNA sequence data. Our results corroborated three independent dispersal events for the origins of Central American Eleutherodactylus: (1) an ancestor of Craugastor entered northern Central America from South American in the early Paleocene, (2) an ancestor of the subgenus Syrrhophus entered northern Central America from the Caribbean at the end of the Eocene, and (3) a wave of independent dispersal events from South America coincided with formation of the Isthmus of Panama during the Pliocene. We elevate the subgenus Craugastor to the genus rank.

The mitochondrial genome organization of the rice frog, Fejervarya limnocharis (Amphibia: Anura): a new gene order in the vertebrate mtDNA

The mitochondrial DNA of the rice frog, Fejervarya limnocharis (Amphibia, Anura), was obtained using long-and-accurate polymerase chain reaction (LA-PCR) combining with subcloning method. The complete nucleotide sequence (17,717 bp) of mitochondrial genome was determined subsequently. This mitochondrial genome is characterized by four distinctive features: the translocation of ND5 gene, a cluster of rearranged tRNA genes (tRNAThr, tRNAPro, tRNA(Leu CUN)), a tandem duplication of tRNAMet gene, and eight large 89-bp tandem repeats in the control region, as well as three short noncoding regions containing two repeated motifs existing in the gene cluster of ND5/tRNAThr/tRNAPro/tRNALeu/tRNAPhe. The tandem duplication of gene regions followed by deletions of supernumerary genes can be invoked to explain the shuffling of tRNAMet and a cluster of tRNA and ND5 genes, as observed in this study. Both ND5 gene translocation and tandem duplication of tRNAMet were first observed in the vertebrate mitochondrial genomes.

Sequencing and comparing whole mitochondrial genomes of animals

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, which can be especially powerful. We describe here the protocols commonly used for physically isolating mitochondrial DNA (mtDNA), for amplifying these by polymerase chain reaction (PCR) or rolling circle amplification (RCA), for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences with determining and comparing complete mitochondrial DNA sequences.

Phylogenetic relationships among amphisbaenian reptiles based on complete mitochondrial genomic sequences

Complete mitochondrial genomic sequences are reported from 12 members in the four families of the reptile group Amphisbaenia. Analysis of 11,946 aligned nucleotide positions (5797 informative) produces a robust phylogenetic hypothesis. The family Rhineuridae is basal and Bipedidae is the sister taxon to the Amphisbaenidae plus Trogonophidae. Amphisbaenian reptiles are surprisingly old, predating the breakup of Pangaea 200 million years before present, because successive basal taxa (Rhineuridae and Bipedidae) are situated in tectonic regions of Laurasia and nested taxa (Amphisbaenidae and Trogonophidae) are found in Gondwanan regions. Thorough sampling within the Bipedidae shows that it is not tectonic movement of Baja California away from the Mexican mainland that is primary in isolating Bipes species, but rather that primary vicariance occurred between northern and southern groups. Amphisbaenian families show parallel reduction in number of limbs and Bipes species exhibit parallel reduction in number of digits. A measure is developed for comparing the phylogenetic information content of various genes. A synapomorphic trait defining the Bipedidae is a shift from the typical vertebrate mitochondrial gene arrangement to the derived state of trnE and nad6. In addition, a tandem duplication of trnT and trnP is observed in Bipes biporus with a pattern of pseudogene formation that varies among populations. The first case of convergent rearrangement of the mitochondrial genome among animals demonstrated by complete genomic sequences is reported. Relative to most vertebrates, the Rhineuridae has the block nad6, trnE switched in order with the block cob, trnT, trnP, as they are in birds.

Molecular phylogenetics and mitochondrial genomic evolution in the chamaeleonidae (Reptilia, Squamata)

A phylogenetic hypothesis for the lizard family Chamaeleonidae is generated from 1503 aligned base positions (883 parsimony-informative) of mitochondrial DNA for specimens representing 59 species (57 ingroup and two outgroup). Sequences are reported for a genomic segment encoding eight transfer RNAs, NADH dehydrogenase component 2 (ND2), and portions of NADH dehydrogenase component 1 (ND1) and cytochrome c oxidase subunit 1 (COI). Newly reported genomic rearrangements and duplications support the hypothesis that mitochondrial gene order and content are destabilized by phylogenetic loss of a functional origin for light-strand replication between the genes encoding tRNA(Asn) and tRNA(Cys). A novel gene order characterizes all sampled Brookesia except B. nasus. Brookesia nasus, the apparent sister taxon of a clade formed by all other Brookesia, has the ancestral gene order but contains a large tandem duplication. An apparently noncoding 220 base pair insertion between the genes encoding ND2 and tRNA(Trp) is reported for Bradypodion tavetanum. Phylogenetic analysis identifies nine clades whose ancestral lineages diverged early in chamaeleonid evolutionary history: (1) Brookesia (possibly excluding B. nasus), (2) Chamaeleo subgenus Chamaeleo (excluding C. namaquensis), (3) Chamaeleo subgenus Trioceros, (4) viviparous Bradypodion, (5) oviparous Bradypodion, (6) genus Furcifer (except F. balteatus), and (7-9) three distinct clades of Calumma. Chamaeleo namaquensis, Brookesia nasus, Furcifer balteatus, Rhampholeon brevicaudatus, and R. spectrum represent ancient lineages dating to approximately the same time. Multiple independent losses and a possible secondary gain of horns are inferred for Trioceros. Viviparity has at least two separate origins in chameleons, one in Bradypodion and

A molecular phylogenetic study of ecological diversification in the Australian lizard genus Ctenophorus

We present phylogenetic analyses of the lizard genus Ctenophorus using 1,639 aligned positions of mitochondrial DNA sequences containing 799 parsimony-informative characters for samples of 22 species of Ctenophorus and 12 additional Australian agamid genera. Sequences from three protein-coding genes (ND1, ND2, and COI) and eight intervening tRNA genes are examined using both parsimony and maximum-likelihood analyses. Species of Ctenophorus form a monophyletic group with Rankinia adelaidensis, which we suggest placing in Ctenophorus. Ecological differentiation among species of Ctenophorus is most evident in the kinds of habitats used for shelter. Phylogenetic analyses suggest that the ancestral condition is to use burrows for shelter, and that habits of sheltering in rocks and shrubs/hummock grasses represent separately derived conditions. Ctenophorus appears to have undergone extensive cladogenesis approximately 10-12 million years ago, with all three major ecological modes being established at that time.

Mitochondrial DNA Evidence and Evolution in Varanoidea (Squamata)

Varanoidea is a monophyletic group of anguimorph lizards, comprising the New World helodermatids, the Bornean earless monitor Lanthanotus borneensis, and the Old World monitors (Varanus). I use mitochondrial DNA sequences and extensive taxonomic sampling to test alternative hypotheses of varanoid relationships. The most parsimonious hypothesis confirms the monophyly of Varanoidea (Heloderma, Lanthanotus, and Varanus) and Varanus, as well as the sister-taxon relationship of Varanus and Lanthanotus. The relationships among Varanus species differ in several respects from previous hypotheses. Three major lineages are recognized within Varanus: an African clade basal to the rest of the group, an Indo-Asian clade, and an Indo-Australian clade. Within the last lineage, the endemic Australian dwarf monitors (Odatria) form a clade sister to the large Australian monitors (the gouldii group). Tests of the effects of rate heterogeneity and homoplasy demonstrate that putative process partitions of data are largely congruent with one another and contribute positive support to the overall hypothesis.

Molecular phylogenetics of western North American frogs of the Rana boylii species group

Phylogenetic relationships among frogs of the genus Rana from western North America are investigated using 2013 aligned bases of mitochondrial DNA sequence from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNA(Ile), tRNA(Gln), tRNA(Met), ND2, tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr), and COI (subunit I of cytochrome c oxidase), plus the origin for light-strand replication (O(L)) between the tRNA(Asn) and tRNA(Cys) genes. The aligned sequences contain 401 phylogenetically informative characters. A well-resolved phylogenetic hypothesis in which the Rana boylii species group (R. aurora, R. boylii, R. cascadae, R. muscosa, and R. pretiosa) is monophyletic is obtained. Molecular sequence divergence suggests that the R. boylii species group is approximately 8 million years old. The traditional hypothesis showing monophyly of the yellow-legged frogs (R. boylii and R. muscosa) is statistically rejected in favor of a hypothesis in which R. aurora, R. cascadae, and R. muscosa form a clade. Reanalyses of published nuclear ribosomal DNA restriction-site data and allozymic data support a monophyletic R. boylii group, but do not effectively resolve relationships among species within this group. Eight populations of R. muscosa form two major clades separated by a biogeographic break in the Sierra Nevada of California. This biogeographic break is broadly concordant with breaks found in four other amphibian and reptilian taxa. The two major clades within R. muscosa are estimated to have diverged approximately 2.2 million years before present. Each of these major clades contains two subgroups showing approximately 1.5 million years divergence, implicating climatic effects of Pleistocene glaciation in vicariance. The four distinct subgroups of R. muscosa separated by at least 1.4 million years of evolutionary divergence are suggested as potential units for conservation.

Bias in the introduction of variation as an orienting factor in evolution

According to New Synthesis doctrine, the direction of evolution is determined by selection and not by "internal causes" that act by way of propensities of variation. This doctrine rests on the theoretical claim that because mutation rates are small in comparison to selection coefficients, mutation is powerless to overcome opposing selection. Using a simple population-genetic model, this claim is shown to depend on assuming the prior availability of variation, so that mutation may act only as a "pressure" on the frequencies of existing alleles, and not as the evolutionary process that introduces novelty. As shown here, mutational bias in the introduction of novelty can strongly influence the course of evolution, even when mutation rates are small in comparison to selection coefficients. Recognizing this mode of causation provides a distinct mechanistic basis for an "internalist" approach to determining the contribution of mutational and developmental factors to evolutionary phenomena such as homoplasy, parallelism, and directionality.

Molecular phylogenetics and historical biogeography among salamandrids of the "true" salamander clade: rapid branching of numerous highly divergent lineages in Mertensiella luschani associated with the rise of Anatolia

Phylogenetic relationships among salamandrids of the "true" salamander clade are investigated using 2019 aligned base positions (713 parsimony informative) of 20 mitochondrial DNA sequences from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNA(Ile), tRNA(Gln), tRNA(Met), ND2, tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr), and COI (subunit I of cytochrome c oxidase), plus the origin for light-strand replication (O(L)) between the tRNA(Asn) and the tRNA(Cys) genes. Parsimony analysis produces a robust phylogenetic estimate for the relationships of the major groups of "true" salamanders. Strong support is provided for the sister taxon relationship of Chioglossa and Mertensiella caucasica and for the placement of Salamandra and Mertensiella luschani as sister taxa. These relationships suggest two vicariant events between Europe and Anatolia caused by the formation of seaways in the Mediterranean Basin. Molecular divergence indicates an Early Miocene separation of Chioglossa and M. caucasica and a Late Miocene separation of Salamandra and M. luschani. The traditional phylogenetic hypothesis of a monophyletic Mertensiella is statistically rejected, indicating that southwestern and northeastern Anatolian populations have separate historical biogeographic origins. Therefore, we recommend placement of M. luschani in the genus Salamandra. Within M. luschani, six highly divergent lineages showing 7.6 to 10.1% pairwise sequence divergence are identified. Tests using four-taxon subsamples suggest that these lineages diverged nearly simultaneously in the Late Miocene, approximately 6 to 8 million years ago, when extensive uplifting of Anatolia occurred in response to the Arabian collision.

Mitochondrial sequence evolution in spiders: intraspecific variation in tRNAs lacking the TPsiC Arm

Analyses of mitochondrial DNA sequences from three species of Habronattus jumping spiders (Chelicerata: Arachnida: Araneae) reveal unusual inferred tRNA secondary structures and gene arrangements, providing new information on tRNA evolution within chelicerate arthropods. Sequences from the protein-coding genes NADH dehydrogenase subunit 1 (ND1), cytochrome oxidase subunit I (COI), and subunit II (COII) were obtained, along with tRNA, tRNA, and large-subunit ribosomal RNA (16S) sequences; these revealed several peculiar features. First, inferred secondary structures of tRNA and, likely, tRNA, lack the TPsiC arm and the variable arm and therefore do not form standard cloverleaf structures. In place of these arms is a 5-6-nt T arm-variable loop (TV) replacement loop such as that originally described from nematode mitochondrial tRNAs. Intraspecific variation occurs in the acceptor stem sequences in both tRNAs. Second, while the proposed secondary structure of the 3' end of 16S is similar to that reported for insects, the sequence at the 5' end is extremely divergent, and the entire gene is truncated about 300 nt with respect to Drosophila yakuba. Third, initiation codons appear to consist of ATY (ATT and ATC) and TTG for ND1 and COII, respectively. Finally, Habronattus shares the same ND1-tRNA-16S gene arrangement as insects and crustaceans, thus illustrating variation in a tRNA gene arrangement previously proposed as a character distinguishing chelicerates from insects and crustaceans.

Complete sequence of the mitochondrial DNA of the primitive opisthobranch gastropod Pupa strigosa: systematic implication of the genome organization

The complete sequence (14,189 bp) of the mitochondrial DNA of the opisthobranch gastropod Pupa strigosa was determined. The genome contains 13 protein, 2 rRNA, and 22 tRNA genes typical of metazoan mtDNA. The Pupa mitochondrial genome is highly compact and shows the following unusual features, like pulmonate land snails: (1) extremely small genome size, (2) absence of lengthy noncoding regions (with the largest intergenic spacer being only 46 nt), (3) size reduction of encoded genes, and (4) many overlapping genes. Several tRNA genes exhibit bizarre secondary structures with reduced T or D stems, and many tRNA genes have unstable acceptor stems that might be corrected by posttranscriptional RNA editing. The Pupa mitochondrial gene arrangement is almost identical to those of pulmonate land snails but is radically divergent from those of the prosobranch gastropod Littorina saxatilis and other molluscs. Our finding that the unique gene arrangement and highly compact genome organization are shared between opisthobranch and pulmonate gastropods strongly suggests their close phylogenetic affinity.