Mitochondrial whims: metabolic rate,longevity and the rate of molecular evolution

The evolutionary rate of mitochondrial DNA (mtDNA) is highly variable across lineages in animals, and particularly in mammals. This variation has been interpreted as reflecting variations in metabolic rate: mitochondrial respiratory activity would tend to generate mutagenic agents, thus increasing the mutation rate. Here we review recent evidence suggesting that a direct, mechanical effect of species metabolic rate on mtDNA evolutionary rate is unlikely. We suggest that natural selection could act to reduce the (somatic) mtDNA mutation rate in long-lived species, in agreement with the mitochondrial theory of ageing.     

Diversity and phylogeny of mitochondrial cytochrome B lineages from six morphospecies of avian Haemoproteus (Haemosporida: Haemoproteidae)

Species of Haemoproteus (Haemosporida: Haemoproteidae), avian haemosporidians, have traditionally been described based on morphology of their gametocytes and on limited experimental information on their vertebrate host specificity. We investigated to what extent the morphological species are represented by monophyletic groups based on DNA sequence data using 2 different fragment lengths of the cytochrome b (cyt. b) gene. Phylogenetic reconstructions of obtained cyt. b lineages from 6 morphospecies of Haemoproteus showed that all lineages formed monophyletic clusters matching the morphospecies. Comparing our data with a recently published study showed that this is not always the case; the morphospecies H. belopolskyi consists of 2 distinct clusters of lineages that apparently have converged in morphology. However, the overall broad congruence between the molecular and morphological clustering of lineages will facilitate the integration of the knowledge obtained by traditional and molecular parasitology. Mean between morphospecies variation was 10-fold higher than the within species variation (5.5% vs. 0.54%), suggesting that Haemoproteus lineages with a genetic differentiation >5% are expected to be morphologically differentiated in most cases. When investigate the utility of 2 different fragment sizes of the cyt. b gene, the partial, 479-bp, cyt. b protocol picked up all mitochondrial (mt)DNA lineages that are found when using the full cyt. b gene, 1073 bp, suggesting that this protocol is sufficient for identification of most mtDNA lineages. All of the mtDNA lineages were associated with unique alleles when amplification was possible at a nuclear locus, strengthening the hypothesis that the designation of lineages based on mtDNA is largely genome-wide representative. We, therefore, propose the use of a cyt. b fragment of this length as a standard gene fragment for a DNA bar-coding system for avian Haemoproteus species.     

A mitochondrial phylogeography of Brachidontes variabilis (Bivalvia: Mytilidae) reveals three cryptic species

This study examined genetic variation across the range of Brachidontes variabilis to produce a molecular phylogeography. Neighbour joining (NJ), minimum evolution (ME) and maximum parsimony (MP) trees based on partial mitochondrial DNA sequences of 16S-rDNA and cytochrome oxidase (COI) genes revealed three monophyletic clades: (1) Brachidontes pharaonis s.l. from the Mediterranean Sea and the Red Sea; (2) B. variabilis from the Indian Ocean; (3) B. variabilis from the western Pacific Ocean. Although the three clades have never been differentiated by malacologists employing conventional morphological keys, they should be ascribed to the taxonomic rank of species. The nucleotide divergences between Brachidontes lineages (between 10.3% and 23.2%) were substantially higher than the divergence between congeneric Mytilus species (2.3–6.7%) and corresponded to interspecific divergences found in other bivalvia, indicating that they should be considered three different species. Analysis of the 16S-rDNA sequences revealed heteroplasmy, indicating dual uniparental inheritance (DUI) of mtDNA in the species of Brachidontes collected in the Indian Ocean, but not in the species in the Pacific nor in the species in the Red Sea and the Mediterranean Sea. When we employed the conventional estimate of the rate of mitochondrial sequence divergence (2% per million years), the divergence times for the three monophyletic lineages were 6–11 Myr for the Indian Ocean and Pacific Ocean Brachidontes sp. and 6.5–9 Myr for the Red Sea and Indian Ocean Brachidontes sp . Thus, these species diverged from one another during the Miocene (23.8–5.3 Myr). We infer that a common ancestor of the three Brachidontes species probably had an Indo-Pacific distribution and that vicariance events, linked to Pleistocene glaciations first and then to the opening of the Red Sea, produced three monophyletic lineages.     

Inverse relationship between longevity and evolutionary rate of mitochondrial proteins in mammals and birds

Recently, an unexpected, positive correlation between the rate of evolution of mitochondrial proteins and longevity was reported. Here we re-analyze this relationship in various mammalian lineages using a bayesian phylogenetic analysis of amino-acid sequences, allowing for variable evolutionary rates across sites and species. A negative relationship between protein evolutionary rate and species longevity is reported for all oxidative phosphorylation complexes. A detailed analysis of the cytochrome b in 528 mammals reinforced this result, which contradicts previous publications. Reconducting the analysis in birds yielded similar results. We explain the discrepancy between this and previous reports by our improved taxon sampling and more appropriate methodology: unlike distance-based methods, the tree-based bayesian approach can take into account the high variation of substitution rate across amino-acid sites, and the resulting multiple substitution events. We discuss how our analysis contradicts Rottenberg’s rationale, but does not dismiss his proposal of a longevity-dependent selective pressure on mitochondrial mutation rate in mammals and birds. This is because his interpretation invokes adaptation as the single evolutionary force at work, disregarding the effects of mutation, genetic drift, and purifying selection.     

Mitochondrial DNA evolution at a turtle's pace: evidence for low genetic variability and reduced microevolutionary rate in the Testudines.

Evidence is compiled suggesting a slowdown in mean microevolutionary rate for turtle mitochondrial DNA (mtDNA). Within each of six species or species complexes of Testudines, representing six genera and three taxonomic families, sequence divergence estimates derived from restriction assays are consistently lower than expectations based on either (a) the dates of particular geographic barriers with which significant mtDNA genetic clades appear associated or (b) the magnitudes of sequence divergence between mtDNA clades in nonturtle species that otherwise exhibit striking phylogeographic concordance with the genetic partitions in turtles. Magnitudes of the inferred rate slowdowns average eightfold relative to the "conventional" mtDNA clock calibration of 2%/Myr sequence divergence between higher animal lineages. Reasons for the postulated deceleration remain unknown, but two intriguing correlates are (a) the exceptionally long generation length most turtles and (b) turtles' low metabolic rate. Both factors have been suspected of influencing evolutionary rates in the DNA sequences of some other vertebrate groups. Uncertainities about the dates of cladogenetic events in these Testudines leave room for alternatives to the slowdown interpretation, but consistency in the direction of the inferred pattern, across several turtle species and evolutionary settings, suggests the need for caution in acceptance of a universal mtDNA-clock calibration for higher animals.     

High mutation rates in the mitochondrial genomes of Daphnia pulex

Despite the great utility of mitochondrial DNA (mtDNA) sequence data in population genetics and phylogenetics, key parameters describing the process of mitochondrial mutation (e.g., the rate and spectrum of mutational change) are based on few direct estimates. Furthermore, the variation in the mtDNA mutation process within species or between lineages with contrasting reproductive strategies remains poorly understood. In this study, we directly estimate the mtDNA mutation rate and spectrum using Daphnia pulex mutation-accumulation (MA) lines derived from sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages. The nearly complete mitochondrial genome sequences of 82 sexual and 47 asexual MA lines reveal high mtDNA mutation rate of 1.37 × 10(-7) and 1.73 × 10(-7) per nucleotide per generation, respectively. The Daphnia mtDNA mutation rate is among the highest in eukaryotes, and its spectrum is dominated by insertions and deletions (70%), largely due to the presence of mutational hotspots at homopolymeric nucleotide stretches. Maximum likelihood estimates of the Daphnia mitochondrial effective population size reveal that between five and ten copies of mitochondrial genomes are transmitted per female per generation. Comparison between sexual and asexual lineages reveals no statistically different mutation rates and highly similar mutation spectra.     

Mitochondrial DNA rearrangements are correlated with a delayed amplification of the mobile intron (plDNA) in a long-lived mutant of Podospora anserina

A new long-lived mutant of Podospora anserina has been isolated and characterized. Its longevity is maternally inherited as revealed by reciprocal crosses. A molecular analysis resulted in the identification of an amplified DNA species (designated pAL2-1) with homology to mitochondrial DNA (mtDNA). The presence of this DNA species is correlated with mtDNA rearrangements and a delayed amplification of the mobile intron (plDNA).     

Mitochondrial diversity of Opsariichthys bidens (Teleostei, Cyprinidae) in three Chinese drainages

We describe the phylogeographic structure of 28 Chinese populations of the cyprinid Opsariichthys bidens across three main Chinese river drainages. Our study is based on the phylogenetic analysis of the complete mitochondrial cytochrome b gene (1140 bp). We combined this analysis with population processes inferred from nested clade analysis (NCA) and mismatch distributions. Both analyses showed that Chinese O. bidens consists of five mtDNA lineages (Opsariichthys 1-5) with high genetic divergence among them. Molecular divergences (TrN+G) higher than 20% among the Opsariichthys 1-5 mtDNA lineages suggest a taxonomic underestimation at the species level. About 92% of the genetic variance among samples was explained by differences among Opsariichthys mtDNA lineages. Drainage-restricted haplotypes with high frequencies and moderate nucleotide diversity show that Opsariichthys populations have evolved independently. NCA results were congruent with the phylogeny, and unimodal mismatch distributions with negative Tajima's D values suggest population expansions in some Opsariichthys lineages. The phylogeographic structure of the Opsariichthys 1-5 mtDNA lineages appears to be related to their long-term interruption of gene flow (theta(ST)>0.97). Our results suggested that fragmentation of ancestral ranges might have caused Opsariichthys diversification in Chinese waters. However, current distribution of common haplotypes across the Yangtze and Pearl drainages suggests a recent river connection that could have favoured gene flow across drainages. Overall, the results indicated that the richness of current Asian widespread species might have been underestimated, and that the cyprinid populations of O. bidens in the Yangtze, Pearl and Hai He drainages may correspond to five species.     

Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.     

A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.     

Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics

Mitochondrial DNA sequences are often used to construct molecular phylogenetic trees among closely related animals. In order to examine the usefulness of mtDNA sequences for deep-branch phylogenetics, genes in previously reported mtDNA sequences were analyzed among several animals that diverged 20–600 million years ago. Unambiguous alignment was achieved for stem-forming regions of mitochondrial tRNA genes by virtue of their conservative secondary structures. Sequences derived from stem parts of the mitochondrial tRNA genes appeared to accumulate much variation linearly for a long period of time: nearly 100 Myr for transition differences and more than 350 Myr for transversion differences. This characteristic could be attributed, in part, to the structural variability of mitochondrial tRNAs, which have fewer restrictions on their tertiary structure than do nonmitochondrial tRNAs. The tRNA sequence data served to reconstruct a well-established phylogeny of the animals with 100% bootstrap probabilities by both maximum parsimony and neighbor joining methods. By contrast, mitochondrial protein genes coding for cytochrome b and cytochrome oxidase subunit I did not reconstruct the established phylogeny or did so only weakly, although a variety of fractions of the protein gene sequences were subjected to tree-building. This discouraging phylogenetic performance of mitochondrial protein genes, especially with respect to branches originating over 300 Myr ago, was not simply due to high randomness in the data. It may have been due to the relative susceptibility of the protein genes to natural selection as compared with the stem parts of mitochondrial tRNA genes. On the basis of these results, it is proposed that mitochondrial tRNA genes may be useful in resolving deep branches in animal phylogenies with divergences that occurred some hundreds of Myr ago. For this purpose, we designed a set of primers with which mtDNA fragments encompassing clustered tRNA genes were successfully amplified from various vertebrates by the polymerase chain reaction.Abbreviations AA stem amino acid-acceptor stem - AC stem anticodon stem - COI cytochrome oxidase subunit I - cytb cytochrome b - D stem dihydrouridine stem - MP maximum parsimony - mtDNA mitochondrial DNA - Myr million years - NJ neighbor joining - PCR polymerase chain reaction - Ti transition - T stem tC stem - Tv transversion Correspondence to: Y. Kumazawa     

Molecular phylogenetics and biogeography of transisthmian and amphi-Atlantic needlefishes (Belonidae: Strongylura and Tylosurus): perspectives on New World marine speciation

Phylogenetic relationships among New World and eastern Atlantic species in the belonid genera Strongylura and Tylosurus were hypothesized using 3689bp of nucleotide sequence; including the entire mitochondrial (mtDNA) ATP synthase 6 and 8 genes; partial cytochrome b; 12S and 16S ribosomal genes; and introns and exons, 2 and 3 of the nuclear-encoded creatine kinase B gene. Concordant mtDNA and nuclear genealogies permitted well-supported inference of species relationships within Strongylura and Tylosurus, and of the chronology of diversification in the two genera. Our phylogenetic hypothesis permitted an assessment of Rosen's [Syst. Zool. 24 (1975) 431] model of species diversification across the eastern Atlantic to eastern Pacific marine biogeographic track. The spatial predictions of the Rosen model were generally supported, but not the temporal predictions. Furthermore, long branches leading to terminal Belonidae indicated that many species have persisted for millions of years or that nucleotide substitution rates were elevated for some clades. Though heterogeneity of nucleotide substitution rate was indicated across some belonid lineages, molecular clock estimates were used to hypothesize biogeographic scenarios for Strongylura across the eastern Pacific and Atlantic region. Furthermore, use of a molecular clock indicated; that early diversification among contemporary Strongylura may have been initiated by changes in Atlantic Ocean circulation precipitated by closure of the Tethys Sea; and provided approximate dates for the isolation of the freshwater species on the American continents.     

Deep sympatric mitochondrial divergence without reproductive isolation in the common redstart Phoenicurus phoenicurus

Mitochondrial DNA usually shows low sequence variation within and high sequence divergence among species, which makes it a useful marker for phylogenetic inference and DNA barcoding. A previous study on the common redstart (Phoenicurus phoenicurus) revealed two very different mtDNA haplogroups (5% K2P distance). This divergence is comparable to that among many sister species; however, both haplogroups coexist and interbreed in Europe today. Herein, we describe the phylogeographic pattern of these lineages and test hypotheses for how such high diversity in mtDNA has evolved. We found no evidence for mitochondrial pseudogenes confirming that both haplotypes are of mitochondrial origin. When testing for possible reproductive barriers, we found no evidence for lineage‐specific assortative mating and no difference in sperm morphology, indicating that they are not examples of cryptic species, nor likely to reflect the early stages of speciation. A gene tree based on a short fragment of cytochrome c oxidase subunit 1 from the common redstart and 10 other Phoenicurus species, showed no introgression from any of the extant congenerics. However, introgression from an extinct congeneric cannot be excluded. Sequences from two nuclear introns did not show a similar differentiation into two distinct groups. Mismatch distributions indicated that the lineages have undergone similar demographic changes. Taken together, these results confirm that deeply divergent mitochondrial lineages can coexist in biological species. Sympatric mtDNA divergences are relatively rare in birds, but the fact that they occur argues against the use of threshold mtDNA divergences in species delineation.     

Phylogeography of the lacertid lizard, Psammodromus algirus, in Iberia and across the Strait of Gibraltar

Aim  To determine genetic substructuring within the lacertid lizard Psammodromus algirus . To compare levels of variation across a geological barrier, the Strait of Gibraltar, and to compare this against the known age of the barrier using a molecular clock hypothesis. To compare the effect of the barrier within this species with previously published data from other organisms.
Location  The Iberian Peninsula and North Africa.
Methods  Partial sequences from the mitochondrial cytochrome b , 12S rRNA and 16S rRNA genes were obtained from 101 specimens belonging to the subfamily Gallotiinae and used in this study. The data set was aligned using C lustal X and phylogenetic trees produced using both maximum-parsimony and maximum-likelihood methods. Maximum likelihood estimates of divergence times for the combined data set (12S + 16S + cytochrome b ) were obtained after discovery of lineage rate constancy across the tree using a likelihood ratio test.
Results  Psammodromus algirus contains divergent eastern and western mtDNA clades within the Iberian Peninsula. The western clade has northern and southern lineages in Iberia and one in North Africa. This phylogeographical pattern indicates that the lizard invaded North Africa after the opening of the Strait, presumably by natural rafting.
Main conclusions  As in several other species, current patterns of genetic diversity within P. algirus are not directly related to the opening of the Strait of Gibraltar. Widespread sampling on both sides of the barrier is necessary to determine its effect on species in this area accurately.     

Evolutionary history of the genus Rhamdia (Teleostei: Pimelodidae) in Central America

We constructed phylogenetic hypotheses for Mesoamerican Rhamdia, the only genus of primary freshwater fish represented by sympatric species across Central America. Phylogenetic relationships were inferred from analysis of 1990 base pairs (bp) of mitochondrial DNA (mtDNA), represented by the complete nucleotide sequences of the cytochrome b (cyt b) and the ATP synthase 8 and 6 (ATPase 8/6) genes. We sequenced 120 individuals from 53 drainages to provide a comprehensive geographic picture of Central American Rhamdia systematics and phylogeography. Phylogeographic analysis distinguished multiple Rhamdia mtDNA lineages, and the geographic congruence across evolutionarily independent Rhamdia clades indicated that vicariance has played a strong role in the Mesoamerican diversification of this genus. Phylogenetic analyses of species-level relationships provide strong support for the monophyly of a trans-Andean clade of three evolutionarily equivalent Rhamdia taxa: R. guatemalensis, R. laticauda, and R. cinerascens. Application of fish-based mitochondrial DNA clocks ticking at 1.3-1.5% sequence divergence per million years (Ma), suggests that the split between cis- and trans-Andean Rhamdia extends back about 8 Ma, and the three distinct trans-Andean Rhamdia clades split about 6 Ma ago. Thus the mtDNA divergence observed between cis- and trans-Andean Rhamdia species is too low to support an ancient colonization of Central America in the Late Cretaceous or Paleocene as had been hypothesized in one colonization model for Mesoamerican fishes. Rather the mtDNA data indicate that Rhamdia most likely colonized Central America in the late Miocene or Pliocene, promoting a strong role for the Isthmus of Panamá in the Mesoamerican expansion of this genus. Basal polytomies suggest that both the R. laticauda and R. guatemalensis clades spread rapidly across the Central American landscape, but differences in the average mtDNA genetic distances among clades comprising the two species, indicate that the R. laticauda spread and diversified across Mesoamerica about 1 million years before R. guatemalensis.     

Mitochondrial rate variation among lineages of passerine birds

The order Passeriformes comprises the majority of extant avian species. Analyses of molecular data have provided important insights into the evolution of this diverse order. However, molecular estimates of the evolutionary and demographic timescales of passerine species have been hindered by a lack of reliable calibrations. This has led to a reliance on the application of standard substitution rates to mitochondrial DNA data, particularly rates estimated from analyses of the gene encoding cytochrome b (CYTB). To investigate patterns of rate variation across passerine lineages, we used a Bayesian phylogenetic approach to analyse the protein‐coding genes of 183 mitochondrial genomes. We found that the most commonly used mitochondrial marker, CYTB, has low variation in rates across passerine lineages. This lends support to its widespread use as a molecular clock in birds. However, we also found that the patterns of among‐lineage rate variation in CYTB are only weakly related to the evolutionary rate of the mitochondrial genome as a whole. Our analyses confirmed the presence of mutational saturation at third codon positions across the protein‐coding genes of the mitochondrial genome, reinforcing the view that these sites should be excluded in studies of deep passerine relationships. The results of our analyses have provided information that will be useful for molecular‐clock studies of passerine evolution.     

Molecular evolutionary relationships in the avian genus Anthus (Pipits: Motacillidae)

Nucleotide sequences for 1035 bp of the mitochondrial cytochrome b gene were used to determine the molecular evolutionary relationships of species in the cosmopolitan avian genus Anthus. Phylogenetic analysis of these mtDNA sequences supported four major clades within the genus: (1) the small-bodied African pipits, (2) a largely Palearctic clade, (3) a largely South American clade, and (4) an African-Eurasian-Australian clade. Anthus hellmayri, A. correndera, and A. rubescens are shown to be paraphyletic. The possibility of paraphyly within A. similis is instead inferred to be the discovery of a new species and supported by reference to the museum voucher specimen. Sequence divergence suggests a Pliocene/Miocene origin for the genus. Although Anthus cytochrome b is found not to be behaving in a clocklike fashion across all taxa, speciation during the Pleistocene epoch can be reasonably inferred for the 66% of sister pairs that are diverging in a clocklike manner. Base compositions at each codon position are similar to those found across a growing number of avian lineages. The resulting phylogenetic hypothesis is compared to previous hypotheses of Anthus relationships, all of which deal with relationships of a particular species or a particular species complex; roughly half of these previous hypotheses are supported.     

Variability of nucleotide sequences of the mitochondrial DNA cytochrome c gene in dolly varden and taranetz char

Nucleotide sequence of the 307-bp fragment of the mitochondrial DNA cytochrome b gene was determined in representatives of the three species of the Salvelinus genus, specifically, dolly varden char (S. malma), taranetz char (S. taranetzi), and white-spotted char (S. leucomaenis). These results pointed to a high level of mitochondrial DNA (mtDNA) divergence between white-spotted char and dolly varden char, on the one hand, and taranetz char, on the other (the mean d value was 5.45%). However, the divergence between the dolly varden char and taranetz char was only 0.81%, which is comparable with the level of intraspecific divergence in the dolly varden char (d = 0.87%). It was shown that the dolly varden char mitochondrial gene pool contained DNA lineages differing from the main mtDNA pool at least in the taranetz char-specific mitochondrial lineages. One of these dolly varden char mtDNA lineages was characterized by the presence of the restriction endonuclease MspI-D variant of the cytochrome b gene. This lineage was widely distributed in the Chukotka populations but it was not detected in the Yana River (Okhotsk sea) populations. These findings suggest that dolly varden char has a more ancient evolutionary lineage, diverging from the common ancestor earlier than did taranetz char.     

Pleistocene effects on the European freshwater fish fauna: double origin of the cobitid genus Sabanejewia in the Danube basin (Osteichthyes: Cobitidae)

Biogeographical hypotheses of European freshwater fishes were inferred using phylogeographic analysis of the complete cytochrome b and ATP synthase 8 and 6 mitochondrial genes (1982bp). To test the relative importance of drainage origin versus Pleistocene glaciations in the origin of primary freshwater fishes in Europe, we reconstructed the phylogenetic relationships of the genus Sabanejewia which is distributed in European waters. The phylogenetic relationships recovered for the genus Sabanejewia (n=75) provide support for the monophyly of six main evolutionary mtDNA lineages: Sabanejewia larvata, Sabanejewia romanica, Sabanejewia aurata/Sabanejewia caucasica, Sabanejewia kubanica, Sabanejewia baltica, and the Danubian-Balkanian complex. The Caucasian-Caspian mtDNA lineages, S. kubanica, S. aurata/S. caucasica, and the Northern European S. baltica represents the sister group of the Danubian-Balkanian complex mtDNAclade, supporting a Caucasian-Northern European origin of most of mtDNA lineages of the Central European freshwater fish fauna. The mtDNA divergence observed between the Danubian Sabanejewia species is too dissimilar to support their contemporary origin. Rather, the mtDNA data suggest that the Danubian Sabanejewia lineages most likely have a double origin, indicating that the European Sabanejewia lineages have experienced different historical processes for the following reasons. First, the origin of the S. larvata and S. romanica mtDNA clades predates the origin of the Danubian-Balkanian complex, and our results showed that the completion of the Alps and the origin of the Danube drainage seem to have promoted the speciation of the earliest Sabanejewia clades in the Miocene. Second, small genetic distances and the geographical pattern found within the Danubian-Balkanian complex clade indicate that the lineages included in this clade spread recently across the Danube and Greek river drainages. The inclusion of the S. balcanica species within all mtDNA lineages suggests that cyclical cold periods during the Pleistocene glaciations have favoured its rapid expansion and genetic homogenisation across Central European and Greek waters.     

Functional respiratory chain analyses in murid xenomitochondrial cybrids expose coevolutionary constraints of cytochrome b and nuclear subunits of complex III

The large number of extant Muridae species provides the opportunity of investigating functional limits of nuclear/mitochondrial respiratory chain (RC) subunit interactions by introducing mitochondrial genomes from progressively more divergent species into Mus musculus domesticus mtDNA-less (rho0) cells. We created a panel of such xenomitochondrial cybrids, using as mitochondrial donors cells from six murid species with divergence from M. m. domesticus estimated at 2 to 12 Myr before present. Species used were Mus spretus, Mus caroli, Mus dunni, Mus pahari, Otomys irroratus, and Rattus norvegicus. Parsimony analysis of partial mtDNA sequences showed agreement with previous molecular phylogenies, with the exception that Otomys did not nest within the murinae as suggested by some recent nuclear gene analyses. Cellular production of lactate, a sensitive indicator of decreased respiratory chain ATP production, correlated with divergence. Functional characterization of the chimeric RC complexes in isolated mitochondria using enzymological analyses demonstrated varying decreases in activities of complexes I, III, and IV, which have subunits encoded in both mitochondrial and nuclear genomes. Complex III showed a striking decline in electron transfer function in the most divergent xenocybrids, being greatly reduced in the Rattus xenocybrid and virtually absent in the Otomys xenocybrid. This suggests that nuclear subunits interacting with cytochrome b face the greatest constraints in the coevolution of murid RC subunits. We sequenced the cytochrome b gene from the species used to identify potential amino acid substitutions involved in such interactions. The greater sensitivity of complex III to xenocybrid dysfunction may result from the encoding of redox center apoproteins in both nuclear and mitochondrial genomes, a unique feature of this RC complex.