Evolution in subdivided plant populations: concepts, recent advances and future directions.

SUMMARY: Research into the evolution of subdivided plant populations has long involved the study of phenotypic variation across plant geographic ranges and the genetic details underlying that variation. Genetic polymorphism at different marker loci has also allowed us to infer the long- and short-term histories of gene flow within and among populations, including range expansions and colonization-extinction dynamics. However, the advent of affordable genome-wide sequences for large numbers of individuals is opening up new possibilities for the study of subdivided populations. In this review, we consider what the new tools and technologies may allow us to do. In particular, we encourage researchers to look beyond the description of variation and to use genomic tools to address new hypotheses, or old ones afresh. Because subdivided plant populations are complex structures, we caution researchers away from adopting simplistic interpretations of their data, and to consider the patterns they observe in terms of the population genetic processes that have given rise to them; here, the genealogical framework of the coalescent will continue to be conceptually and analytically useful.

Characterizing variation in mycorrhiza effect among diverse plant varieties.

Exploitation of arbuscular mycorrhizal fungi may be an important approach for development of reduced-input agriculture. We discuss the use of linear models to analyze variation in mycorrhiza response among diverse plant varieties in order to assess the value of mycorrhizas. Our approach allows elimination of variation linked to differences in plant performance in the absence of mycorrhizas and the selection of plant lines that might harbor genetic variation of use to improve the mycorrhizal symbiosis in agriculture. We illustrate our approach by applying it to previously published and to novel data. We suggest that in dealing with a relative trait such as mycorrhiza effect, the choice of measure used to quantify the trait greatly affects interpretation. In the plant populations under consideration, we find evidence for a greater potential to increase mycorrhiza benefit than previously suggested.

Comparative population genetic structure in a plant-pollinator/seed predator system.

Comparative analyses of spatial genetic structure of populations of plants and the insects they interact with provide an indication of how gene flow, natural selection and genetic drift may jointly influence the distribution of genetic variation and potential for local co-adaptation for interacting species. Here, we analysed the spatial scale of genetic structure within and among nine populations of an interacting species pair, the white campion Silene latifolia and the moth Hadena bicruris, along a latitudinal gradient across Northern/Central Europe. This dioecious, short-lived perennial plant inhabits patchy, often disturbed environments. The moth H. bicruris acts both as its pollinator and specialist seed predator that reproduces by laying eggs in S. latifolia flowers. We used nine microsatellite markers for S. latifolia and eight newly developed markers for H. bicruris. We found high levels of inbreeding in most populations of both plant and pollinator/seed predator. Among populations, significant genetic structure was observed for S. latifolia but not for its pollinator/seed predator, suggesting that despite migration among populations of H. bicruris, pollen is not, or only rarely, carried over between populations, thus maintaining genetic structure among plant populations. There was a weak positive correlation between genetic distances of S. latifolia and H. bicruris. These results indicate that while significant structure of S. latifolia populations creates the potential for differentiation at traits relevant for the interaction with the pollinator/seed predator, substantial gene flow in H. bicruris may counteract this process in at least some populations.

Variation in the intensity of inbreeding depression among successive life-cycle stages and generations in gynodioecious Silene vulgaris (Caryophyllaceae).

Inbreeding depression is one of the hypotheses explaining the maintenance of females within gynodioecious plant populations. However, the measurement of fitness components in selfed and outcrossed progeny depends on life-cycle stage and the history of inbreeding. Comparative data indicate that strong inbreeding depression is more likely to occur at later life-cycle stages. We used hermaphrodite individuals of Silene vulgaris originating from three populations located in different valleys in the Swiss Alps to investigate the effect of two generations of self- and cross-fertilization on fitness components among successive stages of the life cycle in a glasshouse experiment. We detected significant inbreeding depression for most life-cycle stages including: the number of viable and aborted seeds per fruit, probability of germination, above ground biomass, probability of flowering, number of flowers per plant, flower size and pollen viability. Overall, the intensity of inbreeding depression increased among successive stages of the life cycle and cumulative inbreeding depression was significantly stronger in the first generation (delta approximately 0.5) compared with the second generation (delta approximately 0.35). We found no evidence for synergistic epistasis in our experiment. Our finding of more intense inbreeding depression during later stages of the life cycle may help to explain the maintenance of females in gynodioecious populations of S. vulgaris because purging of genetic load is less likely to occur.

Identifying genetic signatures of selection in a non-model species, alpine gentian (Gentiana nivalis L.), using a landscape genetic approach

It is generally accepted that most plant populations are locally adapted. Yet, understanding how environmental forces give rise to adaptive genetic variation is a challenge in conservation genetics and crucial to the preservation of species under rapidly changing climatic conditions. Environmental variation, phylogeographic history, and population demographic processes all contribute to spatially structured genetic variation, however few current models attempt to separate these confounding effects. To illustrate the benefits of using a spatially-explicit model for identifying potentially adaptive loci, we compared outlier locus detection methods with a recently-developed landscape genetic approach. We analyzed 157 loci from samples of the alpine herb Gentiana nivalis collected across the European Alps. Principle coordinates of neighbor matrices (PCNM), eigenvectors that quantify multi-scale spatial variation present in a data set, were incorporated into a landscape genetic approach relating AFLP frequencies with 23 environmental variables. Four major findings emerged. 1) Fifteen loci were significantly correlated with at least one predictor variable (R (adj) (2) > 0.5). 2) Models including PCNM variables identified eight more potentially adaptive loci than models run without spatial variables. 3) When compared to outlier detection methods, the landscape genetic approach detected four of the same loci plus 11 additional loci. 4) Temperature, precipitation, and solar radiation were the three major environmental factors driving potentially adaptive genetic variation in G. nivalis. Techniques presented in this paper offer an efficient method for identifying potentially adaptive genetic variation and associated environmental forces of selection, providing an important step forward for the conservation of non-model species under global change.

Forecasting changes in population genetic structure of alpine plants in response to global warming.

Species range shifts in response to climate and land use change are commonly forecasted with species distribution models based on species occurrence or abundance data. Although appealing, these models ignore the genetic structure of species, and the fact that different populations might respond in different ways because of adaptation to their environment. Here, we introduced ancestry distribution models, that is, statistical models of the spatial distribution of ancestry proportions, for forecasting intra-specific changes based on genetic admixture instead of species occurrence data. Using multi-locus genotypes and extensive geographic coverage of distribution data across the European Alps, we applied this approach to 20 alpine plant species considering a global increase in temperature from 0.25 to 4 °C. We forecasted the magnitudes of displacement of contact zones between plant populations potentially adapted to warmer environments and other populations. While a global trend of movement in a north-east direction was predicted, the magnitude of displacement was species-specific. For a temperature increase of 2 °C, contact zones were predicted to move by 92 km on average (minimum of 5 km, maximum of 212 km) and by 188 km for an increase of 4 °C (minimum of 11 km, maximum of 393 km). Intra-specific turnover-measuring the extent of change in global population genetic structure-was generally found to be moderate for 2 °C of temperature warming. For 4 °C of warming, however, the models indicated substantial intra-specific turnover for ten species. These results illustrate that, in spite of unavoidable simplifications, ancestry distribution models open new perspectives to forecast population genetic changes within species and complement more traditional distribution-based approaches.

Genetic structure of a linear population of Beta vulgaris ssp maritima (sea beet) revealed by isozyme and RFLP analysis

Knowledge of the genetic structure of plant populations is necessary for the understanding of the dynamics of major ecological processes. It also has applications in conservation biology and risk assessment for genetically modified crops. This paper reports the genetic structure of a linear population of sea beet, Beta vulgaris ssp. maritima (the wild relative of sugar beet), on Furzey Island, Poole Harbour. The relative spatial positions of the plants were accurately mapped and the plants were scored for variation at isozyme and RFLP loci. Structure was analysed by repeated subdivision of the population to find the average size of a randomly mating group. Estimates of F-ST between randomly mating units were then made, and gave patterns consistent with the structure of the population being determined largely by founder effects. The implications of these results for the monitoring of transgene spread in wild sea beet populations are discussed.

Reduced genetic diversity, increased isolation and multiple introductions of invasive giant hogweed in the western Swiss Alps.

Abstract The giant hogweed (Heracleum mantegazzianum) has successfully invaded 19 European countries as well as parts of North America. It has become a problematic species due to its ability to displace native flora and to cause public health hazards. Applying population genetics to species invasion can help reconstruct invasion history and may promote more efficient management practice. We thus analysed levels of genetic variation and population genetic structure of H. mantegazzianum in an invaded area of the western Swiss Alps as well as in its native range (the Caucasus), using eight nuclear microsatellite loci together with plastid DNA markers and sequences. On both nuclear and plastid genomes, native populations exhibited significantly higher levels of genetic diversity compared to invasive populations, confirming an important founder event during the invasion process. Invasive populations were also significantly more differentiated than native populations. Bayesian clustering analysis identified five clusters in the native range that corresponded to geographically and ecologically separated groups. In the invaded range, 10 clusters occurred. Unlike native populations, invasive clusters were characterized by a mosaic pattern in the landscape, possibly caused by anthropogenic dispersal of the species via roads and direct collection for ornamental purposes. Lastly, our analyses revealed four main divergent groups in the western Swiss Alps, likely as a consequence of multiple independent establishments of H. mantegazzianum.

Jasmonic acid and the THO/TREX complex participate in phosphate starvation response

All plants are typically confronted to simultaneous biotic and abiotic stress throughout their life cycle. Low inorganic phosphate (Pi) is the most common nutrient deficiency limiting plant growth in natural and agricultural ecosystems while insect herbivory accounts for major losses in plant productivity and impacts on ecological and evolutionary changes in plant populations. Here we report that plants experiencing Pi deficiency induce the jasmonic acid (JA) pathway and enhance their defence against insect herbivory. The phol mutant is impaired in the translocation of Pi from roots to shoots and shows the typical symptoms associated with Pi deficiency, including high anthocyanin and poor shoot growth. These phol shoot phenotypes were significantly attenuated by blocking the JA biosynthesis or signalling pathways. Wounded phol leaves hyper-accumulated JA in comparison to wild type, leading to increased resistance against the generalist herbivore Spodoptera littoralis. Pi deficiency also triggered enhanced resistance to herbivory in wild-type Arabidopsis as well as tomato and tobacco, revealing that the link between Pi deficiency and JA-mediated herbivory resistance is conserved in a diversity of plants, including crops. We performed a phol suppressor screen to identify new components involved in the adaptation of plants to Pi deficiency. We report that the THO RNA TRANSCRIPTION AND EXPORT (THO/TREX) complex is a crucial component involved in modulating the Pi- deficiency response. Knockout mutants of at least three members of the THO/TREX complex, including TEX1, HPR1, and TH06, can suppress the phol shoot phenotype. Grafting experiments showed that loss of function of TEX1 only in the root was sufficient to suppress the reduced shoot growth phenotype of phol while maintaining low Pi contents. This indicates that TEX1 is involved in a long distance root-to-shoot signalling component of the Pi-deficiency response. We identified a small MYB-like transcription factor, RAD LIKE 3 (RL3), as a potential downstream target of the THO/TREX complex. RL3 expression is induced in phol mutants but attenuated in phol-7 texl-4 double mutants. Identical to knockout mutants of the THO/TREX complex, rl3 mutants can suppress the phol shoot phenotypes. Interestingly, RL3 is induced during Pi deficiency and is described in the literature as likely being mobile. It is therefore a promising new candidate involved in the root-to-shoot Pi-deficiency signalling response. Finally, we report that PHOl and its homologue PH01:H3 are involved in the co-regulation of Pi and zinc (Zn) homeostasis. PH01;H3 is up-regulated in response to Zn deficiency and, like PHOl, is expressed in the root vascular cylinder and localizes to the Golgi when expressed transiently in tobacco cells. The phol;h3 mutant accumulates more Pi as compared to wild-type when grown in Zn-deficient medium, but this increase is abolished in the phol phol;h3 double mutant. These results suggest that PH01;H3 restricts the PHOl-mediated root-to-shoot Pi transfer in responsé to Zn deficiency. Résumé Au cours de leur cycle de vie, toutes les plantes sont généralement confrontées à divers stress biotiques et abiotiques. La carence nutritionnelle la plus fréquente, limitant la croissance des plantes dans les écosystèmes naturels et agricoles, est la faible teneur en phosphate inorganique (Pi). Au niveau des stress biotiques, les insectes herbivores sont responsables de pertes majeures de rendement et ont un impact considérable sur les changements écologiques et évolutifs dans les populations des plantes. Au cours de ce travail, nous avons mis en évidence que les plantes en situation de carence en Pi induisent la voie de l'acide jasmonique (JA) et augmentent leur défense contre les insectes herbivores. Le mutant phol est déficient dans le transport du phosphate des racines aux feuilles et démontre les symptômes typiques associés à la carence, tels que la forte concentration en anthocyane et une faible croissance foliaire. Ces phénotypes du mutant phol sont significativement atténués lors d'un blocage de la voie de la biosynthèse ou des voies de signalisation du JA. La blessure des feuilles induit une hyper-accumulation de JA chez phol, résultant en une augmentation de la résistance contre l'herbivore généraliste Spodoptera littoralis. Outre Arabidopsis, la carence en Pi induit une résistance accrue aux insectes herbivores aussi chez la tomate et le tabac. Cette découverte révèle que le lien entre la carence en Pi et la résistance aux insectes herbivores via le JA est conservé dans différentes espèces végétales, y compris les plantes de grandes cultures. Nous avons effectué un crible du suppresseur de phol afin d'identifier de nouveaux acteurs impliqués dans l'adaptation de la plante à la carence en Pi. Nous rapportons que le complexe nommé THO RNA TRANSCRIPTION AND EXPORT (THO/TREX) est un élément crucial participant à la réponse des feuilles à la carence en Pi. Les mutations d'au moins trois des membres que composent le complexe THO/TREX, incluant TEX1, HPR1 et 77/06, peuvent supprimer le phénotype de phol. Des expériences de greffes ont montré que la perte de fonction de TEX1, seulement dans la racine, est suffisante pour supprimer le phénotype de la croissance réduite des parties aériennes observé chez le mutant phol, tout en maintenant de faibles teneurs en Pi foliaire. Ceci indique que TEX1 est impliqué dans la signalisation longue distance entre les racines et les parties aériennes. Nous avons identifié un petit facteur de transcription proche de la famille des MYB, RAD LIKE 3 (RL3), comme une cible potentielle en aval du complexe THO / TREX. L'expression du gène RL3 est induite dans le mutant phol mais atténuée dans le double mutant phol-7 texl-4. Exactement comme les plantes mutées d'un des membres du complexe THO/TREX, le mutant rl3 peut supprimer le phénotype foliaire de phol. RL3 est induit au cours d'une carence en Pi et est décrit dans la littérature comme étant potentiellement mobile. Par conséquent, il serait un nouveau candidat potentiellement impliqué dans la réponse longue distance entre les racines et les parties aériennes lors d'un déficit en Pi. Enfin, nous reportons que PHOl et son homologue PHOl: H3 sont impliqués dans la co- régulation de l'homéostasie du Pi et du zinc (Zn). PHOl; H3 est sur-exprimé en réponse au déficit en Zn et, comme PHOl, est exprimé dans les tissus vasculaires des racines et se localise dans l'appareil de Golgi lorsqu'il est exprimé de manière transitoire dans des cellules de tabac. Le mutant phol; h3 accumule plus de Pi par rapport aux plantes sauvages lorsqu'il est cultivé sur un milieu déficient en Zn, mais cette augmentation en Pi est abolie dans le double mutant phol phol; h3. Ces résultats suggèrent qu'en réponse à une carence en Zn, PHOl; H3 limite l'action de PHOl et diminue le transfert du Pi des racines aux parties aériennes.

Phosphate Deficiency Induces the Jasmonate Pathway and Enhances Resistance to Insect Herbivory.

During their life cycle, plants are typically confronted by simultaneous biotic and abiotic stresses. Low inorganic phosphate (Pi) is one of the most common nutrient deficiencies limiting plant growth in natural and agricultural ecosystems, while insect herbivory accounts for major losses in plant productivity and impacts ecological and evolutionary changes in plant populations. Here, we report that plants experiencing Pi deficiency induce the jasmonic acid (JA) pathway and enhance their defense against insect herbivory. Pi-deficient Arabidopsis (Arabidopsis thaliana) showed enhanced synthesis of JA and the bioactive conjugate JA-isoleucine, as well as activation of the JA signaling pathway, in both shoots and roots of wild-type plants and in shoots of the Pi-deficient mutant pho1 The kinetics of the induction of the JA signaling pathway by Pi deficiency was influenced by PHOSPHATE STARVATION RESPONSE1, the main transcription factor regulating the expression of Pi starvation-induced genes. Phenotypes of the pho1 mutant typically associated with Pi deficiency, such as high shoot anthocyanin levels and poor shoot growth, were significantly attenuated by blocking the JA biosynthesis or signaling pathway. Wounded pho1 leaves hyperaccumulated JA/JA-isoleucine in comparison with the wild type. The pho1 mutant also showed an increased resistance against the generalist herbivore Spodoptera littoralis that was attenuated in JA biosynthesis and signaling mutants. Pi deficiency also triggered increased resistance to S. littoralis in wild-type Arabidopsis as well as tomato (Solanum lycopersicum) and Nicotiana benthamiana, revealing that the link between Pi deficiency and enhanced herbivory resistance is conserved in a diversity of plants, including crops.

Plasticity in sex allocation in the plant Mercurialis annua is greater for hermaphrodites sampled from dimorphic than from monomorphic populations.

Plants are notoriously variable in gender, ranging in sex allocation from purely male through hermaphrodite to purely female. This variation can have both a genetic and an adaptive plastic component. In gynodioecious species, where females co-occur with hermaphrodites, hermaphrodites tend to shift their allocation towards greater maleness when growing under low-resource conditions, either as a result of hermaphrodites shifting away from an expensive female function, or because of enhanced siring advantages in the presence of females. Similarly, in the androdioecious plant Mercurialis annua, where hermaphrodites co-exist with males, hermaphrodites also tend to enhance their relative male allocation under low-resource conditions. Here, we ask whether this response differs between hermaphrodites that have been evolving in the presence of males, in a situation analogous to that supposed for gynodioecious populations, vs. those that have been evolving in their absence. We grew hermaphrodites of M. annua from populations in which males were either present or absent under different levels of nutrient availability and compared their reaction norms. We found that, overall, hermaphrodites from populations with males tended to be more female than those from populations lacking males. Importantly, hermaphrodites' investment in pollen and seed production was more plastic when they came from populations with males than without them, reducing their pollen production at low resource availability and increasing their seed production at high resource availability. These results are consistent with the hypothesis that plasticity in sex allocation is enhanced in hermaphrodites that have likely been exposed to variation in mating opportunities due to fluctuations in the frequency of co-occurring males.

Disjunct populations of European vascular plant species keep the same climatic niches

Aim Previous research on how climatic niches vary across species ranges has focused on a limited number of species, mostly invasive, and has not, to date, been very conclusive. Here we assess the degree of niche conservatism between distant populations of native alpine plant species that have been separated for thousands of years. Location European Alps and Fennoscandia. Methods Of the studied pool of 888 terrestrial vascular plant species occurring in both the Alps and Fennoscandia, we used two complementary approaches to test and quantify climatic-niche shifts for 31 species having strictly disjunct populations and 358 species having either a contiguous or a patchy distribution with distant populations. First, we used species distribution modelling to test for a region effect on each species' climatic niche. Second, we quantified niche overlap and shifts in niche width (i.e. ecological amplitude) and position (i.e. ecological optimum) within a bi-dimensional climatic space. Results Only one species (3%) of the 31 species with strictly disjunct populations and 58 species (16%) of the 358 species with distant populations showed a region effect on their climatic niche. Niche overlap was higher for species with strictly disjunct populations than for species with distant populations and highest for arctic-alpine species. Climatic niches were, on average, wider and located towards warmer and wetter conditions in the Alps. Main conclusion Climatic niches seem to be generally conserved between populations that are separated between the Alps and Fennoscandia and have probably been so for 10,000-15,000 years. Therefore, the basic assumption of species distribution models that a species' climatic niche is constant in space and time - at least on time scales 104 years or less - seems to be largely valid for arctic-alpine plants.

Cadmium hyperaccumulation and genetic differentiation of Thlaspi caerulescens populations

Although the knowledge on heavy metal hyperaccumulation mechanisms is increasing, the genetic basis of cadmium (Cd) hyperaccurnulation remains to be elucidated. Thlaspi caerulescens is an attractive model since Cd accumulation polymorphism observed in this species suggests genetic differences between populations with low versus high Cd hyperaccumulation capacities. In our study, a methodology is proposed to analyse at a regional scale the genetic differentiation of T. caerulescens natural populations in relation to Cd hyperaccumulation capacity while controlling for different environmental, soil, plant parameters and geographic origins of populations. Twenty-two populations were characterised with AFLP markers and cpDNA polymorphism. Over all loci, a partial Mantel test showed no significant genetic structure with regard to the Cd hyperaccumulation capacity. Nevertheless, when comparing the marker variation to a neutral model, seven AFLP fragments (9% of markers) were identified as presenting particularly high genetic differentiation between populations with low and high Cd hyperaccurnulation capacity. Using simulations, the number of outlier loci was showed to be significantly higher than expected at random. These loci presented a genetic structure linked to Cd hyperaccumulation capacity independently of the geography, environment, soil parameters and Zn, Pb, Fe and Cu concentrations in plants. Using a canonical correspondence analysis, we identified three of them as particularly related to the Cd hyperaccumutation capacity. This study demonstrates that populations with low and high hyperaccurnulation capacities can be significantly distinguished based on molecular data. Further investigations with candidate genes and mapped markers may allow identification and characterization of genomic regions linked to factors involved in Cd hyperaccumulation.

Genetic diversity and differentiation processes in the ploidy series of Olea europaea L.: a multiscale approach from subspecies to insular populations.

Geographical isolation and polyploidization are central concepts in plant evolution. The hierarchical organization of archipelagos in this study provides a framework for testing the evolutionary consequences for polyploid taxa and populations occurring in isolation. Using amplified fragment length polymorphism and simple sequence repeat markers, we determined the genetic diversity and differentiation patterns at three levels of geographical isolation in Olea europaea: mainland-archipelagos, islands within an archipelago, and populations within an island. At the subspecies scale, the hexaploid ssp. maroccana (southwest Morocco) exhibited higher genetic diversity than the insular counterparts. In contrast, the tetraploid ssp. cerasiformis (Madeira) displayed values similar to those obtained for the diploid ssp. guanchica (Canary Islands). Geographical isolation was associated with a high genetic differentiation at this scale. In the Canarian archipelago, the stepping-stone model of differentiation suggested in a previous study was partially supported. Within the western lineage, an east-to-west differentiation pattern was confirmed. Conversely, the easternmost populations were more related to the mainland ssp. europaea than to the western guanchica lineage. Genetic diversity across the Canarian archipelago was significantly correlated with the date of the last volcanic activity in the area/island where each population occurs. At the island scale, this pattern was not confirmed in older islands (Tenerife and Madeira), where populations were genetically homogeneous. In contrast, founder effects resulted in low genetic diversity and marked genetic differentiation among populations of the youngest island, La Palma.