Plant molecular phylogeography in China and adjacent regions: Tracing the genetic imprints of Quaternary climate and environmental change in the world's most diverse temperate flora

The Sino-Japanese Floristic Region (SJFR) of East Asia harbors the most diverse of the world's temperate flora, and was the most important glacial refuge for its Tertiary representatives ('relics') throughout Quaternary ice-age cycles. A steadily increasing number of phylogeographic studies in the SJFR of mainland China and adjacent areas, including the Qinghai-Tibetan-Plateau (QTP) and Sino-Himalayan region, have documented the population histories of temperate plant species in these regions. Here we review this current literature that challenges the oft-stated view of the SJFR as a glacial sanctuary for temperate plants, instead revealing profound effects of Quaternary changes in climate, topography, and/or sea level on the current genetic structure of such organisms. There are three recurrent phylogeographic scenarios identified by different case studies that broadly agree with longstanding biogeographic or palaeo-ecological hypotheses: (i) postglacial re-colonization of the QTP from (south-)eastern glacial refugia; (ii) population isolation and endemic species formation in Southwest China due to tectonic shifts and river course dynamics; and (iii) long-term isolation and species survival in multiple localized refugia of (warm-)temperate deciduous forest habitats in subtropical (Central/East/South) China. However, in four additional instances, phylogeographic findings seem to conflict with a priori predictions raised by palaeo-data, suggesting instead: (iv) glacial in situ survival of some hardy alpine herbs and forest trees on the QTP platform itself; (v) long-term refugial isolation of (warm-)temperate evergreen taxa in subtropical China; (vi) 'cryptic' glacial survival of (cool-)temperate deciduous forest trees in North China; and (vii) unexpectedly deep (Late Tertiary/early-to-mid Pleistocene) allopatric-vicariant differentiation of disjunct lineages in the East China-Japan-Korea region due to past sea transgressions. We discuss these and other consequences of the main phylogeographic findings in light of palaeo-environmental evidence, emphasize notable gaps in our knowledge, and outline future research prospects for disentangling the evolution and biogeographic history of the region's extremely diverse temperate flora.

Differential cycles of range contraction and expansion in European high mountain plants during the Late Quaternary: insights from Pritzelago alpina (L.) O. Kuntze (Brassicaceae)

Nuclear DNA sequence variation of the internal transcribed spacer (ITS) and amplified fragment length polymorphisms (AFLPs) were used to illuminate the evolutionary history of Pritzelago alpina, a herbaceous perennial of (sub)alpine to nival habitats of the European high mountains. Maximum likelihood analysis of ITS sequences of P. alpina, Hornungia petraea and Hymenolobus procumbens (the 'Pritzelago alliance') resolved P. alpina and H. petraea as sister taxa. ITS divergence estimates support an origin for P. alpina in the Late Tertiary, while intraspecific diversification started in the Late Quaternary (0.4-0.9 million years ago). AFLP analysis of 76 individuals of P. alpina, representing 24 localities across its entire west-east distribution, identified four mountain lineages in Cantabria, the Pyrenees, (south-) western Alps, and northeastern Alps/Tatras/Carpathians. In an analysis of molecular variance (amova), 14.3% of the total variation derived from this separation. However, relationships among these lineages remained unresolved in neighbour-joining and principal co-ordinates analyses, suggesting a population history of near simultaneous vicariance. Comparison with our previous ITS/AFLP study of Anthyllis montana (Fabaceae) indicates that the two co-distributed but altitudinally differentiated plant species exhibit temporally concordant but spatially discordant patterns of genetic variation. Moreover, levels of AFLP divergence were significantly lower in P. alpina than in the submediterranean, lower-elevation A. montana. Together, these data are consistent with a 'displacement refugia model', which predicts that European mountain plant species associated with lower- and upper-elevation habitats had a different cycle of range contraction into (long-term) glacial and (short-term) interglacial refugia, respectively.

Late Quaternary distributional stasis in the submediterranean mountain plant Anthyllis montana L. (Fabaceae) inferred from ITS sequences and amplified fragment length polymorphism markers

Anthyllis montana is a submediterranean, herbaceous plant of the southern and central European mountains. The internal transcribed spacer (ITS) regions of nuclear ribosomal DNA were sequenced from multiple accessions of the species and several closely related taxa. In addition, amplified fragment length polymorphism (AFLP) was analysed from 71 individuals of A. montana collected in 20 localities, mainly in the Pyrenees, Alps, Italian Peninsula and Balkans. Our ITS phylogeny showed a sequential branching pattern in A. montana, implying a western Mediterranean origin followed by an eastward migration. ITS clock calibrations suggest that speciation of A. montana took place at the Pliocene-Pleistocene boundary, while intraspecific divergence dates to Late Quaternary times (i.e. 0.7 million years ago). The AFLP analyses revealed a major genetic (west/east) subdivision within A. montana, probably caused by the massive glaciation of the Alps during this latter time period. The present-day absence of A. montana from vast parts of the Alps, which appear ecologically suitable for the species, together with the finding of evenly distributed AFLP variability within each of the two western and eastern lineages identified, is taken as evidence for a largely static Late Quaternary history without large-scale migration. High levels of AFLP variation observed among populations, together with weak or absent patterns of isolation by distance, seem to be in accord with long-term population insularization and distributional stasis. However, recent small-scale migration and a narrow hybrid zone between western and eastern lineages need to be postulated to explain the intermediate genetic composition of individuals from the Maritime Alps, a well-known suture-zone for other plant and animal species.

Evolutionary processes in a continental island system: molecular phylogeography of the Aegean Nigella arvensis alliance (Ranunculaceae) inferred from chloroplast DNA

Continental shelf island systems, created by rising sea levels, provide a premier setting for studying the effects of past fragmentation, dispersal, and genetic drift on taxon diversification. We used phylogeographical (nested clade) and population genetic analyses to elucidate the relative roles of these processes in the evolutionary history of the Aegean Nigella arvensis alliance (= 'coenospecies'). We surveyed chloroplast DNA (cpDNA) variation in 455 individuals from 47 populations (nine taxa) of the alliance throughout its core range in the Aegean Archipelago and surrounding mainland areas of Greece and Turkey. The study revealed the presence of three major lineages, with largely nonoverlapping distributions in the Western, Central, and Eastern Aegean. There is evidence supporting the idea that these major lineages evolved in situ from a widespread (pan-Aegean) ancestral stock as a result of multiple fragmentation events, possibly due to the influence of post-Messinian sea flooding, Pleistocene eustatic changes and corresponding climate fluctuations. Over-sea dispersal and founder events appear to have played a rather insignificant role in the group's history. Rather, all analytical approaches identified the alliance as an organism group with poor seed dispersal capabilities and a susceptibility to genetic drift. In particular, we inferred that the observed level of cpDNA differentiation between Kikladian island populations of Nigella degenii largely reflects population history, (viz. Holocene island fragmentation) and genetic drift in the near absence of seed flow since their time of common ancestry. Overall, our cpDNA data for the N. arvensis alliance in general, and N. degenii in particular, indicate that historical events were important in determining the phylogeographical patterns seen, and that genetic drift has historically been relatively more influential on population structure than has cytoplasmic gene flow.

Molecular data and ecological niche modelling reveal a highly dynamic evolutionary history of the East Asian Tertiary relict Cercidiphyllum (Cercidiphyllaceae)

East Asia's temperate deciduous forests served as sanctuary for Tertiary relict trees, but their ages and response to past climate change remain largely unknown. To address this issue, we elucidated the evolutionary and population demographic history of Cercdiphyllum, comprising species in China/Japan (Cercdiphyllum japonicum) and central Japan (Cercdiphyllum magnificum). Fifty-three populations were genotyped using chloroplast and ribosomal DNA sequences and microsatellite loci to assess molecular structure and diversity in relation to past (Last Glacial Maximum) and present distributions based on ecological niche modelling. Late Tertiary climate cooling was reflected in a relatively recent speciation event, dated at the Mio-/Pliocene boundary. During glacials, the warm-temperate C. japonicum experienced massive habitat losses in some areas (north-central China/north Japan) but increases in others (southwest/-east China, East China Sea landbridge, south Japan). In China, the Sichuan Basin and/or the middle-Yangtze were source areas of postglacial northward recolonization; in Japan, this may have been facilitated through introgressive hybridization with the cool-temperate C. magnificum. Our findings challenge the notion of relative evolutionary and demographic stability of Tertiary relict trees, and may serve as a guideline for assessing the impact of Neogene climate change on the evolution and distribution of East Asian temperate plants.

Quaternary diversification in European alpine plants: pattern and process

Molecular clock approaches applied previously to European alpine plants suggest that Primula sect. Auricula, Gentiana sect. Ciminalis and Soldanella diversified at the beginning of the Quaternary or well within this period, whereas Globularia had already started diversifying in the (Late-)Tertiary. In the first part of this paper we present evidence that, in contrast to Globularia and Soldanella, the branching patterns of the molecular internal transcribed spacer phylogenies of both Primula and Gentiana are incompatible with a constant-rates birth-death model. In both of these last two taxa, speciation probably decreased through Quaternary times, perhaps because of some niche-filling process and/or a decrease in specific range size. In the second part, we apply nonlinear regression analyses to the lineage-through-time plots of P. sect. Auricula to test a range of capacity-dependent models of diversification, and the effect of Quaternary climatic oscillations on diversification and extinction. At least for one major clade of sect. Auricula there is firm evidence that both diversification and extinction are a function of temperature. Intriguingly, temperature appears to be correlated positively with extinction, but negatively with diversification. This suggests that diversification did not take place, as previously assumed, in geographical isolation in high-altitude interglacial refugia, but rather at low altitudes in geographically isolated glacial refugia.

Molecular phylogeography and ecological niche modelling of a widespread herbaceous climber, Tetrastigma hemsleyanum (Vitaceae): insights into Plio-Pleistocene range dynamics of evergreen forest in subtropical China

Warm-temperate evergreen (WTE) forest represents the typical vegetation type of subtropical China, but how its component species responded to past environmental change remains largely unknown. Here, we reconstruct the evolutionary history of Tetrastigma hemsleyanum, an herbaceous climber restricted to the WTE forest. Twenty populations were genotyped using chloroplast DNA sequences and nuclear microsatellite loci to assess population structure and diversity, supplemented by phylogenetic dating, ancestral area reconstructions and ecological niche modeling (ENM) of the species distributions during the Last Glacial Maximum (LGM) and at present. Lineages in Southwest vs Central-South-East China diverged through climate/tectonic-induced vicariance of an ancestral southern range during the early Pliocene. Long-term stability in the Southwest contrasts with latitudinal range shifts in the Central-South-East region during the early-to-mid-Pleistocene. Genetic and ENM data strongly suggest refugial persistence in situ at the LGM. Pre-Quaternary environmental changes appear to have had a persistent influence on the population genetic structure of this subtropical WTE forest species. Our findings suggest relative demographic stability of this biome in China over the last glacial-interglacial cycle, in contrast with palaeobiome reconstructions showing that this forest biome retreated to areas of today's tropical South China during the LGM.

Plant speciation in continental island floras as exemplified by Nigella in the Aegean Archipelago

Continental shelf island systems, created by rising sea levels, provide a premier setting for studying the effects of geographical isolation on non-adaptive radiation and allopatric speciation brought about by genetic drift. The Aegean Archipelago forms a highly fragmented complex of mostly continental shelf islands that have become disconnected from each other and the mainland in relatively recent geological times (ca <5.2 Ma). These ecologically fairly homogenous islands thus provide a suitable biogeographic context for assessing the relative influences of past range fragmentation, colonization, gene flow and drift on taxon diversification. Indeed, recent molecular biogeographic studies on the Aegean Nigella arvensis complex, combining phylogenetic, phylogeographic and population level approaches, exemplify the importance of allopatry and genetic drift coupled with restricted gene flow in driving plant speciation in this continental archipelago at different temporal and spatial scales. While the recent (Late Pleistocene) radiation of Aegean Nigella, as well as possible instances of incipient speciation (in the Cyclades), is shown to be strongly conditioned by (palaeo)geographic factors (including changes in sea level), shifts in breeding system (selfing) and associated isolating mechanisms have also contributed to this radiation. By contrast, founder event speciation has probably played only a minor role, perhaps reflecting a migratory situation typical for continental archipelagos characterized by niche pre-emption because of a long established resident flora. Overall, surveys of neutral molecular markers in Aegean Nigella have so far revealed population genetic processes that conform remarkably well to predictions raised by genetic drift theory. The challenge is now to gain more direct insights into the relative importance of the role of genetic drift, as opposed to natural selection, in the phenotypic and reproductive divergence among these Aegean plant species.

Evolution of East Asia's Arcto-Tertiary relict Euptelea (Eupteleaceae) shaped by Late Neogene vicariance and Quaternary climate change

BACKGROUND

The evolutionary origin and historical demography of extant Arcto-Tertiary forest species in East Asia is still poorly understood. Here, we reconstructed the evolutionary and population demographic history of the two extant Euptelea species in China (E. pleiosperma) and Japan (E. polyandra). Chloroplast/nuclear DNA sequences and microsatellite loci were obtained from 36 Euptelea populations to explore molecular structure and diversity in relation to past and present distributions based on ecological niche modelling (ENM). Time-calibrated phylogenetic/phylogeographic inferences and niche-identity tests were used to infer the historical process of lineage formation.

RESULTS

Euptelea pleiosperma diverged from E. polyandra around the Late Miocene and experienced significant ecological differentiation. A near-simultaneous diversification of six phylogroups occurred during the mid-to-late Pliocene, in response to the abrupt uplift of the eastern Tibetan Plateau and an increasingly cooler and drier climate. Populations of E. pleiosperma seem to have been mostly stationary through the last glacial cycles, while those of E. polyandra reflect more recent climate-induced cycles of range contraction and expansion.

CONCLUSIONS

Our results illustrate how Late Neogene climatic/tectonic changes promoted speciation and lineage diversification in East Asia's Tertiary relict flora. They also demonstrate for the first time a greater variation in such species' responses to glacial cycles in Japan when compared to congeners in China.

Late Quaternary history of Hippophaë rhamnoides L. (Elaeagnaceae) inferred from chalcone synthase intron (Chsi) sequences and chloroplast DNA variation

Fossil pollen records indicate that Hippophaë rhamnoides (Sea Buckthorn) was widespread on late- and early postglacial raw soils throughout much of central and northern Europe, but that Early Holocene reforestation restricted populations to northern coastal habitats, or along mountain streams in the Alps, Pyrenees, and Carpathians. We used sequence variation at the nuclear chalcone synthase intron (Chsi), in conjunction with chloroplast DNA-restriction fragment length polymorphism data, to investigate the intraspecific phylogeny, phylogeographic structure, and expansion demographic history of this dioecious and wind-pollinated shrub at its range-wide scale in Europe and Asia Minor. Four major Chsi phylogroups of unresolved relationships were identified with estimated divergences approximately 172,000 years ago. Large-scale phylogeographic structures of nuclear and cytoplasmic markers were congruent in identifying (i) southeastern Europe as the most likely source of colonization into central Europe and Scandinavia, and (ii) the area just north of the Alps as a contact zone between populations from the Alps and the east/central European-Scandinavian lineage. Coalescence-based analyses (i.e. nested clade analysis and mismatch distributions) of Chsi variation were able to detect at least four major episodes of population growth, all within about the last 40,000 years. In particular, these analyses identified a nearly synchronized timing of population expansions in various parts of the species' range in central-eastern Europe/Asia Minor, most likely correlating with the Younger Dryas Stadial ( approximately 13,000-11,600 years ago). It remains to be established whether the phylogeographic history of H. rhamnoides, and particularly its rapid response to the rapid environmental changes of the Younger Dryas cold snap, is unique to the species, or whether it is shared with other cold-tolerant shrub (or grassland) species known from late-glacial raw soils in Europe.

Patterns of genetic and phenotypic variation in Iris haynei and I. atrofusca (Iris sect. Oncocyclus = the royal irises) along an ecogeographical gradient in Israel and the West Bank

Iris haynei and I. atrofusca are two closely related narrow endemics distributed vicariously along an ecogeographical north-south gradient in Israel and the West Bank. To obtain baseline information of the taxonomic status, conservation and population history of these taxa, we investigated patterns of phenotypic variation and the partitioning of genetic variation within and among populations using dominant random amplified polymorphic DNA (RAPD) markers. Multivariate (principal components analysis) and taxonomic distance analyses based on morphometric traits from eight populations revealed no unambiguous separation into two distinct groups. Results of genetic analyses for nine populations differed only slightly when either allele- or marker-based approaches were employed. Mean within-population diversity was high (0.258 for Nei's expected heterozygosity), but there was no significant relationship between genetic diversity and either population size or latitude. Although the range-wide estimate of GST ( approximately 0.20) revealed relatively high differentiation among populations this value was inflated because of a small, but significant, component of molecular variance among regions viz. taxa ( approximately 5%). Limited long-distance dispersal capabilities in conjunction with a linearized habitat distribution are proposed to contribute to the approximate isolation by distance pattern observed. It also appears that extant populations are currently deviating from equilibrium conditions because of primary divergence of a formerly more widespread ancestral population. Given the absence of deep genetic and phenotypic subdivision among northern (I. haynei) vs. central/southern (I. atrofusca) populations, we argue for a revision of their species status. Nonetheless, we recommend conservation attention to these geographically differentiated segments as separate management units, which can be seen as an instructive example of incipient species formation.

A strong 'filter' effect of the East China Sea land bridge for East Asia's temperate plant species: inferences from molecular phylogeography and ecological niche modelling of Platycrater arguta (Hydrangeaceae)

Background

In East Asia, an increasing number of studies on temperate forest tree species find evidence for migration and gene exchange across the East China Sea (ECS) land bridge up until the last glacial maximum (LGM). However, it is less clear when and how lineages diverged in this region, whether in full isolation or in the face of post-divergence gene flow. Here, we investigate the effects of Quaternary changes in climate and sea level on the evolutionary and demographic history of Platycrater arguta, a rare temperate understorey shrub with disjunct distributions in East China (var. sinensis) and South Japan (var. arguta). Molecular data were obtained from 14 P. arguta populations to infer current patterns of molecular structure and diversity in relation to past (Last Interglacial and Last Glacial Maximum) and present distributions based on ecological niche modelling (ENM). A coalescent-based isolation-with-migration (IM) model was used to estimate lineage divergence times and population demographic parameters.

Results

Combining information from nuclear/chloroplast sequence data with nuclear microsatellites, our IM analyses identify the two varieties as genetically distinct units that evolved in strict allopatry since the mid-Pleistocene, c. 0.89 (0.51–1.2) Ma. Together with Bayesian Skyeline Plots, our data further suggest that both lineages experienced post-divergence demographic growth, followed by refugial isolation, divergence, and in the case of var. arguta post-glacial admixture. However, past species distribution modelling indicates that the species’ overall distribution has not greatly changed over the last glacial cycles.

Conclusions

Our findings highlight the important influence of ancient sea-level changes on the diversification of East Asia’s temperate flora. Implicitly, they challenge the notion of general temperate forest expansion across the ECS land bridge, demonstrating instead its ‘filter’ effect owing to an unsuitable environment for certain species and their biological (e.g., recruitment) properties.

Molecular phylogeography of East Asian Kirengeshoma (Hydrangeaceae) in relation to quaternary climate change and landbridge configurations

Kirengeshoma comprises two species inhabiting warm temperate-deciduous forests in East China/South Japan (Kirengeshoma palmata) and South Korea (Kirengeshoma koreana). A survey of chloroplast (cp) DNA and inter-simple sequence repeats (ISSRs) variation in Kirengeshoma was carried out to determine the population history of a plant taxon around the East China Sea (ECS). CpDNA and ISSRs revealed lower genetic divergence between China and Japan relative to the other contrasts, in line with intrageneric classification. Molecular dating suggests that K. koreana diverged at the Plio-Pleistocene boundary from the Japanese populations, whereas the latter migrated into China during the early-to-mid Pleistocene via the ECS basin. Vicariant segregation of Chinese and Japanese populations likely occurred during the mid-Pleistocene. Mismatch distributions and neutrality tests indicated that Chinese populations expanded their range during the Late Pleistocene, probably during a cold period, whereas those from Japan showed no significant population growth. We conclude that the current distribution and differentiation of components of presently isolated warm temperate-deciduous forests in China, Japan and Korea likely resulted from a combination of relatively ancient vicariant and immigration events, and those from Japan were particularly sensitive to range fragmentation and long-term refugial isolation throughout the Late Pleistocene.

Long-distance dispersal vs vicariance: the origin and genetic diversity of alpine plants in the Spanish Sierra Nevada

Here, we investigated the origin and genetic diversity of four alpine plant species co-occurring in the Spanish Sierra Nevada and other high mountains in south-western Europe by analysis of amplified fragment length polymorphisms (AFLPs). In Kernera saxatilis, Silene rupestris and Gentiana alpina we found intraspecific phylogroups corresponding to mountain regions as predicted by the vicariance hypothesis. Moreover, genetic distances between Sierra Nevada and Pyrenees populations were always higher than those between populations from the Pyrenees and the south-western Alps/Massif Central. This suggests successive disruption of gene exchange between mountain ranges as postglacial climatic warming proceeded from south to north. In Papaver alpinum, our data indicate that a central Pyrenean population arose via long-distance dispersal from the Sierra Nevada, and that vicariant separation events between the Sierra Nevada and the Pyrenees and between the Pyrenees and the south-western Alps occurred simultaneously. Overall, Sierra Nevada populations of all species investigated here preserve unexpectedly high (or not exceptionally reduced) genetic diversity. This testifies to the important influence of long-term isolation, i.e. vicariance, on genetic diversity through fostering the accumulation of new mutations and/or the fixation of ancestral ones.

An AFLP clock for the absolute dating of shallow-time evolutionary history based on the intraspecific divergence of southwestern European alpine plant species

The dating of recent events in the history of organisms needs divergence rates based on molecular fingerprint markers. Here, we used amplified fragment length polymorphisms (AFLPs) of three distantly related alpine plant species co-occurring in the Spanish Sierra Nevada, the Pyrenees and the southwestern Alps/Massif Central to establish divergence rates. Within each of these species (Gentiana alpina, Kernera saxatilis and Silene rupestris), we found that the degree of AFLP divergence (D(N72)) between mountain phylogroups was significantly correlated with their time of divergence (as inferred from palaeoclimatic/palynological data), indicating constant AFLP divergence rates. As these rates did not differ significantly among species, a regression analysis based on the pooled data was utilized to generate a general AFLP rate. The application of this latter rate to AFLP data from other herbaceous plant species (Minuartia biflora: Schönswetter et al. 2006; Nigella degenii: Comes et al. 2008) resulted in a plausible timing of the recolonization of the Svalbard Islands and the separation of populations from the Alps and Scandinavia (Minuartia), and of island population separation in the Aegean Archipelago (Nigella). Furthermore, the AFLP mutation rate obtained in our study is of the same magnitude as AFLP mutation rates published previously. The temporal limits of our AFLP rate, which is based on intraspecific vicariance events at shallow (i.e. late glacial/Early Holocene) time scales, remains to be tested.

Plastome organization, genome-based phylogeny and evolution of plastid genes in Podophylloideae (Berberidaceae)

Species of Podophylloideae (Berberidaceae, Ranunculales) are of great pharmacogenetic importance and represent the classic biogeographic disjunction between eastern Asia (EA; 10 ssp.) and eastern North America (ENA; 2 ssp.). However, previous molecular studies of this group suffered from low phylogenetic resolution and/or insufficient marker variability. This study is the first to report whole-plastome sequence data for all 12 species of Podophylloideae (14 individuals) and a close relative, Achlys triphylla. These 15 plastomes proved highly similar in overall size (156,240-157,370 bp), structure, gene order and content, also when compared to other Ranunculales, but also revealed some structural variations caused by the expansion or contraction of the inverted repeats (IRs) into or out of adjacent single-copy regions. Our phylogenomic analysis, based on 63 plastome-derived protein-coding genes (CDS), supported the monophyly of Podophylloideae and its two major genera (EA: Dysosma, EA/ENA: Diphylleia), with Podophyllum peltatum L. (ENA) being more closely related to Diphylleia than to the group's earliest diverging species, Sinopodophyllum hexandrum (EA). Furthermore, within this subfamily/dataset, matK was identified as the fastest evolving gene, which proved to be under positive selection especially in more recently derived, lower-elevation lineages of Dysosma, possibly reflecting an adaptive response to novel environmental (i.e. subtropical compared to higher-elevation/alpine) conditions. Finally, several highly variable noncoding regions were identified in the plastomes of Podophylloideae and Ranunculales. These highly variable loci should be the best choices for future phylogenetic, phylogeographic, and population-level genetic studies. Overall, our results demonstrate the power of plastid phylogenomics to improve phylogenetic resolution, and contribute to a better understanding of plastid gene evolution in Podophylloideae.

Contributions of historical and contemporary geographic and environmental factors to phylogeographic structure in a Tertiary relict species, Emmenopterys henryi (Rubiaceae)

Examining how historical and contemporary geographic and environmental factors contribute to genetic divergence at different evolutionary scales is a central yet largely unexplored question in ecology and evolution. Here, we examine this key question by investigating how environmental and geographic factors across different epochs have driven genetic divergence at deeper (phylogeographic) and shallower (landscape genetic) evolutionary scales in the Chinese Tertiary relict tree Emmenopterys henryi. We found that geography played a predominant role at all levels – phylogeographic clades are broadly geographically structured, the deepest levels of divergence are associated with major geological or pre-Quaternary climatic events, and isolation by distance (IBD) primarily explained population genetic structure. However, environmental factors are clearly also important – climatic fluctuations since the Last Interglacial (LIG) have likely contributed to phylogeographic structure, and the population genetic structure (in our AFLP dataset) was partly explained by isolation by environment (IBE), which may have resulted from natural selection in environments with divergent climates. Thus, historical and contemporary geography and historical and contemporary environments have all shaped patterns of genetic structure in E. henryi, and, in fact, changes in the landscape through time have also been critical factors.

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.

Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end

BACKGROUND

The transition from outcrossing to selfing has long been portrayed as an 'evolutionary dead end' because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a 'rostellum') co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing.

RESULTS

The eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.4-3.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. ≤ 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing.

CONCLUSIONS

The frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the 'rostellum') is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales.