Resolving Recent Plant Radiations: Power and Robustness of Genotyping-by-Sequencing.

Disentangling species boundaries and phylogenetic relationships within recent evolutionary radiations is a challenge due to the poor morphological differentiation and low genetic divergence between species, frequently accompanied by phenotypic convergence, interspecific gene flow and incomplete lineage sorting. Here we employed a genotyping-by-sequencing (GBS) approach, in combination with morphometric analyses, to investigate a small western Mediterranean clade in the flowering plant genus Linaria that radiated in the Quaternary. After confirming the morphological and genetic distinctness of eight species, we evaluated the relative performances of concatenation and coalescent methods to resolve phylogenetic relationships. Specifically, we focused on assessing the robustness of both approaches to variations in the parameter used to estimate sequence homology (clustering threshold). Concatenation analyses suffered from strong systematic bias, as revealed by the high statistical support for multiple alternative topologies depending on clustering threshold values. By contrast, topologies produced by two coalescent-based methods (NJ$_{\mathrm{st}}$, SVDquartets) were robust to variations in the clustering threshold. Reticulate evolution may partly explain incongruences between NJ$_{\mathrm{st}}$, SVDquartets and concatenated trees. Integration of morphometric and coalescent-based phylogenetic results revealed (i) extensive morphological divergence associated with recent splits between geographically close or sympatric sister species and (ii) morphological convergence in geographically disjunct species. These patterns are particularly true for floral traits related to pollinator specialization, including nectar spur length, tube width and corolla color, suggesting pollinator-driven diversification. Given its relatively simple and inexpensive implementation, GBS is a promising technique for the phylogenetic and systematic study of recent radiations, but care must be taken to evaluate the robustness of results to variation of data assembly parameters.

The divergence history of the perennial plant Linaria cavanillesii confirms a recent loss of self-incompatibility.

Many angiosperms prevent inbreeding through a self-incompatibility (SI) system, but the loss of SI has been frequent in their evolutionary history. The loss of SI may often lead to an increase in the selfing rate, with the purging of inbreeding depression and the ultimate evolution of a selfing syndrome, where plants have smaller flowers with reduced pollen and nectar production. In this study, we used approximate Bayesian computation (ABC) to estimate the timing of divergence between populations of the plant Linaria cavanillesii that differ in SI status and in which SI is associated with low inbreeding depression but not with a transition to full selfing or a selfing syndrome. Our analysis suggests that the mixed-mating self-compatible (SC) population may have begun to diverge from the SI populations around 2810 generation ago, a period perhaps too short for the evolution of a selfing syndrome. We conjecture that the SC population of L. cavanillesii is at an intermediate stage of transition between outcrossing and selfing.

Past and future demographic dynamics of alpine species: limited genetic consequences despite dramatic range contraction in a plant from the Spanish Sierra Nevada.

Anthropogenic global climate change is expected to cause severe range contractions among alpine plants. Alpine areas in the Mediterranean region are of special concern because of the high abundance of endemic species with narrow ranges. This study combined species distribution models, population structure analyses and Bayesian skyline plots to trace the past and future distribution and diversity of Linaria glacialis, an endangered narrow endemic species that inhabits summits of Sierra Nevada (Spain). The results showed that: (i) the habitat of this alpine-Mediterranean species in Sierra Nevada suffered little changes during glacial and interglacial stages of late Quaternary; (ii) climatic oscillations in the last millennium (Medieval Warm Period and Little Ice Age) moderately affected the demographic trends of L. glacialis; (iii) future warming conditions will cause severe range contractions; and (iv) genetic diversity will not diminish at the same pace as the distribution range. As a consequence of the low population structure of this species, genetic impoverishment in the alpine zones of Sierra Nevada should be limited during range contraction. We conclude that maintenance of large effective population sizes via high mutation rates and high levels of gene flow may promote the resilience of alpine plant species when confronted with global warming.

Chemotaxonomy of iridoids in Linaria vulgaris.

The phytochemical analysis of the extracts of Linaria vulgaris, has allowed to underline an iridoidic pattern similar to that of the other Linaria plants, with the presence of antirrinoside, antirride, 6-beta-idrossiantirride, 10-beta-glucosilaucubina and a new iridoidic compound, whose structure was demonstrated to be that of 4-carboxy-boonein.

Coalescent simulations reveal hybridization and incomplete lineage sorting in Mediterranean Linaria.

We examined the phylogenetic history of Linaria with special emphasis on the Mediterranean sect. Supinae (44 species). We revealed extensive highly supported incongruence among two nuclear (ITS, AGT1) and two plastid regions (rpl32-trnL(UAG), trnS-trnG). Coalescent simulations, a hybrid detection test and species tree inference in *BEAST revealed that incomplete lineage sorting and hybridization may both be responsible for the incongruent pattern observed. Additionally, we present a multilabelled *BEAST species tree as an alternative approach that allows the possibility of observing multiple placements in the species tree for the same taxa. That permitted the incorporation of processes such as hybridization within the tree while not violating the assumptions of the *BEAST model. This methodology is presented as a functional tool to disclose the evolutionary history of species complexes that have experienced both hybridization and incomplete lineage sorting. The drastic climatic events that have occurred in the Mediterranean since the late Miocene, including the Quaternary-type climatic oscillations, may have made both processes highly recurrent in the Mediterranean flora.

Historical isolation versus recent long-distance connections between Europe and Africa in bifid toadflaxes (Linaria sect. Versicolores).

BACKGROUND: Due to its complex, dynamic and well-known paleogeography, the Mediterranean region provides an ideal framework to study the colonization history of plant lineages. The genus Linaria has its diversity centre in the Mediterranean region, both in Europe and Africa. The last land connection between both continental plates occurred during the Messinian Salinity Crisis, in the late Miocene (5.96 to 5.33 Ma). METHODOLOGY/PRINCIPAL FINDINGS: We analyzed the colonization history of Linaria sect. Versicolores (bifid toadflaxes), which includes c. 22 species distributed across the Mediterranean, including Europe and Africa. Two cpDNA regions (rpl32-trnL(UAG) and trnK-matK) were sequenced from 66 samples of Linaria. We conducted phylogenetic, dating, biogeographic and phylogeographic analyses to reconstruct colonization patterns in space and time. Four major clades were found: two of them exclusively contain Iberian samples, while the other two include northern African samples together with some European samples. The bifid toadflaxes have been split in African and European clades since the late Miocene, and most lineage and speciation differentiation occurred during the Pliocene and Quaternary. We have strongly inferred four events of post-Messinian colonization following long-distance dispersal from northern Africa to the Iberian Peninsula, Sicily and Greece. CONCLUSIONS/SIGNIFICANCE: The current distribution of Linaria sect. Versicolores lineages is explained by both ancient isolation between African and European populations and recent events of long-distance dispersal over sea barriers. This result provides new evidence for the biogeographic complexity of the Mediterranean region.