Parallel rate heterogeneity in chloroplast and mitochondrial genomes of Brazil nut trees (Lecythidaceae) is consistent with lineage effects

We investigated whether relative rates of divergence were correlated between the mitochondrial and chloroplast genomes as expected under lineage effects or were genome specific as expected with locus-specific effects. Five mitochondrial noncoding regions (nad1B_C, nad4exon1_2, nad7exon2_3, nad7exon3_4, and rps14-cob) for 21 samples from Lecythidaceae were sequenced. Three chloroplast regions (rpl20-5'rps12, trnS-trnG, and psbA-trnH) were sequenced to expand the taxa in an existing data set. Absolute rates of nucleotide and insertion and deletion (indel) changes were 13 times faster in the chloroplast genome than in the mitochondrial genome. Similar indel length frequency distributions for both organelles suggested that common mechanisms were responsible for generating indels. Molecular clock tests applied to phylogenetic trees estimated from mitochondrial and chloroplast sequences revealed global rate heterogeneity of nucleotide substitution. Maximum likelihood and Tajima's 1D relative rate tests show that Lecythis zabucajo exhibited a rate acceleration for both the mitochondrial and chloroplast sequences. Whereas Eschweilera romeu-cardosoi showed a significant rate slowdown for chloroplast sequences, the mitochondrial sequences for 3 Eschweilera taxa showed evidence for a rate slowdown only when compared with L. zabucajo. Significant rate heterogeneity was also observed for indel changes in the mitochondrial genome but not for the chloroplast. The lack of mitochondrial nucleotide changes for some taxa as well as chloroplast indel homoplasy may have limited the power of relative rate tests to detect rate variation. Relative ratio tests consistently indicated rate proportionality among branch lengths between the mitochondrial and chloroplast phylogenetic trees. The relative ratio tests showed that taxa possessing rate heterogeneity had parallel relative divergence rates in both mitochondrial and chloroplast sequences as expected under lineage effects. A neutral replication-dependent model of rate heterogeneity for both nucleotide and indel changes provides a simple explanation for common patterns of rate heterogeneity across the 2 organelle genomes in Lecythidaceae. The lineage effects observed here were uncoupled from annual/perennial habit because all the species from this study are perennial.     

Molecular evolution of insertions and deletion in the chloroplast genome of silene

Insertions, deletions, and inversions in the chloroplast genome of higher plants have been shown to be extremely useful for resolving phylogenetic relationships both between closely related taxa and among more basal lineages. Introns and intergenic spacers from the chloroplast genome are now increasingly used for phylogenetic and population genetic studies of populations from a single species, and it is therefore interesting to know whether indels can provide useful data and hence increase the power of intraspecific studies. Here, we show that indels in three cpDNA intergenic spacers and one cpDNA intron for two species of Silene evolve at slightly higher rates than base pair substitutions. Repeat indels appear to have the highest rate of evolution and are thus more prone to homoplasy. We show that coded indel data have high information content for phylogenetic analysis, and indels thus provide useful information to infer phylogenetic relationships at the intraspecific level.     

Evolution of the mitochondrial rps3 intron in perennial and annual angiosperms and homology to nad5 intron 1.

The plant mitochondrial rps3 intron was analyzed for substitution and indel rate variation among 15 monocot and dicot angiosperms from 10 genera, including perennial and annual taxa. Overall, the intron sequence was very conserved among angiosperms. Based on length polymorphism, 10 different alleles were identified among the 10 genera. These allelic differences were mainly attributable to large indels. An insertion of 133 nucleotides, observed in the Alnus intron was partially or completely absent in the other lineages of the family Betulaceae. This insertion was located within domain IV of the secondary-structure model of this group IIA intron. A mobile element of 47 nucleotides that showed homology to sequences located in rice rps3 intron and in intergenic plant mitochondrial genomes was found within this insertion. Both substitution and indel rates were low among the Betulaceae sequences, but substitution rates were increasingly larger than indel rates in comparisons involving more distantly related taxa. From a secondary-structure model, regions involved in helical structures were shown to be well preserved from indels as compared to substitutions, but compensatory changes were not observed among the angiosperm sequences analyzed. Using approximate divergence times based on the fossil record, substitution and indel rate heterogeneity was observed between different pairs of annual and perennial taxa. In particular, the annual petunia and primrose evolved more than 15 and 10 times faster, for substitution and indel rates respectively, than the perennial birch and alder. This is the first demonstration of an evolutionary rate difference between perennial and annual forms in noncoding DNA, lending support to neutral causes such as the generation time, population size, and speciation rate effects to explain such rate heterogeneity. Surprisingly, the sequence from the rps3 intron had a high identity with the sequence of intron 1 from the angiosperm mitochondrial nad5 gene, suggesting a common origin of these two group IIA introns.     

The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis

Chloroplast DNA sequences are a primary source of data for plant molecular systematic studies. A few key papers have provided the molecular systematics community with universal primer pairs for noncoding regions that have dominated the field, namely trnL-trnF and trnK/matK. These two regions have provided adequate information to resolve species relationships in some taxa, but often provide little resolution at low taxonomic levels. To obtain better phylogenetic resolution, sequence data from these regions are often coupled with other sequence data. Choosing an appropriate cpDNA region for phylogenetic investigation is difficult because of the scarcity of information about the tempo of evolutionary rates among different noncoding cpDNA regions. The focus of this investigation was to determine whether there is any predictable rate heterogeneity among 21 noncoding cpDNA regions identified as phylogenetically useful at low levels. To test for rate heterogeneity among the different cpDNA regions, we used three species from each of 10 groups representing eight major phylogenetic lineages of phanerogams. The results of this study clearly show that a survey using as few as three representative taxa can be predictive of the amount of phylogenetic information offered by a cpDNA region and that rate heterogeneity exists among noncoding cpDNA regions.     

Pattern and rate of indel evolution inferred from whole chloroplast intergenic regions in sugarcane, maize and rice.

Microstructural changes such as insertions and deletions (=indels) are a major driving force in the evolution of non-coding DNA sequences. To better understand the mechanisms by which indel mutations arise, as well as the molecular evolution of non-coding regions, the number and pattern of indels and nucleotide substitutions were compared in the whole chloroplast genomes. Comparisons were made for a total of over 38 kb non-coding DNA sequences from 126 intergenic regions in two data sets representing species with different divergence times: sugarcane and maize and Oryza sativa var. indica and japonica. The main findings of this study are: (i) Approximately half of all indels are single nucleotide indels. This observation agrees with previous studies in various organisms. (ii) The distribution and number of indels was different between two data sets, and different patterns were observed for tandem repeat and non-repeat indels. (iii) Distribution pattern of tandem repeat indels showed statistically significant bias towards A/T-rich. (iv) The rate of indel mutation was estimated to be approximately 0.8 +/- 0.04 x 10(-9) per site per year, which was similar to previous estimates in other organisms. (v) The frequencies of nucleotide substitutions and indels were significantly lower in inverted repeat (IR).     

Sequence variation at cpDNA regions of watermelon and related wild species: implications for the evolution of Citrullus haplotypes

Sequencing analysis of one coding and four noncoding cpDNA regions was conducted to infer biogeographic and evolutionary relationships in the genus Citrullus. Eighteen taxa from diverse geographical areas were included. A low number of parsimony informative characters (1.1%) was observed at the ～4 kb section of cpDNA. Variability within Citrullus was detected primarily at noncoding regions of high A + T content. Substitution rates varied from 0-0.48% for ndhF with A + T content of 68.4% to 0.39- 1.69% for the intergenic region of atpA with A + T content of 82.8%, mainly resulting in indels and transversions. Indels at several regions acted as valuable parsimony informative markers. Citrullus lanatus var. lanatus, the cultivated watermelon, and C. ecirrhosus and C. rehmii from Namibia, lacked molecular variability. The genus Citrullus is supported monophyletically and shows two main clades, one of which contains C. colocynthis. In the other clade, C. rehmii is sister to a clade containing C. ecirrhosus and C. lanatus. Two clades were recovered within C. lanatus, consisting of domesticated watermelon and wild citron, var. citroides. Five haplotypes within C. colocynthis were used to deduce colonization routes of the species. Biogeographic patterns point to separate colonization events into Africa and the Far East.     

The pattern of insertion/deletion polymorphism in Arabidopsis thaliana

Little is known about variation of nucleotide insertion/deletions (indels) within species. In Arabidopsis thaliana, we investigated indel polymorphism patterns between two genome sequences and among 96 accessions at 1215 loci. Our study identified patterns in the variation of indel density, size, GC content and distribution, and a correlation between indels and substitutions. We found that the GC content in indel sequences was lower than that in non-indel sequences and that indels typically occur in regions with lower GC content. Patterns of indel frequency distribution among populations were more consistent with neutral expectation than substitution patterns. We also found that the local level of substitutions is positively correlated with indel density and negatively correlated with their distance to the closed indel, suggesting that indels play an important role in nucleotide variation.     

The evolution of the atpbeta-rbcL intergenic spacer in the epacrids (Ericales) and its systematic and evolutionary implications

Sequence data from the noncoding region separating the plastid genes atpbeta and rbcL were gathered for 27 epacrid taxa, representing all previously recognized infrafamilial groups, and four outgroup taxa (Ericaceae), to address several persistent phylogenetic questions in the group. Parsimony analyses were conducted on these data, as well as on a complementary rbcL sequence dataset assembled from the literature and the combined dataset. The atpbeta-rbcL spacer was notable for the high frequency of insertion-deletion mutations (indels); their distributions were coded as binary characters and included as a adjunct matrix in some of the analyses. The phylogenetic patterns derived from the spacer and rbcL data and the combined analyses, both including and excluding the indel data, concur in resolving seven major lineages corresponding to the tribes of Crayn et al. (1998, Aust. J. Bot. 46, 187-200), viz. Prionoteae, Archerieae, Oligarrheneae, Cosmelieae, Richeeae, Epacrideae, and Styphelieae. The relationships of the tribes and within Styphelieae, however, are not convincingly resolved. Minor conflicts in the positions of some taxa between the spacer and the rbcL trees are poorly supported. Among epacrids, the spacer region provided more cladistically informative characters than rbcL and resulted in trees with lower homoplasy. Further, the spacer data, when analyzed alone and when combined with rbcL, resolved several clades that could not be retrieved on rbcL data alone and provided increased support for many other relationships. The evolution of a putative three-base inversion associated with a hairpin secondary structure in the spacer region is discussed in the light of the inferred phylogeny.     

Genetic differentiation of relictual populations of Alsophila spinulosa in southern China inferred from cpDNA trnL-F noncoding sequences

The genetic differentiation and phylogeographical pattern of 11 relictual populations of Alsophila spinulosa distributed across Hainan, Guangdong, and Guangxi in southern China were inferred from sequence variations of trnL-F noncoding regions of chloroplast DNA (cpDNA). The length of trnL-F noncoding sequences varied from 863 to 940 bp. The A + T content was 62.23-63.36%. Sequences were neutral in terms of evolution (Tajima's criterion D=-0.62417, P>0.10 and Fu and Li's test D*=-1.45455, P>0.10; F*=-1.32798, P>0.10). Thirty-four haplotypes were identified based on nucleotide variation. Relatively high levels of haplotype diversity (h=0.929) and nucleotide diversity (Dij=0.022263) were detected in A. spinulosa, probably associated with its long evolutionary history which allowed the accumulation of genetic variation within lineages. Both the minimum spanning network and the strict consensus tree of the most parsimonious trees generated for haplotypes demonstrated that the investigated populations of A. spinulosa were subdivided into two geographical groups: Hainan and Guangdong-Guangxi. An analysis of molecular variance (AMOVA) indicated that most of the genetic variation (87.48%, P<0.001) was partitioned among regions. Spatial structure measurements revealed that population genetic structure was not related to geographical distance. This research suggests that blocked gene flow by Qiongzhou strait and an inbreeding system might result in the geographical subdivision between Hainan and Guangdong-Guangxi (F(ST)=0.92, Nm=0.09). Within each region, the "star like" pattern of phylogeography of haplotypes implied a population expansion process during evolutionary history. Gene genealogies together with coalescent theory were useful tools for uncovering the phylogeography of A. spinulosa.     

Molecular phylogenetics of the lizard genus Microlophus (squamata:tropiduridae): aligning and retrieving indel signal from nuclear introns

We use a multigene data set (the mitochondrial locus and nine nuclear gene regions) to test phylogenetic relationships in the South American "lava lizards" (genus Microlophus) and describe a strategy for aligning noncoding sequences that accounts for differences in tempo and class of mutational events. We focus on seven nuclear introns that vary in size and frequency of multibase length mutations (i.e., indels) and present a manual alignment strategy that incorporates insertions and deletions (indels) for each intron. Our method is based on mechanistic explanations of intron evolution that does not require a guide tree. We also use a progressive alignment algorithm (Probabilistic Alignment Kit; PRANK) and distinguishes insertions from deletions and avoids the "gapcost" conundrum. We describe an approach to selecting a guide tree purged of ambiguously aligned regions and use this to refine PRANK performance. We show that although manual alignment is successful in finding repeat motifs and the most obvious indels, some regions can only be subjectively aligned, and there are limits to the size and complexity of a data matrix for which this approach can be taken. PRANK alignments identified more parsimony-informative indels while simultaneously increasing nucleotide identity in conserved sequence blocks flanking the indel regions. When comparing manual and PRANK with two widely used methods (CLUSTAL, MUSCLE) for the alignment of the most length-variable intron, only PRANK recovered a tree congruent at deeper nodes with the combined data tree inferred from all nuclear gene regions. We take this concordance as an objective function of alignment quality and present a strongly supported phylogenetic hypothesis for Microlophus relationships. From this hypothesis we show that (1) a coded indel data partition derived from the PRANK alignment contributed significantly to nodal support and (2) the indel data set permitted detection of significant conflict between mitochondrial and nuclear data partitions, which we hypothesize arose from secondary contact of distantly related taxa, followed by hybridization and mtDNA introgression.     

The hierarchical spatial distribution of chloroplast DNA polymorphism across the introduced range of Silene vulgaris

Silene vulgaris was introduced into North America sometime prior to 1800. In order to document the population structure that has developed since that time, collections were made from 56 local populations distributed among 9 geographical regions in eastern North America. Individual plants were characterized for chloroplast DNA (cpDNA) haplotype by restriction fragment size analysis of four noncoding regions of cpDNA amplified by polymerase chain reaction. A total of 19 cpDNA haplotypes were detected using this method. The overall gene diversity of 0.85 is quite similar to the diversity detected in these same regions of cpDNA in a previously published sample of S. vulgaris taken from across much of Europe. The spatial distribution of the North American cpDNA diversity was quantified by hierarchical F-statistics that partitioned the genetic variance into variation among local populations within regions, and variation among regions. The average FST among populations within regions was 0.66 and the FST among regions was 0.09. The among-region variation was due to both differences among regions in the frequency of two most common haplotypes, and to the presence of a number of region-specific haplotypes. In order to test for isolation by distance at the regional level, FST values were calculated for all possible pairs of regions, and regressed against the geographical distance between those regions. There was no evidence for isolation by distance. It is suggested that the local population structure is generated by recent extinction/colonization dynamics, and that the among-region structure reflects demographic events associated with range expansion following introduction to North America.     

基于叶绿体DNA变异的湖北梨属种质系统进化及遗传多样性分析

利用trnL intron、trnL-trnF、trnS-psbC和accD-psa I等4个叶绿体DNA片段对来自湖北省的88份梨属种质资源进行系统进化和遗传多样性分析。结果表明,4个cpDNA片段共检测到变异位点11个,其中单一突变位点6个,插入/缺失(Indel)位点5个。acc D-psa I多态性最高,其变异位点数、核苷酸多态性和单倍型多样性均为最高。供试梨种质的核苷酸多样性和单倍型多样性分别为0.00112和0.769;Tajima's D检验值在P0.10水平上均不显著,表明所检测的4个区域以及合并后的片段均遵循中性进化模型;4个序列合并共检测到叶绿体单倍型10个,其中兴山梨种质中检测到的单倍型最多,荆门其次;Hap2和Hap5是2个主要单倍型,分别占总样本数的31.82%和30.68%;中介邻接网络图显示东方梨和西洋梨独立进化,而较为原始的稀有单倍型Hap8和Hap9均位于荆门,暗示该地区可能为砂梨的起源中心或多样性中心之一。     

Phylogeography of Japanese water crowfoot based on chloroplast DNA haplotypes

Water crowfoot, Ranunculus subgenus Batrachium, are submerged macrophytes in the Ranunculaceae. We aimed to infer the Quaternary history of these aquatic macrophytes in the Japanese archipelago. We studied 212 individuals of three perennial and one annual species from 46 populations covering the entire geographic range and found eight haplotypes based on approximately 1800 bp of four spacers in chloroplast DNA. The relationships among haplotypes were resolved using maximum parsimony and parsimony network analyses. To identify the zones of clear genetic boundaries, Monmonier's algorithm was used. The eight haplotypes were distinguished from adjacent haplotypes by one substitution or indel. Each of the 46 populations was fixed for a single haplotype. Inconsistency between cpDNA haplotypes and Batrachium taxa was found, except in one annual species. The distribution of the haplotypes in perennial species was highly geographically structured. An abrupt genetic change was detected between the Tohoku region and more southerly regions. Since the perennial Batrachium are cold-adapted, this genetic differentiation may be due to historical changes in their distributions caused by Quaternary climatic oscillations: interglacial retreats in colder refugia and glacial range expansions. The single haplotype composition of each population may have been shaped by founder effects during colonization and/or by genetic drift. The interglacial refugial populations must have been small enough to deplete haplotype diversity. Inconsistency between cpDNA haplotypes and Batrachium taxa may be due to incomplete lineage sorting of ancestral polymorphic haplotypes.     

SHORT INDELS ARE SUBJECT TO INSERTION‐BIASED GENE CONVERSION

Recombination between homologous loci is accompanied by formation of heteroduplexes. Repairing mismatches in heteroduplexes often leads to single nucleotide substitutions in a process known as gene conversion. Gene conversion was shown to be GC‐biased in different organisms; that is, a W(A or T)→S(G or C) substitution is more likely in this process than a S→W substitution. Here, we show that the insertion/deletion ratio for short noncoding indels that reach fixation between species is positively correlated with the recombination rate in Drosophila melanogaster, Homo sapiens, and Saccharomyces cerevisiae. This correlation is both due to an increase of the fixation rate of insertions and decrease of the fixation rate of deletions in regions of high recombination. Whole‐genome data on indel polymorphism and divergence in D. melanogaster rule out mutation biases and selection as the cause of this trend, pointing to insertion‐biased gene conversion as the most likely explanation. The bias toward insertions is the strongest for single‐nucleotide indels, and decreases with indel length. In regions of high recombination rate this bias leads to an up to ～5‐fold excess of fixed short insertions over deletions, and substantially affects the evolution of DNA segments.     

Comparative sequencing of a microsatellite locus reveals size homoplasy within and between european oak species (<Emphasis Type="Italic">Quercus</Emphasis> spp.)

Microsatellites (simple sequence repeats [SSRs]) are highly variable molecular markers that are a rich and readily assayed source of variation for population genetic studies. Cross-amplification between closely related species is possible when there are no (or few) sequence differences in the primer binding sites. The occurrence of nonhomologous fragments of the same size (size homoplasy) is a contraint of microsatellites. Size homoplasy can be caused by insertions/deletions (indels) in SSR flanking regions. We found that size variation in locus ssrQZAG9 is due to different repeat numbers of the SSR motifs but also to indels in SSR flanking regions. Indels were found within species belonging to sectionsRobur andCerris of genusQuercus and also between species of the 2 sections. In sectionRobur (Quercis robur L.,Quercus petraea [Matt.] Liebl.,Quercus pubescens Willd.), we detected rare alleles with an indel of 57 bp or 62 bp followed by a smaller indel of 12 bp in the SSR flanking regions. These alleles show a size range overlapping with that of alleles amplified inQuercus cerris L. (sectionCerris). Multiple alignments with sequences of sectionRobur revealed the same SSR repeat motif but multiple indels in SSR flanking regions inQ. cerris. We discuss the effects of size homoplasy of SSR loci for the study of interspecific gene flow and on estimates of population differentiation.     

Chloroplast DNA variation and geographical structure of the Aristolochia kaempferi group (Aristolochiaceae)

The present study documents cpDNA variation in the Aristolochia kaempferi group (Aristolochiaceae), which consists of one Chinese and all Japanese and Taiwanese species of the subgenus Siphisia. In a phylogenetic analysis based on the nucleotide sequences of the matK gene, and the atpB-rbcL and trnS-trnG intergenic spacer regions, 38 haplotypes were recognized in the A. kaempferi group and as many as 24 within A. kaempferi. This is the most haplotypes reported for a single species to date. Although six highly significant major clades were identified in the phylogenetic analysis, they were not congruent with previous classifications. This might be attributed to the specific speciation process, such as convergent evolution, incomplete lineage sorting, and/or reticulate evolution. The six major clades had a clear geographical distribution pattern and were significantly associated with geographical distribution of haplotypes in a nested clade analysis and AMOVA. The results allow us to deduce a scenario in which multiple contractions and expansions of the geographical ranges brought about by Quaternary climatic oscillations affected the patterns of genetic diversity. The present geographic patterns of haplotype distribution within the A. kaempferi group can be explained by the last postglacial range expansion from different refugia, and the boundaries may be suture zones.     

Phylogeography of the Arctic-Alpine Saxifraga oppositifolia (Saxifragaceae) and some related taxa based on cpDNA and ITS sequence variation

Saxifraga oppositifolia (Saxifragaceae) is an important model system for the evolution of Arctic-Alpine plant species. Sequences of the psbA-trnH intergenic spacer of chloroplast DNA and of the internal transcribed spacer region, ITS1-5.8S-ITS2, of the nuclear ribosomal DNA were used to investigate the intraspecific evolution and phylogeography of this species. Samples from nearly the species' entire circumpolar distribution were included in the analysis as well as samples from the closely related taxa S. smalliana from Alaska and S. blepharophylla and S. biflora from the Alps (S. aizoides served as outgroup). These latter taxa showed a low number of parsimony informative characters, in both cpDNA and ITS sequences, which separated them from S. oppositifolia. Two main cpDNA haplotypes were detected within S. oppositifolia, one with a Eurasian distribution and one with an East Asian-North American distribution. This confirmed the existence of two cpDNA lineages with different geographical distributions in this species, which had previously been reported based on a RFLP analysis. The ITS phylogeny was not useful with respect to the intraspecific evolution and phylogeography of S. oppositifolia, because it showed a largely unresolved topology with low statistical support. The cpDNA sequence analysis, however, also suggested a putative long-distance dispersal event. All investigated taxa had cpDNA haplotypes that were congruent with their geographical origin rather than their phylogeny. This could point to putative lineage sorting in S. oppositifolia and related taxa.     

Molecular population genetics of sequence length diversity in the Adh region of Drosophila pseudoobscura

Positive and negative selection on indel variation may explain the correlation between intron length and recombination levels in natural populations of Drosophila. A nucleotide sequence analysis of the 3.5 kilobase sequence of the alcohol dehydrogenase (Adh) region from 139 Drosophila pseudoobscura strains and one D. miranda strain was used to determine whether positive or negative selection acts on indel variation in a gene that experiences high levels of recombination. A total of 30 deletion and 36 insertion polymorphisms were segregating within D. pseudoobscura populations and no indels were fixed between D. pseudoobscura and its two sibling species D. miranda and D. persimilis. The ratio of Tajima's D to its theoretical minimum value (D(min)) was proposed as a metric to assess the heterogeneity in D among D. pseudoobscura loci when the number of segregating sites differs among loci. The magnitude of the D/D(min) ratio was found to increase as the rate of population expansion increases, allowing one to assess which loci have an excess of rare variants due to population expansion versus purifying selection. D. pseudoobscura populations appear to have had modest increases in size accounting for some of the observed excess of rare variants. The D/D(min) ratio rejected a neutral model for deletion polymorphisms. Linkage disequilibrium among pairs of indels was greater than between pairs of segregating nucleotides. These results suggest that purifying selection removes deletion variation from intron sequences, but not insertion polymorphisms. Genome rearrangement and size-dependent intron evolution are proposed as mechanisms that limit runaway intron expansion.     

The linkage disequilibrium between chloroplast DNA and mitochondrial DNA haplotypes in Beta vulgaris ssp. maritima (L.): the usefulness of both genomes for population genetic studies

The structure and evolution of the plant mitochondrial genome may allow recurrent appearance of the same mitochondrial variants in different populations. Whether the same mitochondrial variant is distributed by migration or appears recurrently by mutation (creating homoplasy) in different populations is an important question with regard to the use of these markers for population genetic analyses. The genetic association observed between chloroplasts and mitochondria (i.e. two maternally inherited cytoplasmic genomes) may indicate whether or not homoplasy occurs in the mitochondrial genome. Four-hundred and fourteen individuals sampled in wild populations of beets from France and Spain were screened for their mitochondrial and chloroplast polymorphisms. Mitochondrial DNA (mtDNA) polymorphism was investigated with restriction fragment length polymorphism (RFLP) and chloroplast DNA (cpDNA) polymorphism was investigated with polymerase chain reaction PCR-RFLP, using universal primers for the amplification. Twenty and 13 variants for mtDNA and cpDNA were observed, respectively. Most exhibited a widespread geographical distribution. As a very strong linkage disequilibrium was estimated between mtDNA and cpDNA haplotypes, a high rate of recurrent mutation was excluded for the mitochondrial genome of beets. Identical mitochondrial variants found in populations of different regions probably occurred as a result of migration. We concluded from this study that mtDNA is a tool as valuable as cpDNA when a maternal marker is needed for population genetics analyses in beet on a large regional scale.