Spatial analysis of nuclear and cytoplasmic DNA diversity in wild sea beet (Beta vulgaris ssp. maritima) populations: do marine currents shape the genetic structure?

Genetic diversity and genome-wide association study of 13 agronomic traits in 977 Beta vulgaris L. germplasms

BACKGROUND

Sugar beet (Beta vulgaris L.) is an economically essential sugar crop worldwide. Its agronomic traits are highly diverse and phenotypically plastic, influencing taproot yield and quality. The National Beet Medium-term Gene Bank in China maintains more than 1700 beet germplasms with diverse countries of origin. However, it lacks detailed genetic background associated with morphological variability and diversity.

RESULTS

Here, a comprehensive genome-wide association study (GWAS) of 13 agronomic traits was conducted in a panel of 977 sugar beet accessions. Almost all phenotypic traits exhibited wide genetic diversity and high coefficient of variation (CV). A total of 170,750 high-quality single-nucleotide polymorphisms (SNPs) were obtained using the genotyping-by-sequencing (GBS). Neighbour-joining phylogenetic analysis, principal component analysis, population structure and kinship showed no obvious relationships among these genotypes based on subgroups or regional sources. GWAS was carried out using a mixed linear model, and 159 significant associations were detected for these traits. Within the 25 kb linkage disequilibrium decay of the associated markers, NRT1/PTR FAMILY 6.3 (BVRB_5g097760); nudix hydrolase 15 (BVRB_8g182070) and TRANSPORT INHIBITOR RESPONSE 1 (BVRB_8g181550); transcription factor MYB77 (BVRB_2g023500); and ethylene-responsive transcription factor ERF014 (BVRB_1g000090) were predicted to be strongly associated with the taproot traits of root groove depth (RGD); root shape (RS); crown size (CS); and flesh colour (FC), respectively. For the aboveground traits, UDP-glycosyltransferase 79B6 (BVRB_9g223780) and NAC domain-containing protein 7 (BVRB_5g097990); F-box protein At1g10780 (BVRB_6g140760); phosphate transporter PHO1 (BVRB_3g048660); F-box protein CPR1 (BVRB_8g181140); and transcription factor MYB77 (BVRB_2g023500) and alcohol acyltransferase 9 (BVRB_2g023460) might be associated with the hypocotyl colour (HC); plant type (PT); petiole length (PL); cotyledon size (C); and fascicled leaf type (FLT) of sugar beet, respectively. AP-2 complex subunit mu (BVRB_5g106130), trihelix transcription factor ASIL2 (BVRB_2g041790) and late embryogenesis abundant protein 18 (BVRB_5g106150) might be involved in pollen quantity (PQ) variation. The candidate genes extensively participated in hormone response, nitrogen and phosphorus transportation, secondary metabolism, fertilization and embryo maturation.

CONCLUSIONS

The genetic basis of agronomical traits is complicated in heterozygous diploid sugar beet. The putative valuable genes found in this study will help further elucidate the molecular mechanism of each phenotypic trait for beet breeding.

Geographic isolation and long-distance gene flow influence the genetic structure of the blue fan palm Brahea armata (Arecaceae)

The formation of the Baja California Peninsula (BCP) has impacted the microevolutionary dynamics of different species in ways that depend on biological traits such as dispersal capacity. Plants with relatively low levels of vagility have exhibited high genetic divergence between the BCP and Continental mainland. Brahea armata (Arecaceae) is a palm species inhabiting the northern part of the BCP and Sonora; its distribution occurs in isolated oases of vegetation. We aimed to evaluate the influence of the formation of the BCP on the genetic structure of B. armata using nuclear microsatellites and chloroplast markers (cpDNA) to compare patterns of genetic diversity and structure with previous published studies. Because gene flow through seeds is usually more limited compared to pollen flow, we expect to find stronger genetic structure at (cpDNA) than at nuclear markers. Moreover, larger genetic structure might also be explained by the smaller effective population size of cpDNA. We analyzed six microsatellite markers and two cpDNA regions. The main results indicated high levels of genetic differentiation among isolated populations located in the BCP, while low genetic differentiation was found between southern populations of the BCP and Sonora, suggesting long distance gene flow. In contrast, chloroplast markers indicated high levels of genetic structure between BCP and Sonora populations, suggesting asymmetrical gene flow between pollen (measured by nuclear microsatellites) and seed (cpDNA markers). This study provides valuable information on genetic diversity of B. armata that can be relevant for conservation and management; and develops microsatellites markers that can be transferred to other Brahea species.

Effects of contemporary shifts of range margins on patterns of genetic structure and mating system in two coastal plant species

Species' geographical ranges are often restricted due to niche limitation resulting in geographical isolation and reduced population size at range margins. Under the "abundant center" paradigm, static marginal populations are thus expected to show higher genetic differentiation and lower genetic diversity than core populations. Low mate availability may also drive shifts toward higher propensity for selfing in geographically marginal populations. However, these predictions remain to be validated for contemporary range shifts occurring under current environmental change. This study is devoted to bridging this gap and assesses the spatial patterns of genetic structure and mating system across the geographical range of two coastal plant species characterized by contrasting contemporary range dynamics: the receding myrmecochorous Dune pansy (Viola tricolor subsp. curtisii) and the widespread expanding hydrochorous Rock samphire (Crithmum maritimum). Both species exhibited high propensity for selfing, with indications of inbreeding depression acting at early life stages. In Dune pansy, a biogeographical break was observed between core and marginal populations, with trailing-edge populations showing higher levels of genetic differentiation, reduced genetic diversity, and higher levels of selfing estimated through progeny arrays. In contrast, genetic structuring was weak in Rock samphire and no clear spatial trends were observed in genetic diversity nor in mating system, likely the result of efficient long-distance seed dispersal by sea-surface currents. Our study highlights that key species differences in life-history traits related to dispersal and/or mate limitation modify the expectations of genetic diversity loss and mating system shift in contemporary range-expanding populations, as compared with historical core populations.

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.

Insights into the genetic relationships among plants of Beta section Beta using SNP markers

Using a much higher number of SNP markers and larger sample sizes than all the previous studies, we characterized the genetic relationships among wild and cultivated plants of section Beta. We analyzed the genetic variation of Beta section Beta, which includes wild taxa (Beta macrocarpa, B. patula, B. vulgaris subsp. adanensis and B. vulgaris subsp. maritima) and cultivars (fodder beet, sugar beet, garden beet, leaf beet, and swiss chards), using 9724 single nucleotide polymorphism markers. The analyses conducted at the individual level without a priori groups confirmed the strong differentiation of B. macrocarpa and B. vulgaris subsp. adanensis from the other taxa. B. vulgaris subsp. maritima showed a complex genetic structure partly following a geographical pattern, which confounded the differences between this taxon and the cultivated varieties. Cultivated varieties were structured into three main groups: garden beets, fodder and sugar beets, and leaf beets and swiss chards. The genetic structure described here will be helpful to correctly estimate linkage disequilibrium and to test for statistical associations between genetic markers and environmental variables.

Taxonomic, spatial and adaptive genetic variation of Beta section Beta

KEY MESSAGE

The genetic variation of Beta section Beta is structured into four taxonomic and spatial clusters. There are significant associations between molecular markers and environmental variables.

ABSTRACT

We investigated the genetic diversity of Beta section Beta, which includes the wild and cultivated relatives of the sugar beet. The taxa included in the study were: Beta vulgaris subsp. maritima, B. vulgaris subsp. adanensis, B. macrocarpa, B. patula and B. vulgaris subsp. vulgaris (garden beet, leaf beet and swiss chards). We collected 1264 accessions originating from the entire distribution area of these taxa and genotyped them for 4436 DArT markers (DArTs). We showed that the genetic variation of these accessions is structured into four taxonomic and spatial clusters: (1) samples of Beta macrocarpa, (2) samples of Beta vulgaris subsp. adanensis, (3) Mediterranean and Asian samples and (4) Atlantic and Northern European samples. These last two clusters were mainly composed of samples of Beta vulgaris subsp. maritima. We investigated in deeper detail the genetic structure of B. vulgaris subsp. maritima, which constituted the majority (80%) of the wild samples. This subspecies exhibited a clinal genetic variation from South-East to North-West. We detected some markers significantly associated to environmental variables in B. vulgaris subsp. maritima. These associations are interpreted as results of natural selection. The variable most often involved in the associations was annual mean temperature. Therefore, these markers can be useful for the development of frost-tolerant winter beets and drought-tolerant rain-fed beets.

Genetic Diversity and Physiological Performance of Portuguese Wild Beet (Beta vulgaris spp. maritima) from Three Contrasting Habitats

The establishment of stress resilient sugar beets (Beta vulgaris spp. vulgaris) is an important breeding goal since this cash crop is susceptible to drought and salinity. The genetic diversity in cultivated sugar beets is low and the beet wild relatives are useful genetic resources for tolerance traits. Three wild beet populations (Beta vulgaris spp. maritima) from contrasting environments, Vaiamonte (VMT, dry inland hill), Comporta (CMP, marsh) and Oeiras (OEI, coastland), and one commercial sugar beet (Isella variety, SB), are compared. At the genetic level, the use of six microsatellite allowed to detect a total of seventy six alleles. It was observed that CMP population has the highest value concerning the effective number of alleles and of expected heterozygosity. By contrast, sugar beet has the lowest values for all the parameters considered. Loci analysis with STRUCTURE allows defining three genetic clusters, the sea beet (OEI and CMP), the inland ruderal beet (VMT) and the sugar beet (SB). A screening test for progressive drought and salinity effects demonstrated that: all populations were able to recover from severe stress; drought impact was higher than that from salinity; the impact on biomass (total, shoot, root) was population specific. The distinct strategies were also visible at physiological level. We evaluated the physiological responses of the populations under drought and salt stress, namely at initial stress stages, late stress stages, and early stress recovery. Multivariate analysis showed that the physiological performance can be used to discriminate between genotypes, with a strong contribution of leaf temperature and leaf osmotic adjustment. However, the separation achieved and the groups formed are dependent on the stress type, stress intensity and duration. Each of the wild beet populations evaluated is very rich in genetic terms (allelic richness) and exhibited physiological plasticity, i.e., the capacity to physiologically adjust to changing environments. These characteristics emphasize the importance of the wild beet ecotypes for beet improvement programs. Two striking ecotypes are VMT, which is the best to cope with drought and salinity, and CMP which has the highest root to shoot ratio. These genotypes can supply breeding programs with distinct goals.

Isolation by resistance across a complex coral reef seascape

A detailed understanding of the genetic structure of populations and an accurate interpretation of processes driving contemporary patterns of gene flow are fundamental to successful spatial conservation management. The field of seascape genetics seeks to incorporate environmental variables and processes into analyses of population genetic data to improve our understanding of forces driving genetic divergence in the marine environment. Information about barriers to gene flow (such as ocean currents) is used to define a resistance surface to predict the spatial genetic structure of populations and explain deviations from the widely applied isolation-by-distance model. The majority of seascape approaches to date have been applied to linear coastal systems or at large spatial scales (more than 250 km), with very few applied to complex systems at regional spatial scales (less than 100 km). Here, we apply a seascape genetics approach to a peripheral population of the broadcast-spawning coral Acropora spicifera across the Houtman Abrolhos Islands, a high-latitude complex coral reef system off the central coast of Western Australia. We coupled population genetic data from a panel of microsatellite DNA markers with a biophysical dispersal model to test whether oceanographic processes could explain patterns of genetic divergence. We identified significant variation in allele frequencies over distances of less than 10 km, with significant differentiation occurring between adjacent sites but not between the most geographically distant ones. Recruitment probabilities between sites based on simulated larval dispersal were projected into a measure of resistance to connectivity that was significantly correlated with patterns of genetic divergence, demonstrating that patterns of spatial genetic structure are a function of restrictions to gene flow imposed by oceanographic currents. This study advances our understanding of the role of larval dispersal on the fine-scale genetic structure of coral populations across a complex island system and applies a methodological framework that can be tailored to suit a variety of marine organisms with a range of life-history characteristics.

Mitochondrial genome and plant taxonomy

The lability in size, structure, and sequence content of mitochondrial genome (mtDNA) across plant species has sharply limited its use in taxonomic studies. However, due to the new opportunities offered by the availability of complete mtDNA sequence in plant species and the subsequent development of universal primers, the number of mtDNA-based molecular studies has recently increased. Historically, universal primers have enabled to characterize mtDNA polymorphism mainly by the RFLP technique. This methodology has been progressively replaced by Sanger DNA sequencing, which actually provides the full phylogenetic information content of a DNA fragment (single nucleotide, insertion/deletion, and single sequence repeat length polymorphism). This chapter presents a sequencing working protocol to be routinely used in mtDNA-based phylogenetic studies.

A seascape genetic analysis reveals strong biogeographical structuring driven by contrasting processes in the polyploid saltmarsh species Puccinellia maritima and Triglochin maritima

Little is known about the processes shaping population structure in saltmarshes. It is expected that the sea should act as a powerful agent of dispersal. Yet, in contrast, import of external propagules into a saltmarsh is thought to be small. To determine the level of connectivity between saltmarsh ecosystems at a macro-geographical scale, we characterized and compared the population structure of two polyploid saltmarsh species, Puccinellia maritima and Triglochin maritima based on a seascape genetics approach. A discriminant analysis of principal components highlighted a genetic structure for both species arranged according to a regional pattern. Subsequent analysis based on isolation-by-distance and isolation-by-resistance frameworks indicated a strong role of coastal sediment transport processes in delimiting regional structure in P. maritima, while additional overland propagule dispersal was indicated for T. maritima. The identification and comparison of regional genetic structure and likely determining factors presented here allows us to understand the biogeographical units along the UK coast, between which barriers to connectivity occur not only at the species level but at the ecosystem scale. This information is valuable in plant conservation and community ecology and in the management and restoration of saltmarsh ecosystems.

Population genetic structure and demographic history of Atrina pectinata based on mitochondrial DNA and microsatellite markers

The pen shell, Atrina pectinata, is one of the commercial bivalves in East Asia and thought to be recently affected by anthropogenic pressure (habitat destruction and/or fishing pressure). Information on its population genetic structure is crucial for the conservation of A. pectinata. Considering its long pelagic larval duration and iteroparity with high fecundity, the genetic structure for A. pectinata could be expected to be weak at a fine scale. However, the unusual oceanography in the coasts of China and Korea suggests potential for restricted dispersal of pelagic larvae and geographical differentiation. In addition, environmental changes associated with Pleistocene sea level fluctuations on the East China Sea continental shelf may also have strongly influenced historical population demography and genetic diversity of marine organisms. Here, partial sequences of the mitochondrial Cytochrome c oxidase subunit I (COI) gene and seven microsatellite loci were used to estimate population genetic structure and demographic history of seven samples from Northern China coast and one sample from North Korea coast. Despite high levels of genetic diversity within samples, there was no genetic differentiation among samples from Northern China coast and low but significant genetic differentiation between some of the Chinese samples and the North Korean sample. A late Pleistocene population expansion, probably after the Last Glacial Maximum, was also demonstrated for A. pectinata samples. No recent genetic bottleneck was detected in any of the eight samples. We concluded that both historical recolonization (through population range expansion and demographic expansion in the late Pleistocene) and current gene flow (through larval dispersal) were responsible for the weak level of genetic structure detected in A. pectinata.

Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system, influence of marine currents, and footprints of postglacial recolonization routes

Understanding the factors that contribute to population genetic divergence across a species' range is a long-standing goal in evolutionary biology and ecological genetics. We examined the relative importance of historical and ecological features in shaping the present-day spatial patterns of genetic structure in two related plant species, Beta vulgaris subsp. maritima and Beta macrocarpa. Using nuclear and mitochondrial markers, we surveyed 93 populations from Brittany (France) to Morocco – the southern limit of their species' range distribution. Whereas B. macrocarpa showed a genotypic structure and a high level of genetic differentiation indicative of selfing, the population genetic structure of B. vulgaris subsp. maritima was consistent with an outcrossing mating system. We further showed (1) a strong geographic clustering in coastal B. vulgaris subsp. maritima populations that highlighted the influence of marine currents in shaping different lineages and (2) a peculiar genetic structure of inland B. vulgaris subsp. maritima populations that could indicate the admixture of distinct evolutionary lineages and recent expansions associated with anthropogenic disturbances. Spatial patterns of nuclear diversity and differentiation also supported a stepwise recolonization of Europe from Atlantic-Mediterranean refugia after the last glacial period, with leading-edge expansions. However, cytoplasmic diversity was not impacted by postglacial recolonization: stochastic long-distance seed dispersal mediated by major oceanic currents may mitigate the common patterns of reduced cytoplasmic diversity observed for edge populations. Overall, the patterns we documented here challenge the general view of reduced genetic diversity at the edge of a species' range distribution and provide clues for understanding how life-history and major geographic features interact to shape the distribution of genetic diversity.

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.

DNA fingerprinting in botany: past, present, future

Almost three decades ago Alec Jeffreys published his seminal Nature papers on the use of minisatellite probes for DNA fingerprinting of humans (Jeffreys and colleagues Nature 1985, 314:67-73 and Nature 1985, 316:76-79). The new technology was soon adopted for many other organisms including plants, and when Hilde Nybom, Kurt Weising and Alec Jeffreys first met at the very First International Conference on DNA Fingerprinting in Berne, Switzerland, in 1990, everybody was enthusiastic about the novel method that allowed us for the first time to discriminate between humans, animals, plants and fungi on the individual level using DNA markers. A newsletter coined "Fingerprint News" was launched, T-shirts were sold, and the proceedings of the Berne conference filled a first book on "DNA fingerprinting: approaches and applications". Four more conferences were about to follow, one on each continent, and Alec Jeffreys of course was invited to all of them. Since these early days, methodologies have undergone a rapid evolution and diversification. A multitude of techniques have been developed, optimized, and eventually abandoned when novel and more efficient and/or more reliable methods appeared. Despite some overlap between the lifetimes of the different technologies, three phases can be defined that coincide with major technological advances. Whereas the first phase of DNA fingerprinting ("the past") was dominated by restriction fragment analysis in conjunction with Southern blot hybridization, the advent of the PCR in the late 1980s gave way to the development of PCR-based single- or multi-locus profiling techniques in the second phase. Given that many routine applications of plant DNA fingerprinting still rely on PCR-based markers, we here refer to these methods as "DNA fingerprinting in the present", and include numerous examples in the present review. The beginning of the third phase actually dates back to 2005, when several novel, highly parallel DNA sequencing strategies were developed that increased the throughput over current Sanger sequencing technology 1000-fold and more. High-speed DNA sequencing was soon also exploited for DNA fingerprinting in plants, either in terms of facilitated marker development, or directly in the sense of "genotyping-by-sequencing". Whereas these novel approaches are applied at an ever increasing rate also in non-model species, they are still far from routine, and we therefore treat them here as "DNA fingerprinting in the future".

Local differentiation amidst extensive allele sharing in Oryza nivara and O. rufipogon

Genetic variation patterns within and between species may change along geographic gradients and at different spatial scales. This was revealed by microsatellite data at 29 loci obtained from 119 accessions of three Oryza series Sativae species in Asia Pacific: Oryza nivara Sharma and Shastry, O. rufipogon Griff., and O. meridionalis Ng. Genetic similarities between O. nivara and O. rufipogon across their distribution are evident in the clustering and ordination results and in the large proportion of shared alleles between these taxa. However, local-level species separation is recognized by Bayesian clustering and neighbor-joining analyses. At the regional scale, the two species seem more differentiated in South Asia than in Southeast Asia as revealed by F ST analysis. The presence of strong gene flow barriers in smaller spatial units is also suggested in the analysis of molecular variance (AMOVA) results where 64% of the genetic variation is contained among populations (as compared to 26% within populations and 10% among species). Oryza nivara (H E = 0.67) exhibits slightly lower diversity and greater population differentiation than O. rufipogon (H E = 0.70). Bayesian inference identified four, and at a finer structural level eight, genetically distinct population groups that correspond to geographic populations within the three taxa. Oryza meridionalis and the Nepalese O. nivara seemed diverged from all the population groups of the series, whereas the Australasian O. rufipogon appeared distinct from the rest of the species.

Erodium maritimum (Geraniaceae), a species with an uneven and fragmented distribution along the Western Mediterranean and European Atlantic coasts, has a weak genetic structure

Erodium maritimum L. is an annual species presenting heterogeneous, sometimes very small, and distant populations, distributed along a discontinuous coastal strip of the European Atlantic and the central and western Mediterranean basin. The aim of this study is to investigate genetic variation and geographic structure changes across its large distribution. Fourteen populations of E. maritimum were studied using AFLP fingerprints, together with their population sizes, reproductive systems and flower visitors. AFLP markers revealed the genetic structure of the species to be weak. Many individuals from one population clustered together with those of other populations, showing a high degree of genetic admixture. Despite having a self-compatible reproductive system, populations (especially the largest ones) showed high levels of genetic polymorphism, and the majority of genetic variation was contained within populations. The low genetic structure suggests high levels of gene flow, which might be explained through the dispersability of the species' fruits. Finally, recommendations are provided for management strategies to facilitate the conservation of this endangered species.

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.

Polymorphic minisatellites in the mitochondrial DNAs of Oryza and Brassica

Polymorphic analyses of angiosperm mitochondrial DNA are rare in comparison with chloroplast DNA, because few target sequences in angiosperm mitochondrial DNA are known. Minisatellites, a tandem array of repeated sequences with a repeat unit of 10 to ~100 bp, are popular target sequences of animal mitochondria, but Beta vulgaris is the only known angiosperm species for which such an analysis has been conducted. From this lack of information, it was uncertain as to whether polymorphic minisatellites existed in other angiosperm species. Ten plant mitochondrial DNAs were found to contain minisatellite-like repeated sequences, most of which were located in intergenic regions but a few occurred in gene coding and intronic regions. Oryza and Brassica accessions were selected as models for the investigation of minisatellite polymorphism because substantial systematic information existed. PCR analysis of 42 Oryza accessions revealed length polymorphisms in four of the five minisatellites. The mitochondrial haplotypes of the 16 Oryza accessions with chromosomal complement (genome) types of CC, BBCC and CCDD were identical but were clearly distinguished from BB-genome accessions, a result consistent with the notion that the cytoplasmic donor parent of the amphidiploid species might be the CC-genome species. Twenty-nine accessions of six major cultivated species of Brassica were classified into five mitochondrial haplotypes based on two polymorphic minisatellites out of six loci. The haplotypes of Brassica juncea and Brassica carinata accessions were identical to Brassica rapa and Brassica nigra accessions, respectively. The haplotypes of Brassica napus accessions were heterogeneous and unique, results that were consistent with previous studies.

Metapopulation structure and fine-scaled genetic structuring in crop-wild hybrid weed beets

This study explores the microspatial and temporal genetic variation in crop-wild hybrid weed beets that emerged from the seed bank in a cultivated field surveyed over two successive years. We demonstrate the occurrence of demes highly genetically differentiated, kin-structured, characterized by moderate effective population sizes, differing in propensity for selfing, and arising from nonrandom genetic subsets of the seed bank. Only one deme identified in the first survey year significantly contributed to the weed beets that emerged in the second year. Spatial structuring appears to be primarily due to gravity seed dispersal and limited pollen flow among weed beet demes. Within each genetic cluster identified by Bayesian assignments and multivariate analyses, F(IS) estimates and level of biparental inbreeding--revealed by progeny analyses--dropped to non-significant values. This suggests that random mating occurs at the scale of genetically distinct demes over a very short scale. Our results highlight the need to carefully depict genetic discontinuities in weed species, when attempting to describe their local genetic neighborhoods within which genetic drift and selective processes occur.

Landscape genetics of plants

Landscape genetics is the amalgamation of landscape ecology and population genetics to help with understanding microevolutionary processes such as gene flow and adaptation. In this review, we examine why landscape genetics of plants lags behind that of animals, both in number of studies and consideration of landscape elements. The classical landscape distance/resistance approach to study gene flow is challenging in plants, whereas boundary detection and the assessment of contemporary gene flow are more feasible. By contrast, the new field of landscape genetics of adaptive genetic variation, establishing the relationship between adaptive genomic regions and environmental factors in natural populations, is prominent in plant studies. Landscape genetics is ideally suited to study processes such as migration and adaptation under global change.

Pollen limitation of female reproductive success at fine spatial scale in a gynodioecious and wind-pollinated species, Beta vulgaris ssp. maritima

In sexually polymorphic plants, the spatial distribution of sexes is usually not random. Local variation in phenotype frequencies is expected to affect individual fitness of the different phenotypes. For gynodioecious species, with co-occurrence of hermaphrodites and females, if sexual phenotypes are structured in space and pollen flow is spatially restricted, local pollen availability should vary among patches. Female fitness may thus be low when hermaphrodites are locally rare. To test this hypothesis, we analysed how the reproductive output of females varied among patches within two natural study sites of the gynodioecious wind-pollinated Beta vulgaris ssp. maritima. Plants growing in female-biased areas and experiencing pollen limitation were found to have low fruit and seed sets but did not reallocate resources towards better offspring. Our results show that fine-scale processes influence individual fitness and the evolution of sex ratio in sexually polymorphic plants.

Taking the chaos out of genetic patchiness: seascape genetics reveals ecological and oceanographic drivers of genetic patterns in three temperate reef species

Marine species frequently show weak and/or complex genetic structuring that is commonly dismissed as 'chaotic' genetic patchiness and ecologically uninformative. Here, using three datasets that individually feature weak chaotic patchiness, we demonstrate that combining inferences across species and incorporating environmental data can greatly improve the predictive value of marine population genetics studies on small spatial scales. Significant correlations in genetic patterns of microsatellite markers among three species, kelp bass Paralabrax clathratus, Kellet's whelk Kelletia kelletii and California spiny lobster Panulirus interruptus, in the Southern California Bight suggest that slight differences in diversity and pairwise differentiation across sampling sites are not simply noise or chaotic patchiness, but are ecologically meaningful. To test whether interspecies correlations potentially result from shared environmental drivers of genetic patterns, we assembled data on kelp bed size, sea surface temperature and estimates of site-to-site migration probability derived from a high resolution multi-year ocean circulation model. These data served as predictor variables in linear models of genetic diversity and linear mixed models of genetic differentiation that were assessed with information-theoretic model selection. Kelp was the most informative predictor of genetics for all three species, but ocean circulation also played a minor role for kelp bass. The shared patterns suggest a single spatial marine management strategy may effectively protect genetic diversity of multiple species. This study demonstrates the power of environmental and ecological data to shed light on weak genetic patterns and highlights the need for future focus on a mechanistic understanding of the links between oceanography, ecology and genetic structure.

Potential for evolutionary change in the seasonal timing of germination in sea beet (Beta vulgaris ssp. maritima) mediated by seed dormancy

In sea beet (Beta vulgaris ssp. maritima), germination occurs in autumn or spring and is mediated by dormancy which can be released by cold or dry periods. Environmental change such as current climate change may require evolutionary response in seasonal timing. Here, we explore the potential for such evolutionary change. Seed dormancy was studied in a composite population based on seeds from all over the species range in France together with several generations of reciprocal crosses. We found high, repeatable variability for dormancy rate among individuals under greenhouse conditions and confirmed its relevance for germination phenology in the field. Our data fitted best with an exclusively maternal determination of the dormancy phenotype. Narrow-sense heritability, h(2) approximately 0.5 in the composite population and approximately 0.4 in the original local populations, was such that rapid evolutionary change in the relative proportions of autumn and spring germination may be possible.

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.

Populations of weedy crop-wild hybrid beets show contrasting variation in mating system and population genetic structure

Reproductive traits are key parameters for the evolution of invasiveness in weedy crop–wild hybrids. In Beta vulgaris, cultivated beets hybridize with their wild relatives in the seed production areas, giving rise to crop–wild hybrid weed beets. We investigated the genetic structure, the variation in first-year flowering and the variation in mating system among weed beet populations occurring within sugar beet production fields. No spatial genetic structure was found for first-year populations composed of F1 crop–wild hybrid beets. In contrast, populations composed of backcrossed weed beets emerging from the seed bank showed a strong isolation-by-distance pattern. Whereas gametophytic self-incompatibility prevents selfing in wild beet populations, all studied weed beet populations had a mixed-mating system, plausibly because of the introgression of the crop-derived Sf gene that disrupts self-incompatibility. No significant relationship between outcrossing rate and local weed beet density was found, suggesting no trends for a shift in the mating system because of environmental effects. We further reveal that increased invasiveness of weed beets may stem from positive selection on first-year flowering induction depending on the B gene inherited from the wild. Finally, we discuss the practical and applied consequences of our findings for crop-weed management.

Fine-scale geographical structure of genetic diversity in inland wild beet populations

Introgression arising from crop-to-wild gene flow provides novel sources of genetic variation in plant species complexes. Hybridization within the Beta vulgaris species complex is of immediate concern; crop lineages (B. vulgaris ssp. vulgaris) hybridize easily with their wild relatives (B. vulgaris ssp. maritima) thereby threatening wild beet gene diversity with genetic swamping. Hybridization 'hotspots' occur in European seed production areas because inland ruderal wild beets occur and reproduce in sympatry with cultivated beets. We studied gene flow occurring between seed-producing cultivars and ruderal wild B. vulgaris in southwestern France to determine whether feral beets, arising from unharvested cultivated seed, represent an opportunity for crop-to-wild gene flow. We surveyed 42 inland ruderal beet populations located near seed production fields for nucleo-cytoplasmic variation and used a cytoplasmic marker diagnostic of cultivated lines. Occurrence of cultivated-type cytoplasm within ruderal populations clearly reflected events of crop seed escape. However, we found no genetic signatures of nuclear cultivated gene introgression, which suggests past introgression of cultivated cytoplasm into a wild nuclear background through seed escape rather than recent direct pollen flow. Overall, patterns of genetic structure suggested that inland ruderal wild beet populations act as a metapopulation, with founding events involving a few sib groups, followed by low rates of seed or pollen gene flow after populations are established. Altogether, our results indicate that a long-lived seed bank plays a key role in maintaining cultivated-type cytoplasm in the wild and highlight the need for careful management of seed production areas where wild and cultivated relatives co-occur.

Fine-scale population genetic structure of Zhikong scallop (Chlamys farreri): do local marine currents drive geographical differentiation?

Marine scallops, with extended planktonic larval stages which can potentially disperse over large distances when advected by marine currents, are expected to possess low geographical differentiation. However, the sessile lifestyle as adult tends to form discrete "sea beds" with unique population dynamics and structure. The narrow distribution of Zhikong scallop (Chlamys farreri), its long planktonic larval stage, and the extremely hydrographic complexity in its distribution range provide an interesting case to elucidate the impact of marine currents on geographical differentiation for marine bivalves at a fine geographical scale. In this study, we analyzed genetic variation at nine microsatellite DNA loci in six locations throughout the distribution of Zhikong scallop in the Northern China. Very high genetic diversity was present in all six populations. Two populations sampled from the same marine gyre had no detectable genetic differentiation (F (ST) = 0.0013); however, the remaining four populations collected from different marine gyres or separated by strong marine currents showed low but significant genetic differentiation (F (ST) range 0.0184-0.0602). Genetic differentiation was further analyzed using the Monmonier algorithm to identify genetic barriers and using the assignment test conducted by software GeneClass2 to ascertain population membership of individuals. The genetic barriers fitting the orientation of marine gyres/currents were clearly identified, and the individual assignment analysis indicated that 95.6% of specimens were correctly allocated to one of the six populations sampled. The results support the hypothesis that significant population structure is present in Zhikong scallop at a fine geographical scale, and marine currents can be responsible for the genetic differentiation.

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.

Long distance pollen-mediated gene flow at a landscape level: the weed beet as a case study

Gene flow is a crucial parameter that can affect the organization of genetic diversity in plant species. It has important implications in terms of conservation of genetic resources and of gene exchanges between crop to wild relatives and within crop species complex. In the Beta vulgaris complex, hybridization between crop and wild beets in seed production areas is well documented and the role of the ensuing hybrids, weed beets, as bridges towards wild forms in sugar beet production areas have been shown. Indeed, in contrast to cultivated beets that are bi-annual, weed beets can bolt, flower and reproduce in the same crop season. Nonetheless, the extent of pollen gene dispersal through weedy lineages remains unknown. In this study, the focus is directed towards weed-to-weed gene flow, and we report the results of a pollen-dispersal analysis within an agricultural landscape composed of five sugar beet fields with different levels of infestation by weed beets. Our results, based on paternity analysis of 3240 progenies from 135 maternal plants using 10 microsatellite loci, clearly demonstrate that even if weedy plants are mostly pollinated by individuals from the same field, some mating events occur between weed beets situated several kilometres apart (up to 9.6 km), with rates of interfield-detected paternities ranging from 11.3% to 17.5%. Moreover, we show that pollen flow appears to be more restricted when individuals are aggregated as most mating events occurred only for short-distance classes. The best-fit dispersal curves were fat-tailed geometric functions for populations exhibiting low densities of weed beets and thin-tailed Weibull function for fields with weed beet high densities. Thus, weed beet populations characterized by low density with geographically isolated individuals may be difficult to detect but are likely to act as pollen traps for pollen emitted by close and remote fields. Hence, it appears evident that interfield pollen-mediated gene flow between weed beets is almost unavoidable and could contribute to the diffusion of (trans)genes in the agricultural landscape.

Mitochondrial DNA phylogeny of cultivated and wild beets: relationships among cytoplasmic male-sterility-inducing and nonsterilizing cytoplasms

Cytoplasmic male sterility (CMS), the maternally inherited failure to produce functional pollen, has been used in the breeding of sugar beet (Beta vulgaris ssp. vulgaris). At least three different sources of CMS can be distinguished from one another as well as from normal fertile cytoplasm by polymorphisms in their mitochondrial genomes. Here we analyzed 50 accessions of cultivated and wild beets to investigate the phylogenetic relationships among male-sterility-inducing and normal cytoplasms. The haplotypes were characterized by the nucleotide sequence of the mitochondrial cox2-cox1 spacer region and mitochondrial minisatellite loci. The results indicated that (1) a normal cytoplasm line, cv. TK81-O, was situated at the major core node of the haplotype network, and (2) the three sterilizing cytoplasms in question derived independently from the core haplotype. The evolutionary pathway was investigated by physical mapping study of the mitochondrial genome of a wild beet (B. vulgaris ssp. orientalis) accession BGRC56777 which shared the same mitochondrial haplotype with TK81-O, but was not identical to TK81-O for the RFLP profiles of mitochondrial DNA. Interestingly, three sets of inverted repeated sequences appeared to have been involved in a series of recombination events during the course of evolution between the BGRC56777 and the TK81-O mitochondrial genomes.