The aim of the present study was to evaluate the current distribution of suitable niches for the invasive orchid species, Epipactis helleborine, and to estimate the possibility of its further expansion. Moreover, niche modeling tools were used to explain its rapid expansion in North America and to test the niche conservatism of the species. The maximum entropy method was used to create models of the suitable niche distribution. A database of E. helleborine localities was prepared based on the examination of herbarium specimens, information from electronic databases as well as data gathered during field works. The differences between the niches occupied by native and invasive populations were evaluated using the niche overlap and niche identity test indexes. Moreover, the coverage of the most suitable habitats for the species was measured for three future scenarios as well as for the present time model. Populations of E. helleborine occupy North American west coast habitats very similar to those preferred by native, Eurasian populations, while the expansion in the east coast is related to the niche shift. The created models of suitable niche distribution indicate that the species does not realize its potential niche in the native range. The total surface of the habitats potentially available for E. helleborine will decrease in all climate change scenarios created for 2080.
Pleistocene climatic oscillations have played a major role in structuring present-day biodiversity. The southern Mediterranean peninsulas have long been recognized as major glacial refugia, from where Northern Europe was post-glacially colonized. However, recent studies have unravelled numerous additional refugia also in northern regions. We investigated the phylogeographic pattern of the widespread Western Palaearctic lizard Podarcis muralis, using a range-wide multilocus approach, to evaluate whether it is concordant with a recent expansion from southern glacial refugia or alternatively from a combination of Mediterranean and northern refugia.
We analyzed DNA sequences of two mitochondrial (cytb and nd4) and three nuclear (acm4, mc1r, and pdc) gene fragments in individuals from 52 localities across the species range, using phylogenetic and phylogeographic methods. The complex phylogeographic pattern observed, with 23 reciprocally monophyletic allo- parapatric lineages having a Pleistocene divergence, suggests a scenario of long-term isolation in multiple ice-age refugia across the species distribution range. Multiple lineages were identified within the three Mediterranean peninsulas – Iberia, Italy and the Balkans - where the highest genetic diversity was observed. Such an unprecedented phylogeographic pattern - here called “refugia within all refugia” – compasses the classical scenario of multiple southern refugia. However, unlike the southern refugia model, various distinct lineages were also found in northern regions, suggesting that additional refugia in France, Northern Italy, Eastern Alps and Central Balkans allowed the long-term persistence of this species throughout Pleistocene glaciations.
The phylogeography of Podarcis muralis provides a paradigm of temperate species survival in Mediterranean and extra-Mediterranean glacial refugia. Such refugia acted as independent biogeographic compartments for the long-term persistence of this species, for the differentiation of its genetic lineages, and for the short-distance post-glacial re-colonization of neighbouring areas. This finding echoes previous findings from recent phylogeographic studies on species from temperate ecoregions, thus suggesting the need for a reappraisal of the role of northern refugia for glacial persistence and post-glacial assembly of Holarctic biota.
Podarcis muralis; Phylogeography; Western Palaearctic; Glacial refugia; Mediterranean peninsulas; Genetic diversity; Temperate species
The Alpine Region, constituting the Alps and the Dinaric Alps, has played a major role in the formation of current patterns of biodiversity either as a contact zone of postglacial expanding lineages or as the origin of genetic diversity. In our study, we tested these hypotheses for two widespread, sympatric microgastropod taxa – Carychium minimum O.F. Müller, 1774 and Carychium tridentatum (Risso, 1826) (Gastropoda, Eupulmonata, Carychiidae) – by using COI sequence data and species potential distribution models analyzed in a statistical phylogeographical framework. Additionally, we examined disjunct transatlantic populations of those taxa from the Azores and North America. In general, both Carychium taxa demonstrate a genetic structure composed of several differentiated haplotype lineages most likely resulting from allopatric diversification in isolated refugial areas during the Pleistocene glacial periods. However, the genetic structure of Carychium minimum is more pronounced, which can be attributed to ecological constraints relating to habitat proximity to permanent bodies of water. For most of the Carychium lineages, the broader Alpine Region was identified as the likely origin of genetic diversity. Several lineages are endemic to the broader Alpine Region whereas a single lineage per species underwent a postglacial expansion to (re)colonize previously unsuitable habitats, e.g. in Northern Europe. The source populations of those expanding lineages can be traced back to the Eastern and Western Alps. Consequently, we identify the Alpine Region as a significant ‘hot-spot’ for the formation of genetic diversity within European Carychium lineages. Passive dispersal via anthropogenic means best explains the presence of transatlantic European Carychium populations on the Azores and in North America. We conclude that passive (anthropogenic) transport could mislead the interpretation of observed phylogeographical patterns in general.
The Iberian Peninsula is recognized as an important refugial area for species survival and diversification during the climatic cycles of the Quaternary. Recent phylogeographic studies have revealed Iberia as a complex of multiple refugia. However, most of these studies have focused either on species with narrow distributions within the region or species groups that, although widely distributed, generally have a genetic structure that relates to pre-Quaternary cladogenetic events. In this study we undertake a detailed phylogeographic analysis of the lizard species, Lacerta lepida, whose distribution encompasses the entire Iberian Peninsula. We attempt to identify refugial areas, recolonization routes, zones of secondary contact and date demographic events within this species.
Results support the existence of 6 evolutionary lineages (phylogroups) with a strong association between genetic variation and geography, suggesting a history of allopatric divergence in different refugia. Diversification within phylogroups is concordant with the onset of the Pleistocene climatic oscillations. The southern regions of several phylogroups show a high incidence of ancestral alleles in contrast with high incidence of recently derived alleles in northern regions. All phylogroups show signs of recent demographic and spatial expansions. We have further identified several zones of secondary contact, with divergent mitochondrial haplotypes occurring in narrow zones of sympatry.
The concordant patterns of spatial and demographic expansions detected within phylogroups, together with the high incidence of ancestral haplotypes in southern regions of several phylogroups, suggests a pattern of contraction of populations into southern refugia during adverse climatic conditions from which subsequent northern expansions occurred. This study supports the emergent pattern of multiple refugia within Iberia but adds to it by identifying a pattern of refugia coincident with the southern distribution limits of individual evolutionary lineages. These areas are important in terms of long-term species persistence and therefore important areas for conservation.
Species with strict ecological requirements may provide new insights into the forces that shaped the geographic variation of genetic diversity. The Pyrenean desman, Galemys pyrenaicus, is a small semi-aquatic mammal that inhabits clean streams of the northern half of the Iberian Peninsula and is endangered in most of its geographic range, but its genetic structure is currently unknown. While the stringent ecological demands derived from its aquatic habitat might have caused a partition of the genetic diversity among river basins, Pleistocene glaciations would have generated a genetic pattern related to glacial refugia.
To study the relative importance of historical and ecological factors in the genetic structure of G. pyrenaicus, we used mitochondrial and intronic sequences of specimens covering most of the species range. We show, first, that the Pyrenean desman has very low levels of genetic diversity compared to other mammals. In addition, phylogenetic and dating analyses of the mitochondrial sequences reveal a strong phylogeographic structure of a Middle Pleistocene origin, suggesting that the main lineages arose during periods of glacial isolation. Furthermore, both the spatial distribution of nuclear and mitochondrial diversity and the results of species distribution modeling suggest the existence of a major glacial refugium in the northwestern part of the Iberian Peninsula. Finally, the main mitochondrial lineages show a striking parapatric distribution without any apparent exchange of mitochondrial haplotypes between the lineages that came into secondary contact (although with certain permeability to nuclear genes), indicating incomplete mixing after the post-glacial recolonization. On the other hand, when we analyzed the partition of the genetic diversity among river basins, the Pyrenean desman showed a lower than expected genetic differentiation among main rivers.
The analysis of mitochondrial and intronic markers in G. pyrenaicus showed the predominant effects of Pleistocene glaciations on the genetic structure of this species, while the distribution of the genetic diversity was not greatly influenced by the main river systems. These results and, particularly, the discovery of a marked phylogeographic structure, may have important implications for the conservation of the Pyrenean desman.
Conservation genetics; Introns; Mammals; Mitochondrial genes; Nuclear genes; Pyrenean desman; Niche modeling; Iberian Peninsula; Endemism
It is estimated that floral deception has evolved in at least 7500 species of angiosperms, of which two thirds are orchids. Epipactis veratrifolia (Orchidaceae) is a model system of aphid mimicry as aphidophagous hoverflies lay eggs on false brood sites on their flowers. To understand the evolutionary ecology of floral deception, we investigated the pollination biology of E. veratrifolia across 10 populations in the Eastern Himalayas. We reconstructed the phylogeny of Epipactis and mapped the known pollination systems of previously studied species onto the tree.
Some inflorescences of E. veratrifolia were so infested with aphids while they were still in bud that the some larvae of hoverflies developed to the third instar while flower buds opened. This indicated that adult female hoverflies were partly rewarded for oviposition. Although flowers failed to secrete nectar, they mimicked both alarm pheromones and aphid coloring of to attract female hoverflies as their exclusive pollinators. Phylogenetic mapping indicate that pollination by aphidophagous hoverflies is likely an ancestral condition in the genus Epipactis. We suggest that the biological interaction of aphid (prey), orchid (primary producer) and hoverfly (predator) may represent an intermediate stage between mutualism and deception in the evolution of pollination-by-deceit in E. veratrifolia.
Our analyses indicate that this intermediate stage may be used as a model system to interpret the origin of oviposition (brood site) mimicry in Epipactis. We propose the hypothesis that some deceptive pollination systems evolved directly from earlier (partly) mutualistic systems that maintained the fidelity of the original pollinator(s) even though rewards (nectar/ brood site) were lost.
Anther cap; Aphids; Floral mimicry; Hoverflies; Intermediate stage; Pollinator
Population genetic theory holds that oceanic island populations are expected to have lower levels of genetic variation than their mainland counterparts, due to founder effect after island colonization from the continent. Cistus monspeliensis (Cistaceae) is distributed in both the Canary Islands and the Mediterranean region. Numerous phylogenetic results obtained in the last years allow performing further phylogeographic analyses in Cistus.
We analyzed sequences from multiple plastid DNA regions in 47 populations of Cistus monspeliensis from the Canary Islands (21 populations) and the Mediterranean basin (26 populations). The time-calibrated phylogeny and phylogeographic analyses yielded the following results: (1) a single, ancestral haplotype is distributed across the Mediterranean, whereas 10 haplotypes in the Canary Islands; (2) four haplotype lineages are present in the Canarian Islands; (3) multiple colonization events across the archipelago are inferred; (4) the earliest split of intraspecific lineages occurred in the Early to Middle Pleistocene (<930,000 years BP).
The contrasting pattern of cpDNA variation is best explained by genetic bottlenecks in the Mediterranean during Quaternary glaciations, while the Canarian archipelago acted as a refugium of high levels of genetic diversity. Active colonization across the Canarian islands is supported not only by the distribution of C. monspeliensis in five of the seven islands, but also by our phylogeographic reconstruction in which unrelated haplotypes are present on the same island. Widespread distribution of thermophilous habitats on every island, as those found throughout the Mediterranean, has likely been responsible for the successful colonization of C. monspeliensis, despite the absence of a long-distance dispersal mechanism. This is the first example of a plant species with higher genetic variation among oceanic island populations than among those of the continent.
Present day distributions of Palearctic taxa in northern latitudes mainly result from populations having survived in local patches during the Late Pleistocene and/or from recolonizing populations from southern temperate refugia. If well-studied Mediterranean and eastern European refugia are widely accepted, some recent biogeographical assumptions still remain unclear, such as the occurrence of multiple glacial refugia in Iberia and cryptic refugia in northern Europe during the last glaciations. The Lusitanian snail Elona quimperiana has a remarkably disjunct distribution, limited to northwestern France (Brittany), northwestern Spain and the Basque Country. By describing the phylogeographical structure of this species across its entire range, the present study attempts to identify refugia and subsequent recolonization routes.
Results based on 16S and COI gene sequences showed that the low genetic diversity observed in the Brittany populations should be associated with a recent demographic expansion. By contrast, populations from Spain exhibit several differentiated lineages and are characterized by demographic equilibrium, while the Basque populations are the only ones harboring typical distinct haplotypes. The center of the star-like networks of both gene sequences is occupied by a common ancestral-like haplotype found in Brittany and Spain, which might have originated from the middle of Northern Spain (i.e. Asturias, eastern Lugo and western Cantabria). Estimates of the divergence time between the Spain-Brittany and Basque lineages strongly suggest that E. quimperiana survived the Pleistocene glaciations in distinct refugia on the Iberian Peninsula, one of which is situated in Picos de Europa, and the other in the Basque Country. The occurrence of a northern refugium in France cannot be rejected as of yet.
Present results confirm the Iberian origin of the land snail E. quimperiana and strongly support the emerging phylogeographic hypothesis of multiple refugia in Iberia during the last glaciations. The scenario of a spatial expansion of E. quimperiana from an Iberian refuge located in Asturias to northern areas provides the most probable explanation for the present distribution of this land snail. By harboring distinct haplotypes, the Basque Country populations appear to be of great importance in terms of potential adaptation, long term persistence and hence, the conservation of E. quimperiana.
A central aim of island biogeography is to understand the colonization history of insular species using current distributions, fossil records and genetic diversity. Here, we analyze five plastid DNA regions of the endangered Juniperus brevifolia, which is endemic to the Azores archipelago.
The phylogeny of the section Juniperus and the phylogeographic analyses of J. brevifolia based on the coalescence theory of allele (plastid) diversity suggest that: (1) a single introduction event likely occurred from Europe; (2) genetic diversification and inter-island dispersal postdated the emergence of the oldest island (Santa Maria, 8.12 Ma); (3) the genetic differentiation found in populations on the islands with higher age and smaller distance to the continent is significantly higher than that on the younger, more remote ones; (4) the high number of haplotypes observed (16), and the widespread distribution of the most frequent and ancestral ones across the archipelago, are indicating early diversification, demographic expansion, and recurrent dispersal. In contrast, restriction of six of the seven derived haplotypes to single islands is construed as reflecting significant isolation time prior to colonization.
Our phylogeographic reconstruction points to the sequence of island emergence as the key factor to explain the distribution of plastid DNA variation. The reproductive traits of this juniper species (anemophily, ornithochory, multi-seeded cones), together with its broad ecological range, appear to be largely responsible for recurrent inter-island colonization of ancestral haplotypes. In contrast, certain delay in colonization of new haplotypes may reflect intraspecific habitat competition on islands where this juniper was already present.
Phylogeographic patterns of temperate species from the Mediterranean peninsulas have been investigated intensively. Nevertheless, as more phylogeographies become available, either unique patterns or new lines of concordance continue to emerge, providing new insights on the evolution of regional biotas. Here, we investigated the phylogeography and evolutionary history of the Italian crested newt, Triturus carnifex, through phylogenetic, molecular dating and population structure analyses of two mitochondrial gene fragments (ND2 and ND4; overall 1273 bp). We found three main mtDNA lineages having parapatric distribution and estimated divergence times between Late Pliocene and Early Pleistocene. One lineage (S) was widespread south of the northern Apennine chain and was further geographically structured into five sublineages, likely of Middle Pleistocene origin. The second lineage (C) was widespread throughout the Padano–Venetian plain and did not show a clear phylogeographic structure. The third lineage (N) was observed in only two populations located on western Croatia/Slovenia. Results of analysis of molecular variance suggested that partitioning populations according to the geographic distribution of these lineages and sublineages explains 76% of the observed genetic variation. The phylogeographic structure observed within T. carnifex and divergence time estimates among its lineages, suggest that responses to Pleistocene environmental changes in this single species have been as diverse as those found previously among several codistributed temperate species combined. Consistent with the landscape heterogeneity, physiographic features, and palaeogeographical evolution of its distribution range, these responses encompass multiple refugia along the Apennine chain, lowland refugia in large peri-coastal plains, and a ‘cryptic’ northern refugium.
Floral nectar of animal-pollinated plants is commonly infested with microorganisms, yet little is known about the microorganisms inhabiting the floral nectar of orchids. In this study, we investigated microbial communities occurring in the floral nectar of seven Epipactis (Orchidaceae) species. Culturable bacteria and yeasts were isolated and identified by partially sequencing the small subunit (SSU) ribosomal RNA (rRNA) gene and the D1/D2 domains of the large subunit (LSU) rRNA gene, respectively. Using three different culture media, we found that bacteria were common inhabitants of the floral nectar of Epipactis. The most widely distributed bacterial operational taxonomic units (OTUs) in nectar of Epipactis were representatives of the family of Enterobacteriaceae, with an unspecified Enterobacteriaceae bacterium as the most common. In contrast to previous studies investigating microbial communities in floral nectar, very few yeast species (mainly of the genus Cryptococcus) were observed, and most of them occurred in very low densities. Total OTU richness (i.e., the number of bacterial and yeast OTUs per orchid species) varied between 4 and 20. Cluster analysis revealed that microbial communities of allogamous species differed from those of autogamous and facultatively autogamous species. This study extends previous efforts to identify microbial communities in floral nectar and indicates that the floral nectar of the orchids investigated mainly contained bacterial communities with moderate phylogenetic diversity.
Bacteria; floral nectar; microbial communities; orchids; yeasts.
The origin and colonisation history after the Quaternary ice ages remain largely unresolved for many plant lineages, mainly owing to a lack of fine-scale studies. Here, we present a molecular phylogeny and a phylogeographic analysis of Antirrhinum, an important model system in plant biology, in the Pyrenees range. Our goal was to reconstruct the evolutionary and colonisation history of four taxa endemic to this region (A. majus subsp. majus, A. majus. subsp. striatum, A. molle, and A. sempervirens) by using a dense sampling strategy, with a total of 452 individuals from 99 populations whose collective distribution spans nearly the entirety of the Pyrenees and adjacent mountains.
Phylogenetic and phylogeographic analyses of the sequences of two plastid (trnS-trnG and trnK-matK) regions revealed the following: (i) historical relationship between the Pyrenees and Iberia (but not with the Alps); (ii) the long persistence of populations in the Pyrenees, at least since the Late Pleistocene; (iii) three different colonisation histories for populations from the Western, Central, and Eastern Pyrenees; (iv) the deep phylogeographic separation of the eastern and western populations; and (v) the colonisation of southern France from the Eastern Pyrenees.
The present study underlines the enormous influence of the glacial history of the mountain ranges on the current configuration of intra- and inter-specific genetic diversity in Antirrhinum, as well as the importance of periglacial areas for the survival of species during glacial periods of the Quaternary.
Antirrhinum; Phylogeny; Phylogeography; Pyrenees; Quaternary
Accurate identification of Arctic plant species is critical for understanding potential climate-induced changes in their diversity and distributions. To facilitate rapid identification we generated DNA barcodes for the core plastid barcode loci (rbcL and matK) for 490 vascular plant species, representing nearly half of the Canadian Arctic flora and 93% of the flora of the Canadian Arctic Archipelago. Sequence recovery was higher for rbcL than matK (93% and 81%), and rbcL was easier to recover than matK from herbarium specimens (92% and 77%). Distance-based and sequence-similarity analyses of combined rbcL + matK data discriminate 97% of genera, 56% of species, and 7% of infraspecific taxa. There is a significant negative correlation between the number of species sampled per genus and the percent species resolution per genus. We characterize barcode variation in detail in the ten largest genera sampled (Carex, Draba, Festuca, Pedicularis, Poa, Potentilla, Puccinellia, Ranunculus, Salix, and Saxifraga) in the context of their phylogenetic relationships and taxonomy. Discrimination with the core barcode loci in these genera ranges from 0% in Salix to 85% in Carex. Haplotype variation in multiple genera does not correspond to species boundaries, including Taraxacum, in which the distribution of plastid haplotypes among Arctic species is consistent with plastid variation documented in non-Arctic species. Introgression of Poa glauca plastid DNA into multiple individuals of P. hartzii is problematic for identification of these species with DNA barcodes. Of three supplementary barcode loci (psbA–trnH, psbK–psbI, atpF–atpH) collected for a subset of Poa and Puccinellia species, only atpF–atpH improved discrimination in Puccinellia, compared with rbcL and matK. Variation in matK in Vaccinium uliginosum and rbcL in Saxifraga oppositifolia corresponds to variation in other loci used to characterize the phylogeographic histories of these Arctic-alpine species.
Despite its key location between the rest of the continent and Europe, research on the phylogeography of north African species remains very limited compared to European and North American taxa. The Mediterranean land mollusc Cornu aspersum (= Helix aspersa) is part of the few species widely sampled in north Africa for biogeographical analysis. It then provides an excellent biological model to understand phylogeographical patterns across the Mediterranean basin, and to evaluate hypotheses of population differentiation. We investigated here the phylogeography of this land snail to reassess the evolutionary scenario we previously considered for explaining its scattered distribution in the western Mediterranean, and to help to resolve the question of the direction of its range expansion (from north Africa to Europe or vice versa). By analysing simultaneously individuals from 73 sites sampled in its putative native range, the present work provides the first broad-scale screening of mitochondrial variation (cyt b and 16S rRNA genes) of C. aspersum.
Phylogeographical structure mirrored previous patterns inferred from anatomy and nuclear data, since all haplotypes could be ascribed to a B (West) or a C (East) lineage. Alternative migration models tested confirmed that C. aspersum most likely spread from north Africa to Europe. In addition to Kabylia in Algeria, which would have been successively a centre of dispersal and a zone of secondary contacts, we identified an area in Galicia where genetically distinct west and east type populations would have regained contact.
Vicariant and dispersal processes are reviewed and discussed in the light of signatures left in the geographical distribution of the genetic variation. In referring to Mediterranean taxa which show similar phylogeographical patterns, we proposed a parsimonious scenario to account for the "east-west" genetic splitting and the northward expansion of the western (B) clade which roughly involves (i) the dispersal of ancestral (eastern) types through Oligocene terranes in the Western Mediterranean (ii) the Tell Atlas orogenesis as gene flow barrier between future west and east populations, (iii) the impact of recurrent climatic fluctuations from mid-Pliocene to the last ice age, (iv) the loss of the eastern lineage during Pleistocene northwards expansion phases.
Background and Aims
Anatolia is a biologically diverse, but phylogeographically under-explored region. It is described as either a centre of origin and long-term Pleistocene refugium, or as a centre for genetic amalgamation, fed from distinct neighbouring refugia. These contrasting hypotheses are tested through a global phylogeographic analysis of the arctic–alpine herb, Arabis alpina.
Herbarium and field collections were used to sample comprehensively the entire global range, with special focus on Anatolia and Levant. Sequence variation in the chloroplast DNA trnL-trnF region was examined in 483 accessions. A haplotype genealogy was constructed and phylogeographic methods, demographic analysis and divergence time estimations were used to identify the centres of diversity and to infer colonization history.
Fifty-seven haplotypes were recovered, belonging to three haplogroups with non-overlapping distributions in (1) North America/Europe/northern Africa, (2) the Caucuses/Iranian Plateau/Arabian Peninsula and (3) Ethiopia–eastern Africa. All haplogroups occur within Anatolia, and all intermediate haplotypes linking the three haplogroups are endemic to central Anatolia and Levant, where haplotypic and nucleotide diversities exceeded all other regions. The local pattern of haplotype distribution strongly resembles the global pattern, and the haplotypes began to diverge approx. 2·7 Mya, coinciding with the climate cooling of the early Middle Pleistocene.
The phylogeographic structure of Arabis alpina is consistent with Anatolia being the cradle of origin for global genetic diversification. The highly structured landscape in combination with the Pleistocene climate fluctuations has created a network of mountain refugia and the accumulation of spatially arranged genotypes. This local Pleistocene population history has subsequently left a genetic imprint at the global scale, through four range expansions from the Anatolian diversity centre into Europe, the Near East, Arabia and Africa. Hence this study also illustrates the importance of sampling and scaling effects when translating global from local diversity patterns during phylogeographic analyses.
Anatolia; centre of origin; Pleistocene glaciations; chloroplast trnL-F; divergence times; alpine plants; Arabis alpina
Ligia isopods are widely distributed in the Pacific rocky intertidal shores from central California to central Mexico, including the Gulf of California. Yet, their biological characteristics restrict them to complete their life cycles in a very narrow range of the rocky intertidal supralittoral. Herein, we examine phylogeographic patterns of Ligia isopods from 122 localities between central California and central Mexico. We expect to find high levels of allopatric diversity. In addition, we expect the phylogeographic patterns to show signatures of past vicariant events that occurred in this geologically dynamic region.
We sequenced two mitochondrial genes (Cytochrome Oxidase I and 16S ribosomal DNA). We conducted Maximum Likelihood and Bayesian phylogenetic analyses. We found many divergent clades that, in general, group according to geography. Some of the most striking features of the Ligia phylogeographic pattern include: (1) deep mid-peninsular phylogeographic breaks on the Pacific and Gulf sides of Baja peninsula; (2) within the Gulf lineages, the northern peninsula is most closely related to the northern mainland, while the southern peninsula is most closely related to the central-southern mainland; and, (3) the southernmost portion of the peninsula (Cape Region) is most closely related to the southernmost portion of mainland.
Our results shed light on the phylogenetic relationships of Ligia populations in the study area. This study probably represents the finest-scale phylogeographic examination for any organism to date in this region. Presence of highly divergent lineages suggests multiple Ligia species exist in this region. The phylogeographic patterns of Ligia in the Gulf of California and Baja peninsula are incongruent with a widely accepted vicariant scenario among phylogeographers, but consistent with aspects of alternative geological hypotheses and phylo- and biogeographic patterns of several other taxa. Our findings contribute to the ongoing debate regarding the geological origin of this important biogeographic region.
There is increasing evidence that most parapatric cryptic/sister taxa are reproductively compatible across their areas of contact. Consequently, the biological species concept, which assumes absence of interbreeding, is becoming a not so effective criterion in evolutionary ecology. Nevertheless, the few parapatric sister taxa showing complete reproductive barriers represent interesting models to study speciation processes and the evolution of reproductive isolation. In this study, we examined contact populations in northwestern Italy of two butterfly species, Zerynthia polyxena and Z. cassandra, characterized by different genitalic morphotypes. We studied levels of divergence among 21 populations distributed from Sicily to France using three genetic markers (the mitochondrial COI and ND1 genes and the nuclear wingless gene) and genitalic geometric morphometrics. Moreover, we performed species distribution modelling to estimate different climatic requirements of Z. polyxena and Z. cassandra. We projected climatic data into glacial maximum scenarios in order to verify if and to which extent glacial cycles could have contributed to speciation processes. Genetic and morphometric analyses identified two main groups. All specimens showed a concordant pattern of diversification, including those individuals sampled in the contact area. Haplotype distribution and climatic models showed that during glacial maxima both species experienced a strong range contraction and presumably remained separated into different microrefugia in southern France, in the Italian Peninsula and on the islands of Elba and Sicily. Long term separation was probably favoured by reduced dispersal ability and high phylopatry, while genitalic diversification probably favoured interbreeding avoidance. Conversely, the aposematic wing pattern remained almost identical. We compared our results with those obtained in other species and concluded that Z. polyxena and Z. cassandra represent a valuable model in the study of speciation.
Background and Aims
Laeliinae are a neotropical orchid subtribe with approx. 1500 species in 50 genera. In this study, an attempt is made to assess generic alliances based on molecular phylogenetic analysis of DNA sequence data.
Six DNA datasets were gathered: plastid trnL intron, trnL-F spacer, matK gene and trnK introns upstream and dowstream from matK and nuclear ITS rDNA. Data were analysed with maximum parsimony (MP) and Bayesian analysis with mixed models (BA).
Although relationships between Laeliinae and outgroups are well supported, within the subtribe sequence variation is low considering the broad taxonomic range covered. Localized incongruence between the ITS and plastid trees was found. A combined tree followed the ITS trees more closely, but the levels of support obtained with MP were low. The Bayesian analysis recovered more well-supported nodes. The trees from combined MP and BA allowed eight generic alliances to be recognized within Laeliinae, all of which show trends in morphological characters but lack unambiguous synapomorphies.
By using combined plastid and nuclear DNA data in conjunction with mixed-models Bayesian inference, it is possible to delimit smaller groups within Laeliinae and discuss general patterns of pollination and hybridization compatibility. Furthermore, these small groups can now be used for further detailed studies to explain morphological evolution and diversification patterns within the subtribe.
Laeliinae; Orchidaceae; ITS; trnL intron; trnL-F spacer; matK
Climatic and sea-level fluctuations throughout the last Pleistocene glacial cycle (~130-0 ka) profoundly influenced present-day distributions and genetic diversity of Northern Hemisphere biotas by forcing range contractions in many species during the glacial advance and allowing expansion following glacial retreat ('expansion-contraction’ model). Evidence for such range dynamics and refugia in the unglaciated Gulf-Atlantic Coastal Plain stems largely from terrestrial species, and aquatic species Pleistocene responses remain relatively uninvestigated. Heterandria formosa, a wide-ranging regional endemic, presents an ideal system to test the expansion-contraction model within this biota. By integrating ecological niche modeling and phylogeography, we infer the Pleistocene history of this livebearing fish (Poeciliidae) and test for several predicted distributional and genetic effects of the last glaciation.
Paleoclimatic models predicted range contraction to a single southwest Florida peninsula refugium during the Last Glacial Maximum, followed by northward expansion. We inferred spatial-population subdivision into four groups that reflect genetic barriers outside this refuge. Several other features of the genetic data were consistent with predictions derived from an expansion-contraction model: limited intraspecific divergence (e.g. mean mtDNA p-distance = 0.66%); a pattern of mtDNA diversity (mean Hd = 0.934; mean π = 0.007) consistent with rapid, recent population expansion; a lack of mtDNA isolation-by-distance; and clinal variation in allozyme diversity with higher diversity at lower latitudes near the predicted refugium. Statistical tests of mismatch distributions and coalescent simulations of the gene tree lent greater support to a scenario of post-glacial expansion and diversification from a single refugium than to any other model examined (e.g. multiple-refugia scenarios).
Congruent results from diverse data indicate H. formosa fits the classic Pleistocene expansion-contraction model, even as the genetic data suggest additional ecological influences on population structure. While evidence for Plio-Pleistocene Gulf Coast vicariance is well described for many freshwater species presently codistributed with H. formosa, this species demography and diversification departs notably from this pattern. Species-specific expansion-contraction dynamics may therefore have figured more prominently in shaping Coastal Plain evolutionary history than previously thought. Our findings bolster growing appreciation for the complexity of phylogeographical structuring within North America’s southern refugia, including responses of Coastal Plain freshwater biota to Pleistocene climatic fluctuations.
Species that are widespread throughout historically glaciated and currently non-glaciated areas provide excellent opportunities to investigate the role of Pleistocene climatic change on the distribution of North American biodiversity. Many studies indicate that northern animal populations exhibit low levels of genetic diversity over geographically widespread areas whereas southern populations exhibit relatively high levels. Recently, paleoclimatic data have been combined with niche-based distribution modeling to locate possible refugia during the Last Glacial Maximum. Using phylogeographic, population, and paleoclimatic data, we show that the distribution and mitochondrial data for the millipede genus Narceus are consistent with classical examples of Pleistocene refugia and subsequent post-glacial population expansion seen in other organismal groups.
The phylogeographic structure of Narceus reveals a complex evolutionary history with signatures of multiple refugia in southeastern North America followed by two major northern expansions. Evidence for refugial populations were found in the southern Appalachian Mountains and in the coastal plain. The northern expansions appear to have radiated from two separate refugia, one from the Gulf Coastal Plain area and the other from the mid-Atlantic coastal region. Distributional models of Narceus during the Last Glacial Maximum show a dramatic reduction from the current distribution, with suitable ecological zones concentrated along the Gulf and Atlantic coastal plain. We found a strong correlation between these zones of ecological suitability inferred from our paleo-model with levels of genetic diversity derived from phylogenetic and population estimates of genetic structuring.
The signature of climatic change, during and after the Pleistocene, on the distribution of the millipede genus Narceus is evident in the genetic data presented. Niche-based historical distribution modeling strengthens the conclusions drawn from the genetic data and proves useful in identifying probable refugia. Such interdisciplinary biogeographic studies provide a comprehensive approach to understanding these processes that generate and maintain biodiversity as well as the framework necessary to explore questions regarding evolutionary diversification of taxa.
Background and Aims
Four species of Pinus subsection Australes occur in the Caribbean Basin: P. caribaea, P. cubensis, P. maestrensis and P. occidentalis. This study analyses the phylogeography of these species to assess possible colonization events from Central America to the islands and subsequent population expansions during glacial periods driven by both drier climate and larger emerged land areas.
Allele size data were obtained for plastid microsatellites for 314 individuals from 24 populations, covering the distribution range of subsection Australes in the Caribbean Basin.
In total, 113 plastid haplotypes were identified. The highest genetic diversity was found in populations of P. caribaea. Overall, Caribbean Basin populations fit the isolation by distance model. Significant phylogeographical structure was found (RST = 0·671 > permuted RST = 0·101; P < 0·0001). The haplotype network and a Bayesian analysis of population structure (BAPS) indicated different Central American origins for P. caribaea var. bahamensis and P. caribaea var. caribaea plastids, with Central America populations in northern and south-eastern groups. Sudden expansion times for BAPS clusters were close to three glacial maxima.
Central America contains ancestral plastid haplotypes. Population expansion has played a major role in the distribution of genetic diversity in P. caribaea var. hondurensis. Two colonization events gave rise to the P. caribaea var. bahamensis and P. caribaea var. caribaea lineages. Plastid variation in the eastern species (P. cubensis, P. maestrensis and P. occidentalis) evolved independently from that in P. caribaea var. caribaea. Incomplete lineage sorting between P. cubensis and P. maestrensis is apparent. Inferred expansion times for P. caribaea var. bahamensis and for the eastern lineages correspond to glacial maxima, whereas those for P. caribaea var. hondurensis correspond to the beginning of the temperature decrease that led to Marine Isotope Stage 8.
Caribbean Basin; historical demography; phylogeography; Pinus; plastid microsatellites; population genetics; subsection Australes
Understanding the evolutionary origin and the phylogeographic patterns of asexual taxa can shed light on the origin and maintenance of sexual reproduction. We assessed the geographic origin, genetic diversity, and phylogeographic history of obligate parthenogen diploid Artemia parthenogenetica populations, a widespread halophilic crustacean.
We analysed a partial sequence of the Cytochrome c Oxidase Subunit I mitochondrial gene from an extensive set of localities (including Eurasia, Africa, and Australia), and examined their phylogeographic patterns and the phylogenetic relationships of diploid A. parthenogenetica and its closest sexual relatives. Populations displayed an extremely low level of mitochondrial genetic diversity, with one widespread haplotype shared by over 79% of individuals analysed. Phylogenetic and phylogeographic analyses indicated a multiple and recent evolutionary origin of diploid A. parthenogenetica, and strongly suggested that the geographic origin of parthenogenesis in Artemia was in Central Asia. Our results indicate that the maternal sexual ancestors of diploid A. parthenogenetica were an undescribed species from Kazakhstan and A. urmiana.
We found evidence for multiple origin of parthenogenesis in Central Asia. Our results indicated that, shortly after its origin, diploid A. parthenogenetica populations underwent a rapid range expansion from Central Asia towards the Mediterranean region, and probably to the rest of its current geographic distribution. This contrasts with the restricted geographic distribution, strong genetic structure, and regional endemism of sexual Artemia lineages and other passively dispersed sexual continental aquatic invertebrates. We hypothesize that diploid parthenogens might have reached their current distribution in historical times, with a range expansion possibly facilitated by an increased availability of suitable habitat provided by anthropogenic activities, such as the spread of solar saltworks, aided by their natural dispersal vectors (i.e., waterbirds).
Through the study of the phylogeographic structure and demographic history of the common goby, Pomatoschistus microps, the influence of Quaternary climatic changes on the evolutionary history of coastal and marine fishes is investigated. Because of its sedentary life cycle in Mediterranean lagoons, it is also a good model to study more specifically if the formation of lagoons during the Holocene had an impact on population structure and demography. Mitochondrial sequences of Northeastern Atlantic and Western Mediterranean specimens were used for phylogenetic reconstructions as well as divergence time estimates, demographic history and population structure analyses. Pomatoschistus microps was a highly supported monophyletic clade including four lineages. It may have appeared 77,000 yr ago, and the divergence of its lineages likely occured shortly thereafter (between 61,000 and 54,000 yr). Most lineages had polytomic topologies, low nucleotide diversity and demographic analyses providing evidence of population expansion. Each lineage was characterized by a large number of private haplotypes. Most haplotypes found in Mediterranean localities were endemic, and one was dominant. Complex reticulated relationships connecting North European, Atlantic and Mediterranean haplotypes were observed. Moderate to high population structure was underlined. Contrary to previous published studies, no significant differentiation was observed between Atlantic and Mediterranean populations, indicating that the Gibraltar Strait is not a phylogeographic break for P. microps. Indeed, molecular dating combined with the tree topologies, phylogeographic and demographic analyses as well as high haplotype diversity underline a recent and rapid population divergence during the last glacial. However, population structure indicates that differentiation is an ongoing process. From an ancestral population trapped in the Atlantic, this goby colonized first northern Europe and later the Mediterranean Sea. Shared haplotypes could have dispersed in the western Mediterranean basin before the lagoon formation, while most private haplotypes, evidencing a recent isolation, probably diverged in lagoons after their closure.
Background and Aims
Hymenaea stigonocarpa (Fabaceae: Caesalpinioideae) is an endemic tree from the Brazilian cerrado (savanna vegetation), a biome classified as a hotspot for conservation priority. This study investigates the phylogeographic structure of H. stigonocarpa, in order to understand the processes that have led to its current spatial genetic pattern.
The polymorphism level and spatial distribution of variants of the plastid non-coding region between the genes psbC and trnS were investigated in 175 individuals from 17 populations, covering the greater part of the total distribution of the species. Molecular diversity indices were calculated and intra-specific relationships were inferred by the construction of haplotype networks using the median-joining method. Genetic differentiation among populations and main geographical groups was evaluated using spatial analysis of molecular variance (SAMOVA).
Twenty-three different haplotypes were identified. The level of differentiation among the populations analysed was relatively high (FST = 0·692). Phylogeographic analyses showed a clear association between the haplotype network and geographic distribution of populations, revealing three main geographical groups: western, central and eastern. SAMOVA corroborated this finding, indicating that most of the variation can be attributed to differences among these three groups (58·8 %), with little difference among populations within groups (FSC = 0·252).
The subdivision of the geographic distribution of H. stigonocarpa populations into three genetically differentiated groups can be associated with Quaternary climatic changes. The data suggest that during glacial times H. stigonocarpa populations became extinct in most parts of the southern present-day cerrado area. Milder climatic conditions in the north and eastern portions of the cerrado resulted in maintenance of populations in these regions. Thus it is inferred that the most southern part of the present-day cerrado was re-colonized by different lineages from northern parts of this biome, after postglacial climate amelioration.
Biogeography; cerrado; genetic structure; Quaternary climate changes; Fabaceae; Leguminosae; Hymenaea stigonocarpa; neotropical savannas; Pleistocene; phylogeography; psbC-trnS
The evolution and current distribution of the Sino-Tibetan flora have been greatly affected by historical geological events, such as the uplift of the Qinghai-Tibetan Plateau (QTP), and Quaternary climatic oscillations. Rhodiola kirilowii, a perennial herb with its distribution ranging from the southeastern QTP and the Hengduan Mountains (HM) to adjacent northern China and central Asia, provides an excellent model to examine and disentangle the effect of both geological orogeny and climatic oscillation on the evolutionary history of species with such distribution patterns. We here conducted a phylogeographic study using sequences of two chloroplast fragments (trnL-F and trnS-G) and internal transcribed spacers in 29 populations of R. kirilowii. A total of 25 plastid haplotypes and 12 ITS ribotypes were found. Molecular clock estimation revealed deep divergence between the central Asian populations and other populations from the HM and northern China; this split occurred ca. 2.84 million year ago. The majority of populations from the mountains of northern China were dominated by a single haplotype or ribotype, while populations of the HM harbored both high genetic diversity and high haplotype diversity. This distribution pattern indicates that HM was either a diversification center or a refugium for R. kirilowii during the Quaternary climatic oscillations. The present distribution of this species on mountains in northern China may have resulted from a rapid glacial population expansion from the HM. This expansion was confirmed by the mismatch distribution analysis and negative Tajima's D and Fu's FS values, and was dated to ca. 168 thousand years ago. High genetic diversity and population differentiation in both plastid and ITS sequences were revealed; these imply restricted gene flow between populations. A distinct isolation-by-distance pattern was suggested by the Mantel test. Our results show that in old lineages, populations may harbour divergent genetic forms that are sufficient to maintain or even increase overall genetic diversity despite fragmentation and low within-population variation.