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

Mitochondrial phylogeny and distribution of cytoplasmic male sterility-associated genes in Beta vulgaris

Cytoplasmic male sterility (CMS) is a mitochondrial-encoded trait that confers reproductive defects in males but not in females or any vegetative function. Why CMS is so often found in plants should be investigated from the viewpoint of mitochondrial phylogeny. Beta vulgaris, including the wild subspecies maritima and cultivated subspecies vulgaris (e.g., sugar beet), is known to be mitochondrially polymorphic, from which multiple CMS mitochondria have been found, but their evolutionary relationship has been obscure. We first refined the B. vulgaris reference mitochondrial genome to conduct a more accurate phylogenetic study. We identified mitochondrial single-nucleotide polymorphic sites from 600 B. vulgaris accessions. Principal component analysis, hierarchical clustering analysis, and creation of a phylogenetic tree consistently suggested that B. vulgaris mitochondria can be classified into several groups whose geographical distribution tends to be biased toward either the Atlantic or Mediterranean coasts. We examined the distribution of CMS-associated mitochondrial genes from Owen, E- and G-type CMS mitochondria. About one-third of cultivated beets had Owen-type CMS, which reflects the prevalence of using Owen-type CMS in hybrid breeding. Occurrence frequencies for each of the three CMS genes in wild beet were less than 4%. CMS genes were tightly associated with specific mitochondrial groups that are phylogenetically distinct, suggesting their independent origin. However, homologous sequences of the Owen type CMS gene occurred in several different mitochondrial groups, for which an intricate explanation is necessary. Whereas the origin of cultivated beet had been presumed to be Greece, we found an absence of Owen-type mitochondria in Greek accessions.

Utility of the Mitochondrial Genome in Plant Taxonomic Studies

Size, structure, and sequence content lability of plant mitochondrial genome (mtDNA) across species has sharply limited its use in taxonomic studies. Historically, mtDNA variation has been first investigated with RFLPs, while the development of universal primers then allowed studying sequence polymorphisms within short genomic regions (<3 kb). The recent advent of NGS technologies now offers new opportunities by greatly facilitating the assembly of longer mtDNA regions, and even full mitogenomes. Phylogenetic works aiming at comparing signals from different genomic compartments (i.e., nucleus, chloroplast, and mitochondria) have been developed on a few plant lineages, and have been shown especially relevant in groups with contrasted inheritance of organelle genomes. This chapter first reviews the main characteristics of mtDNA and the application offered in taxonomic studies. It then presents tips for best sequencing protocol based on NGS data to be routinely used in mtDNA-based phylogenetic studies.

Beyond balancing selection: frequent mitochondrial recombination contributes to high-female frequencies in gynodioecious Lobelia siphilitica (Campanulaceae)

Gynodioecy is a sexual system in which females and hermaphrodites co-occur. In most gynodioecious angiosperms, sex is determined by an interaction between mitochondrial male-sterility genes (CMS) that arise via recombination and nuclear restorer alleles that evolve to suppress them. In theory, gynodioecy occurs when multiple CMS types are maintained at equilibrium frequencies by balancing selection. However, some gynodioecious populations contain very high frequencies of females. High female frequencies are not expected under balancing selection, but could be explained by the repeated introduction of novel CMS types. To test for balancing selection and/or the repeated introduction of novel CMS, we characterised cytoplasmic haplotypes from 61 populations of Lobelia siphilitica that vary widely in female frequency. We confirmed that mitotype diversity and female frequency were positively correlated across populations, consistent with balancing selection. However, while low-female populations hosted mostly common mitotypes, high-female populations and female plants hosted mostly rare, recombinant mitotypes likely to carry novel CMS types. Our results suggest that balancing selection maintains established CMS types across this species, but extreme female frequencies result from frequent invasion by novel CMS types. We conclude that balancing selection alone cannot account for extreme population sex-ratio variation within a gynodioecious species.

Heteroplasmy and Patterns of Cytonuclear Linkage Disequilibrium in Wild Carrot

Organellar genomes are considered to be strictly uniparentally-inherited. Uniparental inheritance allows for cytonuclear coevolution and the development of highly coordinated cytonuclear interactions. Yet, instances of biparental inheritance have been documented across eukaryotes. Biparental inheritance in otherwise uniparentally-inherited organelles is termed leakage (maternal or paternal) and allows for the presence of multiple variants of the same organellar genome within an individual, called heteroplasmy. It is unclear what, if any, evolutionary consequences are placed on nuclear and/or organellar genomes due to heteroplasmy. One way of accessing cytonuclear interactions and potential coevolution is through calculating cytonuclear linkage disequilibrium (cnLD), or the non-random association of alleles between nuclear and organellar genomes. Patterns of cnLD can indicate positive or negative cytonuclear selection, coevolution between the nuclear and organellar genomes, non-traditional organellar inheritance, or instances of ancestral heteroplasmy. In plants, cytonuclear interactions have been shown to play a role in cytoplasmic male sterility which occurs in gynodioecious species and is associated with leakage. We used the gynodioecious species, Daucus carota L. spp. carota, or wild carrot, to investigate cnLD. We genotyped a total of 265 individuals from two regions of the USA at 15 nuclear microsatellites, the mitochondrial genes cox1 and atp9, and an intergenic region between trnS and trnG (StoG) in the plastid genome to calculate nuclear–nuclear LD (nucLD), cnLD, and organellar LD (i.e., within the mtDNA and between mtDNA and ptDNA) within the two regions. We were further able to identify cox1 and StoG heteroplasmy and calculate some of the same LD measures within heteroplasmic and homoplasmic (non-heteroplasmic) datasets. We used a Z-transformation test to demonstrate that heteroplasmic individuals display significantly higher levels of cnLD within both regions. In spite of this, within and between organellar LD is low to moderate. Given these patterns of LD in two regions of the USA in which gene flow has been shown to occur between crop and wild carrot, we suggest that heteroplasmy is an evolutionary mechanism which permits the maintenance of cnLD while also acting to disrupt organellar LD.

Chloroplastic and nuclear diversity of wild beets at a large geographical scale: Insights into the evolutionary history of the Beta section

Historical demographic processes and mating systems are believed to be major factors in the shaping of the intraspecies genetic diversity of plants. Among Caryophyllales, the Beta section of the genus Beta, within the Amaranthaceae/Chenopodiaceae alliance, is an interesting study model with species and subspecies (Beta macrocarpa, Beta patula, Beta vulgaris maritima and B.v. adanensis) differing in geographical distribution and mating system. In addition, one of the species, B. macrocarpa, mainly diploid, varies in its level of ploidy with a tetraploid cytotype described in the Canary Islands and in Portugal. In this study, we analyzed the nucleotide diversity of chloroplastic and nuclear sequences on a representative sampling of species and subspecies of the Beta section (except B. patula). Our objectives were (1) to assess their genetic relationships through phylogenetic and multivariate analyses, (2) relate their genetic diversity to their mating system, and (3) reconsider the ploidy status and the origin of the Canarian Beta macrocarpa.

Prospects on the evolutionary mitogenomics of plants: A case study on the olive family (Oleaceae)

The mitogenome is rarely used to reconstruct the evolutionary history of plants, contrary to nuclear and plastid markers. Here, we evaluate the usefulness of mitochondrial DNA for molecular evolutionary studies in Oleaceae, in which cases of cytoplasmic male sterility (CMS) and of potentially contrasted organelle inheritance are known. We compare the diversity and the evolution of mitochondrial and chloroplast genomes by focusing on the olive complex and related genera. Using high-throughput techniques, we reconstructed complete mitogenomes (ca. 0.7 Mb) and plastomes (ca. 156 kb) for six olive accessions and one Chionanthus. A highly variable organization of mitogenomes was observed at the species level. In olive, two specific chimeric genes were identified in the mitogenome of lineage E3 and may be involved in CMS. Plastid-derived regions (mtpt) were observed in all reconstructed mitogenomes. Through phylogenetic reconstruction, we demonstrate that multiple integrations of mtpt regions have occurred in Oleaceae, but mtpt regions shared by all members of the olive complex derive from a common ancestor. We then assembled 52 conserved mitochondrial gene regions and complete plastomes of ten additional accessions belonging to tribes Oleeae, Fontanesieae and Forsythieae. Phylogenetic congruence between topologies based on mitochondrial regions and plastomes suggests a strong disequilibrium linkage between both organellar genomes. Finally, while phylogenetic reconstruction based on plastomes fails to resolve the evolutionary history of maternal olive lineages in the Mediterranean area, their phylogenetic relationships were successfully resolved with complete mitogenomes. Overall, our study demonstrates the great potential of using mitochondrial DNA in plant phylogeographic and metagenomic studies.

Bidirectional natural hybridization between sympatric Ligularia vellerea and L. subspicata

Natural hybridization has been regarded as a crucial pathway of speciation and provides the raw materials for the evolution of biodiversity. The interspecific natural hybridization of the genus Ligularia Cass. is universal and has been considered to be an important factor driving the high diversity of Ligularia species in the Hengduan Mountains, China. Although the natural hybridization between L. vellerea and L. subspicata was reported previously, the direction of hybridization was uncertain due to the limitation of sampling. Thus, in this study, we sampled more individuals and increased two fragments of chloroplast DNA on the basis of the previous study to further verify the natural hybridization between L. vellerea and L. subspicata and confirm the direction of hybridization. Based on DNA sequences (atpB-rbcL, trnL-rpl32, trnQ-5'rps16, and nuclear ribosomal internal transcribed spacer region) data, we concluded that putative hybrids were primary products of hybridization between L. vellerea and L. subspicata and the hybridization was bidirectional. Moreover, sympatric L. tongolensis was not apparently involved in the hybridization. Surprisingly, some pure L. subspicata individuals showed the disaccordance between morphology and DNA data, which might indicate that introgression occurs between L. vellerea and L. subspicata.

Understanding invasion history and predicting invasive niches using genetic sequencing technology in Australia: case studies from Cucurbitaceae and Boraginaceae

Part of the challenge in dealing with invasive plant species is that they seldom represent a uniform, static entity. Often, an accurate understanding of the history of plant introduction and knowledge of the real levels of genetic diversity present in species and populations of importance is lacking. Currently, the role of genetic diversity in promoting the successful establishment of invasive plants is not well defined. Genetic profiling of invasive plants should enhance our understanding of the dynamics of colonization in the invaded range. Recent advances in DNA sequencing technology have greatly facilitated the rapid and complete assessment of plant population genetics. Here, we apply our current understanding of the genetics and ecophysiology of plant invasions to recent work on Australian plant invaders from the Cucurbitaceae and Boraginaceae. The Cucurbitaceae study showed that both prickly paddy melon (Cucumis myriocarpus) and camel melon (Citrullus lanatus) were represented by only a single genotype in Australia, implying that each was probably introduced as a single introduction event. In contrast, a third invasive melon, Citrullus colocynthis, possessed a moderate level of genetic diversity in Australia and was potentially introduced to the continent at least twice. The Boraginaceae study demonstrated the value of comparing two similar congeneric species; one, Echium plantagineum, is highly invasive and genetically diverse, whereas the other, Echium vulgare, exhibits less genetic diversity and occupies a more limited ecological niche. Sequence analysis provided precise identification of invasive plant species, as well as information on genetic diversity and phylogeographic history. Improved sequencing technologies will continue to allow greater resolution of genetic relationships among invasive plant populations, thereby potentially improving our ability to predict the impact of these relationships upon future spread and better manage invaders possessing potentially diverse biotypes and exhibiting diverse breeding systems, life histories and invasion histories.

Association between Chloroplast and Mitochondrial DNA sequences in Chinese Prunus genotypes (Prunus persica, Prunus domestica, and Prunus avium)

BACKGROUND

The nuclear DNA is conventionally used to assess the diversity and relatedness among different species, but variations at the DNA genome level has also been used to study the relationship among different organisms. In most species, mitochondrial and chloroplast genomes are inherited maternally; therefore it is anticipated that organelle DNA remains completely associated. Many research studies were conducted simultaneously on organelle genome. The objectives of this study was to analyze the genetic relationship between chloroplast and mitochondrial DNA in three Chinese Prunus genotypes viz., Prunus persica, Prunus domestica, and Prunus avium.

RESULTS

We investigated the genetic diversity of Prunus genotypes using simple sequence repeat (SSR) markers relevant to the chloroplast and mitochondria. Most of the genotypes were genetically similar as revealed by phylogenetic analysis. The Y2 Wu Xing (Cherry) and L2 Hong Xin Li (Plum) genotypes have a high similarity index (0.89), followed by Zi Ye Li (0.85), whereas; L1 Tai Yang Li (plum) has the lowest genetic similarity (0.35). In case of cpSSR, Hong Tao (Peach) and L1 Tai Yang Li (Plum) genotypes demonstrated similarity index of 0.85 and Huang Tao has the lowest similarity index of 0.50. The mtSSR nucleotide sequence analysis revealed that each genotype has similar amplicon length (509 bp) except M5Y1 i.e., 505 bp with CCB256 primer; while in case of NAD6 primer, all genotypes showed different sizes. The MEHO (Peach), MEY1 (Cherry), MEL2 (Plum) and MEL1 (Plum) have 586 bps; while MEY2 (Cherry), MEZI (Plum) and MEHU (Peach) have 585, 584 and 566 bp, respectively. The CCB256 primer showed highly conserved sequences and minute single polymorphic nucleotides with no deletion or mutation. The cpSSR (ARCP511) microsatellites showed the harmonious amplicon length. The CZI (Plum), CHO (Peach) and CL1 (Plum) showed 182 bp; whileCHU (Peach), CY2 (Cherry), CL2 (Plum) and CY1 (Cherry) showed 181 bp amplicon lengths.

CONCLUSIONS

These results demonstrated high conservation in chloroplast and mitochondrial genome among Prunus species during the evolutionary process. These findings are valuable to study the organelle DNA diversity in different species and genotypes of Prunus to provide in depth insight in to the mitochondrial and chloroplast genomes.

Genetic diversity and linkage disequilibrium analysis in elite sugar beet breeding lines and wild beet accessions

Linkage disequilibrium decay in sugar beet is strongly affected by the breeding history, and varies extensively between and along chromosomes, allowing identification of known and unknown signatures of selection. Genetic diversity and linkage disequilibrium (LD) patterns were investigated in 233 elite sugar beet breeding lines and 91 wild beet accessions, using 454 single nucleotide polymorphisms (SNPs) and 418 SNPs, respectively. Principal coordinate analysis suggested the existence of three groups of germplasm, corresponding to the wild beets, the seed parent and the pollen parent breeding pool. LD was investigated in each of these groups, with and without correction for genetic relatedness. Without correction for genetic relatedness, in the pollen as well as the seed parent pool, LD persisted beyond 50 centiMorgan (cM) on four (2, 3, 4 and 5) and three chromosomes (2, 4 and 6), respectively; after correction for genetic relatedness, LD decayed after <6 cM on all chromosomes in both pools. In the wild beet accessions, there was a strong LD decay: on average LD disappeared after 1 cM when LD was calculated with a correction for genetic relatedness. Persistence of LD was not only observed between distant SNPs on the same chromosome, but also between SNPs on different chromosomes. Regions on chromosomes 3 and 4 that harbor disease resistance and monogermy loci showed strong genetic differentiation between the pollen and seed parent pools. Other regions, on chromosomes 8 and 9, for which no a priori information was available with respect to their contribution to the phenotype, still contributed to clustering of lines in the elite breeding material.

Phylogeographic analysis of a temperate-deciduous forest restricted plant (Bupleurum longiradiatum Turcz.) reveals two refuge areas in China with subsequent refugial isolation promoting speciation

This study investigates the influence of climate-induced oscillations and complicated geological conditions on the evolutionary processes responsible for species formation in presently fragmented temperate forest habitats, located in continental East Asia. In addition to this, we also investigate the heavily debated issue of whether temperate forests migrated southwards during such glacial periods or, alternatively, whether there existed refugia within north China, enabling localized survival of temperate forests within this region. In order to achieve these, we surveyed the phylogeography of Bupleurum longiradiatum Turcz. (a herbaceous plant solely confined to temperate forests) constructed from sequence variation in three chloroplast (cp) DNA fragments: trnL-trnF, psbA-trnH and rps16. Our analyses show high genetic diversity within species (h(T)=0.948) and pronounced genetic differentiation among groups (yellow and purple flowers) with a significant phylogeographical pattern (N(ST)>G(ST), P<0.05). Forty-three haplotypes were identified and clustered into two lineages (the purple-flowered lineage and the yellow-flowered lineage). Two corresponding refuge areas, one in Qinling and its adjacent regions and one in the Changbai Mountains/eastern China, were revealed across the entire distribution ranges of Bupleurum longiradiatum. These results provide evidence for the hypothesis that independent refugia were maintained across the range of temperate forests in northern China during the last glacial maximum or earlier cold periods. Bupleurum longiradiatum var. porphyranthum formed a single taxon based on molecular data. This specific formation process suggests that the historical vicariance factors, i.e. climate-induced eco-geographic isolation through the biotic displacement of temperate-deciduous forest habitats, enhanced the divergence of the yellow and purple flower lineages at different spatial-temporal scales and over glacial and interglacial periods. Additionally, geological conditions that restricted gene flow might also be responsible for the observed high genetic and geographic differentiation. A nested clade analysis (NCA) revealed that allopatric fragmentation was a major factor responsible for the phylogeographic pattern observed, and also supported a role for historical vicariance factors. Our results therefore support the inference that Quaternary refugial isolation promoted allopatric speciation of temperate plants in East Asia. This may help to explain the existence of high diversity and endemism of plant species in East Asia.

Disentangling the causes of heterogeneity in male fecundity in gynodioecious Beta vulgaris ssp. maritima

Variation among individuals in reproductive success is advocated as a major process driving evolution of sexual polymorphisms in plants, such as gynodioecy where females and hermaphrodites coexist. In gynodioecious Beta vulgaris ssp. maritima, sex determination involves cytoplasmic male sterility (CMS) genes and nuclear restorers of male fertility. Both restored CMS and non-CMS hermaphrodites co-occur. Genotype-specific differences in male fitness are theoretically expected to explain the maintenance of cytonuclear polymorphism. Using genotypic information on seedlings and flowering plants within two metapopulations, we investigated whether male fecundity was influenced by ecological, phenotypic and genetic factors, while taking into account the shape and scale of pollen dispersal. Along with spatially restricted pollen flow, we showed that male fecundity was affected by flowering synchrony, investment in reproduction, pollen production and cytoplasmic identity of potential fathers. Siring success of non-CMS hermaphrodites was higher than that of restored CMS hermaphrodites. However, the magnitude of the difference in fecundity depended on the likelihood of carrying restorer alleles for non-CMS hermaphrodites. Our results suggest the occurrence of a cost of silent restorers, a condition supported by scarce empirical evidence, but theoretically required to maintain a stable sexual polymorphism in gynodioecious species.

Gynodioecy in structured populations: understanding fine-scale sex ratio variation in Beta vulgaris ssp. maritima

Natural selection, random processes and gene flow are known to generate sex ratio variations among sexually polymorphic plant populations. In gynodioecious species, in which hermaphrodites and females coexist, the relative effect of these processes on the maintenance of sex polymorphism is still up for debate. The aim of this study was to document sex ratio and cytonuclear genetic variation at a very local scale in wind-pollinated gynodioecious Beta vulgaris ssp. maritima and attempt to elucidate which processes explained the observed variation. The study sites were characterized by geographically distinct patches of individuals and appeared to be dynamic entities, with recurrent establishment of distinct haplotypes through independent founder events. Along with substantial variation in sex ratio and unexpectedly low gene flow within study sites, our results showed a high genetic differentiation among a mosaic of genetically distinct demes, with isolation by distance or abrupt genetic discontinuities taking place within a few tens of metres. Overall, random founder events with restricted gene flow could be primary determinants of sex structure, by promoting the clumping of sex-determining genes. Such high levels of sex structure provide a landscape for differential selection acting on sex-determining genes, which could modify the conditions of maintenance of gynodioecy in structured populations.

Sex-specific fitness variation in gynodioecious Beta vulgaris ssp. maritima: do empirical observations fit theoretical predictions?

In gynodioecious species, in which hermaphroditic and female plants co-occur, the maintenance of sexual polymorphism relies on the genetic determination of sex and on the relative fitness of the different phenotypes. Flower production, components of male fitness (pollen quantity and pollen quality) and female fitness (fruit and seed set) were measured in gynodioecious Beta vulgaris spp. maritima, in which sex is determined by interactions between cytoplasmic male sterility (CMS) genes and nuclear restorers of male fertility. The results suggested that (i) female had a marginal advantage over hermaphrodites in terms of flower production only, (ii) restored CMS hermaphrodites (carrying both CMS genes and nuclear restorers) suffered a slight decrease in fruit production compared to non-CMS hermaphrodites and (iii) restored CMS hermaphrodites were poor pollen producers compared to non-CMS hermaphrodites, probably as a consequence of complex determination of restoration. These observations potentially have important consequences for the conditions of maintenance of sexual polymorphism in B. vulgaris and are discussed in the light of existing theory on evolutionary dynamics of gynodioecy.

Structural and content diversity of mitochondrial genome in beet: a comparative genomic analysis

Despite their monophyletic origin, mitochondrial (mt) genomes of plants and animals have developed contrasted evolutionary paths over time. Animal mt genomes are generally small, compact, and exhibit high mutation rates, whereas plant mt genomes exhibit low mutation rates, little compactness, larger sizes, and highly rearranged structures. We present the (nearly) whole sequences of five new mt genomes in the Beta genus: four from Beta vulgaris and one from B. macrocarpa, a sister species belonging to the same Beta section. We pooled our results with two previously sequenced genomes of B. vulgaris and studied genome diversity at the species level with an emphasis on cytoplasmic male-sterilizing (CMS) genomes. We showed that, contrary to what was previously assumed, all three CMS genomes belong to a single sterile lineage. In addition, the CMSs seem to have undergone an acceleration of the rates of substitution and rearrangement. This study suggests that male sterility emergence might have been favored by faster rates of evolution, unless CMS itself caused faster evolution.

Cytoplasmic phylogeny and evidence of cyto-nuclear co-adaptation in Arabidopsis thaliana

In recent years Arabidopsis thaliana has become a model species for genomic variability and adaptation studies. Although impressive quantities of data have been gathered on the nuclear genomic diversity of this species, little has been published regarding its cytoplasmic diversity. We analyzed the diversity of plastid (pt) and mitochondrial (mt) genomes among 95 accessions, covering most Arabidopsis geographic origins. Four intergenic regions of the pt genome were sequenced, and a total of 68 polymorphisms and 65 pt haplotypes were identified. Several strategies were developed to identify mt polymorphisms among a subset of 14 accessions. Fifteen polymorphisms were further developed as PCR-based markers and used to analyze the whole set of 95 accessions. Using statistical parsimony, we built pt and mt phylogenetic networks of haplotype groups. To root the pt network, the pt intergenic regions of two related Arabidopsis species, Arabidopsis lyrata and Arabidopsis arenosa, were also sequenced. The mt and pt phylogenies are highly congruent and could be combined into a single cytoplasmic phylogeny. To estimate whether co-adaptation between nuclear and cytoplasmic genomes exists in A. thaliana, we tested the germination capacity in challenging conditions of 27 pairs of reciprocal F(2) families. We found that the cytoplasm donor had a significant effect on the germination capacity of some F(2) families.

Effects of fine-scale genetic structure on male mating success in gynodioecious Beta vulgaris ssp. maritima

Plant mating systems are known to influence population genetic structure because pollen and seed dispersal are often spatially restricted. However, the reciprocal outcomes of population structure on the dynamics of polymorphic mating systems have received little attention. In gynodioecious sea beet (Beta vulgaris ssp. maritima), three sexual types co-occur: females carrying a cytoplasmic male sterility (CMS) gene, hermaphrodites carrying a non-CMS cytoplasm and restored hermaphrodites that carry CMS genes and nuclear restorer alleles. This study investigated the effects of fine-scale genetic structure on male reproductive success of the two hermaphroditic forms. Our study population was strongly structured and characterized by contrasting local sex-ratios. Pollen flow was constrained over short distances and depended on local plant density. Interestingly, restored hermaphrodites sired significantly more seedlings than non-CMS hermaphrodites, despite the previous observation that the former produce pollen of lower quality than the latter. This result was explained by the higher frequency of females in the local vicinity of restored (CMS) hermaphrodites as compared to non-CMS hermaphrodites. Population structure thus strongly influences individual fitness and may locally counteract the expected effects of selection, suggesting that understanding fine scale population processes is central to predicting the evolution of gender polymorphism in angiosperms.

Phylogenetic analysis of mitochondrial substitution rate variation in the angiosperm tribe Sileneae

BACKGROUND

Recent phylogenetic studies have revealed that the mitochondrial genome of the angiosperm Silene noctiflora (Caryophyllaceae) has experienced a massive mutation-driven acceleration in substitution rate, placing it among the fastest evolving eukaryotic genomes ever identified. To date, it appears that other species within Silene have maintained more typical substitution rates, suggesting that the acceleration in S. noctiflora is a recent and isolated evolutionary event. This assessment, however, is based on a very limited sampling of taxa within this diverse genus.

RESULTS

We analyzed the substitution rates in 4 mitochondrial genes (atp1, atp9, cox3 and nad9) across a broad sample of 74 species within Silene and related genera in the tribe Sileneae. We found that S. noctiflora shares its history of elevated mitochondrial substitution rate with the closely related species S. turkestanica. Another section of the genus (Conoimorpha) has experienced an acceleration of comparable magnitude. The phylogenetic data remain ambiguous as to whether the accelerations in these two clades represent independent evolutionary events or a single ancestral change. Rate variation among genes was equally dramatic. Most of the genus exhibited elevated rates for atp9 such that the average tree-wide substitution rate for this gene approached the values for the fastest evolving branches in the other three genes. In addition, some species exhibited major accelerations in atp1 and/or cox3 with no correlated change in other genes. Rates of non-synonymous substitution did not increase proportionally with synonymous rates but instead remained low and relatively invariant.

CONCLUSION

The patterns of phylogenetic divergence within Sileneae suggest enormous variability in plant mitochondrial mutation rates and reveal a complex interaction of gene and species effects. The variation in rates across genomic and phylogenetic scales raises questions about the mechanisms responsible for the evolution of mutation rates in plant mitochondrial genomes.

Nuclear and cytoplasmic genetic diversity in weed beet and sugar beet accessions compared to wild relatives: new insights into the genetic relationships within the Beta vulgaris complex species

Hybridization between cultivated species and their wild relatives is now widely considered to be common. In the Beta vulgaris complex, the sugar beet seed multiplication areas have been the scene of inadvertent pollination of sugar beet seed bearers by wild ruderal pollen donors, generating a weedy form of beet which infests sugar beet fields in European countries. Up to now, investigations of evolutionary dynamics of genetic diversity within the B. vulgaris complex were addressed using few genetical markers and few accessions. In this study, we tackled this issue using a panel of complementary markers: five nuclear microsatellite loci, four mitochondrial minisatellite loci and one chloroplastic PCR-RFLP marker. We sampled 1,640 individuals that illustrate the actual distribution of inland ruderal beets of South Western France, weed beets and wild sea beets of northern France as well as the diversity of 35 contemporary European diploid cultivars. Nuclear genetic diversity in weed beets appeared to be as high as those of ruderal beets and sea beets, whereas the narrowness of cultivar accessions was confirmed. This genetic bottleneck in cultivars is even more important in the cytoplasmic genome as only one haplotype was found among all sugar beet cultivars. The large majority of weed beet populations also presented this unique cytoplasmic haplotype, as expected owing to their maternal cultivated origin. Nonetheless, various cytoplasmic haplotypes were found within three populations of weed beets, implying wild-to-weed seed flows. Finally, our findings gave new insights into the genetical relationships between the components of the B. vulgaris complex: (1) we found a very strong genetic divergence between wild sea beet and other relatives, which was unexpected given the recent evolutionary history and the full cross-compatibility of all taxa and (2) we definitely confirmed that the classification into cultivated, wild, ruderal and weed forms according to their geographical location, phenotype or their domesticated status is clearly in accordance with genetic clustering despite the very recent domestication process of sugar beet.

Hybrid male sterility in Mimulus (Phrymaceae) is associated with a geographically restricted mitochondrial rearrangement

Cytoplasmic male sterility (CMS) and nuclear fertility restoration (Rf) involves intergenomic coevolution. Although male-sterile phenotypes are rarely expressed in natural populations of angiosperms, CMS genes are thought to be common. The evolutionary dynamics of CMS/Rf systems are poorly understood, leaving gaps in our understanding of mechanisms and consequences of cytonuclear interactions. We characterized the molecular basis and geographic distribution of a CMS gene in Mimulus guttatus. We used outcrossing M. guttatus (with CMS and Rf) to self-fertilizing M. nasutus (lacking CMS and Rf) to generate hybrids segregating for CMS. Mitochondrial transcripts containing an essential gene (nad6) were perfectly associated with male sterility. The CMS mitotype was completely absent in M. nasutus, present in all genotypes collected from the original collection site, but in only two individuals from 34 other M. guttatus populations. This pattern suggests that the CMS likely originated at a single locality, spread to fixation within the population, but has not spread to other populations, indicating possible ecological or genetic constraints on dispersal of this CMS mitotype between populations. Extreme localization may be characteristic of CMS in hermaphroditic species, in contrast to geographically widespread mitotypes commonly found in gynodioecious species, and could directly contribute to hybrid incompatibilities in nature.

Organellar genome analysis of rye (Secale cereale) representing diverse geographic regions

Rye (Secale cereale) is an important diploid (2n = 14, RR) crop species of the Triticeae and a better understanding of its organellar genome variation can aid in its improvement. Previous genetic analyses of rye focused on the nuclear genome. In the present study, the objective was to investigate the organellar genome diversity and relationships of 96 accessions representing diverse geographic regions using chloroplast (cp) and mitochondrial (mt) DNA PCR-RFLPs. Seven cpDNA and 4 mtDNA coding and noncoding regions were amplified using universal cpDNA and mtDNA primer pairs. Each amplified fragment was digested with 13 different restriction enzymes. mtDNA analysis indicated that the number of polymorphic loci (20) was low and genetic differentiation (GST) was 0.60, excluding the outgroups (hexaploid wheat, Triticum aestivum, 2n = 6x = 42, AABBDD; triticale, xTriticosecale Wittmack, 2n = 6x = 42, AABBRR). cpDNA analysis revealed a low level of polymorphism (40%) among the accessions, and GST was 0.39. Of the 96 genotypes studied, 70 could not be differentiated using cpDNA PCR-RFLPs even though they are from different geographic regions. This is most likely due to germplasm exchange, indicating that genotypes might have a common genetic background. Two cpDNA and 3 mtDNA fragments were significantly correlated to the site of germplasm collection. However, there was no clear trend. These results indicate that the level of organellar polymorphism is low among the cultivated rye genotypes. The cpDNA and mtDNA PCR-RFLP markers used in the present study could be used as molecular markers in rye genetics and breeding programs.

Inheritance of chloroplast DNA is not strictly maternal in Silene vulgaris (Caryophyllaceae): evidence from experimental crosses and natural populations

Chloroplast DNA (cpDNA) is maternally inherited in the majority, but not all, of angiosperm species. The mode of inheritance of cpDNA is a critical determinant of its molecular evolution and of its population genetic structure. Here, we present the results of investigations of the inheritance of cpDNA in Silene vulgaris, a plant used in a variety of studies in which cpDNA is an important component. PCR/RFLP markers were used to compare mother and offspring cpDNA genotypes sampled from two natural populations, and mother, father, and offspring genotypes obtained from controlled greenhouse crosses. Ten of 215 offspring cpDNA genotypes studied in the controlled crosses and three of 156 offspring from natural populations did not match that of the mother, demonstrating rare nonmaternal inheritance. That the chloroplast genome is occasionally transmitted through pollen is discussed in the context of using S. vulgaris cpDNA as a marker in studies of seed dispersal and when considering the joint evolution of the chloroplast and mitochondrial genomes.

Nonneutral evolution of organelle genes in Silene vulgaris

Knowledge of mitochondrial gene evolution in angiosperms has taken a dramatic shift within the past decade, from universal slow rates of nucleotide change to a growing realization of high variation in rates among lineages. Additionally, evidence of paternal inheritance of plant mitochondria and recombination among mitochondrial genomes within heteroplasmic individuals has led to speculation about the potential for independent evolution of organellar genes. We report intraspecific mitochondrial and chloroplast sequence variation in a cosmopolitan sample of 42 Silene vulgaris individuals. There was remarkably high variation in two mitochondrial genes (atp1, atp9) and additional variation within a third gene (cob). Tests for patterns of nonneutral evolution were significant for atp1 and atp9, indicative of the maintenance of balanced polymorphisms. Two chloroplast genes (matK, ndhF) possessed less, but still high, variation and no divergence from neutral expectations. Phylogenetic patterns of organelle genes in both the chloroplast and mitochondria were incongruent, indicating the potential for independent evolutionary trajectories. Evidence indicated reassociation among cytoplasmic genomes and recombination between mitochondrial genes and within atp1, implying transient heteroplasmy in ancestral lineages. Although the mechanisms for long-term maintenance of mitochondrial polymorphism are currently unknown, frequency-dependent selection on linked cytoplasmic male sterility genes is a potential candidate.

Emergence of gynodioecy in wild beet (Beta vulgaris ssp. maritima L.): a genealogical approach using chloroplastic nucleotide sequences

Gynodioecy is a breeding system where both hermaphroditic and female individuals coexist within plant populations. This dimorphism is the result of a genomic interaction between maternally inherited cytoplasmic male sterility (CMS) genes and bi-parentally inherited nuclear male fertility restorers. As opposed to other gynodioecious species, where every cytoplasm seems to be associated with male sterility, wild beet Beta vulgaris ssp. maritima exhibits a minority of sterilizing cytoplasms among numerous non-sterilizing ones. Many studies on population genetics have explored the molecular diversity of different CMS cytoplasms, but questions remain concerning their evolutionary dynamics. In this paper we report one of the first investigations on phylogenetic relationships between CMS and non-CMS lineages. We investigated the phylogenetic relationships between 35 individuals exhibiting different mitochondrial haplotypes. Relying on the high linkage disequilibrium between chloroplastic and mitochondrial genomes, we chose to analyse the nucleotide sequence diversity of three chloroplastic fragments (trnK intron, trnD-trnT and trnL-trnF intergenic spacers). Nucleotide diversity appeared to be low, suggesting a recent bottleneck during the evolutionary history of B. vulgaris ssp. maritima. Statistical parsimony analyses revealed a star-like genealogy and showed that sterilizing haplotypes all belong to different lineages derived from an ancestral non-sterilizing cytoplasm. These results suggest a rapid evolution of male sterility in this taxon. The emergence of gynodioecy in wild beet is confronted with theoretical expectations, describing either gynodioecy dynamics as the maintenance of CMS factors through balancing selection or as a constant turnover of new CMSs.

Linkage disequilibrium and association studies in higher plants: present status and future prospects

During the last two decades, DNA-based molecular markers have been extensively utilized for a variety of studies in both plant and animal systems. One of the major uses of these markers is the construction of genome-wide molecular maps and the genetic analysis of simple and complex traits. However, these studies are generally based on linkage analysis in mapping populations, thus placing serious limitations in using molecular markers for genetic analysis in a variety of plant systems. Therefore, alternative approaches have been suggested, and one of these approaches makes use of linkage disequilibrium (LD)-based association analysis. Although this approach of association analysis has already been used for studies on genetics of complex traits (including different diseases) in humans, its use in plants has just started. In the present review, we first define and distinguish between LD and association mapping, and then briefly describe various measures of LD and the two methods of its depiction. We then give a list of different factors that affect LD without discussing them, and also discuss the current issues of LD research in plants. Later, we also describe the various uses of LD in plant genomics research and summarize the present status of LD research in different plant genomes. In the end, we discuss briefly the future prospects of LD research in plants, and give a list of softwares that are useful in LD research, which is available as electronic supplementary material (ESM).

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.

INVITED REVIEW: Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations

Plants offer excellent models to investigate how gene flow shapes the organization of genetic diversity. Their three genomes can have different modes of transmission and will hence experience varying levels of gene flow. We have compiled studies of genetic structure based on chloroplast DNA (cpDNA), mitochondrial DNA (mtDNA) and nuclear markers in seed plants. Based on a data set of 183 species belonging to 103 genera and 52 families, we show that the precision of estimates of genetic differentiation (G(ST)) used to infer gene flow is mostly constrained by the sampling of populations. Mode of inheritance appears to have a major effect on G(ST). Maternally inherited genomes experience considerably more subdivision (median value of 0.67) than paternally or biparentally inherited genomes (approximately 0.10). G(ST) at cpDNA and mtDNA markers covary narrowly when both genomes are maternally inherited, whereas G(ST) at paternally and biparentally inherited markers also covary positively but more loosely and G(ST) at maternally inherited markers are largely independent of values based on nuclear markers. A model-based gross estimate suggests that, at the rangewide scale, historical levels of pollen flow are generally at least an order of magnitude larger than levels of seed flow (median of the pollen-to-seed migration ratio: 17) and that pollen and seed gene flow vary independently across species. Finally, we show that measures of subdivision that take into account the degree of similarity between haplotypes (N(ST) or R(ST)) make better use of the information inherent in haplotype data than standard measures based on allele frequencies only.

Tracing back seed and pollen flow within the crop-wild Beta vulgaris complex: genetic distinctiveness vs. hot spots of hybridization over a regional scale

Hybrids between transgenic crops and wild relatives have been documented successfully in a wide range of cultivated species, having implications on conservation and biosafety management. Nonetheless, the magnitude and frequency of hybridization in the wild is still an open question, in particular when considering several populations at the landscape level. The Beta vulgaris complex provides an excellent biological model to tackle this issue. Weed beets contaminating sugar beet fields are expected to act as a relay between wild populations and crops and from crops-to-crops. In one major European sugar beet production area, nine wild populations and 12 weed populations were genetically characterized using cytoplasmic markers specific to the cultivated lines and nuclear microsatellite loci. A tremendous overall genetic differentiation between neighbouring wild and weed populations was depicted. However, genetic admixture analyses at the individual level revealed clear evidence for gene flow between wild and weed populations. In particular, one wild population displayed a high magnitude of nuclear genetic admixture, reinforced by direct seed flow as evidenced by cytoplasmic markers. Altogether, weed beets were shown to act as relay for gene flow between crops to wild populations and crops to crops by pollen and seeds at a landscape level.

Comparative analysis of chloroplast DNA variability in wild and cultivated Citrullus species

PCR amplification and restriction site analysis of chloroplast (cp) DNA regions was used to detect inter- and intraspecific differences in the genus Citrullus. More than 55 C. lanatus and 15 C. colocynthis accessions collected from diverse geographical areas, C. ecirrhosus and C. rehmii were used. Most of the cpDNA variation within Citrullus was the result of large indels and transitions and transversions. Indels at the ndhA, trnS- trnfM and trnC- trnD regions and several substitutions at restriction enzyme sites can be used to separate C. colocynthis from the other Citrullus species. A nucleotide substitution at a restriction enzyme site at the 3' flanking region of ndhF provided a diagnostic haplotype for C. lanatus var. lanatus, the cultivated watermelon. Similarly, a nucleotide substitution at an intergenic spacer region of the trnC- trnD region resulted in a diagnostic haplotype for citron, C. lanatus var. citroides. Several C. lanatus var. citroides accessions showed the var. lanatus haplotype. C. rehmii showed almost the same haplotype as C. lanatus var. citroides with the exception of a unique insertion at a cpSSR site. Since C. ecirrhosus lacks the derived diagnostic nucleotide substitutions of C. lanatus, it is probably the progenitor of the cultivated watermelon. Intraspecific haplotypes detected within C. colocynthis were associated with geographic origin.

Cytoplasmic composition in Pinus densata and population establishment of the diploid hybrid pine

Sequence and restriction site analyses of the paternally inherited chloroplast rbcL gene and maternally inherited mitochondrial nad1 fragments from the same set of populations and individuals were used to investigate cytoplasmic composition and population establishment of Pinus densata, a diploid pine that originated through hybridization between P. tabuliformis and P. yunnanensis. Two variable sites and three chlorotypes (TT, TC and GC) were detected on the rbcL gene of the three pines. P. densata harboured the three chlorotypes, two of which (TT, GC) were characteristic of the parental species, respectively. The third chlorotype (TC) was distributed extensively in seven of the 10 P. densata populations analysed, and might represent a mutation type or have been derived from an extinct parent. The distribution of chlorotypes, together with that of mitotypes, indicated that significant founder effect and backcross happened during the population establishment of the hybrid pine. P. tabuliformis and P. yunnanensis had acted as both mother and father donors, i.e. bi-directional gene flow existed between the two parental species in the past. Population differentiation of P. densata is high, as detected from the cytoplasmic genomes: GST = 0.533 for cpDNA and GST = 0.905 for mtDNA. The differences in cytoplasmic composition among the hybrid populations suggest that the local populations have undergone different evolutionary histories.

Evidence for the existence of some dissociation in an otherwise strong linkage disequilibrium between mitochondrial and chloroplastic genomes in Cyclobalanopsis glauca

Variations in mitochondrial DNA in Cyclobalanopsis glauca (Thunb. ex Murray) Oerst. were studied in 140 trees from 32 populations collected from within the tree's natural range. By sequencing two mitochondrial DNA intron fragments (nad4/3-nad4/4r and nad7/2-nad7/3r), we revealed a total of 1788 bp and five polymorphic sites which allowed us to distinguish six mitotypes. The mitochondrial DNA markers provided replicated data to support population phylogeographical scenarios suggested previously using chloroplastic DNA markers. The gene genealogical tree of mitochondrial DNA was partially congruent with the chloroplastic DNA tree owing to the slower mutation rate and different mutational direction. Significant linkage disequilibrium existed between the two organellar genomes. Further paring analyses between fragments synthesized using different primers, accompanied by exclusion of polymorphic sites, showed that the random association could be attributed specifically to one of the polymorphic sites of the petG-trnP fragment of the chloroplastic genome, and the three polymorphic sites of the nad4/3-nad4/4r fragment of the mitochondrial genome. The former was inferred to derive from paternal leakage, and the latter from recurrent mutation. These polymorphic sites were also responsible for uncoupling of the combined gene tree of mitotype and chlorotype. In conclusion, specific fragments found in this study contribute to the incomplete congruence of the two organellar lineages that otherwise associate well phylogeographically.

Inferences from mitochondrial DNA patterns on the domestication history of alfalfa (Medicago sativa)

The evolutionary history of a wild-cultivated complex is the outcome of different factors, among which are the domestication of the cultivated form, its geographical and demographic expansion, as well as gene flow between natural and cultivated populations. To clarify this history for alfalfa (Medicago sativa), we adopted a phylogeographical approach based on a sample of natural and cultivated populations of the M. sativa species complex. This sample was characterized for mtDNA variation through restriction fragment length polymorphisms. Twenty-two mitotypes were identified in the whole data set (155 individuals). In the wild pool from the presumed area of origin (Near East to Central Asia), the diversity was high but the absence of geographical differentiation hinders a more precise location of the centre of domestication. Within the cultivated alfalfa, the geographical structure suggests strongly the existence of at least two independent routes of dissemination of alfalfa from its centre of origin. In particular, original mitotypes detected in some regions suggest that wild populations not sampled in this study, and located outside the presumed area of domestication, may have contributed to the diversity of the cultivated pool. Finally, analysis of the mtDNA variation shows that natural populations collected in Spain constitute an endemic wild pool now introgressed partially by cultivated alfalfa.

A MODEL FOR CHLOROPLAST CAPTURE

Chloroplast capture, the introgression of a chloroplast from one species into another, has been frequently suggested as the explanation for inconsistencies between gene trees based on nuclear and cytoplasmic markers in plants. We use a genetic model to determine the conditions for capture to occur, and we find that they are somewhat more general than those given in earlier verbal arguments. Chloroplast capture can occur if cytoplasm substitution provides an advantage in seed production. This can happen through reallocation to the female function when cytonuclear incompatibilities cause partial male sterility, but also under more general conditions. Capture is promoted by nuclear incompatibilities between the two genomes (or a low heterosis in F1 hybrids) and by partial selfing when hybridization causes a decrease in the selfing rate and inbreeding depression is strong. We discuss empirical predictions that can be used to test this mechanism.

Environmental implications of gene flow from sugar beet to wild beet - current status and future research needs

Gene flow via seed or pollen is a basic biological process in plant evolution. The ecological and genetic consequences of gene flow depend on the amount and direction of gene flow as well as on the fitness of hybrids. The assessment of potential risks of transgenic plants should take into account the fact that conventional crops can often cross with wild plants. The precautionary approach in risk management of genetically modified plants (GMPs) may make it necessary to monitor significant wild and weed populations that might be affected by transgene escape. Gene flow is hard to control in wind-pollinated plants like beet (Beta vulgaris). In addition, wild beet populations potentially can undergo evolutionary changes which might expand their geographical distribution. Unintended products of cultivated beets pollinated by wild beets are weed beets that bolt and flower during their first year of planting. Weed beets cause yield losses and can delay harvest. Wild beets are important plant genetic resources and the preservation of wild beet diversity in Europe has been considered in biosafety research. We present here the methodology and research approaches that can be used for monitoring the geographical distribution and diversity of Beta populations. It has recently been shown that a century of gene flow from Beta vulgaris ssp. vulgaris has not altered the genetic diversity of wild Beta vulgaris L. ssp. maritima (L.) Arcang. in the Italian sugar beet seed production area. Future research should focus on the potential evolution of transgenic wild beet populations in comparison to these baseline data. Two monitoring models are presented describing how endpoints can be measured: (1) "Pre-post" crop commercialization against today's baseline and (2) "Parallel" to crop commercialization against GMP free reference areas/ populations. Model 2 has the advantage of taking ongoing changes in genetic diversity and population dynamics into account. Model 1 is more applicable if gene flow is so strong that most areas/populations contain GMPs. Important traits that may change the ecology of populations are genes that confer tolerance to biotic and abiotic stress. An assessment of environmental effects can realistically only be based on endpoints and consequences of gene introgression, which may include economic values of biodiversity in littoral and other ecosystems containing wild beet. In general, there is still a great need to harmonize worldwide monitoring systems by the development of appropriate methods to evaluate the environmental impact of introgressed transgenes.

Phylogeography of maritime pine inferred with organelle markers having contrasted inheritance

Range-wide variation of maritime pine was studied at maternally inherited and paternally inherited markers (mitochondrial DNA and chloroplast DNA). While chloroplast DNA exhibits the highest diversity, phylogeographic inferences from this marker are blurred by homoplasy and extensive pollen flow. In contrast, the only three mitochondrial haplotypes found provide a clear picture of nonoverlapping areas colonized from different refugia, with no single population having a mixed composition (GST = 1). Comparison of the genetic structure inferred from both organelle genomes allows the investigation of differential seed and pollen dispersal, pointing to pollen, but not seed, dispersal across the Strait of Gibraltar (from Morocco into Iberia). A comparison with already available genetic information, especially that of one of the maritime pine's most threatening insect pests, the bast scale Matsucoccus feytaudi, further completes the picture.

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.

Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure

We describe a method for comparing nuclear and organelle population differentiation (F(ST)) in seed plants to test the hypothesis that pollen and seed gene flow rates are equal. Wright's infinite island model is used, with arbitrary levels of self-fertilization and biparental organelle inheritance. The comparison can also be applied to gene flow in animals. Since effective population sizes are smaller for organelle genomes than for nuclear genomes and organelles are often uniparentally inherited, organelle F(ST) is expected to be higher at equilibrium than nuclear F(ST) even if pollen and seed gene flow rates are equal. To reject the null hypothesis of equal seed and pollen gene flow rates, nuclear and organelle F(ST)'s must differ significantly from their expected values under this hypothesis. Finite island model simulations indicate that infinite island model expectations are not greatly biased by finite numbers of populations (>/=100 subpopulations). The power to distinguish dissimilar rates of pollen and seed gene flow depends on confidence intervals for fixation index estimates, which shrink as more subpopulations and loci are sampled. Using data from the tropical tree Corythophora alta, we rejected the null hypothesis that seed and pollen gene flow rates are equal but cannot reject the alternative hypothesis that pollen gene flow is 200 times greater than seed gene flow.

Mitochondrial DNA diversity, population structure, and gender association in the gynodioecious plant Silene vulgaris

A highly variable mitochondrial DNA (mtDNA) restriction fragment length polymorphism (RFLP) locus is used to assess the population structure of mitochondrial genomes in the gynodioecious plant Silene vulgaris at two spatial scales. Thirteen mtDNA haplotypes were identified within 250 individuals from 18 populations in a 20-km diameter region of western Virginia. The population structure of these mtDNA haplotypes was estimated as thetaST = 0.574 (+/- 0.066 SE) and, surprisingly, genetic differentiation among populations was negatively correlated with geographic distance (Mantel r = -0.246, P < 0.002). Additionally, mtDNA haplotypes were spatially clumped at the scale of meters within one population. Gender in S. vulgaris is determined by an interaction between autosomal male fertility restorers and cytoplasmic male sterility (CMS) factors, and seed fitness is affected by an interaction between gender and population sex ratio; thus, selection acting on gender could influence the distribution of mtDNA RFLP haplotypes. The sex ratio (females:hermaphrodites) varied among mtDNA haplotypes across the entire metapopulation, possibly because the haplotypes were in linkage disequilibrium with different CMS factors. The gender associated with some of the most common haplotypes varied among populations, suggesting that there is also population structure in male fertility restorer genes. In comparison with reports of mtDNA variation from other published studies, we found that S. vulgaris exhibits a large number of mtDNA haplotypes relative to that observed in other species.

Phylogeography of Kandelia candel in East Asiatic mangroves based on nucleotide variation of chloroplast and mitochondrial DNAs

Vivipary with precocious seedlings in mangrove plants was thought to be a hindrance to long-range dispersal. To examine the extent of seedling dispersal across oceans, we investigated the phylogeny and genetic structure among East Asiatic populations of Kandelia candel based on organelle DNAs. In total, three, 28 and seven haplotypes of the chloroplast DNA (cpDNA) atpB-rbcL spacer, cpDNA trnL-trnF spacer, and mitochondrial DNA (mtDNA) internal transcribed spacer (ITS) were identified, respectively, from 202 individuals. Three data sets suggested consistent phylogenies recovering two differentiated lineages corresponding to geographical regions, i.e. northern South-China-Sea + East-China-Sea region and southern South-China-Sea region (Sarawak). Phylogenetically, the Sarawak population was closely related to the Ranong population of western Peninsula Malaysia instead of other South-China-Sea populations, indicating its possible origin from the Indian Ocean Rim. No geographical subdivision was detected within the northern geographical region. An analysis of molecular variance (AMOVA) revealed low levels of genetic differentiation between and within mainland and island populations (phiCT = 0.015, phiSC = 0.037), indicating conspicuous long-distance seedling dispersal across oceans. Significant linkage disequilibrium excluded the possibility of recurrent homoplasious mutations as the major force causing phylogenetic discrepancy between mtDNA and the trnL-trnF spacer within the northern region. Instead, relative ages of alleles contributed to non-random chlorotype-mitotype associations and tree inconsistency. Widespread distribution and random associations (chi2 = 0.822, P = 0.189) of eight hypothetical ancestral cytotypes indicated the panmixis of populations of the northern geographical region as a whole. In contrast, rare and recently evolved alleles were restricted to marginal populations, revealing some preferential directional migration.

Organelle DNA phylogeography of Cycas taitungensis, a relict species in Taiwan

The phylogegraphic pattern of Cycas taitungensis, an endemic species with two remaining populations in Taiwan, was investigated based on genetic variability and phylogeny of the atpB-rbcL noncoding spacer of chloroplast DNA (cpDNA) and the ribosomal DNA (rDNA) internal transcribed spacer (ITS) of mitochondrial DNA (mtDNA). High levels of genetic variation at both organelle loci, due to frequent intramolecular recombination, and low levels of genetic differentiation were detected in the relict gymnosperm. The apportionment of genetic variation within and between populations agreed with a migrant-pool model, which describes a migratory pattern with colonists recruited from a random sample of earlier existing populations. Phylogenies obtained from cpDNA and mtDNA were discordant according to neighbour-joining analyses. In total four chlorotypes (clades I-IV) and five mitotypes (clades A-E) were identified based on minimum spanning networks of each locus. Significant linkage disequilibrium in mitotype-chlorotype associations excluded the possibility of the recurrent homoplasious mutations as the major force causing phylogenetic inconsistency. The most abundant chlorotype I was associated with all mitotypes and the most abundant mitotype C with all chlorotypes; no combinations of rare mitotypes with rare chlorotypes were found. According to nested clade analyses, such nonrandom associations may be ascribed to relative ages among alleles associated with the geological history through which cycads evolved. Nested in networks as interior nodes coupled with wide geographical distribution, the most dominant cytotypes of CI and EI may represent ancestral haplotypes of C. taitungensis with a possible long existence prior to the Pleistocene glacial maximum. In contrast, rare chlorotypes and mitotypes with restricted and patchy distribution may have relatively recent origins. Newly evolved genetic elements of mtDNA, with a low frequency, were likely to be associated with the dominant chlorotype, and vice versa, resulting in the nonrandom mitotype-chlorotype associations. Paraphyly of CI and EI cytotypes, leading to the low level of genetic differentiation between cycad populations, indicated a short period for isolation, which allowed low possibilities of the attainment of coalescence at polymorphic ancestral alleles.

Chloroplast DNA phylogeography of Cunninghamia konishii (Cupressaceae), an endemic conifer of Taiwan

In this study, we investigated the genetic structure and phylogeographic pattern of the genus Cunninghamia, a member of the Cupressaceae restricted to mainland China and Taiwan, based on sequences of the trnD–trnT noncoding spacer of the chloroplast DNA. Maternal inheritance of chloroplasts was determined experimentally. No paternal leakage was detected. Both parsimony and neighbor-joining analyses revealed the polyphyly of Cunninghamia konishii, populations of which were nested in clades of C. lanceolata from mainland China. The nucleotide diversity of chloroplast DNA sequences within C. konishii (0.0118) was higher than that between species (0.0104), which agrees with a previous allozyme investigation. Based on mutational differences between sequences, a minimum spanning network consisting of five clades was constructed. Significant genetic differentiation (ΦST = 0.130, P < 0.001) was detected between the clades based on AMOVA analyses. We infer several possible refugia in the Yunnan, Zhejiang, and Guangdong provinces of south China, all located in the minimum network as interior nodes. We also infer possible migration routes of Cunninghamia populations. The phylogeographic pattern shown in the reconstructed network suggests that the present-day Cunninghamia populations in Taiwan were derived from six different sources in continental Asia via long-distance seed dispersal. A migrant-pool model explains the heterogeneous composition of the organelle DNA in Taiwan's populations and the low differentiation between populations of Taiwan and China (ΦCT = 0.012, P = 0.454). In contrast with the genetic heterogeneity within geographic populations, many local populations have attained coalescence at the trnD–trnT alleles, which has led to significant differentiation at the population level.Key words: AMOVA, coalescence, cpDNA, Cunninghamia konishii, Cunninghamia lanceolata, minimum spanning network, phylogeography.