A proteomics approach to investigate the process of Zn hyperaccumulation in Noccaea caerulescens (J & C. Presl) F.K. Meyer

Zinc (Zn) is an essential trace element in all living organisms, but is toxic in excess. Several plant species are able to accumulate Zn at extraordinarily high concentrations in the leaf epidermis without showing any toxicity symptoms. However, the molecular mechanisms of this phenomenon are still poorly understood. A state‐of‐the‐art quantitative 2D liquid chromatography/tandem mass spectrometry (2D‐LC‐MS/MS) proteomics approach was used to investigate the abundance of proteins involved in Zn hyperaccumulation in leaf epidermal and mesophyll tissues of Noccaea caerulescens. Furthermore, the Zn speciation in planta was analyzed by a size‐exclusion chromatography/inductively coupled plasma mass spectrometer (SEC‐ICP‐MS) method, in order to identify the Zn‐binding ligands and mechanisms responsible for Zn hyperaccumulation. Epidermal cells have an increased capability to cope with the oxidative stress that results from excess Zn, as indicated by a higher abundance of glutathione S‐transferase proteins. A Zn importer of the ZIP family was more abundant in the epidermal tissue than in the mesophyll tissue, but the vacuolar Zn transporter MTP1 was equally distributed. Almost all of the Zn located in the mesophyll was stored as Zn–nicotianamine complexes. In contrast, a much lower proportion of the Zn was found as Zn–nicotianamine complexes in the epidermis. However, these cells have higher concentrations of malate and citrate, and these organic acids are probably responsible for complexation of most epidermal Zn. Here we provide evidence for a cell type‐specific adaptation to excess Zn conditions and an increased ability to transport Zn into the epidermal vacuoles.     

Chloroplast DNA variation and postglacial recolonization of common ash (Fraxinus excelsior L.) in Europe

We used chloroplast polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and chloroplast microsatellites to assess the structure of genetic variation and postglacial history across the entire natural range of the common ash (Fraxinus excelsior L.), a broad-leaved wind-pollinated and wind-dispersed European forest tree. A low level of polymorphism was observed, with only 12 haplotypes at four polymorphic microsatellites in 201 populations, and two PCR-RFLP haplotypes in a subset of 62 populations. The clear geographical pattern displayed by the five most common haplotypes was in agreement with glacial refugia for ash being located in Iberia, Italy, the eastern Alps and the Balkan Peninsula, as had been suggested from fossil pollen data. A low chloroplast DNA mutation rate, a low effective population size in glacial refugia related to ash's life history traits, as well as features of postglacial expansion were put forward to explain the low level of polymorphism. Differentiation among populations was high (GST= 0.89), reflecting poor mixing among recolonizing lineages. Therefore, the responsible factor for the highly homogeneous genetic pattern previously identified at nuclear microsatellites throughout western and central Europe (Heuertz et al. 2004) must have been efficient postglacial pollen flow. Further comparison of variation patterns at both marker systems revealed that nuclear microsatellites identified complex differentiation patterns in south-eastern Europe which remained undetected with chloroplast microsatellites. The results suggest that data from different markers should be combined in order to capture the most important genetic patterns in a species.     

Phytochelatin synthesis is not responsible for Cd tolerance in the Zn/Cd hyperaccumulator Thlaspi caerulescens (J. & C. Presl)

Thlaspi caerulescens (J. & C. Presl, "Prayon") is a heavy-metal hyperaccumulator that accumulates Zn and Cd to high concentrations (40,000 and 4,000 mg kg DW-1 respectively) without phytotoxicity. The mechanism of Cd tolerance has not been characterized but reportedly involves vacuolar sequestration. The role of phytochelatins (PCs) in metal tolerance in T. caerulescens and the related non-accumulator T. arvense was examined. Although PCs were produced by both species in response to Cd, these peptides do not appear to be involved in metal tolerance in the hyperaccumulator. Leaf and root PC levels for both species showed a similar positive correlation with tissue Cd, but total PC levels in the hyperaccumulator were generally lower, despite correspondingly higher metal concentrations. The lack of a role for PCs in the hyperaccumulator's response to metal stress suggests that other mechanisms are responsible Cd tolerance. The lower level of leaf PCs in T. caerulescens also implies that Cd in the shoot is sequestered in a compartment or form that does not elicit a PC response.     

Molecular data reveal isolation by distance and past population expansion for the shea tree (Vitellaria paradoxa C.F. Gaertn) in West Africa

While the genetic structure of many tree species in temperate, American and Asian regions is largely explained by climatic oscillations and subsequent habitat contractions and expansions, little is known about Africa. We investigated the genetic diversity and structure of shea tree (Vitellaria paradoxa,) in Western Africa, an economically important tree species in the Sudano-Sahelian zone. Eleven nuclear microsatellites (nuc) were used to genotype 673 trees selected in 38 populations. They revealed moderate to high within-population diversity: allelic richness ranged from R(nuc) = 3.99 to 5.63. This diversity was evenly distributed across West Africa. Populations were weakly differentiated (F(STnuc) = 0.085; P < 0.0001) and a pattern of isolation by distance was noted. No phylogeographic signal could be detected across the studied sample. Additionally, two chloroplast microsatellite loci, leading to 11 chlorotypes, were used to analyse a sub-set of 370 individuals. Some variation in chloroplast allelic richness among populations could be detected (R(cp) = 0.00 to 4.36), but these differences were not significant. No trend with latitude and longitude were observed. Differentiation was marked (G(STcp) = 0.553; P < 0.0001), but without a significant phylogeographical signal. Population expansion was detected considering the total population using approximate Bayesian computation (nuclear microsatellites) and mismatch distribution (chloroplast microsatellites) methods. This expansion signal and the isolation by distance pattern could be linked to the past climatic conditions in West Africa during the Pleistocene and Holocene which should have been favourable to shea tree development. In addition, human activities through agroforestry and domestication (started 10,000 bp) have probably enhanced gene flow and population expansion.     

Multiple origin of metallicolous populations of the pseudometallophyte Arabidopsis halleri (Brassicaceae) in central Europe: the cpDNA testimony

The population structure of the pseudo-metallophyte herb, Arabidopsis halleri, was studied using Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) on chloroplast DNA (cpDNA). The history of metallicolous (M) populations showing increased zinc tolerance was investigated. Eight primer-enzyme combinations out of 72 tested were applied to a total of 625 individuals from 28 widespread populations, 14 of them being M. Eleven distinct chlorotypes were found: five were common to nonmetallicolous (NM) and M populations, whereas six were only observed in one edaphic type (five in NM and one in M). No difference in chlorotype diversity between edaphic types was detected. Computed on the basis of chlorotype frequencies, the level of population differentiation was high but remained the same when taking into account levels of molecular divergence between chlorotypes. Isolation by distance was largely responsible for population differentiation. Geographically isolated groups of M populations were more genetically related to their closest NM populations than to each other. Our results suggest that M populations have been founded separately from distinct NM populations without suffering founding events and that the evolution towards increased tolerance observed in the distinct M population groups occurred independently.     

Influence of hydrogeographic history and hybridization on the distribution of genetic variation in the pupfishes Cyprinodon atrorus and C. bifasciatus

The evolutionary importance of hybridization in animals has been subject of much debate. In this study, we examined the influence of hydrogeographic history and hybridization on the present distribution of nuclear and mitochondrial DNA variation in two pupfish species, Cyprinodon atrorus and Cyprinodon bifasciatus. Results presented here indicate that there has been limited introgression of nuclear genes; however, mtDNA introgression has been substantial, with complete replacement of the C. bifasciatus mitochondrial genome by that of C. atrorus. Subsequent to this replacement, there has been diversification of mitochondrial haplotypes along major geographic regions in the basin. Evidence was also found that mitochondrial replacement follows a predictable, cyclical pattern in this system, with isolation and diversification followed by re-contact and replacement of C. bifasciatus mitochondrial haplotypes by those of C. atrorus. This pattern is best explained by a combination of a numeric bias towards C. atrorus and mating site selection rather than selection for C. atrorus mitochondrial genome. These results demonstrate the important role hybridization can play in evolution.     

Anthropocene invasion of an ecosystem engineer: resolving the history of Corophium volutator (Amphipoda: Corophiidae) in the North Atlantic

Resolving the natural histories of species is important for the interpretation of ecological patterns, as it provides evolutionary context for the interactions between organisms and their environment. Despite playing an integral role on the intertidal mudflats of the North Atlantic as an abundant food source for predators and as an ecosystem engineer that alters the soft sediment environment, no previous studies have provided empirical evidence to determine the biogeographical origin of the amphipod Corophium volutator. To resolve its status as introduced or indigenous in Europe and North America, we analyzed sequence data for two mitochondrial loci and two nuclear markers, aiming to determine whether the present range of C. volutator is the result of unresolved taxonomy, persistence in glacial refugia, natural trans‐Atlantic dispersal, or human‐mediated introduction. Our results demonstrate a reduced genetic diversity in North American populations that is a subsample of diversity in European populations, with coalescent analyses of mitochondrial and nuclear DNA supporting different models of multiple introductions from Europe to the Bay of Fundy and Gulf of Maine in North America. These results suggest that C. volutator was introduced to North America prior to the first surveys of local biota in the 20th Century, which has broad implications for interpretations of community and ecosystem interactions in the North Atlantic intertidal. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 115 , 288–304.     

Chloroplast DNA phylogeography of European ashes, Fraxinus sp. (Oleaceae): roles of hybridization and life history traits

We investigated range-wide phylogeographic variation in three European ash species (Fraxinus sp., Oleaceae). Chloroplast DNA (cpDNA) microsatellites were typed in the thermophilous Fraxinus angustifolia and Fraxinus ornus and the observed haplotypes and the geographic distribution of diversity were compared to cpDNA data previously obtained in the more cold-tolerant Fraxinus excelsior. We found wide-ranging haplotype sharing between the phylogenetically close F. angustifolia and F. excelsior, suggesting hybridization (i) in common glacial refuges in the Iberian Peninsula, northern Italy, the eastern and/or Dinaric Alps and the Balkan Peninsula, and/or (ii) during postglacial recolonization. The data allowed us to propose additional glacial refuges for F. angustifolia in southern Italy and in Turkey, and populations from the latter region were particularly polymorphic. There was evidence for refuge areas in Italy, the Balkan Peninsula and Turkey for F. ornus, which did not share any single chloroplast haplotype with the other species. In both F. angustifolia and F. ornus, cpDNA diversity (h(S) = 0.027 and h(S) = 0.009, respectively) was lower and fixation levels (G(ST) = 0.964 and G(ST) = 0.983, respectively) higher than in sympatric F. excelsior (h(S) = 0.096, G(ST) = 0.870). These diversity patterns could be due to temperature tolerance or the demographic history.     

Glacial in situ survival in the Western Alps and polytopic autopolyploidy in Biscutella laevigata L. (Brassicaceae)

Past climatic changes and especially the ice ages have had a great impact on both the distribution and the genetic composition of plant populations, but whether they promoted speciation is still controversial. The autopolyploid complex Biscutella laevigata is a classical example of polyploidy linked to glaciations and is an interesting model to explore migration and speciation driven by climate changes in a complex alpine landscape. Diploid taxa survived the last glacial maximum in several never-glaciated areas and autotetraploids are clearly dominant in the central parts of the Alps; however, previous range-wide studies failed to identify their diploid ancestor(s). This study highlights the phylogeographical relationships of maternal lineages in the Western Alps and investigates the polyploidy process using plastid DNA sequences (trnS-trnG and trnK-intron) combined with plastid DNA length polymorphism markers, which were transferable among Brassicaceae species. Twenty-one distinct plastid DNA haplotypes were distinguished in 67 populations densely sampled in the Western Alps and main lineages were identified by a median-joining network. The external Alps harboured high levels of genetic diversity, while the Central Alps contained only a subset of haplotypes due to postglacial recolonization. Several haplotypes were restricted to local peripheral refugia and evidence of in situ survival in central nunataks was detected by the presence of highly differentiated haplotypes swamped by frequent ones. As hierarchical genetic structure pointed to an independent evolution of the species in different biogeographical districts, and since tetraploids displayed haplotypes belonging to different lineages restricted to either the northern or the southern parts of the Alpine chain, polytopic autopolyploidy was also apparent in the Western Alps.     

Postglacial recolonization in a cold climate specialist in western Europe: patterns of genetic diversity in the adder (Vipera berus) support the central–marginal hypothesis

Understanding the impact of postglacial recolonization on genetic diversity is essential in explaining current patterns of genetic variation. The central–marginal hypothesis (CMH) predicts a reduction in genetic diversity from the core of the distribution to peripheral populations, as well as reduced connectivity between peripheral populations. While the CMH has received considerable empirical support, its broad applicability is still debated and alternative hypotheses predict different spatial patterns of genetic diversity. Using microsatellite markers, we analysed the genetic diversity of the adder (Vipera berus) in western Europe to reconstruct postglacial recolonization. Approximate Bayesian Computation (ABC) analyses suggested a postglacial recolonization from two routes: a western route from the Atlantic Coast up to Belgium and a central route from the Massif Central to the Alps. This cold‐adapted species likely used two isolated glacial refugia in southern France, in permafrost‐free areas during the last glacial maximum. Adder populations further from putative glacial refugia had lower genetic diversity and reduced connectivity; therefore, our results support the predictions of the CMH. Our study also illustrates the utility of highly variable nuclear markers, such as microsatellites, and ABC to test competing recolonization hypotheses.     

Species-level phylogeographical history of Myricaria plants in the mountain ranges of western China and the origin of M. laxiflora in the Three Gorges mountain region

Myricaria species in China occur mostly in the major high-altitude mountain areas in and around the Qinghai-Tibetan Plateau. The one major exception to this is M. laxiflora which is restricted to the Three Gorges mountain region. In this study, we investigate species-level phylogeographical patterns of Myricaria species in western China and the origin of M. laxiflora . The results show that most chloroplast haplotypes are species-specific, except for one haplotype which is shared by three widespread species. Higher haplotype diversity within the Qinghai-Tibetan Plateau region supports the hypothesis that the Himalayas are the centre of origin for Myricaria . The phylogeny of Myricaria was geographically structured, and an estimated Bayesian chronology suggested the main divergence events occurred during the Late Pliocene and Early Pleistocene (~1.46–2.30 million years ago). The overall phylogeographical pattern was characterized by vicariance events and regional demographical expansion, reflecting a major influence of geological and climatic events on the evolution of Myricaria species. Our data suggest that M. laxiflora has an ancient origin, but has experienced recent population expansion through the Three Gorges Valley. The origin of M. laxiflora was estimated to be during the Early Pleistocene but its demographical expansion was more recent at about 0.015 million years ago. This highlights the unique phylogeographical history of the Three Gorges mountain region, and the deep imprint of the watercourse connections of the Yangtze River Valley on the phylogeographical structure of the species in this region.     

Bayesian inference of a complex invasion history revealed by nuclear and chloroplast genetic diversity in the colonizing plant,Silene latifolia

Species invading new ranges are subject to a series of demographic events that can strongly shape genetic diversity. Describing this demographic history is important for understanding where invasive species come from and how they spread, and is critical to testing hypotheses of postinvasion adaptation. Here, we analyse nuclear and chloroplast genetic diversity to study the invasion history of the widespread colonizing weed, Silene latifolia (Caryophyllaceae). Bayesian clustering and PCA revealed strong population structure in the native range of Europe, and although genotypes from multiple native sources were present in the introduced range of North America, the spatial distribution of genetic variance was dramatically reorganized. Using approximate Bayesian computation (ABC), we compared support for different invasion scenarios, including the number and size of independent introduction events and the amount of admixture occurring between sources of introduced genotypes. Our results supported independent introductions into eastern and western North America, with the latter forming a bridgehead for a secondary invasion into the Great Lakes region of central North America. Despite small estimated founder population sizes, the duration of the demographic bottleneck after the initial introduction appeared extremely short‐lived. This pattern of repeated colonization and rapid expansion has effectively eroded the strong population structure and cytonuclear associations present in Europe, but has retained overall high genetic diversity since invasion. Our results highlight the flexibility of the ABC approach for constructing a narrative of the demographic history of species invasions and provide baseline for future studies of evolutionary changes in introduced S. latifolia populations.     

Random amplified polymorphic DNA (RAPD) and quantitative trait analyses across a major phylogeographical break in the Mediterranean ragwort Senecio gallicus Vill. (Asteraceae)

Random amplified polymorphic DNA (RAPD) and quantitative trait variation of the widespread and ephemeral Senecio gallicus were surveyed in 11 populations sampled from the Iberian Peninsula and southern France. The aim of the study was to compare population relationships and levels of geographical differentiation with chloroplast (cp) DNA and allozyme variation assessed previously in the same populations. Employing multivariate statistics, a moderate level of intraspecific differentiation was observed among populations from Iberian coastal and inland regions for both RAPDs and quantitative traits. However, RAPDs provided greater resolution in identifying additional population structure within the hypothesized, Pleistocene refugial source area of the species in coastal Iberia. A major part of the geographical subdivision in RAPD and quantitative traits was concordant with the coastal vs. inland divergence as previously inferred from cpDNA haplotype frequencies, but strongly contrasted with the geographical uniformity of the species for allozymes. This concordance across various nuclear and cytoplasmic markers (RAPDs/quantitative traits, cpDNA) suggests that geographical uniformity for allozymes is more attributable to low rates of evolution and/or small genome sampling rather than high rates of pollen dispersal, slow rates of nuclear lineage sorting, or indirect balancing selection. The present study underscores the value of using additional classes of nuclear markers for narrowing the numbers of competing causal hypotheses about intraspecific cpDNA-allozyme discrepancies and their underlying evolutionary processes.     

Phylogeography and postglacial dispersion of the chub (Leuciscus cephalus) in Europe

A phylogeographic analysis of mitochondrial DNA variation was performed in order to test the hypothesis of a postglacial recolonization of mid- and north-European rivers from a Danubian refuge. Over 345 chub specimens from European rivers covering most of the species' native range were investigated using 600 bp of the cytochrome b gene. Chub in European rivers belong to four highly divergent mitochondrial groups (lineages) differing by mean divergence estimates from 5.2% to 7.89%. These four lineages have a largely allopatric distribution, implying four geographical sets: two Mediterranean, and two north-European sets. This pattern provided strong evidence for: (i) the eradication of this species from most of Europe during maximum ice extent; (ii) its survival in four refugia (Adriatic side of the Balkans, eastern Greece (Aegean Rivers), southern tributaries of the Danube, and periphery of Black and Caspian Seas); (iii) a differential postglacial recolonization of mid- and northern Europe from the last two refugia only; (iv) the occurrence of this recolonization in two steps for the Danubian (western) lineage that entered western Europe (Rhine-Rhone-Loire drainages) during the Riss-Würm interglacial period and survived the last glaciation there before colonizing Garonne, UK and German drainages up to the Elbe during the Holocene; and (v) the occurrence of this recolonization in a single step for the Ponto-Caspian (eastern) lineage that entered the Baltic area as far as the Oder in the Holocene. Both lineages came into contact in the River Elbe without evident mixing.     

Testing the role of ecology and life history in structuring genetic variation across a landscape: a trait‐based phylogeographic approach

Hypotheses to explain phylogeographic structure traditionally invoke geographic features, but often fail to provide a general explanation for spatial patterns of genetic variation. Organisms' intrinsic characteristics might play more important roles than landscape features in determining phylogeographic structure. We developed a novel comparative approach to explore the role of ecological and life‐history variables in determining spatial genetic variation and tested it on frog communities in Panama. We quantified spatial genetic variation within 31 anuran species based on mitochondrial DNA sequences, for which hierarchical approximate Bayesian computation analyses rejected simultaneous divergence over a common landscape. Regressing ecological variables, on genetic divergence allowed us to test the importance of individual variables revealing that body size, current landscape resistance, geographic range, biogeographic origin and reproductive mode were significant predictors of spatial genetic variation. Our results support the idea that phylogeographic structure represents the outcome of an interaction between organisms and their environment, and suggest a conceptual integration we refer to as trait‐based phylogeography.     

Phylogeography of Sophora davidii (Leguminosae) across the ‘Tanaka‐Kaiyong Line’, an important phytogeographic boundary in Southwest China

The ‘Tanaka‐Kaiyong Line’ (TKL) is a major phytogeographic boundary in Southwest China, separating East Asia's Sino‐Himalayan and Sino‐Japanese Floras. However, little is known about the importance of this boundary in promoting intraspecific phylogeographic subdivision and divergence. Using chloroplast (cpDNA) and nuclear‐intron (nDNA) sequence data, we reconstructed the population history of Sophora davidii, a drought‐tolerant riparian shrub widely distributed on either side of the TKL. Specifically, we aimed at testing two long‐standing explanations for possible vicariant events across the TKL: (i) Late Pliocene (c. 3 Ma) geological uplift of the eastern Qinghai‐Tibetan Plateau (QTP) or (ii) a sharp environmental gradient associated with the establishment of different monsoon regimes on either side of the TKL during the (Late) Pleistocene. Our genealogical analyses detected a major west–east split in cpDNA, geographically largely consistent with the TKL, and dated to c. 1.28 Ma (95% HPD: 0.21–2.96 Ma), hence postdating the latest phase of eastern QTP uplift. Furthermore, integrating cpDNA phylogeographic patterns with mismatch analyses, we found multiple refugial isolation and long‐term demographic stability of populations in the west (Hengduan Mountain Range) compared with extensive range expansions in the east, possibly during the last glacial period(s) and followed by differentiation into regional sublineages (southeast: Yunnan‐Guizhou Plateau vs. northeast: Qinling Mts./Loess Plateau). Although nuclear differentiation was less marked, the geographical pattern of nDNA haplotypes provided some further indication of the species' eastward expansion, possibly from source populations located just east of the TKL (lower Jinshajiang region). Overall, the present data reject the geological (tectonic) explanation for the TKL and, instead, provide supportive evidence for its role as a climatically driven barrier to present‐day plant dispersal. In addition, our study highlights changing temperatures and vegetation types during the last glacial period(s), along with aspects of regional topography, to be important determinants of the glacial eastward expansion of S. davidii. In consequence, our study lends support to a ‘glacial out‐of‐Hengduan Mts’. hypothesis for the xerophytic‐riparian flora of Southwest China, which in turn is inconsistent with the traditional view of the TKL as a ‘classical’ vicariant‐biogeographic boundary.