Divergent evolution of life span associated with mitochondrial DNA evolution

Mitochondria play a key role in ageing. The pursuit of genes that regulate variation in life span and ageing have shown that several nuclear‐encoded mitochondrial genes are important. However, the role of mitochondrial encoded genes (mtDNA) is more controversial and our appreciation of the role of mtDNA for the evolution of life span is limited. We use replicated lines of seed beetles that have been artificially selected for long or short life for >190 generations, now showing dramatic phenotypic differences, to test for a possible role of mtDNA in the divergent evolution of ageing and life span. We show that these divergent selection regimes led to the evolution of significantly different mtDNA haplotype frequencies. Selection for a long life and late reproduction generated positive selection for one specific haplotype, which was fixed in most such lines. In contrast, selection for reproduction early in life led to both positive selection as well as negative frequency‐dependent selection on two different haplotypes, which were both present in all such lines. Our findings suggest that the evolution of life span was in part mediated by mtDNA, providing support for the emerging general tenet that adaptive evolution of life‐history syndromes may involve mtDNA.     

An ancient mitochondrial polymorphism in Adalis bipunctata linked to a sex-ratio-distorting bacterium

Sex-ratio-distorting microbes are common parasites of arthropods. Although the reasons they have invaded and spread though populations are well understood, their subsequent dynamics within those populations are virtually unknown. We have found that different strains of a male-killing Rickettsia bacterium infecting the beetle Adalia bipunctata are associated with distinct mitochondrial haplotypes, which is expected as both the mitochondria and the bacteria are maternally transmitted. These mitochondrial haplotypes shared a common ancestor >2 million years ago, and their overall diversity is significantly greater than expected under neutrality from comparisons with a nuclear gene. Furthermore, a variety of statistical tests show strong deviations from neutrality in mitochondrial but not in nuclear genes. We therefore conclude that natural selection is probably maintaining a polymorphism of different Rickettsia strains in this species. Despite the age of the different mitochondrial haplotypes, there is very little genetic diversity within them. Furthermore, there is considerable variation in mitochondrial haplotype and bacterial strain frequency between populations, despite it being thought that this species has fairly low levels of population structure. We conclude that the fitness of these male killers may be negatively frequency dependent or different strains may be favored in different populations. These hypotheses await experimental confirmation.     

Within‐population genetic effects of mtDNA on metabolic rate in Drosophila subobscura

A growing body of research supports the view that within‐species sequence variation in the mitochondrial genome (mtDNA) is functional, in the sense that it has important phenotypic effects. However, most of this empirical foundation is based on comparisons across populations, and few studies have addressed the functional significance of mtDNA polymorphism within populations. Here, using mitonuclear introgression lines, we assess differences in whole‐organism metabolic rate of adult Drosophila subobscura fruit flies carrying either of three different sympatric mtDNA haplotypes. We document sizeable, up to 20%, differences in metabolic rate across these mtDNA haplotypes. Further, these mtDNA effects are to some extent sex specific. We found no significant nuclear or mitonuclear genetic effects on metabolic rate, consistent with a low degree of linkage disequilibrium between mitochondrial and nuclear genes within populations. The fact that mtDNA haplotype variation within a natural population affects metabolic rate, which is a key physiological trait with important effects on life‐history traits, adds weight to the emergent view that mtDNA haplotype variation is under natural selection and it revitalizes the question as to what processes act to maintain functional mtDNA polymorphism within populations.     

Disequilibrium Pattern Analysis. I. Theory

We have developed a method, disequilibrium pattern analysis, for examining the disequilibrium distribution of the entire array of two locus multiallelic haplotypes in a population. It is shown that a selected haplotype will produce a distinct pattern of linkage disequilibrium values for all generations while the selection is acting. This pattern will also presumably be maintained for many generations after the selection event, until the disequilibrium pattern is eventually broken down by genetic drift and recombination. Related haplotypes, sharing an allele with a selected haplotype, assume a value of linkage disequilibrium proportional to the frequency of the unshared allele and have a single negative value of the normalized linkage disequilibrium. The analysis assumes zero linkage disequilibrium for all allelic combinations initially. The same basic results continue to apply if the selection involves a new mutant, the occurrence of which creates linkage disequilibrium for some haplotypes. The disequilibrium pattern predicted under selection is robust with respect to the influence of migration and random genetic drift. This method is applicable to population data having linked polymorphic loci including that determined from protein or DNA sequencing.     

Genetic similarity of disjunct populations of the giant sea bass Stereolepis gigas

Sequence variation in nuclear and mitochondrial genes of the giant sea bass Stereolepis gigas collected from the Pacific coast and the northern Sea of Cortez was examined. Restriction fragment length polymorphism analysis and direct sequencing showed extremely low mtDNA sequence diversity (13 closely related haplotypes with no evidence of geographical population subdivision). The mitochondrial haplotype mismatch distribution is consistent with a population expansion following the Last Pleistocene glaciation. Differences in single nucleotide polymorphism frequencies between Pacific and Sea of Cortez populations were detected at two of four nuclear loci, which may reflect natural selection or genetic drift in populations with low effective numbers of males. Although Pacific coast and Sea of Cortez populations of giant sea bass do not exhibit the mitochondrial phylogenetic break characteristic of many species with disjunct Pacific and Gulf populations, the possibility of genetic differentiation at nuclear loci suggests that a cautious approach to broodstock selection for captive breeding and restoration programmes be exercised.     

Evidence for Non-Neutrality of Mitochondrial DNA Haplotypes in Drosophila Pseudoobscura

Mitochondrial DNA (mtDNA) haplotypes usually are assumed to be neutral, unselected markers of evolving female lineages. This assumption was tested by monitoring haplotype frequencies in 12 experimental populations of Drosophila pseudoobscura which were polymorphic for mtDNA haplotypes. Populations were maintained for at least 10 generations, and in one case for 32 generations, while tests of mtDNA selective neutrality were conducted. In an initial population, formed from a mixture of two strains with different mitochondrial haplotypes, the frequency of the Bogota haplotype increased 46% in 3 generations, reaching an apparent equilibrium frequency of 82% after 32 generations. Perturbation of this equilibrium by addition of the less common haplotype resulted in a rapid, dramatic increase in frequency of the second haplotype, and a return to essentially the same equilibrium frequency as before perturbation. This behavior is not consistent with mtDNA neutrality, nor is the equilibrium consistent with a simple model of constant selection on the haploid mtDNAs. Replicate cage experiments with mtDNA haplotypes did not always generate the same result as the initial cage. Several lines of evidence, including manipulations of the nuclear genome, support the idea that both nuclear and mitochondrial genomes are involved in the dramatic mtDNA frequency changes. In another experiment, strong female viability selection was implicated via mtDNA frequency changes. Although the causes of the dramatic mtDNA frequency changes in our populations are not obvious, it is clear that Drosophila mitochondrial haplotypes are not always simply neutral markers. Our findings are relevant to the introduction of a novel mtDNA variant from one species or one population into another. Such introductions could be strongly favored by selection, even if it is sporadic.     

Back to basics: using colour polymorphisms to study evolutionary processes

Here, I suggest that colour polymorphic study systems have been underutilized to answer general questions about evolutionary processes, such as morph frequency dynamics between generations and population divergence in morph frequencies. Colour polymorphisms can be used to study fundamental evolutionary processes like frequency‐dependent selection, gene flow, recombination and correlational selection for adaptive character combinations. However, many previous studies of colour polymorphism often suffer from weak connections to population genetic theory. I argue that too much focus has been directed towards noticeable visual traits (colour) at the expense of understanding the evolutionary processes shaping genetic variation and covariation associated with polymorphisms in general. There is thus no need for a specific evolutionary theory for colour polymorphisms beyond the general theory of the maintenance of polymorphisms in spatially or temporally variable environments or through positive or negative frequency‐dependent selection. I outline an integrative research programme incorporating these processes and suggest some fruitful avenues in future investigations of colour polymorphisms.     

THE EFFECT OF POLLEN VERSUS SEED FLOW ON THE MAINTENANCE OF NUCLEAR‐CYTOPLASMIC GYNODIOECY

Gynodioecy, where females co‐occur with hermaphrodites, is a relatively common sexual system in plants that is often the result of a genetic conflict between maternally inherited male sterility genes in the mitochondrial genome and the biparentally inherited male fertility restorer genes in the nucleus. Previous models have shown that nuclear‐cytoplasmic gynodioecy can be maintained under certain conditions by negative frequency‐dependent selection, but the effect of other evolutionary processes such as genetic drift and population subdivision is only partially understood. Here, we investigate the joint effects of frequency‐dependent selection, drift, and migration through either pollen or seeds on the maintenance of nuclear‐cytoplasmic gynodioecy in a subdivided population. We find that the combination of drift and selection causes the loss of gynodioecy under scenarios that would maintain it under the influence of selection alone, and that both seed and, more surprisingly, pollen flow can maintain the polymorphism. In particular, although pollen flow could not avoid the loss of cytoplasmic polymorphism within demes, it allowed the maintenance of nuclear‐cytoplasmic polymorphism at the metapopulation level.     

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium ( r ²   >   0.8) in genomic regions of ∼50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26–113 kb by various definitions, which was larger than that observed in humans (∼10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.     

Evolutionary implications of mitochondrial genetic variation: mitochondrial genetic effects on OXPHOS respiration and mitochondrial quantity change with age and sex in fruit flies

The ancient acquisition of the mitochondrion into the ancestor of modern‐day eukaryotes is thought to have been pivotal in facilitating the evolution of complex life. Mitochondria retain their own diminutive genome, with mitochondrial genes encoding core subunits involved in oxidative phosphorylation. Traditionally, it was assumed that there was little scope for genetic variation to accumulate and be maintained within the mitochondrial genome. However, in the past decade, mitochondrial genetic variation has been routinely tied to the expression of life‐history traits such as fertility, development and longevity. To examine whether these broad‐scale effects on life‐history trait expression might ultimately find their root in mitochondrially mediated effects on core bioenergetic function, we measured the effects of genetic variation across twelve different mitochondrial haplotypes on respiratory capacity and mitochondrial quantity in the fruit fly, Drosophila melanogaster. We used strains of flies that differed only in their mitochondrial haplotype, and tested each sex separately at two different adult ages. Mitochondrial haplotypes affected both respiratory capacity and mitochondrial quantity. However, these effects were highly context‐dependent, with the genetic effects contingent on both the sex and the age of the flies. These sex‐ and age‐specific genetic effects are likely to resonate across the entire organismal life‐history, providing insights into how mitochondrial genetic variation may contribute to sex‐specific trajectories of life‐history evolution.     

Genetic variation in populations of the endangered fish Ladigesocypris ghigii and its implications for conservation

1. The genetic variation of the endangered freshwater fish Ladigesocypris ghigii, endemic to the island of Rhodes (Greece), was investigated for nine populations, originating from seven different stream systems and a reservoir, both at the mtDNA and nuclear level, in order to suggest conservation actions. 2. Both restriction fragment length polymorphism analysis of five segments of mitochondrial DNA (ND‐5/6, COI and 12S‐16S rRNA) amplified by polymerase chain reaction, and random amplified polymorphic DNA analysis, revealed extremely low levels of intra‐population polymorphism. It is highly likely that the low intra‐population variability is the result of successive bottleneck events evident in shrinkage and expansion of the populations year after year, which may have led to a complete loss of several genotypes and haplotypes, and an increased degree of inbreeding. 3. Inter‐population genetic structuring was high, with fixation of haplotypes within six of the nine populations and fixation of alleles within populations originating from different waterbodies. It is probable that all haplotypes and/or alleles found were initially represented in all populations. However, because of the long time of isolation coupled with successive bottleneck and subsequent genetic drift, common mtDNA haplotypes and alleles among the populations may have become rare or extinct through stochastic lineage loss. 4. Although nucleotide divergence among haplotypes was very shallow, half of the haplotypes recorded (three of six), resulted from nucleotide changes on the 12S–16S rRNA segments, which are the most conserved part of the mitochondrial genome. This fact may indicate that the observed genetic variation did not necessarily result only from the retention of ancestral polymorphism, but may have arisen through mutation and complete lineage sorting over a relatively small number of generations, once the populations had become isolated from one another. 5. Our data suggest that two of the L. ghigii populations may be on independent evolutionary trajectories. Considering that each population appears so far well adapted within each site, all populations should be managed and conserved separately.     

Widespread mito‐nuclear discordance with evidence for introgressive hybridization and selective sweeps in Lycaeides

We investigated the extent and potential cause(s) of mitochondrial introgression within the polytypic North American Lycaeides species complex (Lepidoptera). By comparing population genetic structure based on mitochondrial DNA (COI and COII) and nuclear DNA (251 polymorphic amplified fragment length polymorphism markers), we detected substantial mito‐nuclear discordance, primarily involving a single mitochondrial haplotype (h01), which is likely due to mitochondrial introgression between differentiated Lycaeides populations and/or species. We detected reduced mitochondrial genetic diversity relative to nuclear genetic diversity in populations where mitochondrial haplotype h01 occurs, suggesting that the spread of this haplotype was facilitated by selection. We found no evidence that haplotype h01 is associated with increased fitness (in terms of survival to eclosion, fresh adult weight, and adult longevity) in a polymorphic Lycaeides melissa population. However, we did find a positive association between mitochondrial haplotype h01 and infection by the endoparasitic bacterium Wolbachia in one out of three lineages tested. Linkage disequilibrium between mitochondrial haplotype h01 and Wolbachia infection status may have resulted in indirect selection favouring the spread of haplotype h01 in at least one lineage of North American Lycaeides. These results illustrate the potential for introgressive hybridization to produce substantial mito‐nuclear discordance and demonstrate that an individual's mitochondrial and nuclear genome may have strikingly different evolutionary histories resulting from non‐neutral processes and intrinsic differences in the inheritance and biology of these genomes.     

Associations among cytoplasmic molecular markers,gender, and components of fitness in Silene vulgaris,a gynodioecious plant

It has been suggested that the dynamics of chloroplast DNA (cpDNA) or mitochondrial DNA (mtDNA) genetic markers used in studies of plant populations could be influenced by natural selection acting elsewhere in the genome. This could be particularly true in gynodioecious plants if cpDNA or mtDNA genetic markers are in gametic disequilibrium with genes responsible for sex expression. In order to investigate this possibility, a natural population of the gynodioecious plant Silene vulgaris was used to study associations among mtDNA haplotype, cpDNA haplotype, sex and some components of fitness through seed. Individuals were sampled for mtDNA and cpDNA haplotype as determined using restriction fragment length polymorphism (RFLP) methods, sex (female or hermaphrodite), fruit number, fruit set, seeds/fruit and seed germination. The sex of surviving germinating seeds was also noted. All individuals in the population fell into one of two cytoplasmic categories, designated haplotypes f and g by a unique electrophoretic signature in both the mtDNA and cpDNA. The subset of the population carrying haplotype g included a significantly higher proportion of females when compared with the sex ratio of the subset carrying the f haplotype. Haplotype g had a significantly higher fitness when measured by fruit number, fruit set and seeds/fruit, whereas haplotype f had significantly higher fitness when measured by seed germination. Offspring of individuals carrying haplotype g included a significantly greater proportion of females when compared with offspring of individuals carrying the f haplotype. Other studies of gynodioecious plants have shown that females generally have higher fitness through seed than hermaphrodites, but in this study not all fitness differences between haplotypes could be predicted from differences in haplotype-specific sex ratio alone. Rather, some differences in haplotype-specific fitness were due to differences in fitness between individuals of the same sex, but carrying different haplotypes. The results are discussed with regard to the potential for hitchhiking selection to influence the dynamics of the noncoding regions used to designate the cpDNA and mtDNA haplotypes.     

Complex mitonuclear interactions and metabolic costs of mating in male seed beetles

The lack of evolutionary response to selection on mitochondrial genes through males predicts the evolution of nuclear genetic influence on male‐specific mitochondrial function, for example by gene duplication and evolution of sex‐specific expression of paralogs involved in metabolic pathways. Intergenomic epistasis may therefore be a prevalent feature of the genetic architecture of male‐specific organismal function. Here, we assess the role of mitonuclear genetic variation for male metabolic phenotypes [metabolic rate and respiratory quotient (RQ)] associated with ejaculate renewal, in the seed beetle Callosobruchus maculatus, by assaying lines with crossed combinations of distinct mitochondrial haplotypes and nuclear lineages. We found a significant increase in metabolic rate following mating relative to virgin males. Moreover, processes associated with ejaculate renewal showed variation in metabolic rate that was affected by mitonuclear interactions. Mitochondrial haplotype influenced mating‐related changes in RQ, but this pattern varied over time. Mitonuclear genotype and the energy spent during ejaculate production affected the weight of the ejaculate, but the strength of this effect varied across mitochondrial haplotypes showing that the genetic architecture of male‐specific reproductive function is complex. Our findings unveil hitherto underappreciated metabolic costs of mating and ejaculate renewal, and provide the first empirical demonstration of mitonuclear epistasis on male reproductive metabolic processes.     

Differences in population structure estimated within maternally‐ and paternally‐inherited forms of mitochondria in Lampsilis siliquoidea (Bivalvia: Unionidae)

Mussels in several orders possess two separate mitochondrial lineages: a standard female‐inherited form and one inherited only through males. This system of doubly uniparental inheritance (DUI) for mitochondrial genes provides an opportunity to compare the population structure of gene‐lineages passed either mother‐to‐daughter or father‐to‐son. In the present study, we contrast variation in the male and female haplotype lineages of the American freshwater mussel species, Lampsilis siliquoidea (sometimes called Lampsilis radiata luteola), throughout the Lake Erie, Ohio River, and upper Mississippi River watersheds, and contrast variation with the sequences obtained for the related species/subspecies Lampsilis radiata radiata from Maine. The genetic markers were fragments of the cytochrome c oxidase subunit I gene (COI), which occurs in both mitochondrial types, F (female) and M (male). High haplotype diversity was found in the two independent lineages, although purifying selection against amino acid change appeared to be stronger in the female than the male lineage. Phylogeographical patterns also varied between mitochondria passing through females and males. The female lineage exhibited more population structure, with the occurrence of private or nearly‐private haplotypes within two streams, and three others showed restricted haplotype distributions. By contrast to the F‐haplotypes, complex phylogenetic structure occurred for M‐haplotypes, yet this phylogenetic variation coincided with almost no geographical pattern within haplotypes. Basically, F‐haplotypes showed isolation, especially above physical barriers, whereas M‐haplotypes did not. A few individuals in the eastern Lake Erie watershed even possessed M‐haplotypes of an Atlantic Slope (L. radiata radiata) origin, although their F‐haplotypes were typical of Midwestern L. siliquoidea. The finding that mussels package sperm as spermatozuegmata, which float downstream, may underlie greater gene mobility in male‐inherited mitochondria. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 229–240.     

Nuclear–mitochondrial epistasis: a gene's eye view of genomic conflict

We use population genetic models to investigate the cooperative and conflicting synergistic fitness effects between genes from the nucleus and the mitochondrion. By varying fitness parameters, we examine the scope for conflict relative to cooperation among genomes and the utility of the “gene's eye view” analytical approach, which is based on the marginal average fitness of specific alleles. Because sexual conflict can maintain polymorphism of mitochondrial haplotypes, we can explore two types of evolutionary conflict (genomic and sexual) with one epistatic model. We find that the nuclear genetic architecture (autosomal, X‐linked, or Z‐linked) and the mating system change the regions of parameter space corresponding to the evolution by sexual and genomic conflict. For all models, regardless of conflict or cooperation, we find that population mean fitness increases monotonically as evolution proceeds. Moreover, we find that the process of gene frequency change with positive, synergistic fitnesses is self‐accelerating, as the success of an allele in one genome or in one sex increases the frequency of the interacting allele upon which its success depends. This results in runaway evolutionary dynamics caused by the positive intergenomic associations generated by selection. An inbreeding mating system tends to further accelerate these runaway dynamics because it maintains favorable host–symbiont or male–female gene combinations. In contrast, where conflict predominates, the success of an allele in one genome or in one sex diminishes the frequency of the corresponding allele in the other, resulting in considerably slower evolutionary dynamics. The rate of change of mean fitness is also much faster with positive, synergistic fitnesses and much slower where conflict is predominant. Consequently, selection rapidly fixes cooperative gene combinations, while leaving behind a slowing evolving residue of conflicting gene combinations at mutation–selection balance. We discuss how an emphasis on marginal fitness averages may obscure the interdependence of allelic fitness across genomes, making the evolutionary trajectories appear independent of one another when they are not.     

Mitochondrial DNA haplotype frequencies in natural and experimental populations of Drosophila subobscura.

The evolution of Drosophila subobscura mitochondrial DNA has been studied in experimental populations, founded with flies from a natural population from Esporles (Majorca, Balearic Islands, Spain). This population, like other European ones, is characterized by the presence of two very common (>96%) mitochondrial haplotypes (called I and II) and rare and endemic haplotypes that appear at very low frequencies. There is no statistical evidence of positive Darwinian selection acting on the mitochondrial DNA variants according to Tajima''s neutrality test. Two experimental populations, with one replicate each, were established with flies having a heterogeneous nuclear genetic background, which was representative of the composition of the natural population. Both populations were started with the two most frequent mitochondrial haplotypes, but at different initial frequencies. After 13 to 16 generations, haplotype II reached fixation in three cages and its frequency was 0.89 by generation 25 in the fourth cage. Random drift can be rejected as the force responsible for the observed changes in haplotype frequencies. There is not only statistical evidence of a linear trend favoring a mtDNA (haploid) fitness effect, but also of a significant nonlinear deviation that could be due to a nuclear component.     

Mitochondrial gene diversity in Cepaea– population structure, history and positive selection

Previous studies on the geographical distributions of both shell polymorphisms and mitochondrial haplotypes in Cepaea land snails have suggested varying contributions of natural selection, random genetic drift and population history to the origin and persistence of variation. We combine previous studies of polymorphism within two species of Cepaea , with new molecular data from two mitochondrial genes. The distributions of mitochondrial variants suggest that elements of population history may have been influential in creating patterns of diversity. However, some patterns of amino acid substitution that are found in a protein coding gene are also consistent with the action of diversifying selection. This suggests the possibility that in addition to stochastic processes such as repeated founder events, local extinctions and random genetic drift within a structured population, adaptive molecular change may have affected mitochondrial diversity.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society, 2006, 87 , 167–184.     

Comparison between mitochondrial and chloroplast DNA variation in the native range of Silene vulgaris

A detailed survey of mitochondrial and chloroplast diversity in eight populations of Silene vulgaris from Central Europe was conducted for comparison with previously published data on diversity from S. vulgaris populations in the introduced range. Mitochondrial DNA (mtDNA) variation around the coxI gene was assessed with Southern blotting/restriction fragment length polymorphism methods. Chloroplast variation was assessed by sequencing the intergenic spacer separating the trnH and psbA genes. Thirty mtDNA haplotypes and 24 chloroplast DNA (cpDNA) haplotypes were found within 86 individuals. The overall genetic diversity h (0.941 for mitochondrial, and 0.893 for chloroplast markers) and within-population diversity were higher than reported in previous population studies of S. vulgaris in the USA and Europe. The frequency of private alleles was surprisingly high - more than 90% for both kinds of markers. Most of our populations were large and located in relatively undisturbed meadows, whereas surveys in Virginia consisted of smaller roadside populations. The slow rate of population turnover in European populations is discussed as a factor responsible for the relatively high diversity of S. vulgaris in undisturbed areas of its native range. Association between mtDNA and cpDNA haplotypes was also demonstrated. Finally, gender and mtDNA haplotype were associated in the Alps populations, where females were very rare.     

Mitochondrial DNA origins of the Latvian clefting population

Latvia has one of the highest prevalence of isolated cleft lip with or without cleft palate (CL/P) in Europe. To clarify the genetic origins of the Latvian cleft population and establish a method for genetic mapping, mitochondrial DNA variation was studied in a population affected with clefting. One-hundred and seven subjects and 351 samples from unrelated healthy volunteers representing four anthropologically, archaeologically and ethno-linguistically different regions of Latvia were selected. The case group showed a higher frequency of haplogroups U4 (p=0.02) and U5 (p=0.0003) than in non-U haplogroups. We hypothesize that U4 and U5 mtDNA haplotype carriers may also carry susceptibility genes for clefts. Future studies will take into consideration these definitions based on mtDNA haplotypes when analyzing genetic variations and their possible contribution to CL/P susceptibility.