Nucleotide substitution rate of mammalian mitochondrial genomes

p53 in mitochondria enhances the accuracy of DNA synthesis

Mitochondrial localization of p53 was observed in stressed and unstressed cells. p53 is involved in DNA repair and apoptosis. It exerts physical and functional interactions with mitochondrial DNA and DNA polymerase gamma (pol gamma). The functional cooperation of p53 and pol gamma during DNA synthesis was examined in the mitochondrial fraction of p53-null H1299 cells, as the source of pol gamma. The results show that p53 may affect the accuracy of DNA synthesis in mitochondria: (1) the excision of a misincorporated nucleotide increases in the presence of (a) recombinant wild-type p53 (wtp53); (b) cytoplasmic fraction of LCC2 cells expressing endogenous wtp53 (but not specifically pre-depleted fraction); (c) cytoplasmic extract of H1299 cells overexpressing wtp53, but not exonuclease-deficient mutant p53-R175H. (2) Mitochondrial extracts of HCT116(p53+/+) cells display higher exonuclease activity compared with that of HCT116(p53-/-) cells. Addition of exogenous p53 complements the HCT116(p53-/-) mitochondrial extract mispair excision. Furthermore, the misincorporation was lower in the mitochondrial fraction of HCT116(p53+/+) cells as compared with that of HCT116(p53-/-) cells. (3) Irradiation-induced mitochondrial translocation of endogenous p53 in HCT116(p53+/+) cells correlates with the enhancement of error-correction activities. Taken together, the data suggest that p53 in mitochondria may be a component of an error-repair pathway and serve as guardian of the mitochondrial genome. The function of p53 in DNA repair and apoptosis is discussed.

Comparative analysis of mitochondrial genomes in Bombina (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.

Characterization of mitochondrial genome of Chinese wild mulberry silkworm, Bomyx mandarina (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Chinese Bombyx mandarina (ChBm) was determined. The circular genome is 15682 bp long, and contains a typical gene complement, order, and arrangement identical to that of Bombyx mori (B. mori) and Japanese Bombyx mandarina (JaBm) except for two additional tRNA-like structures: tRNA( Ser(TGA))-like and tRN( AIle(TAT))-like. All protein-coding sequences are initiated with a typical ATN codon except for the COI gene, which has a 4-bp TTAG putative initiator codon. Eleven of 13 protein-coding genes (PCGs) have a complete termination codon (all TAA), but the remaining two genes terminate with incomplete codons. All tRNAs have the typical clover-leaf structures of mitochondrial tRNAs, with the exception of tRNA( Ser(TGA))-like, with a four stem-and-loop structure. The length of the A+T-rich region of ChBm is 484 bp, shorter than those of JaBm (747 bp) and B. mori (494-499 bp). Phylogenetic analysis among B. mori, ChBm, JaBm, and Antheraea pernyi (Anpe) showed that B. mori is more closely related to ChBm than JaBm. The earliest divergence time estimate for B. mori-ChBm and B. mori-JaBm is about 1.08+/-0.18-1.41+/-0.24 and 1.53+/-0.20-2.01+/-0.26 Mya, respectively. ChBm and JaBm diverged around 1.11+/-0.16-1.45+/-0.21 Mya.

Adaptive evolution and functional redesign of core metabolic proteins in snakes

BACKGROUND

Adaptive evolutionary episodes in core metabolic proteins are uncommon, and are even more rarely linked to major macroevolutionary shifts.

METHODOLOGY/PRINCIPAL FINDINGS

We conducted extensive molecular evolutionary analyses on snake mitochondrial proteins and discovered multiple lines of evidence suggesting that the proteins at the core of aerobic metabolism in snakes have undergone remarkably large episodic bursts of adaptive change. We show that snake mitochondrial proteins experienced unprecedented levels of positive selection, coevolution, convergence, and reversion at functionally critical residues. We examined Cytochrome C oxidase subunit I (COI) in detail, and show that it experienced extensive modification of normally conserved residues involved in proton transport and delivery of electrons and oxygen. Thus, adaptive changes likely altered the flow of protons and other aspects of function in CO, thereby influencing fundamental characteristics of aerobic metabolism. We refer to these processes as "evolutionary redesign" because of the magnitude of the episodic bursts and the degree to which they affected core functional residues.

CONCLUSIONS/SIGNIFICANCE

The evolutionary redesign of snake COI coincided with adaptive bursts in other mitochondrial proteins and substantial changes in mitochondrial genome structure. It also generally coincided with or preceded major shifts in ecological niche and the evolution of extensive physiological adaptations related to lung reduction, large prey consumption, and venom evolution. The parallel timing of these major evolutionary events suggests that evolutionary redesign of metabolic and mitochondrial function may be related to, or underlie, the extreme changes in physiological and metabolic efficiency, flexibility, and innovation observed in snake evolution.

Human mitochondrial transfer RNAs: role of pathogenic mutation in disease

The human mitochondrial genome encodes 13 proteins. All are subunits of the respiratory chain complexes involved in energy metabolism. These proteins are translated by a set of 22 mitochondrial transfer RNAs (tRNAs) that are required for codon reading. Human mitochondrial tRNA genes are hotspots for pathogenic mutations and have attracted interest over the last two decades with the rapid discovery of point mutations associated with a vast array of neuromuscular disorders and diverse clinical phenotypes. In this review, we use a scoring system to determine the pathogenicity of the mutations and summarize the current knowledge of structure-function relationships of these mutant tRNAs. We also provide readers with an overview of a large variety of mechanisms by which mutations may affect the mitochondrial translation machinery and cause disease.

Mitochondrial genome evolution in the social amoebae

Most mitochondria contain a core set of genes required for mitochondrial function, but beyond this base there are variable genomic features. The mitochondrial genome of the model species Dictyostelium discoideum demonstrated that the social amoebae mitochondrial genomes have a size between those of metazoans and plants, but no comparative study of social amoebae mitochondria has been performed. Here, we present a comparative analysis of social amoebae mitochondrial genomes using D. discoideum, Dictyostelium citrinum, Dictyostelium fasciculatum, and Polysphondylium pallidum. The social amoebae mitochondria have similar sizes, AT content, gene content and have a high level of synteny except for one segmental rearrangement and extensive displacement of tRNAs. From the species that contain the rearrangement, it can be concluded that the event occurred late in the evolution of social amoebae. A phylogeny using 36 mitochondrial genes produced a well-supported tree suggesting that the pairs of D. discoideum/D. citrinum and D. fasciculatum/P. pallidum are sister species although the position of the root is not certain. Group I introns and endonucleases are variable in number and location in the social amoebae. Phylogenies of the introns and endonucleases suggest that there have been multiple recent duplications or extinctions and confirm that endonucleases have the ability to insert into new areas. An analysis of dN/dS ratios in mitochondrial genes revealed that among groups of genes, adenosine triphosphate synthase complex genes have the highest ratio, whereas cytochrome oxidase and nicotinamide adenine dinucleotide (NADH) dehydrogenase genes had the lowest ratio. The genetic codes of D. citrinum, P. pallidum, and D. fasciculatum are the universal code although D. fasciculatum does not use the TGA stop codon. In D. fasciculatum, we demonstrate for the first time that a mitochondrial genome without the TGA stop codon still uses the release factor RF2 that recognizes TGA. Theories of how the genetic code can change and why RF2 may be a constraint against switching codes are discussed.

Ancient hybridization and mitochondrial capture between two species of chipmunks

Models that posit speciation in the face of gene flow are replacing classical views that hybridization is rare between animal species. We use a multilocus approach to examine the history of hybridization and gene flow between two species of chipmunks (Tamias ruficaudus and T. amoenus). Previous studies have shown that these species occupy different ecological niches and have distinct genital bone morphologies, yet appear to be incompletely isolated reproductively in multiple areas of sympatry. We compared data from four sequenced nuclear loci and from seven microsatellite loci to published cytochrome b sequences. Interspecific gene flow was primarily restricted to introgression of the T. ruficaudus mitochondrial genome into a sympatric subspecies of T. amoenus, T. a. canicaudus, with the four sequenced nuclear loci showing little to no interspecific allele sharing. Microsatellite data were consistent with high levels of differentiation between the species and also showed no current gene flow between broadly sympatric populations of T. a. canicaudus and T. ruficaudus. Coalescent analyses date the mtDNA introgression event from the mid-Pleistocene to late Pliocene. Overall, these data indicate that introgression has had a minimal impact on the nuclear genomes of T. amoenus and T. ruficaudus despite multiple independent hybridization events. Our findings challenge long-standing assumptions on patterns of reproductive isolation in chipmunks and suggest that there may be other examples of hybridization among the 23 species of Tamias that occur in western North America.

Mitochondrial base excision repair of uracil and AP sites takes place by single-nucleotide insertion and long-patch DNA synthesis

Base excision repair (BER) corrects a variety of small base lesions in DNA. The UNG gene encodes both the nuclear (UNG2) and the mitochondrial (UNG1) forms of the human uracil-DNA glycosylase (UDG). We prepared mitochondrial extracts free of nuclear BER proteins from human cell lines. Using these extracts we show that UNG is the only detectable UDG in mitochondria, and mitochondrial BER (mtBER) of uracil and AP sites occur by both single-nucleotide insertion and long-patch repair DNA synthesis. Importantly, extracts of mitochondria carry out repair of modified AP sites which in nuclei occurs through long-patch BER. Such lesions may be rather prevalent in mitochondrial DNA because of its proximity to the electron transport chain, the primary site of production of reactive oxygen species. Furthermore, mitochondrial extracts remove 5' protruding flaps from DNA which can be formed during long-patch BER, by a "flap endonuclease like" activity, although flap endonuclease (FEN1) is not present in mitochondria. In conclusion, combined short- and long-patch BER activities enable mitochondria to repair a broader range of lesions in mtDNA than previously known.

Ancient DNA analysis of the Japanese sea lion (Zalophus californianus japonicus Peters, 1866): preliminary results using mitochondrial control-region sequences

In this study, we successfully extracted ancient DNA from skeletal remains of the Japanese sea lion-a species that is practically extinct-from archaeological sites and determined a partial sequence of its mitochondrial DNA control region. A molecular phylogenetic tree constructed by the neighbor-joining (NJ) method showed that the sequences from Japanese sea lions clustered together, with a high bootstrap value, and that this cluster was closest to the California sea lion cluster. The distinctly divergent cluster of Japanese sea lions reflected the morphological classification of these animals as a distinct species of the genus Zalophus; however, proximity to the California sea lion cluster simultaneously implied conformation with the traditional classification of these animals as a subspecies of Zalophus californianus. The average amount of nucleotide substitution between the Japanese and California sea lions was 7.02%. The Japanese and California sea lions were estimated to have diverged 2.2 million years ago, i.e., in the late Pliocene Epoch. This is the first report on a genetic analysis of the Japanese sea lion.

Recruitment of mitochondrial tRNA genes as auxiliary variability markers for both intra- and inter-species analysis: The paradigm of brown hare (Lepus europaeus)

We sequenced and analyzed the mitochondrial tRNA(Thr) and tRNA(Pro) genes from brown hare (Lepus europaeus) individuals of different geographic distribution and we investigated the role of various nucleotide substitutions that were detected. We compared these tRNAs with the respective available mitochondrial tRNA genes sequences within Lepus species and among mammals. The mutations that were detected represent specific and conserved polymorphisms that do not seem to affect the structural and functional features that are required for participation of tRNA molecules in mitochondrial protein synthesis. These changes however, possibly reflect on the evolutionary background of the species, which is based on the high intra-genomic variability and the evolutionary dynamic of the mitochondrial DNA. In an attempt to compare the phylogeny that is based on these specific tRNA genes with the phylogeny that is produced from sequencing data of the mitochondrial variable loop, we came up with results that indicate similar phylogeographic clusters. This observation implies that the tRNA mutations that were used for the present study have been well tolerated during evolution and they define an additional genetic and biochemical tag that can be used for such studies. Based on this notion and according to our results, we propose that mitochondrial tRNA genes can be used as valuable auxiliary molecular markers for contemporaneous and linked biochemical and genetic analyses.

mtDNA diversity and genetic lineages of eighteen cattle breeds from Bos taurus and Bos indicus in China

In order to clarify the origin and genetic diversity of indigenous cattle breeds in China, we carried out phylogenetic analysis of representatives of those breeds by employing mitochondrial gene polymorphism. Complete cyt b gene sequences, 1140 bp in length, were determined for a total of 136 individuals from 18 different breeds and these sequences were clustered into two distinct genetic lineages: taurine (Bos taurus) and zebu (Bos indicus). In analysis of the cyt b gene diversity, Chinese cattle showed higher nucleotide (0.00923) and haplotype diversity (0.848) than the reports from other studies, and the animals from the taurine lineage indicated higher nucleotide diversity (0.00330) and haplotype diversity (0.746) than the ones from the zebu lineage (0.00136; 0.661). The zebu mtDNA dominated in the southern breeds (63.3-100%), while the taurine dominated in the northern breeds (81.8-100%). Six cattle breeds from the central area of China exhibited intermediate frequencies of zebu mtDNA (25-71.4%). This polymorphism revealed a declining south-to-north gradient of female zebu introgression and a geographical hybrid zone of Bos taurus and Bos indicus in China.

Comparative mitochondrial genomics of snakes: extraordinary substitution rate dynamics and functionality of the duplicate control region

BACKGROUND

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

RESULTS

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

CONCLUSION

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.

THE AUTOMATION AND EVALUATION OF NESTED CLADE PHYLOGEOGRAPHIC ANALYSIS

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.

Positive and negative selection on the mitochondrial genome

Several recent studies have confirmed that mitochondrial DNA variation and evolution are not consistent with the neutral theory of molecular evolution and might be inappropriate for estimating effective population sizes. Evidence for the action of both positive and negative selection on mitochondrial genes has been put forward, and the complex genetics of mitochondrial DNA adds to the challenge of resolving this debate. The solution could lie in distinguishing genetic drift from 'genetic draft' and in dissecting the physiology of mitochondrial fitness.

Accelerated molecular evolution in Microtus (Rodentia) as assessed via complete mitochondrial genome sequences

Microtus is one of the most taxonomically diverse mammalian genera, including over 60 extant species. These rodents have evolved rapidly, as the genus originated less than 2 million years ago. If these numbers are taken at face value, then an average of 30 microtine speciation events have occurred every million years. One explanation for the rapid rate of cladogenesis in Microtus could be the karyotypic differentiation exhibited across the genus: diploid numbers range from 17 to 64. Despite the striking chromosomal variability within Microtus, phenotypic variation is unremarkable. To determine whether nucleotide substitution rates are also elevated in voles, we sequenced the entire mitochondrial DNA (mtDNA) genome of the Eurasian sibling vole (Microtus rossiaemeridionalis). We compared this genome to another previously sequenced vole mtDNA genome (Microtus kikuchii) and performed pairwise sequence comparisons with the mtDNA genomes of ten additional mammalian genera. We found that microtine mtDNA genomes are evolving more rapidly than any other mammalian lineage we sampled, as gauged by the rate of nucleotide substitution across the entire mtDNA genome as well as at each individual protein-coding gene. Additionally, we compared substitution rates within the cytochrome b gene to seven other rodent genera and found that Microtus mtDNA is evolving fastest. The root cause of accelerated evolution in Microtus remains uncertain, but merits further investigation.

Repetitive DNA sequences include retrotransposons in genomes of the Glomeromycota

Twenty-five repetitive elements are first described in the genomes of the arbuscular mycorrhizal (AM) fungi Gigaspora margarita, Gig. rosea and Glomus mosseae. Nineteen repetitive DNA sequences isolated by genomic library screening and four by self-priming PCR had no homology to known DNA sequences, except for two Gig. margarita sequences and one Gig. rosea sequence which showed amino acid similarity to retrotransposons. Part of the Gig. rosea sequence was also similar to a DNA transposon. Two other retrotransposon sequences were isolated using PCR targeting of reverse transcriptase and ribonuclease H domains. Evidence is provided for three gypsy-like LTR retrotransposon and two non-LTR retrotransposon sequences in the AM fungal genomes. Four contain stop codons indicating that they cannot be active. Expression of three retrotransposons was not detected in germinating spores or intraradical hyphae of Gig. margarita. Southern blot analyses indicated that these three sequences are dispersed in the genome and that two are methylated. Sequence analysis of different GmarLTR1 copies showed they have undergone mutations by transitions, which may have been induced by cytosine methylation. Transposable elements may have played a major role in shaping genome structure and size during evolution of the Glomeromycota.

Cloning cattle: the methods in the madness

Somatic cell nuclear transfer (SCNT) is much more widely and efficiently practiced in cattle than in any other species, making this arguably the most important mammal cloned to date. While the initial objective behind cattle cloning was commercially driven--in particular to multiply genetically superior animals with desired phenotypic traits and to produce genetically modified animals-researchers have now started to use bovine SCNT as a tool to address diverse questions in developmental and cell biology. In this paper, we review current cattle cloning methodologies and their potential technical or biological pitfalls at any step of the procedure. In doing so, we focus on one methodological parameter, namely donor cell selection. We emphasize the impact of epigenetic and genetic differences between embryonic, germ, and somatic donor cell types on cloning efficiency. Lastly, we discuss adult phenotypes and fitness of cloned cattle and their offspring and illustrate some of the more imminent commercial cattle cloning applications.

Variability of the mitochondrial genome in mammals at the inter-species/intra-species boundary

Genomic variations represent the molecular basis of the biodiversity of living organisms on which selection operates to generate evolution. In eukaryotes, genomic variability can be experienced in both nuclear and organellar, i.e. mitochondrial and plastid (where present), genomes, which can follow completely different evolution pathways, as revealed by comparative genomics analyses. In Metazoa, for which a substantial number of complete genome sequences are available (nuclear, but mainly mitochondrial), we are just starting to grasp the selective pressures operating on some basic features of the genome as a whole. In this brief review, we discuss the variability of the mitochondrial metazoan genome, with particular reference to mitochondrial DNA in mammals. In light of the recent assumption that a small segment of mitochondrial DNA may be used, particularly in Metazoa, as a species marker, some data on mitochondrial gene variability at the inter-species/intra-species boundary are reported. Intra-species variability has been evaluated in four mammalian species, Homo sapiens, Bos taurus, Sus scrofa and Canis familiaris, whereas the relationship between intra- and inter-species variability has been investigated in Bos taurus and Bos indicus.

Comparative phylogeography of eastern chipmunks and white-footed mice in relation to the individualistic nature of species

Palaeoecological studies have demonstrated that ecological communities as a whole did not remain stable throughout the climatic fluctuations of the Quaternary. The result is that long-term associations of species cannot be inferred by contemporary associations in ecological communities. Therefore, the evolutionary significance of any contemporary ecological interactions among species and of the biotic community within which species have evolved also cannot be assumed from contemporary conditions. Comparative phylogeographic data provide a method to identify species within ecological communities that have shared biogeographic histories. We present an example of a long-term association between populations of two mammalian species, eastern chipmunks (Tamias striatus) and white-footed mice (Peromyscus leucopus), which are commonly associated with deciduous forest habitats. The distribution of mitochondrial DNA variation in T. striatus and P. leucopus from previously glaciated regions of the eastern United States support the hypothesis that, in at least part of their range, genetic lineages of the two species have expanded from similar population sources since the Last Glacial Maximum. In addition, the spatial concordance of genetic lineages of T. striatus and P. leucopus with the oak-savannah forest formations of Wisconsin and Illinois, suggest that populations associated with this community colonized the area in association with a set of arboreal species that comprise their deciduous forest habitat.

Five ovine mitochondrial lineages identified from sheep breeds of the near East

Archaeozoological evidence indicates that sheep were first domesticated in the Fertile Crescent. To search for DNA sequence diversity arising from previously undetected domestication events, this survey examined nine breeds of sheep from modern-day Turkey and Israel. A total of 2027 bp of mitochondrial DNA (mtDNA) sequence from 197 sheep revealed a total of 85 haplotypes and a high level of genetic diversity. Six individuals carried three haplotypes, which clustered separately from the known ovine mtDNA lineages A, B, and C. Analysis of genetic distance, mismatch distribution, and comparisons with wild sheep confirmed that these represent two additional mtDNA lineages denoted D and E. The two haplogroup E sequences were found to link the previously identified groups A and C. The single haplogroup D sequence branched with the eastern mouflon (Ovis orientalis), urial (O. vignei), and argali (O. ammon) sheep. High sequence diversity (K = 1.86%, haplogroup D and O. orientalis) indicates that the wild progenitor of this domestic lineage remains unresolved. The identification in this study of evidence for additional domestication events adds to the emerging view that sheep were recruited from wild populations multiple times in the same way as for other livestock species such as goat, cattle, and pig.

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.

Similar compositional biases are caused by very different mutational effects

Compositional replication strand bias, commonly referred to as GC skew, is present in many genomes of prokaryotes, eukaryotes, and viruses. Although cytosine deamination in ssDNA (resulting in C-->T changes on the leading strand) is often invoked as its major cause, the precise contributions of this and other substitution types are currently unknown. It is also unclear if the underlying mutational asymmetries are the same among taxa, are stable over time, or how closely the observed biases are to mutational equilibrium. We analyzed nearly neutral sites of seven taxa each with between three and six complete bacterial genomes, and inferred the substitution spectra of fourfold degenerate positions in nonhighly expressed genes. Using a bootstrap procedure, we extracted compositional biases associated with replication and identified the significant asymmetries. Although all taxa showed an overrepresentation of G relative to C on the leading strand (and imbalances between A and T), widely variable substitution asymmetries are noted. Surprisingly, all substitution types show significant asymmetry in at least one taxon, but none were universally biased in all taxa. Notably, in the two most biased genomes, A-->G, rather than C-->T, shapes the compositional bias. Given the variability in these biases, we propose that the process is multifactorial. Finally, we also find that most genomes are not at compositional equilibrium, and suggest that mutational-based heterotachy is deeply imprinted in the history of biological macromolecules. This shows that similar compositional biases associated with the same essential well-conserved process, replication, do not reflect similar mutational processes in different genomes, and that caution is required in inferring the roles of specific mutational biases on the basis of contemporary patterns of sequence composition.

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.

Balaenoptera omurai is a newly discovered baleen whale that represents an ancient evolutionary lineage

Balaenoptera omurai, formerly classified as a small form of Bryde's whale, was recently reclassified as a new baleen whale species of the family Balaenopteridae. Although researchers have investigated the evolutionary history of Balaenopteridae and their relatives using molecular phylogenetic methods, the taxonomy of the ordinary Bryde's whale (Balaenoptera brydei) and small-form Bryde's whales (Balaenoptera edeni and B. omurai) remains unclear. We have used complete mtDNA sequences and short interspersed repetitive element (SINE) insertion patterns to construct the evolutionary history of both B. omurai and the taxonomically redefined species, B. edeni. The combined results demonstrate that B. omurai forms a monophyletic lineage with B. musculus, B. brydei, B. edeni and B. borealis and that B. omurai and B. musculus successively diverged from their common ancestor. In addition, we also showed that B. edeni constitutes a sister taxon to B. brydei. Our data suggest that B. omurai evolved as an ancient independent lineage that diverged much earlier than B. borealis, B. brydei and B. edeni, which were previously believed to be closely related to B. omurai.

Fidelity of the human mitochondrial DNA polymerase

We have quantified the fidelity of polymerization of DNA by human mitochondrial DNA polymerase using synthetic DNA oligonucleotides and recombinant holoenzyme and examining each of the possible 16-base pair combinations. Although the kinetics of incorporation for all correct nucleotides are similar, with an average Kd of 0.8 microM and an average k(pol) of 37 s(-1), the kinetics of misincorporation vary widely. The ground state binding Kd of incorrect bases ranges from a low of 25 microM for a dATP:A mispair to a high of 360 microM for a dCTP:T mispair. Similarly, the rates of incorporation of incorrect bases vary from 0.0031 s(-1) for a dCTP:C mispair to 1.16 s(-1) for a dGTP:T mispair. Due to the variability in the kinetic parameters for misincorporation, the estimates of fidelity range from 1 error in 3563 nucleotides for dGTP:T to 1 error in 2.3 x 10(6) nucleotides for dCTP:C. Interestingly, the discrimination against a dGTP:T mismatch is 16.5 times lower than that of a dTTP:G mismatch due to a tighter Kd for ground state binding and a faster rate of incorporation of the dGTP:T mismatch relative to the dTTP:G mismatch. We calculate an average fidelity of 1 error in 440,000 nucleotides.

Mitochondrial diversity in European and Chinese pigs is consistent with population expansions that occurred prior to domestication

Mitochondrial DNA (mtDNA) diversity in European and Asian pigs was assessed using 1536 samples representing 45 European and 21 Chinese breeds. Diagnostic nucleotide differences in the cytochrome b (Cytb) gene between the European and Asian mtDNA variants were determined by pyrosequencing as a rapid screening method. Subsequently, 637bp of the hypervariable control region was sequenced to further characterize mtDNA diversity. All sequences belonged to the D1 and D2 clusters of pig mtDNA originating from ancestral wild boar populations in Europe and Asia, respectively. The average frequency of Asian mtDNA haplotypes was 29% across European breeds, but varied from 0 to 100% within individual breeds. A neighbour-joining (NJ) tree of control region sequences showed that European and Asian haplotypes form distinct clusters consistent with the independent domestication of pigs in Asia and Europe. The Asian haplotypes found in the European pigs were identical or closely related to those found in domestic pigs from Southeast China. The star-like pattern detected by network analysis for both the European and Asian haplotypes was consistent with a previous demographic expansion. Mismatch analysis supported this notion and suggested that the expansion was initiated before domestication.

The complete mitochondrial genome of the whiting, Merlangius merlangus and the haddock, Melanogrammus aeglefinus: a detailed genomic comparison among closely related species of the Gadidae family

We determined the first complete mitochondrial DNA (mtDNA) sequences for the whiting (Merlangius merlangus, family Gadidae, order Gadiformes) and the haddock (Melanogrammus aeglefinus, family Gadidae, order Gadiformes). The entire mitogenomes were amplified and sequenced by primer walking using newly designed specific internal primers. Lengths were 16,569 and 16,585 bases for whiting and haddock respectively, lengths which lie within the range of previously reported gadiform sequences from Atlantic cod (Gadus morhua, 16,696 bases) and walleye pollock (Theragra chalcogramma, 16,570 bases). Gene arrangement in both species conformed to the order seen in most vertebrate mitochondrial genomes. We identified a long intergenic spacer located between the tRNA(Thr) and tRNA(Pro) genes (of 100 and 70 bp long for whiting and haddock, respectively), as previously described for other species of the order Gadiformes. Using nucleotide and amino acid divergence data of four complete gadoid mitogenomes (M. merlangius, M. aeglefinus, G. morhua and T. chalcogramma), we examined in detail the relative mtDNA mutation patterns across genes and among Gadidae species and tested for the performance of each protein-coding, transfer RNA and ribosomal RNA gene in depicting the expected phylogeny among the four species, as compared with the whole genome dataset. This comparison may be particularly useful in phylogenetic analyses of such a diverse fish family, as well as for the understanding of the patterns of nucleotide substitution of the mtDNA at low levels of divergence.

History, ocean channels, and distance determine phylogeographic patterns in three widespread Philippine fruit bats (Pteropodidae)

The comparative phylogeography of widespread, codistributed species provides unique insights into regional biodiversity and diversification patterns. I used partial DNA sequences of the mitochondrial genes ND2 and cyt b to investigate phylogeographic structure in three widespread Philippine fruit bats. Ptenochirus jagori is endemic to the oceanic region of the Philippines and is most abundant in lowland primary forest. Macroglossus minimus and Cynopterus brachyotis are most common in disturbed and open habitats and are not endemic. In all three, genetic differentiation is present at multiple spatial scales and is associated to some degree with Pleistocene landbridge island groups. In P. jagori and C. brachyotis, genetic distance is correlated with geographic distance; in C. brachyotis and M. minimus, it is correlated with the sea-crossing distance between islands. P. jagori has the least overall genetic structure of these three species, whereas C. brachyotis and M. minimus have more geographic association among haplotypes, suggesting that phylogeographic patterns are linked to ecology and habitat preference. However, contrary to expectation, the two widespread, disturbed habitat species have more structure than the endemic species. Mismatch distributions suggest rapid changes in effective population size in C. brachyotis and P. jagori, whereas M. minimus appears to be demographically more stable. Geologic and geographic history are important in structuring variation, and phylogeographic patterns are the result of dynamic long-term processes rather than simply reflecting current conditions.

Choosing the right molecular genetic markers for studying biodiversity: from molecular evolution to practical aspects

The use of molecular genetic markers (MGMs) has become widespread among evolutionary biologists, and the methods of analysis of genetic data improve rapidly, yet an organized framework in which scientists can work is lacking. Elements of molecular evolution are summarized to explain the origin of variation at the DNA level, its measures, and the relationships linking genetic variability to the biological parameters of the studied organisms. MGM are defined by two components: the DNA region(s) screened, and the technique used to reveal its variation. Criteria of choice belong to three categories: (1) the level of variability, (2) the nature of the information (e.g. dominance vs. codominance, ploidy, ... ) which must be determined according to the biological question and (3) some practical criteria which mainly depend on the equipment of the laboratory and experience of the scientist. A three-step procedure is proposed for drawing up MGMs suitable to answer given biological questions, and compiled data are organized to guide the choice at each step: (1) choice, determined by the biological question, of the level of variability and of the criteria of the nature of information, (2) choice of the DNA region and (3) choice of the technique.

Mitochondria and the success of somatic cell nuclear transfer cloning: from nuclear-mitochondrial interactions to mitochondrial complementation and mitochondrial DNA recombination

The overall success of somatic cell nuclear transfer (SCNT) cloning is rather unsatisfactory, both in terms of efficacy and from an animal health and welfare point of view. Most research activities have concentrated on epigenetic reprogramming problems as one major cause of SCNT failure. The present review addresses the limited success of mammalian SCNT from yet another viewpoint, the mitochondrial perspective. Mitochondria have a broad range of critical functions in cellular energy supply, cell signalling and programmed cell death and, thus, affect embryonic and fetal development, suggesting that inadequate or perturbed mitochondrial functions may adversely affect SCNT success. A survey of perinatal clinical data from human subjects with deficient mitochondrial respiratory chain activity has revealed a plethora of phenotypes that have striking similarities with abnormalities commonly encountered in SCNT fetuses and offspring. We discuss the limited experimental data on nuclear-mitochondrial interaction effects in SCNT and explore the potential effects in the context of new findings about the biology of mitochondria. These include mitochondrial fusion/fission, mitochondrial complementation and mitochondrial DNA recombination, processes that are likely to be affected by and impact on SCNT cloning. Furthermore, we indicate pathways that could link epigenetic reprogramming and mitochondria effects in SCNT and address questions and perspectives for future research.

Disentangling reasons for low Y chromosome variation in the greater white-toothed shrew (Crocidura russula)

Y chromosome variation is determined by several confounding factors including mutation rate, effective population size, demography, and selection. Disentangling these factors is essential to better understand the evolutionary properties of the Y chromosome. We analyzed genetic variation on the Y chromosome, X chromosome, and mtDNA of the greater white-toothed shrew, a species with low variance in male reproductive success and limited sex-biased dispersal, which enables us to control to some extent for life-history effects. We also compared ancestral (Moroccan) to derived (European) populations to investigate the role of demographic history in determining Y variation. Recent colonization of Europe by a small number of founders (combined with low mutation rates) is largely responsible for low diversity observed on the European Y and X chromosomes compared to mtDNA. After accounting for mutation rate, copy number, and demography, the Y chromosome still displays a deficit in variation relative to the X in both populations. This is possibly influenced by directional selection, but the slightly higher variance in male reproductive success is also likely to play a role, even though the difference is small compared to that in highly polygynous species. This study illustrates that demography and life-history effects should be scrutinized before inferring strong selective pressure as a reason for low diversity on the Y chromosome.

Evidence of three maternal lineages in Near Eastern sheep supporting multiple domestication events

The variability of mtDNA was analysed in local sheep breeds reared throughout Turkey, for which a fragment of the D-loop region and the complete cytochrome b were sequenced. Phylogenetic analyses performed independently for the D-loop and the Cyt b gene revealed three clearly separated clusters indicating three major maternal lineages, two of which had been previously described as types B and A. The new type, C, was present in all the breeds analysed and showed considerable mtDNA variability. Divergence time was obtained on the basis of Cyt b gene and was estimated to be around 160,000-170,000 years ago for lineages B and A, whereas the divergence of lineage C proved to have occurred earlier (between 450,000 and 750,000 years ago). These times greatly predate domestication and suggest that the origin of modern sheep breeds was more complex than previously thought and that at least three independent sheep domestication events occurred. Our results, together with archaeological information and the current wild sheep populations in the Near East region support the high importance of this area in the sheep domestication process. Finally, the evidence of a third maternal lineage has important implications regarding the history of modern sheep.

Population histories of right whales (Cetacea: Eubalaena) inferred from mitochondrial sequence diversities and divergences of their whale lice (Amphipoda: Cyamus)

Right whales carry large populations of three 'whale lice' (Cyamus ovalis, Cyamus gracilis, Cyamus erraticus) that have no other hosts. We used sequence variation in the mitochondrial COI gene to ask (i) whether cyamid population structures might reveal associations among right whale individuals and subpopulations, (ii) whether the divergences of the three nominally conspecific cyamid species on North Atlantic, North Pacific, and southern right whales (Eubalaena glacialis, Eubalaena japonica, Eubalaena australis) might indicate their times of separation, and (iii) whether the shapes of cyamid gene trees might contain information about changes in the population sizes of right whales. We found high levels of nucleotide diversity but almost no population structure within oceans, indicating large effective population sizes and high rates of transfer between whales and subpopulations. North Atlantic and Southern Ocean populations of all three species are reciprocally monophyletic, and North Pacific C. erraticus is well separated from North Atlantic and southern C. erraticus. Mitochondrial clock calibrations suggest that these divergences occurred around 6 million years ago (Ma), and that the Eubalaena mitochondrial clock is very slow. North Pacific C. ovalis forms a clade inside the southern C. ovalis gene tree, implying that at least one right whale has crossed the equator in the Pacific Ocean within the last 1-2 million years (Myr). Low-frequency polymorphisms are more common than expected under neutrality for populations of constant size, but there is no obvious signal of rapid, interspecifically congruent expansion of the kind that would be expected if North Atlantic or southern right whales had experienced a prolonged population bottleneck within the last 0.5 Myr.

Worldwide structure of mtDNA diversity among Cuvier's beaked whales (Ziphius cavirostris): implications for threatened populations

We present the first description of phylogeographic structure among Cuvier's beaked whales (Ziphius cavirostris) worldwide using mitochondrial DNA (mtDNA) control region sequences obtained from strandings (n = 70), incidental fisheries takes (n = 11), biopsy (n = 1), and whale-meat markets (n = 5). Over a 290-base pair fragment, 23 variable sites defined 33 unique haplotypes among the total of 87 samples. Nucleotide diversity at the control region was relatively low (pi = 1.27%+/- 0.723%) compared to wide-ranging baleen whales, but higher than strongly matrifocal sperm, pilot and killer whales. Phylogenetic reconstruction using maximum likelihood revealed four distinct haplotype groups, each of which displayed strong frequency differences among ocean basins, but no reciprocal monophyly or fixed character differences. Consistent with this phylogeographic pattern, an analysis of molecular variance showed high levels of differentiation among ocean basins (F(ST) = 0.14, Phi ST = 0.42; P < 0.001). Estimated rates of female migration among ocean basins were low (generally < or = 2 individuals per generation). Regional sample sizes were too small to detect subdivisions within oceans except in the North Atlantic, where the Mediterranean Sea (n = 12) was highly differentiated due to the presence of two private haplotypes. One market product purchased in South Korea grouped with other haplotypes found only in the North Atlantic, suggesting a violation of current agreements banning international trade in cetacean species. Together, these results demonstrate a high degree of isolation and low maternal gene flow among oceanic, and in some cases, regional populations of Cuvier's beaked whales. This has important implications for understanding the threats of human impact, including fisheries by-catch, direct hunting, and disturbance or mortality from anthropogenic sound.

The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms

Complete mitochondrial genome sequences for the schistosomes Schistosoma haematobium and Schistosoma. spindale have been characterized. S. haematobium is the causative agent of urinary schistosomiasis in humans and S. spindale uses ruminants as its definitive host; both are transmitted by freshwater snail intermediate hosts. Results confirm a major gene order rearrangement among schistosomes in all traditional Schistosoma species groups other than Schistosoma japonicum; i.e., species groups S. mansoni, S. haematobium, and S. indicum. These data lend support to the 'out of Asia' (East and Southeast Asia) hypothesis for Schistosoma. The gene order change involves translocation of atp6-nad2-trnA and a rearrangement of nad3-nad1 relative to other parasitic flatworm mt genomes so far sequenced. Gene order and tRNA secondary structure changes (loss and acquisition of the DHU and/or TPsiC arms of trnC, trnF, and trnR) between mitochondrial genomes of these and other (digenean and cestode) flatworms were inferred by character mapping onto a phylogeny estimated from nuclear small subunit rRNA gene sequences of these same species, in order to find additional rare genomic changes suitable as synapomorphies. Denser and wider taxon sampling of mt genomes across the Platyhelminthes will validate these putative characters.

Comparative analysis of gender-associated complete mitochondrial genomes in marine mussels (Mytilus spp.)

Marine mussels of the genus Mytilus have an unusual mode of mitochondrial DNA (mtDNA) transmission termed doubly uniparental inheritance (DUI). Female mussels are homoplasmic for the F mitotype, which is inherited maternally, while males are usually heteroplasmic, carrying a mixture of the maternal F mitotype and the paternally inherited M genome. Two classes of M genomes have been observed: "standard" M genomes and "recently masculinized" M genomes. The latter are more similar to F genomes at the sequence level but are transmitted paternally like standard M genomes. In this study we report the complete sequences of two standard male M. edulis and one recently masculinized male M. trossulus mitochondrial genome. A comparative analysis, including the previously sequenced M. edulis F and M. galloprovincialis F and M mtDNAs, reveals that these genomes are identical in gene order, but highly divergent in nucleotide and amino acid sequence. The large amount (>20%) of nucleotide substitutions that fall in coding regions implies that there are several amino acid replacements between the F and M genomes, which likely have an impact on the structural and functional properties of the mitochondrial proteome. Correlation of the divergence rate of different protein-coding genes indicates that mtDNA-encoded proteins of the M genome are still under selective constraints, although less highly than genes of the F genome. The mosaic F/M control region of the masculinized F genome provides evidence for lineage-specific sequences that may be responsible for the different mode of transmission genetics. This analysis shows the value of comparative genomics to better understand the mechanisms of maintenance and segregation of mtDNA sequence variants in mytilid mussels.

Detecting past population bottlenecks using temporal genetic data

Population bottlenecks wield a powerful influence on the evolution of species and populations by reducing the repertoire of responses available for stochastic environmental events. Although modern contractions of wild populations due to human-related impacts have been documented globally, discerning historic bottlenecks for all but the most recent and severe events remains a serious challenge. Genetic samples dating to different points in time may provide a solution in some cases. We conducted serial coalescent simulations to assess the extent to which temporal genetic data are informative regarding population bottlenecks. These simulations demonstrated that the power to reject a constant population size hypothesis using both ancient and modern genetic data is almost always higher than that based solely on modern data. The difference in power between the modern and temporal DNA approaches depends significantly on effective population size and bottleneck intensity and less significantly on sample size. The temporal approach provides more power in cases of genetic recovery (via migration) from a bottleneck than in cases of demographic recovery (via population growth). Choice of genetic region is critical, as mutation rate heavily influences the extent to which temporal sampling yields novel information regarding the demographic history of populations.

New DNA data from a transthyretin nuclear intron suggest an Oligocene to Miocene diversification of living South America opossums (Marsupialia: Didelphidae)

Phylogenetic relationships of 19 species of didelphid marsupials were studied using two nuclear markers, the non-coding transthyretin intron 1 (TTR) and the coding interphotoreceptor retinoid binding protein exon 1 (IRBP), and two mitochondrial genes, the protein-coding cytochrome b (cyt-b) and the structural 12S ribosomal DNA (12S rDNA). Evolutionary dynamics of these four markers were compared to each other, revealing the appropriate properties presented by TTR intron 1 together with its well supported and resolved phylogenetic signal. Nuclear markers supported the monophyly of medium and large-sized opossums Metachirus+(Chironectes, Lutreolina, Didelphis, Philander), and the paraphyly of mouse-sized opossums, with the genera Gracilinanus, Thylamys, and Marmosops as a sister group to medium and large-sized didelphids. Conflicting branching patterns between mitochondrial and nuclear data involved the phylogenetic position of Marmosa-Micoureus-Monodelphis relative to other mouse-sized opossums. Nuclear phylogenetic inferences among genera were confirmed by the presence of synapomorphic indels observed in TTR intron 1. A Bayesian relaxed molecular clock dating of didelphid evolution using nuclear markers estimated their origin in the Middle Eocene (39.8 million years ago), with subsequent diversification during the Oligocene (Deseadan) and Miocene.

Mitochondrial phylogeography of the Eurasian beaver Castor fiber L

Nucleotide variation in an approximately 490 bp fragment of the mitochondrial DNA control region (mtDNA CR) was used to describe the genetic variation and phylogeographical pattern in the Eurasian beaver (Castor fiber) over its entire range. The sampling effort was focused on the relict populations that survived a drastic population bottleneck, caused by overhunting, at the end of the 19th century. A total of 152 individuals grouped into eight populations representing all currently recognized subspecies were studied. Sixteen haplotypes were detected, none of them shared among populations. Intrapopulation sequence variation was very low, most likely a result of the severe bottleneck. Extreme genetic structure could result from human-mediated extinction of intermediate populations, but it could also be an effect of prior substantial structuring of the beaver populations with watersheds of major Eurasian rivers acting as barriers to gene flow. Phylogenetic analysis revealed the presence of two mtDNA lineages: eastern (Poland, Lithuania, Russia and Mongolia) and western (Germany, Norway and France), the former comprising more divergent haplotypes. The low level of sequence divergence of the entire cytochrome b gene among six individuals representing six subspecies suggests differentiation during the last glacial period and existence of multiple glacial refugia. At least two evolutionary significant units (ESU) can be identified, the western and the eastern haplogroup. The individual relict populations should be regarded as management units, the eastern subspecies possibly also as ESUs. Guidelines for future translocations and reintroductions are proposed.

Mitochondrial Sequence Reveals High Levels of Gene Flow Between Breeds of Domestic Sheep from Asia and Europe

Sequence variation present within the mitochondrial genome was used to investigate genetic diversity within sheep breeds from Asia and Europe. Comparison of 2027 bp of sequence from 121 animals revealed 44 phylogenetically informative nucleotide positions and a single insertion/deletion. A total of 57 haplotypes were observed which formed two distinct clades. Type A haplotypes were found in breeds from Asia (India, Indonesia, Mongolia, and Tibet), while type B haplotypes were observed at the highest frequency in breeds sourced from Europe (nine breeds from Austria, Aland, Finland, Spain, and northwestern Russia). The distribution of haplotypes indicates sheep appear to have the weakest population structure and the highest rate of intercontinental dispersal of any domestic animal reported to date. Only 2.7% of the sequence variation observed was partitioned between continents, which is lower than both goat (approximately 10%) and cattle (approximately 50%). Diagnostic restriction fragment length polymorphism polymerase chain reaction (RFLP-PCR) tests which distinguish type A and B haplotypes were used to test an additional 223 animals from 17 breeds of European and Asian origin. A mixture of the two lineages was found in every breed except Suffolk and the Indian Garole, indicating introgression has played a major part during breed development and subsequent selection.

Sequence variation in the tRNA genes of human mitochondrial DNA

Recent analyses have shown that nonsynonymous variation in human mitochondrial DNA (mtDNA) contains nonneutral variants, suggesting the presence of mildly deleterious mutations. Many of the disease-causing mutations in mtDNA occur in the genes encoding the tRNAs. Nucleotide sequence variation in these genes has not been studied in human populations, nor have the structural consequences of nucleotide substitutions in tRNA molecules been examined. We therefore determined the nucleotide sequences of the 22 tRNA genes in the mtDNA of 477 Finns and, also, obtained 435 European sequences from the MitoKor database. No differences in population polymorphism indices were found between the two data sets. We assessed selective constraints against various tRNA domains by comparing allele frequencies between these domains and the synonymous and nonsynonymous sites, respectively. All tRNA domains except the variable loop were more conserved than synonymous sites, and T stem and D stem were more conserved than the respective loops. We also analyzed the energetic consequences of the 96 polymorphisms recovered in the two data sets or in the Mitomap database. The minimum free energy (DeltaG) was calculated using the free energy rules as implemented in mfold version 3.1. The DeltaG's were normally distributed among the 22 wild-type tRNA genes, whereas the 96 polymorphic tRNAs departed significantly from a normal distribution. The largest differences in DeltaG between the wild-type and the polymorphic tRNAs in the Finnish population tended to be in the polymorphisms that were present at low frequencies. Allele frequency distributions and minimum free energy calculations both suggested that some polymorphisms in tRNA genes are nonneutral.