Concordant genetic breaks, identified by combining clustering and tessellation methods, in two co-distributed alpine plant species

Spatial and temporal patterns in the population genomics of the European cockchafer Melolontha melolontha in the Alpine region

The European cockchafer Melolontha melolontha is an agricultural pest in many European countries. Populations have a synchronized 3 or 4 years life cycle, leading to temporally isolated populations. Despite the economic importance and availability of comprehensive historical as well as current records on cockchafer occurrence, population genomic analyses of M. melolontha are missing. For example, the effects of geographic separation caused by the mountainous terrain of the Alps and of temporal isolation on the genomic structure of M. melolontha still remain unknown. To address this gap, we genotyped 475 M. melolontha adults collected during 3 years from 35 sites in a central Alpine region. Subsequent population structure analyses discriminated two main genetic clusters, i.e., the South Tyrol cluster including collections located southeast of the Alpine mountain range, and a northwestern alpine cluster with all the other collections, reflecting distinct evolutionary history and geographic barriers. The "passo di Resia" linking South and North Tyrol represented a regional contact zone of the two genetic clusters, highlighting genomic differentiation between the collections from the northern and southern regions. Although the collections from northwestern Italy were assigned to the northwestern alpine genetic cluster, they displayed evidence of admixture with the South Tyrolean genetic cluster, suggesting shared ancestry. A linear mixed model confirmed that both geographic distance and, to a lower extent, also temporal isolation had a significant effect on the genetic distance among M. melolontha populations. These effects may be attributed to limited dispersal capacity and reproductive isolation resulting from synchronized and non-synchronized swarming flights, respectively. This study contributes to the understanding of the phylogeography of an organism that is recognized as an agricultural problem and provides significant information on the population genomics of insect species with prolonged temporally shifted and locally synchronized life cycles.

Repeated hybridization increased diversity in the door snail complex Charpentieria itala in the Southern Alps

The door snail species complex Charpentieria itala is widely distributed in the Southern Alps and subdivided into several morphologically differentiated subspecies. Thus, it can be used as a model group for understanding migration and differentiation processes in the Southern Alps. We generated genome-wide double digest Restriction Site Associated DNA (ddRAD) sequencing data for 166 specimens from 36 populations of the door snail Charpentieria itala and for 8 specimens of the other three Charpentieria species to reconstruct their evolutionary history and phylogeography. Phylogenetic and structure analyses based on the ddRAD data indicated that the repeated separation of the populations in western and eastern groups by the Garda glacier during the glacials was the process that most strongly shaped the population structure of C. itala. This process may also explain a similar phylogeographic boundary in many other southern Alpine animal and plant species. Our study revealed that some populations that resemble Charpentieria stenzii morphologically and ecologically, the 'stenzioid' subspecies, originated by a hybridization event with Charpentieria stenzii. A further hybridization event between stenzioid populations that survived the glacials in mountain refuges and non-stenzioid populations that probably came into contact with stenzioid populations as a result of climate warming during an interglacial resulted in the origin of a hybrid subspecies that is adapted to intermediate altitudes. Our study demonstrated that the origin of new differentiated taxa by hybridization, is more frequent than previously assumed.

Diversification and independent polyploid origins in the disjunct species Alyssum repens from the Southeastern Alps and the Carpathians

PREMISE

Disjunct distributions have been commonly observed in mountain plant species and have stimulated phylogeographic and phylogenetic research. Here we studied Alyssum repens, a member of the polyploid species complex A. montanum-A. repens, which exhibits SE Alpine-Carpathian disjunctions with a large elevational span and consists of diploid and tetraploid populations. We aimed to investigate the species' genetic and cytotype structure in the context of its distribution patterns, to elucidate the polyploid origins and to propose an appropriate taxonomic treatment.

METHODS

We combined AFLP fingerprinting markers, sequence variation of the highly repetitive ITS region of rDNA and the low-copy DET1 nuclear gene, genome size, and morphometric data.

RESULTS

We identified four geographically structured genetic lineages. One consisted of diploid populations from the foothills of the Southeastern Alps and neighboring regions, and the three others were allopatric montane to alpine groups comprising diploids and tetraploids growing in the Southeastern Carpathians and the Apuseni Mts. in Romania.

CONCLUSIONS

We inferred a vicariance scenario associated with Quaternary climatic oscillations, accompanied by one auto- and two allopolyploidization events most likely involving a northern Balkan relative. Whereas genetic differentiation and allopatric distribution would favor the taxonomic splitting of this species, the genetic lineages largely lack morphological distinguishability, and their ecological, cytotype and genome size divergence is only partial. Even though we probably face here a case of incipient speciation, we propose to maintain the current taxonomic treatment of Alyssum repens as a single, albeit variable, species.

Three in One--Multiple Faunal Elements within an Endangered European Butterfly Species

Ice ages within Europe forced many species to retreat to refugia, of which three major biogeographic basic types can be distinguished: "Mediterranean", "Continental" and "Alpine / Arctic" species. However, this classification often fails to explain the complex phylogeography of European species with a wide range of latitudinal and altitudinal distribution. Hence, we tested for the possibility that all three mentioned faunal elements are represented within one species. Our data was obtained by scoring 1,307 Euphydryas aurinia individuals (46 European locations) for 17 allozyme loci, and sequencing a subset of 492 individuals (21 sites) for a 626 base pairs COI fragment. Genetic diversity indices, F statistics, hierarchical analyses of molecular variance, individual-based clustering, and networks were used to explore the phylogeographic patterns. The COI fragment represented 18 haplotypes showing a strong geographic structure. All but one allozyme loci analysed were polymorphic with a mean F_ST of 0.20, supporting a pronounced among population structure. Interpretation of both genetic marker systems, using several analytical tools, calls for the recognition of twelve genetic groups. These analyses consistently distinguished different groups in Iberia (2), Italy, Provence, Alps (3), Slovenia, Carpathian Basin, the lowlands of West and Central Europe as well as Estonia, often with considerable additional substructures. The genetic data strongly support the hypothesis that E. aurinia survived the last glaciation in Mediterranean, extra-Mediterranean and perialpine refugia. It is thus a rare example of a model organism that combines attributes of faunal elements from all three of these sources. The observed differences between allozymes and mtDNA most likely result from recent introgression of mtDNA into nuclear allozyme groups. Our results indicate discrepancies with the morphologically-based subspecies models, underlining the need to revise the current taxonomy.

Environmental factors contribute to the formation and maintenance of the contact zone observed in deciduous broad-leaved tree species in Japan

Contact zones are defined as areas where populations from different refugia meet during a postglacial expansion and distinct DNA lineages are mixedly distributed. In Japan, contact zones of various plants and animals were reported from the Kinki-Chugoku region. These contact zones appear to be maintained without any drastic topographic barriers such as those observed in the Alps and Pyrenees Mountains. In this study, the mechanisms underlying the formation and/or maintenance of these contact zones were investigated using six deciduous broad-leaved tree species (Carpinus laxiflora, C. tschonoskii, C. japonica, Magnolia obovata, Padus grayana, and Euonymus oxyphyllus). First, the precise location of the contact zones was examined by intensive genetic analysis of the six species. Second, the relationships between the geographic location of the contact zone and various environmental factors, including climate and topography, were investigated by generalized additive models to reveal the mechanisms of the formation and maintenance of the contact zones. As a result, four of the six examined plant species clearly showed a geographically common contact zone in Hyogo Prefecture and its adjacent areas. The results of the generalized additive models indicate that the pattern of low habitat suitability estimated by ecological niche modeling was the most important factor for determining the location of the common contact zone. These results suggest that areas with low habitat suitability in Hyogo Prefecture restrict the migration and gene flow of the four species in this region, and thus, they maintain the pattern of the contact zones. This study suggests that there are major effects of habitat suitability on the formation and maintenance of the contact zones.

Genetic and ecological data reveal species boundaries between viviparous and oviparous lizard lineages

Identification of cryptic species is an essential aim for conservation biologists to avoid premature extinctions of 'unrecognized' species. Integrating different types of data can undoubtedly aid in resolving the issue of species delimitation. We studied here two lineages of the common lizard Zootoca vivipara that display different reproductive mode (the viviparous Z. v. vivipara and the oviparous Z. v. carniolica) and that overlap their distributional ranges in the European Alps. With the purpose of delimiting species' boundaries, we analyzed their ecological, genetic and natural history features. More than 300 samples were collected and analyzed at cytochrome b and 11 microsatellites loci for investigating genetic variation, population structure, individual relatedness and evolutionary histories of the two lineages. Additionally, we compared their ecological niches using eight ecological variables. Genetic data showed contrasting patterns of genetic structure between the two lineages, different demographic dynamics and no hybridization events. Also strong ecological differences (such as temperature) emerged between the two lineages, and niche overlap was limited. Taken together, these results indicate that Z. v. vivipara and Z. v. carniolica should be recognized as two separate species, and particular conservation consideration should be given to the oviparous lineage that tends to live in areas threatened by increasing impact of human activities. However, recent and rapid climate warming might determine an increasing risk for the persistence of the viviparous lineage, being adapted to cold environments.

Mitochondrial and nuclear DNA survey of Zootoca vivipara across the eastern Italian Alps: evolutionary relationships, historical demography and conservation implications

The European common lizard Zootoca vivipara exhibits reproductive bimodality, with populations being either viviparous or oviparous. In the central-eastern Italian Alps oviparous populations (Z. v. carniolica) and viviparous populations (Z. v. vivipara) partly overlap geographically. Studying the evolutionary relationship between these taxa presents an interesting opportunity to gain insight into the evolution of this trait. We aim to: i) test whether Z. v. carniolica, which is endangered, constitutes an ESU (Evolutionary Significant Unity); ii) infer mtDNA divergence time between the Z. v. carniolica clade and all the other Z. vivipara subspecies with the aid of an external calibration point; and iii) describe the phylogeographical and demographic scenarios in the area. To do so we sequenced about 200 individuals for mitochondrial variation; 64 of them were also analysed for three nuclear genes. Furthermore, we analysed the same nuclear markers in 17 individuals from the other oviparous subspecies Z. v. louislantzi and 11 individuals of Z. v. vivipara from widespread geographical origins. The mtDNA and nDNA loci that we examined supported the monophyly of Z. v. carniolica. The mtDNA-based estimate of divergence time between Z. v. carniolica and all the other subspecies indicated a separation at 4.5 Mya (95% CI 6.1–2.6), with about 5% of sequence divergence. Considering that Z. v. carniolica harbours higher genetic diversity, while Z. v. vivipara from central-eastern Alps shows a signature of recent population and spatial expansion, we argue that Z. v. carniolica represents a distinct evolutionary unit, with a presumably long-term evolutionary history of separation. Z. v. carniolica populations, occurring at higher latitudes and altitudes than insofar supposed, live in peat bogs, a seriously threatened habitat: taking into account also its evolutionary distinctness, specific conservation measures should be considered.

Glacial history affected phenotypic differentiation in the alpine plant, Campanula thyrsoides

Numerous widespread Alpine plant species show molecular differentiation among populations from distinct regions. This has been explained as the result of genetic drift during glacial survival in isolated refugia along the border of the European Alps. Since genetic drift may affect molecular markers and phenotypic traits alike, we asked whether phenotypic differentiation mirrors molecular patterns among Alpine plant populations from different regions. Phenotypic traits can be under selection, so we additionally investigated whether part of the phenotypic differentiation can be explained by past selection and/or current adaptation. Using the monocarpic Campanula thyrsoides as our study species, a common garden experiment with plants from 21 populations from four phylogeographic groups located in regions across the Alps and the Jura Mountains was performed to test for differentiation in morphological and phenological traits. Past selection was investigated by comparing phenotypic differentiation among and within regions with molecular differentiation among and within regions. The common garden results indicated regional differentiation among populations for all investigated phenotypic traits, particularly in phenology. Delayed flowering in plants from the South-eastern Alps suggested adaptation to long sub-mediterranean summers and contrasted with earlier flowering of plants experiencing shorter growing seasons in regions with higher elevation to the West. Comparisons between molecular and phenotypic differentiation revealed diversifying selection among regions in height and biomass, which is consistent with adaptation to environmental conditions in glacial refugia. Within regions, past selection acted against strong diversification for most phenotypic traits, causing restricted postglacial adaptation. Evidence consistent with post-glacial adaptation was also given by negative correlation coefficients between several phenotypic traits and elevation of the population's origin. In conclusion, our study suggests that, irrespective of adaptation of plants to their current environment, glacial history can have a strong and long-lasting influence on the phenotypic evolution of Alpine plants.

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.

Nuclear and chloroplast DNA phylogeography reveals vicariance among European populations of the model species for the study of metal tolerance, Arabidopsis halleri (Brassicaceae)

Arabidopsis halleri is a pseudometallophyte involved in numerous molecular studies of the adaptation to anthropogenic metal stress. In order to test the representativeness of genetic accessions commonly used in these studies, we investigated the A. halleri population genetic structure in Europe. Microsatellite and nucleotide polymorphisms from the nuclear and chloroplast genomes, respectively, were used to genotype 65 populations scattered over Europe. The large-scale population structure was characterized by a significant phylogeographic signal between two major genetic units. The localization of the phylogeographic break was assumed to result from vicariance between large populations isolated in southern and central Europe, on either side of ice sheets covering the Alps during the Quaternary ice ages. Genetic isolation was shown to be maintained in western Europe by the high summits of the Alps, whereas admixture was detected in the Carpathians. Considering the phylogeographic literature, our results suggest a distinct phylogeographic pattern for European species occurring in both mountain and lowland habitats. Considering the evolution of metal adaptation in A. halleri, it appears that recent adaptations to anthropogenic metal stress that have occurred within either phylogeographic unit should be regarded as independent events that potentially have involved the evolution of a variety of genetic mechanisms.

Extensive range persistence in peripheral and interior refugia characterizes Pleistocene range dynamics in a widespread Alpine plant species (Senecio carniolicus, Asteraceae)

Recent evidence suggests that survival of arctic-alpine organisms in peripheral or interior glacial refugia are not mutually exclusive and may both be involved in shaping an organism’s Pleistocene history, yet potentially at different time levels. Here, we test this hypothesis in a high-mountain plant (diploid lineage of Senecio carniolicus, Asteraceae) from the Eastern European Alps, in which patterns of morphological variation and current habitat requirements suggest survival in both types of refugia. To this end, we used AFLPs, nuclear and plastid DNA sequences and analysed them, among others, within a graph theoretic framework and using novel Bayesian methods of phylogeographic inference. On the basis of patterns of genetic diversity, occurrence of rare markers, distribution of distinct genetic lineages and patterns of range connectivity both interior refugia in the formerly strongly glaciated central Alps and peripheral refugia along the southern margin of the Alps were identified. The presence of refugia congruently inferred by markers resolving at different time levels suggests that these refugia acted as such throughout several glacial cycles. The high degree of range persistence together with gradual range expansion, which contrasts with the extent of range shifts implied for other Alpine species, is likely responsible for incipient lineage differentiation evident from the genetic data. Replacing a simplistic peripheral vs. interior refugia dualism by more complex models involving both types of refugia and considering different time levels will help identifying common phylogeographic patterns with respect to, for instance, location of refugia and colonization routes and elucidating their underlying genetic and/or ecological causes.

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.

Barriers to gene flow between emerging populations of Watermelon mosaic virus in Southeastern France

Since 1999, "emerging" (EM) strains of Watermelon mosaic virus (WMV) have been detected in cucurbit crops of southeastern France, probably as a result of recent introductions. Population genetic approaches were used to study the structure of EM isolates in southeastern France and to identify factors involved in their spatial distribution. A population clustering method (SAMOVA) and a maximum-difference algorithm (Monmonier's algorithm) were combined to visualize and quantify barriers to gene flow between populations. Both methods yielded similar results and two main barriers were identified. A North/South oriented barrier may be related to physical obstacles to gene flow (Rhône River, presence of an area with few cucurbit crops). Although the barrier was very strong, some "crossing" events were detected. A second barrier, oriented Northwest to Southeast, was not correlated with obvious geographical features. The two methods used here are complementary and confirm the limited spread of WMV-EM isolates. This approach can be useful in epidemiology studies to characterize the structure of viral populations and identify barriers to gene flow.

The Strait of Gibraltar as a major biogeographic barrier in Mediterranean conifers: a comparative phylogeographic survey

The Strait of Gibraltar (SG) is reputed for being both a bridge and a geographic barrier to biological exchanges between Europe and Africa. Major genetic breaks associated with this strait have been identified in various taxa, but it is unknown whether these disjunctions have been produced simultaneously or by independent biogeographic processes. Here, the genetic structure of five conifers distributed on both sides of the SG was investigated using mitochondrial (nad1 b/c, nad5-1, nad5-4 and nad7-1) and chloroplast (Pt1254, Pt15169, Pt30204, Pt36480, Pt71936 and Pt87268) DNA markers. The distribution of genetic variation was partially congruent between types of markers within the same species. Across taxa, there was a significant overlapping between the SG and the genetic breaks detected, especially for the four Tertiary species surveyed (Abies pinsapo complex, Pinus nigra, Pinus pinaster and Taxus baccata). For most of these taxa, the divergence of populations across the SG could date back to long before the Pleistocene glaciations. However, their strongly different cpDNA G(ST) and R(ST) values point out that they have had dissimilar population histories, which might include contrasting amounts of pollen-driven gene flow since their initial establishment in the region. The fifth species, Pinus halepensis, was genetically depauperated and homogenous on both sides of the SG. A further analysis of nuclear DNA sequences with coalescent-based isolation with migration models suggests a Pleistocene divergence of P. halepensis populations across the SG, which is in sharp contrast with the pre-Pleistocene divergence dates obtained for P. pinaster. Altogether, these results indicate that the genetic breaks observed across this putative biogeographical barrier have been produced by independent evolutionary processes related to the biological history of each individual species instead of a common vicariant phenomenon.