Traces of ancient range shifts in a mountain plant group (Androsace halleri complex, Primulaceae)

Is the incidence of survival in interior Pleistocene refugia (nunataks) underestimated? Phylogeography of the high mountain plant Androsace alpina (Primulaceae) in the European Alps revisited

Temperate mountain ranges such as the European Alps have been strongly affected by the Pleistocene glaciations. Glacial advances forced biota into refugia, which were situated either at the periphery of mountain ranges or in their interior. Whereas in the Alps peripheral refugia have been repeatedly and congruently identified, support for the latter scenario, termed "nunatak hypothesis," is still limited and no general pattern is recognizable yet. Here, we test the hypothesis of nunatak survival for species growing in the high alpine to subnival zones on siliceous substrate using the cushion plant Androsace alpina (Primulaceae), endemic to the European Alps, as our model species. To this end, we analyzed AFLP and plastid DNA sequence data obtained from a dense and range-wide sampling. Both AFLPs and plastid sequence data identified the southwestern-most population as the most divergent one. AFLP data did not allow for discrimination of interior and peripheral populations, but rather identified two to three longitudinally separated major gene pools. In contrast, in the eastern half of the Alps several plastid haplotypes of regional or local distribution in interior ranges-the Alpine periphery mostly harbored a widespread haplotype-were indicative for the presence of interior refugia. Together with evidence from other Alpine plant species, this study shows that in the eastern Alps silicicolous species of open habitats in the alpine and subnival zone survived, also or exclusively so, in interior refugia. As the corresponding genetic structure may be lost in mostly nuclear-derived, rapidly homogenizing marker systems such as AFLPs or RAD sequencing tags, markers not prone to homogenization, as is the case for plastid sequences (Sanger-sequenced or extracted from an NGS data set) will continue to be important for detecting older, biogeographically relevant patterns.

Plant evolution in alkaline magnesium-rich soils: A phylogenetic study of the Mediterranean genus Hormathophylla (Cruciferae: Alysseae) based on nuclear and plastid sequences

Habitats with alkaline edaphic substrates are often associated with plant speciation and diversification. The tribe Alysseae, in the family Brassicaceae, epitomizes this evolutionary trend. In this lineage, some genera, like Hormathophylla, can serve as a good case for testing the evolutionary framework. This genus is centered in the western Mediterranean. It grows on different substrates, but mostly on alkaline soils. It has been suggested that diversification in many lineages of the tribe Alysseae and in the genus Hormathophylla is linked to a tolerance for high levels of Mg+2 in xeric environments. In this study, we investigated the controversial phylogenetic placement of Hormathophylla in the tribe, the generic limits and the evolutionary relationships between the species using ribosomal and plastid DNA sequences. We also examined the putative association between the evolution of different ploidy levels, trichome morphology and the type of substrates. Our analyses demonstrated the monophyly of the genus Hormathophylla including all previously described species. Nuclear sequences revealed two lineages that differ in basic chromosome numbers (x = 7 and x = 8 or derived 11, 15) and in their trichome morphology. Contrasting results with plastid genes indicates more complex relationships between these two lineages involving recent hybridization processes. We also found an association between chloroplast haplotypes and substrate, especially in populations growing on dolomites. Finally, our dated phylogeny demonstrates that the origin of the genus took place in the mid-Miocene, during the establishment of temporal land bridges between the Tethys and Paratethys seas, with a later diversification during the upper Pliocene.

Evolution, ecology and systematics of Soldanella (Primulaceae) in the southern Apennines (Italy)

BACKGROUND

The populations of Soldanella (Primulaceae) of the southern Apennines (Italy) are unique within the genus for their distribution and ecology. Their highly fragmented distribution range, with three main metapopulations on some of the highest mountains (Gelbison, Sila and Aspromonte massifs) of the area, poses intriguing questions about their evolutionary history and biogeography, and about the possibility of local endemisms.

AIMS AND METHODS

In order to clarify the phylogeny and biogeography of the three metapopulations of Soldanella in the southern Apennines, attributed to S. calabrella to date, and to identify possible local endemisms, a comparative approach based on the study of molecular, morphological and ecological characteristics of the populations was employed. Specifically, one nuclear (total ITS) and two plastid (rbcL and trnL) markers were used for the phylogenetic analyses, performed through both maximum likelihood and Bayesian techniques. Among the morphological features, the glandular hair and leaf biometric traits were analysed, and the environment in which the populations grew was characterised for altitude, forest canopy composition and soil pH, C, N and organic matter.

RESULTS AND CONCLUSIONS

Our findings demonstrate that the lineage of Soldanella of southern Italy diverged from the Carpathians lineage during the Middle Pleistocene, and underwent an evolutionary radiation during the Late Pleistocene. The populations of the Sila and Aspromonte massifs diverged from the populations of the Gelbison massif around 380000 years ago and are probably undergoing a progressive differentiation due to their isolation. The populations on the Gelbison massif, moreover, have different morphological features from those of the Sila and Aspromonte massifs and a different ecological niche. The molecular, morphological and ecological data clearly demonstrate that the metapopulation of Soldanella on the Gelbison massif belongs to a new taxonomic unit at the species level, which we name Soldanella sacra A. & L. Bellino from the name of the massif on which it was discovered, the "Holy Mountain".

Automation and Evaluation of the SOWH Test with SOWHAT

The Swofford–Olsen–Waddell–Hillis (SOWH) test evaluates statistical support for incongruent phylogenetic topologies. It is commonly applied to determine if the maximum likelihood tree in a phylogenetic analysis is significantly different than an alternative hypothesis. The SOWH test compares the observed difference in log-likelihood between two topologies to a null distribution of differences in log-likelihood generated by parametric resampling. The test is a well-established phylogenetic method for topology testing, but it is sensitive to model misspecification, it is computationally burdensome to perform, and its implementation requires the investigator to make several decisions that each have the potential to affect the outcome of the test. We analyzed the effects of multiple factors using seven data sets to which the SOWH test was previously applied. These factors include a number of sample replicates, likelihood software, the introduction of gaps to simulated data, the use of distinct models of evolution for data simulation and likelihood inference, and a suggested test correction wherein an unresolved “zero-constrained” tree is used to simulate sequence data. To facilitate these analyses and future applications of the SOWH test, we wrote SOWHAT, a program that automates the SOWH test. We find that inadequate bootstrap sampling can change the outcome of the SOWH test. The results also show that using a zero-constrained tree for data simulation can result in a wider null distribution and higher p-values, but does not change the outcome of the SOWH test for most of the data sets tested here. These results will help others implement and evaluate the SOWH test and allow us to provide recommendations for future applications of the SOWH test. SOWHAT is available for download from https://github.com/josephryan/SOWHAT.

Genetic divergence and phylogeographic history of two closely related species (Leucomeris decora and Nouelia insignis) across the 'Tanaka Line' in Southwest China

BACKGROUND

Leucomeris decora and Nouelia insignis (Asteraceae) are narrowly and allopatrically distributed species, separated by the important biogeographic boundary Tanaka Line in Southwest China. Previous morphological, cytogenetic and molecular studies suggested that L. decora is sister to N. insignis. However, it is less clear how the two species diverged, whether in full isolation or occurring gene flow across the Tanaka Line. Here, we performed a molecular study at the population level to characterize genetic differentiation and decipher phylogeographic history in two closely related species based on variation examined in plastid and nuclear DNAs using a coalescent-based approach.

RESULTS

These morphologically distinct species share plastid DNA (cpDNA) haplotypes. In contrast, Bayesian analysis of nuclear DNA (nDNA) uncovered two distinct clusters corresponding to L. decora and N. insignis. Based on the IMa analysis, no strong indication of migration was detected based on both cpDNA and nDNA sequences. The molecular data pointed to a major west-east split in nuclear DNA between the two species corresponding with the Tanaka Line. The coalescent time estimate for all cpDNA haplotypes dated to the Mid-Late Pleistocene. The estimated demographic parameters showed that the population size of L. decora was similar to that of N. insignis and both experienced limited demographic fluctuations recently.

CONCLUSIONS

The study revealed comprehensive species divergence and phylogeographic histories of N. insignis and L. decora divided by the Tanaka Line. The phylogeographic pattern inferred from cpDNA reflected ancestrally shared polymorphisms without post-divergence gene flow between species. The marked genealogical lineage divergence in nDNA provided some indication of Tanaka Line for its role as a barrier to plant dispersal, and lent support to its importance in promoting strong population structure and allopatric divergence.

The intraspecific genetic variability of siliceous and calcareous Gentiana species is shaped by contrasting demographic and re-colonization processes

The Ciminalis section of Gentiana comprises seven species, two of them growing on siliceous substrates (G. alpina and G. acaulis), the other ones being calcareous taxa (G. clusii, G. angustifolia, G. ligustica, G. occidentalis and G. dinarica). A total of 515 individuals from 183 populations over the entire Ciminalis distribution range was analyzed using four chloroplast loci (trnH-psbA, matK, rpoB and rpoC1) and the nuclear ribosomal marker ITS2. The siliceous species display only two chloroplast haplotypes each and are both characterized by patterns of range expansions all over the Alps. Conversely, the calcareous species are on average more diverse (two to 13 haplotypes per species) with strong patterns of local structuring. We suggest that the occurrence of many calcareous refugia at the periphery of the Alps must have led to local adaptation and morphological diversification, and helped preserving intraspecific diversities during the last glaciations for the associated taxa. ITS2 was more efficient in delineating species boundaries than the chloroplast markers for which several haplotypes are shared among species. This might be either due to chloroplast capture among species and/or to recent divergence. Species adapted to the same substrate are generally only distantly related when they co-occur in the same place. For both types of markers, G. clusii is found genetically distant from all other species.

The evolutionary history of Antirrhinum suggests that ancestral phenotype combinations survived repeated hybridizations

The model species Antirrhinum majus (the garden snapdragon) has over 20 close wild relatives that are morphologically diverse and adapted to different Mediterranean environments. Hybrids between Antirrhinum species have been used successfully to identify genes underlying their phenotypic differences, and to infer how selection acts on them. To better understand the genetic basis for this diversity, we have examined the evolutionary relationships between Antirrhinum species and how these relate to geography and patterns of phenotypic variation in the genus as a whole. Large population samples and both plastid and multilocus nuclear genotypes resolved the relationships between many species and provided some support for the traditional taxonomic division of the genus into morphological subsections. Morphometric analysis of plants grown in controlled conditions supported the phenotypic distinction of the two largest subsections, and the involvement of multiple underlying genes. Incongruence between nuclear and plastid genotypes further suggested that several species have arisen after hybridization between subsections, and that some species continue to hybridize. However, all potential hybrids appear to have retained a phenotype similar to one of their ancestors, suggesting that ancestral combinations of characters are maintained by selection at many different loci.

Disjunct distribution of chloroplast DNA haplotypes in the understory perennial Veratrum album ssp. oxysepalum (Melanthiaceae) in Japan as a result of ancient introgression

• The Quaternary climatic changes resulted in range shifts of species, providing chances for hybridization. However, the genetic signatures of such ancient introgression have rarely been reported. To investigate such signatures, we performed a phylogeographical study on the perennial plant Veratrum album ssp. oxysepalum, which may have hybridized long ago with another congeneric species, V. stamineum. • Sequence variations in chloroplast DNA (cpDNA) were examined in 43 populations in Japan and adjacent areas. Phylogenetic analyses of different cpDNA haplotypes were conducted on the basis of cpDNA and nuclear ribosomal internal transcribed spacer (nrITS) variations. • In the Japanese archipelago, two major groups of haplotypes were detected, one of which was distributed in a disjunct pattern. The major haplotype, occupying the central part of the species' distribution, formed a monophyletic group with V. stamineum in phylogenetic trees on the basis of cpDNA variation, although the two species did not form a monophyletic group in phylogenetic trees on the basis of nrITS variation. • Historical hybridization between V. album ssp. oxysepalum and V. stamineum in refugia during the Quaternary climatic oscillations, and the resulting chloroplast capture of V. stamineum by V. album ssp. oxysepalum, are most probably responsible for the disjunct distribution of cpDNA in V. album ssp. oxysepalum.

Bayesian hypothesis testing supports long-distance Pleistocene migrations in a European high mountain plant (Androsace vitaliana, Primulaceae)

Colonization of the south-western European mountain ranges is suggested to have predominantly progressed from the Iberian Peninsula eastwards, but this hypothesis has never been tested in a statistical framework. Here, we test this hypothesis using Androsace vitaliana, a high elevation species with eight mostly allopatric subspecies, which is widely but disjunctly distributed across all major south-western European mountain ranges. To this end, we use plastid and nuclear sequence data as well as fingerprint (amplified fragment length polymorphisms) data and employ Bayesian methods, which allow co-estimation of genealogy and divergence times using explicit demographic models, as well as hypothesis testing via Bayes factors. Irrespective of the ambiguity concerning where A. vitaliana started to diversify -- both the Alps and the mountain ranges of the Iberian Peninsula outside the Pyrenees were possible -- colonization routes were not simply unidirectional, but involved Pleistocene connections between the Alps and mountain ranges of the Iberian Peninsula bypassing the interjacent Pyrenees via long-distance dispersal. In contrast, the species' post-glacial history is shaped by regional gene pool homogenization resulting in the genetic pattern showing good congruence with geographical proximity in agreement with a vicariance model, but only partly supporting current taxonomy.

Biogeographical and evolutionary importance of the European high mountain systems

Europe is characterised by several high mountain systems dominating major parts of its area, and these structures have strongly influenced the evolution of taxa. For species now restricted to these high mountain systems, characteristic biogeographical patterns of differentiation exist. (i) Many local endemics are found in most of the European high mountain systems especially in the Alps and the more geographically peripheral regions of Europe. Populations isolated in these peripheral mountain ranges often have strongly differentiated endemic genetic lineages, which survived and evolved in the vicinity of these mountain areas over long time periods. (ii) Populations of taxa with wide distributions in the Alps often have two or more genetic lineages, which in some cases even have the status of cryptic species. In many cases, these lineages are the results of several centres of glacial survival in the perialpine areas. Similar patterns also apply to the other geographically extended European high mountain systems, especially the Pyrenees and Carpathians. (iii) Populations from adjoining high mountain systems often show similar genetic lineages, a phenomenon best explained by postglacial retreat to these mountains from one single differentiation centre between them. (iv) The populations of a number of species show gradients of genetic diversity from a genetically richer East to a poorer West. This might indicate better glacial survival conditions for this biogeographical group of species in the more eastern parts of Europe.

Reciprocal Pleistocene origin and postglacial range formation of an allopolyploid and its sympatric ancestors (Androsace adfinis group, Primulaceae)

The biogeographic history of polyploids and their lower-ploid ancestors is an important feature to achieve a better understanding of polyploid evolution. This is exemplified here using the ecologically congruent members of the Androsace adfinis group (Primulaceae) endemic to the southwestern European Alps. Employing relative genome size, AFLP fingerprint and chloroplast sequence haplotype data, we show that Androsace brigantiaca is a recent (probably no more than 0.2 million years) allopolyploid derivative of the geographically close A. adfinis and A. puberula, which formed reciprocally in a comparatively restricted area in the southern Southwestern Alps. Bayesian admixture analysis--also of artificial additive AFLP profiles--shows that the nuclear genome of A. brigantiaca is significantly biased towards the puberula-genome irrespective of maternal parentage. Nevertheless, there is no evidence for genetic interaction (hybridization, introgression) of A. brigantiaca with either of its ancestors, including the widely sympatric A. puberula. Sympatry might be facilitated by ecological displacement on a local scale or might be a transitory phase on the way to competitive replacement via, for instance, polyploid superiority.