Comparative phylogeography of unglaciated eastern North America

Cryptic lineage divergence in marine environments: genetic differentiation at multiple spatial and temporal scales in the widespread intertidal goby Gobiosoma bosc

The adaptive radiation of the seven-spined gobies (Gobiidae: Gobiosomatini) represents a classic example of how ecological specialization and larval retention can drive speciation through local adaptation. However, geographically widespread and phenotypically uniform species also do occur within Gobiosomatini. This lack of phenotypic variation across large geographic areas could be due to recent colonization, widespread gene flow, or stabilizing selection acting across environmental gradients. We use a phylogeographic approach to test these alternative hypotheses in the naked goby Gobiosoma bosc, a widespread and phenotypically invariable intertidal fish found along the Atlantic Coast of North America. Using DNA sequence from 218 individuals sampled at 15 localities, we document marked intraspecific genetic structure in mitochondrial and nuclear genes at three main geographic scales: (i) between Gulf of Mexico and Atlantic Coast, (ii) between the west coast of the Florida peninsula and adjacent Gulf of Mexico across the Apalachicola Bay, and (iii) at local scales of a few hundred kilometers. Clades on either side of Florida diverged about 8 million years ago, whereas some populations along the East Cost show divergent phylogroups that have differentiated within the last 200,000 years. The absence of noticeable phenotypic or ecological differentiation among lineages suggests the role of stabilizing selection on ancestral phenotypes, together with isolation in allopatry due to reduced dispersal and restricted gene flow, as the most likely explanation for their divergence. Haplotype phylogenies and spatial patterns of genetic diversity reveal frequent population bottlenecks followed by rapid population growth, particularly along the Gulf of Mexico. The magnitude of the genetic divergence among intraspecific lineages suggests the existence of cryptic species within Gobiosoma and indicates that modes of speciation can vary among lineages within Gobiidae.

Population structure within the one-dimensional range of a coastal plain katydid

Biogeography plays a significant role in species’ dispersal, and in turn population structure, across the landscape. The North American katydid Neoconocephalus melanorhinus belongs to a genus with high mobility. Unlike other members of the genus, N. melanorhinus is a salt marsh specialist restricted to a narrow corridor along the Atlantic and Gulf coasts. In addition, their range crosses at least one known biogeographic barrier and possesses biogeographic characteristics of the stepping-stone as well as the hierarchical island model of dispersal. Using AFLP markers we searched for areas that conform to the predictions of isolation by distance and for areas of non-uniform increases in genetic variance, indicative of isolation by barrier. We found significant genetic differentiation between all twelve sampled sites. Isolation by distance was the predominant pattern of variation across their range. In addition, we saw possible evidence of two biogeographic barriers to gene flow, one at the Atlantic-Gulf divide and the other along the Gulf coast. We also observed a change in body size across the range. Body size, as measured by male hind femur length, correlated closely with latitude, a possible indication of differential selection across the species range.

Dispersal corridors for plant species in the Poyang Lake Basin of southeast China identified by integration of phylogeographic and geospatial data

Measuring the dispersal of wildlife through landscapes is notoriously difficult. Recently, the categorical least cost path algorithm that integrates population genetic data with species distribution models has been applied to reveal population connectivity. In this study, we use this method to identify the possible dispersal corridors of five plant species (Castanopsis tibetana, Schima superba, Cyclocarya paliurus, Sargentodoxa cuneata, Eomecon chionantha) in the Poyang Lake Basin (PLB, largely coinciding with Jiangxi Province), China, in the late Quaternary. The results showed that the strongest population connectivity for the five species occurred in the Wuyi Mountains and the Yu Mountains of the eastern PLB (East Corridor) during the late Quaternary. In the western PLB, populations of the five species were connected by the Luoxiao Mountains and the Jiuling Mountains (West Corridor) but with a lower degree of connectivity. There were some minor connections between the eastern and the western populations across the Gannan Hills. When the corridors of five species were overlaid, the East Corridor and the West Corridor were mostly shared by multiple species. These results indicate that plant species in the PLB could have responded to the Quaternary climate changes by moving along the East Corridor and the West Corridor. Given that dispersal corridors have seldom been considered in the governmental strategies of biodiversity conservation in the PLB, preserving and restoring natural vegetation along these corridors should be prioritized to mitigate the effects of anthropogenic climate change by facilitating migration of plant species and other biota.

Diverse range dynamics and dispersal routes of plants on the Tibetan Plateau during the late Quaternary

Phylogeographical studies have suggested that several plant species on the Tibetan Plateau (TP) underwent recolonization during the Quaternary and may have had distinct range dynamics in response to the last glacial. To further test this hypothesis and locate the possible historical dispersal routes, we selected 20 plant species from different parts of the TP and modeled their geographical distributions over four time periods using species distribution models (SDMs). Furthermore, we applied the least-cost path method together with SDMs and shared haplotypes to estimate their historical dispersal corridors. We identified three general scenarios of species distribution change during the late Quaternary: the 'contraction-expansion' scenario for species in the northeastern TP, the 'expansion-contraction' scenario for species in the southeast and the 'stable' scenario for widespread species. During the Quaternary, we identified that these species were likely to recolonize along the low-elevation valleys, huge mountain ranges and flat plateau platform (e.g. the Yarlung Zangbo Valley and the Himalaya). We inferred that Quaternary cyclic glaciations along with the various topographic and climatic conditions of the TP could have resulted in the diverse patterns of range shift and dispersal of Tibetan plant species. Finally, we believe that this study would provide valuable insights for the conservation of alpine species under future climate change.

Alternative glacial-interglacial refugia demographic hypotheses tested on Cephalocereus columna-trajani (Cactaceae) in the intertropical Mexican drylands

Historic demography changes of plant species adapted to New World arid environments could be consistent with either the Glacial Refugium Hypothesis (GRH), which posits that populations contracted to refuges during the cold-dry glacial and expanded in warm-humid interglacial periods, or with the Interglacial Refugium Hypothesis (IRH), which suggests that populations contracted during interglacials and expanded in glacial times. These contrasting hypotheses are developed in the present study for the giant columnar cactus Cephalocereus columna-trajani in the intertropical Mexican drylands where the effects of Late Quaternary climatic changes on phylogeography of cacti remain largely unknown. In order to determine if the historic demography and phylogeographic structure of the species are consistent with either hypothesis, sequences of the chloroplast regions psbA-trnH and trnT-trnL from 110 individuals from 10 populations comprising the full distribution range of this species were analysed. Standard estimators of genetic diversity and structure were calculated. The historic demography was analysed using a Bayesian approach and the palaeodistribution was derived from ecological niche modelling to determine if, in the arid environments of south-central Mexico, glacial-interglacial cycles drove the genetic divergence and diversification of this species. Results reveal low but statistically significant population differentiation (F_ST = 0.124, P < 0.001), although very clear geographic clusters are not formed. Genetic diversity, haplotype network and Approximate Bayesian Computation (ABC) demographic analyses suggest a population expansion estimated to have taken place in the Last Interglacial (123.04 kya, 95% CI 115.3–130.03). The species palaeodistribution is consistent with the ABC analyses and indicates that the potential area of palaedistribution and climatic suitability were larger during the Last Interglacial and Holocene than in the Last Glacial Maximum. Overall, these results suggest that C. columna-trajani experienced an expansion following the warm conditions of interglacials, in accordance with the GRH.

Distribution of Karyotypes of the Cryptocercus punctulatus Species Complex (Blattodea: Cryptocercidae) in Great Smoky Mountains National Park

During the period between 1999 and 2006, wood-feeding cockroaches in the Cryptocercus punctulatus Scudder species complex were collected throughout Great Smoky Mountains National Park, USA. The chromosome numbers of insects from 59 sites were determined, and phylogenetic analyses were performed based on mitochondrial COII and nuclear ITS2 DNA. The distribution of the three male karyotypes found in the park (2n = 37, 39, and 45) is mapped and discussed in relation to recent disturbances and glacial history. Clades of the three karyotype groups meet near the ridgeline separating North Carolina from Tennessee in the center of the park, suggesting that these may have originated from separate lower elevation refugia after the last glacial maximum. The timing of divergence and a significant correlation between elevation difference and genetic distance in two of the clades supports this hypothesis. The ecological role of the cockroaches in the park is discussed.

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0803-0) contains supplementary material, which is available to authorized users.

Phylogeographic Insights into a Peripheral Refugium: The Importance of Cumulative Effect of Glaciation on the Genetic Structure of Two Endemic Plants

Quaternary glaciations and mostly last glacial maximum have shaped the contemporary distribution of many species in the Alps. However, in the Maritime and Ligurian Alps a more complex picture is suggested by the presence of many Tertiary paleoendemisms and by the divergence time between lineages in one endemic species predating the Late Pleistocene glaciation. The low number of endemic species studied limits the understanding of the processes that took place within this region. We used species distribution models and phylogeographical methods to infer glacial refugia and to reconstruct the phylogeographical pattern of Silene cordifolia All. and Viola argenteria Moraldo & Forneris. The predicted suitable area for last glacial maximum roughly fitted current known distribution. Our results suggest that separation of the major clades predates the last glacial maximum and the following repeated glacial and interglacial periods probably drove differentiations. The complex phylogeographical pattern observed in the study species suggests that both populations and genotypes extinction was minimal during the last glacial maximum, probably due to the low impact of glaciations and to topographic complexity in this area. This study underlines the importance of cumulative effect of previous glacial cycles in shaping the genetic structure of plant species in Maritime and Ligurian Alps, as expected for a Mediterranean mountain region more than for an Alpine region.

Inferences of population structure and demographic history for Taxodium distichum, a coniferous tree in North America, based on amplicon sequencing analysis

PREMISE OF THE STUDY

Studies of natural genetic variation can elucidate the genetic basis of phenotypic variation and the past population structure of species. Our study species, Taxodium distichum, is a unique conifer that inhabits the flood plains and swamps of North America. Morphological and ecological differences in two varieties, T. distichum var. distichum (bald cypress) and T. distichum var. imbricarium (pond cypress), are well known, but little is known about the level of genetic differentiation between the varieties and the demographic history of local populations.

METHODS

We analyzed nucleotide polymorphisms at 47 nuclear loci from 96 individuals collected from the Mississippi River Alluvial Valley (MRAV), and Gulf Coastal populations in Texas, Louisiana, and Florida using high-throughput DNA sequencing. Standard population genetic statistics were calculated, and demographic parameters were estimated using a composite-likelihood approach.

KEY RESULTS

Taxodium distichum in North America can be divided into at least three genetic groups, bald cypress in the MRAV and Texas, bald cypress in Florida, and pond cypress in Florida. The levels of genetic differentiation among the groups were low but significant. Several loci showed the signatures of positive selection, which might be responsible for local adaptation or varietal differentiation.

CONCLUSIONS

Bald cypress was genetically differentiated into two geographical groups, and the boundary was located between the MRAV and Florida. This differentiation could be explained by population expansion from east to west. Despite the overlap of the two varieties' ranges, they were genetically differentiated in Florida. The estimated demographic parameters suggested that pond cypress split from bald cypress during the late Miocene.

Multi-locus sequence data illuminate demographic drivers of Pleistocene speciation in semi-arid southern Australian birds (Cinclosoma spp.)

BACKGROUND

During the Pleistocene, shifts of species distributions and their isolation in disjunct refugia led to varied outcomes in how taxa diversified. Some species diverged, others did not. Here, we begin to address another facet of the role of the Pleistocene in generating today's diversity. We ask which processes contributed to divergence in semi-arid southern Australian birds. We isolated 11 autosomal nuclear loci and one mitochondrial locus from a total of 29 specimens of the sister species pair, Chestnut Quail-thrush Cinclosoma castanotum and Copperback Quail-thrush C. clarum.

RESULTS

A population clustering analysis confirmed the location of the current species boundary as a well-known biogeographical barrier in southern Australia, the Eyrean Barrier. Coalescent-based analyses placed the time of species divergence to the Middle Pleistocene. Gene flow between the species since divergence has been low. The analyses suggest the effective population size of the ancestor was 54 to 178 times smaller than populations since divergence. This contrasts with recent multi-locus studies in some other Australian birds (butcherbirds, ducks) where a lack of phenotypic divergence was accompanied by larger historical population sizes. Post-divergence population size histories of C. clarum and C. castanotum were inferred using the extended Bayesian skyline model. The population size of C. clarum increased substantially during the late Pleistocene and continued to increase through the Last Glacial Maximum and Holocene. The timing of this expansion across its vast range is broadly concordant with that documented in several other Australian birds. In contrast, effective population size of C. castanotum was much more constrained and may reflect its smaller range and more restricted habitat east of the Eyrean Barrier compared with that available to C. clarum to the west.

CONCLUSIONS

Our results contribute to awareness of increased population sizes, following significant contractions, as having been important in shaping diversity in Australian arid and semi-arid zones. Further, we improve knowledge of the role of Pleistocene climatic shifts in areas of the planet that were not glaciated at that time but which still experienced that period's cyclical climatic fluctuations.

Diversification in wild populations of the model organism Anolis carolinensis: A genome-wide phylogeographic investigation

The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.

High migration rates shape the postglacial history of amphi-Atlantic bryophytes

Paleontological evidence and current patterns of angiosperm species richness suggest that European biota experienced more severe bottlenecks than North American ones during the last glacial maximum. How well this pattern fits other plant species is less clear. Bryophytes offer a unique opportunity to contrast the impact of the last glacial maximum in North America and Europe because about 60% of the European bryoflora is shared with North America. Here, we use population genetic analyses based on approximate Bayesian computation on eight amphi-Atlantic species to test the hypothesis that North American populations were less impacted by the last glacial maximum, exhibiting higher levels of genetic diversity than European ones and ultimately serving as a refugium for the postglacial recolonization of Europe. In contrast with this hypothesis, the best-fit demographic model involved similar patterns of population size contractions, comparable levels of genetic diversity and balanced migration rates between European and North American populations. Our results thus suggest that bryophytes have experienced comparable demographic glacial histories on both sides of the Atlantic. Although a weak, but significant genetic structure was systematically recovered between European and North American populations, evidence for migration from and towards both continents suggests that amphi-Atlantic bryophyte population may function as a metapopulation network. Reconstructing the biogeographic history of either North American or European bryophyte populations therefore requires a large, trans-Atlantic geographic framework.

Idiosyncratic responses of evergreen broad-leaved forest constituents in China to the late Quaternary climate changes

Subtropical evergreen broad-leaved forest (EBLF) is one of the most important vegetation types in China. Inferences from palaeo-biome reconstruction (PBR) and phylogeography regarding range shift history of EBLF during the late Quaternary are controversial and should be reconciled. We compared phylogeographic patterns of three EBLF constituents in China, Castanopsis tibetana, Machilus thunbergii and Schima superba. Contrary to a chorus of previous phylogeographic studies and the results of species distribution modelling (SDM) of this study (in situ survival during the LGM), the three species displayed three different phylogeographic patterns that conform to either an in situ survival model or an expansion-contraction model. These results are partially congruent with the inference of PBR that EBLF was absent to the north of 24° N at the LGM. This study suggests that the constituents of EBLF could have responded idiosyncratically to climate changes during the Late Quaternary. The community assemblages of EBLF could have been changing over time, resulting in no palaeo-analogs to modern-day EBLF, which may be the main reason responsible for the failure of PBR to detect the occurrence of EBLF north of 24° N at the LGM.

Evolutionary relationships within the lamioid tribe Synandreae (Lamiaceae) based on multiple low-copy nuclear loci

The subfamily Lamioideae (Lamiaceae) comprises ten tribes, of which only Stachydeae and Synandreae include New World members. Previous studies have investigated the phylogenetic relationships among the members of Synandreae based on plastid and nuclear ribosomal DNA loci. In an effort to re-examine the phylogenetic relationships within Synandreae, the current study incorporates data from four low-copy nuclear loci, PHOT1, PHOT2, COR, and PPR. Our results confirm previous studies based on chloroplast and nuclear ribosomal markers in supporting the monophyly of tribe Synandreae, as well as sister relationships between Brazoria and Warnockia, and between that pair of genera and a monophyletic Physostegia. However, we observe incongruence in the relationships of Macbridea and Synandra. The placement of Synandreae within Lamioideae is poorly resolved and incongruent among different analyses, and the sister group of Synandreae remains enigmatic. Comparison of the colonization and migration patterns corroborates a single colonization of the New World by Synandreae during the Late Miocene/Tortonian age. This is in contrast to the only other lamioid tribe that includes New World members, Stachydeae, which colonized the New World at least twice—during the mid-Miocene and Pliocene. Edaphic conditions and intolerance of soil acidity may be factors that restricted the distribution of most genera of Synandreae to southeastern and south–central North America, whereas polyploidy could have increased the colonizing capability of the more wide-ranging genus, Physostegia.

Toward a paradigm shift in comparative phylogeography driven by trait-based hypotheses

For three decades, comparative phylogeography has conceptually and methodologically relied on the concordance criterion for providing insights into the historical/biogeographic processes driving population genetic structure and divergence. Here we discuss how this emphasis, and the corresponding lack of methods for extracting information about biotic/intrinsic contributions to patterns of genetic variation, may bias our general understanding of the factors driving genetic structure. Specifically, this emphasis has promoted a tendency to attribute discordant phylogeographic patterns to the idiosyncracies of history, as well as an adherence to generic null expectations of concordance with reduced predictive power. We advocate that it is time for a paradigm shift in comparative phylogeography, especially given the limited utility of the concordance criterion as genomic data provide ever-increasing levels of resolution. Instead of adhering to the concordance-discordance dichotomy, comparative phylogeography needs to emphasize the contribution of taxon-specific traits that will determine whether concordance is a meaningful criterion for evaluating hypotheses or may predict discordant phylogeographic structure. Through reference to some case studies we illustrate how refined hypotheses based on taxon-specific traits can provide improved predictive frameworks to forecast species responses to climatic change or biogeographic barriers while gaining unique insights about the taxa themselves and their interactions with their environment. We outline a potential avenue toward a synthetic comparative phylogeographic paradigm that includes addressing some important conceptual and methodological challenges related to study design and application of model-based approaches for evaluating support of trait-based hypotheses under the proposed paradigm.

Chemotypic Variation of Conocephalum salebrosum in the Southeastern Appalachian Range: A Search for Cryptic Plant Biodiversity Around the Tennessee River Valley

The chemotaxonomy of the Conocephalum spp. complex, based on GC-MS analysis of the volatile compositions, has helped to reveal cryptic biodiversity and delineate actual distribution patterns of chemotypes. In the Appalachian Mountains, two samples from eastern central part of the range were previously shown to be C. salebrosum. Additionally, it has recently come to light that stress can alter the volatile composition of C. conicum. This study address a previously unsampled region of the southeastern Appalachians, a region that is a biodiversity epicenter, to determine if more chemotypic diversity remains to be seen for the Conocephalum spp. complex. A common garden experiment was performed, but yielded more of a common stress experiment, and significantly altered the volatile compositions. Wild-collected controls and a meta-analysis of these data and those from previous works suggest that the common garden experiment caused stress and that the liverworts sampled belong to the C. salebrosum clade of of the Conocephalum spp. complex.

Evolutionary History of a Desert Shrub Ephedra przewalskii (Ephedraceae): Allopatric Divergence and Range Shifts in Northwestern China

Based on two chloroplast DNA sequences, psbA-trnH and trnT-trnF, phylogeographical patterns of a desert shrub, Ephedra przewalskii, were examined across most of its geographic range in northwestern China. A total of sixteen haplotypes were detected. There was a common haplotype in each basin, that was haplotype A in Tarim Basin, haplotype G in Junggar Basin, and haplotype M in Qaidam Basin. Genetic variance mainly occurred among populations, geographic regions, and eleven geographic groups subdivided by SAMOVA analysis. E. przewalskii likely had a smaller and more fragmented geographic range during the Last Glacial Maximum, which was determined based on ecological niche modelling. Three groups of E. przewalskii populations were detected to have experience range expansion, and this was based on significant values of Fu’s F_S, Tajima’s D, and unimodel mismatch distributions. The cold and dry climate during the glacial period of the Quaternary is postulated to have been a driver for significant genetic isolation and divergence among populations or groups in E. przewalskii, whereas the warmer and wetter climate during the interglacial period is speculated to have provided favourable conditions for range expansion of the species.

The population genomic basis of geographic differentiation in North American common ragweed (Ambrosia artemisiifolia L.)

Common ragweed (Ambrosia artemisiifolia L.) is an invasive, wind-pollinated plant nearly ubiquitous in disturbed sites in its eastern North American native range and present across growing portions of Europe, Africa, Asia, and Australia. Phenotypic divergence between European and native-range populations has been described as rapid evolution. However, a recent study demonstrated major human-mediated shifts in ragweed genetic structure before introduction to Europe and suggested that native-range genetic structure and local adaptation might fully explain accelerated growth and other invasive characteristics of introduced populations. Genomic differentiation that potentially influenced this structure has not yet been investigated, and it remains unclear whether substantial admixture during historical disturbance of the native range contributed to the development of invasiveness in introduced European ragweed populations. To investigate fine-scale population genetic structure across the species' native range, we characterized diallelic SNP loci via a reduced-representation genotyping-by-sequencing (GBS) approach. We corroborate phylogeographic domains previously discovered using traditional sequencing methods, while demonstrating increased power to resolve weak genetic structure in this highly admixed plant species. By identifying exome polymorphisms underlying genetic differentiation, we suggest that geographic differentiation of this important invasive species has occurred more often within pathways that regulate growth and response to defense and stress, which may be associated with survival in North America's diverse climatic regions.

Phylogeographic analyses and genetic structure illustrate the complex evolutionary history of Phragmites australis in Mexico

PREMISE OF THE STUDY

Genetic data suggest that three lineages of Phragmites australis are found in North America: the Native North American lineage, the Gulf Coast lineage, and the Invasive lineage. In Mexico, P. australis is a common species, but nothing is known about the distribution or ecology of these lineages. We examined the phylogeography of P. australis to analyze the current geographic distribution of genetic variation, demographic history, and dispersal patterns to better understand its evolutionary history in Mexico.

METHODS

We sampled 427 individuals from 28 populations. We used two noncoding regions of chloroplast DNA to estimate the levels of genetic variation and identified the genetic groups across the species' geographical range in Mexico. We compared the genealogical relationships among haplotypes with those previously reported. A hypothesis of demographic expansion was also tested for the Mexican P. australis lineages.

KEY RESULTS

We found 13 new haplotypes native to Mexico that might be undergoing an active process of expansion and diversification. Genealogical analyses provided evidence that two independent lineages of P. australis are present in Mexico. The invasive lineage was not detected with our sampling. Our estimates of population expansions in Mexico ranged from 0.202 to 0.726 mya.

CONCLUSIONS

Phragmites australis is a native species that has been in Mexico for thousands of years. Genetic data suggest that climatic changes during the Pleistocene played an important role in the demographic expansion of the populations that constitute the different genetic groups of P. australis in Mexico.

Phylogeographic concordance factors quantify phylogeographic congruence among co-distributed species in the Sarracenia alata pitcher plant system

Comparative phylogeographic investigations have identified congruent phylogeographic breaks in co-distributed species in nearly every region of the world. The qualitative assessments of phylogeographic patterns traditionally used to identify such breaks, however, are limited because they rely on identifying monophyletic groups across species and do not account for coalescent stochasticity. Only long-standing phylogeographic breaks are likely to be obvious; many species could have had a concerted response to more recent landscape events, yet possess subtle signs of phylogeographic congruence because ancestral polymorphism has not completely sorted. Here, we introduce Phylogeographic Concordance Factors (PCFs), a novel method for quantifying phylogeographic congruence across species. We apply this method to the Sarracenia alata pitcher plant system, a carnivorous plant with a diverse array of commensal organisms. We explore whether a group of ecologically associated arthropods have co-diversified with the host pitcher plant, and identify if there is a positive correlation between ecological interaction and PCFs. Results demonstrate that multiple arthropods share congruent phylogeographic breaks with S. alata, and provide evidence that the level of ecological association can be used to predict the degree of similarity in the phylogeographic pattern. This study outlines an approach for quantifying phylogeographic congruence, a central concept in biogeographic research.

Model-based analysis supports interglacial refugia over long-dispersal events in the diversification of two South American cactus species

Pilosocereus machrisii and P. aurisetus are cactus species within the P. aurisetus complex, a group of eight cacti that are restricted to rocky habitats within the Neotropical savannas of eastern South America. Previous studies have suggested that diversification within this complex was driven by distributional fragmentation, isolation leading to allopatric differentiation, and secondary contact among divergent lineages. These events have been associated with Quaternary climatic cycles, leading to the hypothesis that the xerophytic vegetation patches which presently harbor these populations operate as refugia during the current interglacial. However, owing to limitations of the standard phylogeographic approaches used in these studies, this hypothesis was not explicitly tested. Here we use Approximate Bayesian Computation to refine the previous inferences and test the role of different events in the diversification of two species within P. aurisetus group. We used molecular data from chloroplast DNA and simple sequence repeats loci of P. machrisii and P. aurisetus, the two species with broadest distribution in the complex, in order to test if the diversification in each species was driven mostly by vicariance or by long-dispersal events. We found that both species were affected primarily by vicariance, with a refuge model as the most likely scenario for P. aurisetus and a soft vicariance scenario most probable for P. machrisii. These results emphasize the importance of distributional fragmentation in these species, and add support to the hypothesis of long-term isolation in interglacial refugia previously proposed for the P. aurisetus species complex diversification.

Molecular Phylogeny and Biogeography of the Amphidromous Fish Genus Dormitator Gill 1861 (Teleostei: Eleotridae)

Species of the genus Dormitator, also known as sleepers, are representatives of the amphidromous freshwater fish fauna that inhabit the tropical and subtropical coastal environments of the Americas and Western Africa. Because of the distribution of this genus, it could be hypothesized that the evolutionary patterns in this genus, including a pair of geminate species across the Central American Isthmus, could be explained by vicariance following the break-up of Gondwana. However, the evolutionary history of this group has not been evaluated. We constructed a time-scaled molecular phylogeny of Dormitator using mitochondrial (Cytochrome b) and nuclear (Rhodopsin and β-actin) DNA sequence data to infer and date the cladogenetic events that drove the diversification of the genus and to relate them to the biogeographical history of Central America. Two divergent lineages of Dormitator were recovered: one that included all of the Pacific samples and another that included all of the eastern and western Atlantic samples. In contrast to the Pacific lineage, which showed no phylogeographic structure, the Atlantic lineage was geographically structured into four clades: Cameroon, Gulf of Mexico, West Cuba and Caribbean, showing evidence of potential cryptic species. The separation of the Pacific and Atlantic lineages was estimated to have occurred ~1 million years ago (Mya), whereas the four Atlantic clades showed mean times of divergence between 0.2 and 0.4 Mya. The splitting times of Dormitator between ocean basins are similar to those estimated for other geminate species pairs with shoreline estuarine preferences, which may indicate that the common evolutionary histories of the different clades are the result of isolation events associated with the closure of the Central American Isthmus and the subsequent climatic and oceanographic changes.

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation-interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.

Revision of widespread red squirrels (genus: Tamiasciurus) highlights the complexity of speciation within North American forests

Integration of molecular methods, ecological modeling, and statistical hypothesis testing are increasing our understanding of differentiation within species and phylogenetic relationships among species by revealing environmental connections to evolutionary processes. Within mammals, novel diversity is being discovered and characterized as more complete geographic sampling is coupled with newer multi-disciplinary approaches. North American red squirrels exemplify a forest obligate genus whose species are monitored as indicators of forest ecosystem condition, yet phylogenetic relationships reflecting evolutionary history within this genus remain tentative. Through testing of competing systematic and niche-based divergence hypotheses, we recognize three species, Tamiasciurus douglasii, T. hudsonicus, and T. fremonti. Our data provide evidence of regional differences in evolutionary dynamics and continental gradients of complexity that are important both for future management and for investigating multiple pathways that can lead to the formation of new species.

Comparative phylogeography of black mangroves (Avicennia germinans) and red mangroves (Rhizophora mangle) in Florida: Testing the maritime discontinuity in coastal plants

PREMISE OF THE STUDY

Previous studies of the comparative phylogeography of coastal and marine species in the southeastern United States revealed that phylogenetically diverse taxa share a phylogeographic break at the southern tip of Florida (the maritime discontinuity). These studies have focused nearly exclusively on animals; few coastal plant species in Florida have been analyzed phylogeographically. We investigated phylogeographic patterns of black mangroves (Avicennia germinans) and red mangroves (Rhizophora mangle), two coastal trees that occur on both coasts of the peninsula of Florida.

METHODS

We sampled and genotyped 150 individuals each of A. germinans and R. mangle, using eight microsatellite loci per species. We used observed and expected heterozygosity to quantify genetic diversity in each sampling location and allele frequencies to identify putative phylogeographic breaks and measure gene flow using BayesAss and Migrate-n. We tested the hypothesis that both species would exhibit a phylogeographic break at the southern tip of Florida.

KEY RESULTS

We did not find any significant phylogeographic breaks in either species. Rhizophora mangle exhibits greater genetic structure than A. germinans, contrary to expectations based on propagule dispersal capability. However, directional gene flow from the Gulf to the Atlantic was more pronounced in R. mangle, indicating that the Gulf Stream may affect genetic patterns in R. mangle more than in A. germinans.

CONCLUSIONS

The high dispersal capability of these species may lead to high genetic connectivity between sampling locations and little geographic structure. We also identified several locations that, based on genetic data, should be the focus of conservation efforts.

Post-glacial phylogeography and evolution of a wide-ranging highly-exploited keystone forest tree, eastern white pine (Pinus strobus) in North America: single refugium, multiple routes

Background

Knowledge of the historical distribution and postglacial phylogeography and evolution of a species is important to better understand its current distribution and population structure and potential fate in the future, especially under climate change conditions, and conservation of its genetic resources. We have addressed this issue in a wide-ranging and heavily exploited keystone forest tree species of eastern North America, eastern white pine (Pinus strobus). We examined the range-wide population genetic structure, tested various hypothetical population history and evolutionary scenarios and inferred the location of glacial refugium and post-glacial recolonization routes. Our hypothesis was that eastern white pine survived in a single glacial refugium and expanded through multiple post-glacial recolonization routes.

Results

We studied the range-wide genetic diversity and population structure of 33 eastern white pine populations using 12 nuclear and 3 chloroplast microsatellite DNA markers. We used Approximate Bayesian Computation approach to test various evolutionary scenarios. We observed high levels of genetic diversity, and significant genetic differentiation (F_ST = 0.104) and population structure among eastern white pine populations across its range. A south to north trend of declining genetic diversity existed, consistent with repeated founder effects during post-glaciation migration northwards. We observed broad consensus from nuclear and chloroplast genetic markers supporting the presence of two main post-glacial recolonization routes that originated from a single southern refugium in the mid-Atlantic plain. One route gave rise to populations at the western margin of the species’ range in Minnesota and western Ontario. The second route gave rise to central-eastern populations, which branched into two subgroups: central and eastern. We observed minimal sharing of chloroplast haplotypes between recolonization routes but there was evidence of admixture between the western and west-central populations.

Conclusions

Our study reveals a single southern refugium, two recolonization routes and three genetically distinguishable lineages in eastern white pine that we suggest to be treated as separate Evolutionarily Significant Units. Like many wide-ranging North American species, eastern white pine retains the genetic signatures of post-glacial recolonization and evolution, and its contemporary population genetic structure reflects not just the modern distribution and effects of heavy exploitation but also routes northward from its glacial refugium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0624-1) contains supplementary material, which is available to authorized users.

Phylogeography and evolution of the Red Salamander (Pseudotriton ruber)

Phylogeographic studies frequently result in the elevation of subspecific taxa to species given monophyly, or the synonymy of subspecies that are not monophyletic. However, given limited or incongruent datasets, retention of subspecies can be useful to describe hypothesized incipient species or to illustrate interesting biological phenomena driving morphological diversity. Four subspecific taxa have been used to describe largely allopatric geographic variation within the species Pseudotriton ruber, a plethodontid salamander occupying stream and spring habitats across eastern North America: P. r. vioscai occurs in lowland Coastal Plain habitats, while P. r. ruber, P. r. nitidus, and P. r. schencki occupy upland regions in and around the Appalachian Mountains. Pseudotriton ruber co-occurs through its distribution with the aposematic newt Notophthalmus viridescens, and both species are hypothesized to be part of a Müllerian mimicry complex. In this study, we sequenced regions of two mitochondrial (cytochrome b, NADH dehydrogenase subunit 2) and one single copy nuclear protein-coding gene (pro-opiomelanocortin) from individuals sampled across much of the distribution of P. ruber and then used maximum-likelihood and Bayesian phylogenetic inference to test the monophyly of subspecies, reconstruct biogeographic history, and make inferences about morphological evolution. Phylogeographic hypotheses from mitochondrial and nuclear datasets described structure among populations of P. ruber which separated Coastal Plain and upland Appalachian populations, but subspecies were not monophyletic. Biogeographic reconstruction estimated the ancestor of all populations to have occupied and initially diverged in the Coastal Plain during the Pliocene (∼3.6mya), before one lineage subsequently invaded upland areas of Appalachia. Bold bright coloration of high elevation subspecies P. r. nitidus and P. r. schencki appears to have evolved twice. We hypothesize that the Müllerian mimicry complex with N. viridescens and P. ruber may provide a selective mechanism driving the co-evolution of striking bright and dull morphological variation among populations of both species. While P. ruber subspecies were not consistent with our criteria for diagnosing species (monophyly) and therefore could not be elevated to species, we advocate for the retention of subspecies because they describe hypotheses about an incipient species (P. r. vioscai) and how Müllerian mimicry may shape morphological diversity of species.

Biogeographic barriers drive co-diversification within associated eukaryotes of the Sarracenia alata pitcher plant system

Understanding if the members of an ecological community have co-diversified is a central concern of evolutionary biology, as co-diversification suggests prolonged association and possible coevolution. By sampling associated species from an ecosystem, researchers can better understand how abiotic and biotic factors influence diversification in a region. In particular, studies of co-distributed species that interact ecologically can allow us to disentangle the effect of how historical processes have helped shape community level structure and interactions. Here we investigate the Sarracenia alata pitcher plant system, an ecological community where many species from disparate taxonomic groups live inside the fluid-filled pitcher leaves. Direct sequencing of the eukaryotes present in the pitcher plant fluid enables us to better understand how a host plant can shape and contribute to the genetic structure of its associated inquilines, and to ask whether genetic variation in the taxa are structured in a similar manner to the host plant. We used 454 amplicon-based metagenomics to demonstrate that the pattern of genetic diversity in many, but not all, of the eukaryotic community is similar to that of S. alata, providing evidence that associated eukaryotes share an evolutionary history with the host pitcher plant. Our work provides further evidence that a host plant can influence the evolution of its associated commensals.

Pliocene-Pleistocene lineage diversifications in the Eastern Indigo Snake (Drymarchon couperi) in the Southeastern United States

Indigo Snakes (Drymarchon; with five currently recognized species) occur from northern Argentina, northward to the United States in southern Texas and eastward in disjunct populations in Florida and Georgia. Based on this known allopatry and a difference in supralabial morphology the two United States taxa previously considered as subspecies within D. corais (Boie 1827), the Western Indigo Snake, D. melanurus erebennus (Cope 1860), and Eastern Indigo Snake, D. couperi (Holbrook 1842), are currently recognized as separate species. Drymarchon couperi is a Federally-designated Threatened species by the United States Fish and Wildlife Service under the Endangered Species Act, and currently being incorporated into a translocation program. This, combined with its disjunct distribution makes it a prime candidate for studying speciation and genetic divergence. In this study, we (1) test the hypothesis that D. m. erebennus and D. couperi are distinct lineages by analyzing 2411 base pairs (bp) of two mitochondrial (mtDNA) loci and one single copy nuclear (scnDNA) locus; (2) estimate the timing of speciation using a relaxed phylogenetics method to determine if Milankovitch cycles during the Pleistocene might have had an influence on lineage diversifications; (3) examine historical population demography to determine if identified lineages have undergone population declines, expansions, or remained stable during the most recent Milankovitch cycles; and (4) use this information to assist in an effective and scientifically sound translocation program. Our molecular data support the initial hypothesis that D. melanurus and D. couperi should be recognized as distinct species, but further illustrate that D. couperi is split into two distinct genetic lineages that correspond to historical biogeography and sea level changes in peninsular Florida. These two well-supported genetic lineages (herein termed Atlantic and Gulf lineages) illustrate a common biogeographic distributional break previously identified for other plants and animals, suggesting that these organisms might have shared a common evolutionary history related to historic sea level changes caused by Milankovitch cycles. Our estimated divergence times suggest that the most recent common ancestor (MRCA) between D. melanurus and southeastern United States Drymarchon occurred ca. 5.9Ma (95% HPD=2.5-9.8Ma; during the late Blancan of the Pleistocene through the Hemphillian of the Miocene), whereas the MRCA between the Atlantic and Gulf lineages in the southeastern United States occurred ca. 2.0Ma (95% HPD=0.7-3.7Ma; during the Irvingtonian of the Pleistocene through the Blancan of the Pliocene). During one or more glacial intervals within these times, these two lineages must have become separated and evolved independently. Despite numerous Milankovitch cycles along with associated forming of physical barriers (i.e., sea level fluctuations, high elevation sand ridges, clayey soils, and/or insufficient habitats) since their initial lineage diversification, these two lineages have likely come in and out of contact with each other many times, yet today they still illustrate near discrete geographic distributions. Although the Atlantic and Gulf lineages appear to be cryptic, a thorough study examining morphological characters should be conducted. We believe that our molecular data is crucial and should be incorporated in making conscious decisions in the management of a translocation program. We suggest that source populations for translocations include maintaining the integrity of the known genetic lineages found herein, as well as those coming from the closest areas that currently support sizable Drymarchon populations.

Contact zone dynamics during early stages of speciation in a chorus frog (Pseudacris crucifer)

Characterizing the genetic and behavioural consequences of contact between previously geographically isolated lineages provides insights into the mechanisms underlying diversification and ultimately speciation. The spring peeper (Pseudacris crucifer) is a widespread Nearctic chorus frog with six divergent mitochondrial DNA (mtDNA) lineages, many of which came into secondary contact during the Holocene. We examined genetics, morphology, advertisement calls and female preference for two lineages that began diverging in allopatry in the Pliocene and now overlap in southwestern Ontario, Canada. We found non-coincident clines in mtDNA and nuclear DNA, mirroring directionality of premating isolation barriers. We also found divergence in a range of traits between these two lineages, displacement in male call attributes and female preference for calls of their natal lineage in sympatry. Hybrids were morphologically distinct from both parental lineages, but hybrid male calls were acoustically intermediate. Female hybrids showed asymmetrical preference for Eastern male calls. These results considered together provide evidence of either unidirectional hybridization or selection against hybrids, potentially implying reproductive character displacement. Our work demonstrates the utility of integrated, multi-character approaches to understanding the processes of divergence and the nature of speciation.

The aggregate site frequency spectrum for comparative population genomic inference

Understanding how assemblages of species responded to past climate change is a central goal of comparative phylogeography and comparative population genomics, an endeavour that has increasing potential to integrate with community ecology. New sequencing technology now provides the potential to perform complex demographic inference at unprecedented resolution across assemblages of nonmodel species. To this end, we introduce the aggregate site frequency spectrum (aSFS), an expansion of the site frequency spectrum to use single nucleotide polymorphism (SNP) data sets collected from multiple, co-distributed species for assemblage-level demographic inference. We describe how the aSFS is constructed over an arbitrary number of independent population samples and then demonstrate how the aSFS can differentiate various multispecies demographic histories under a wide range of sampling configurations while allowing effective population sizes and expansion magnitudes to vary independently. We subsequently couple the aSFS with a hierarchical approximate Bayesian computation (hABC) framework to estimate degree of temporal synchronicity in expansion times across taxa, including an empirical demonstration with a data set consisting of five populations of the threespine stickleback (Gasterosteus aculeatus). Corroborating what is generally understood about the recent postglacial origins of these populations, the joint aSFS/hABC analysis strongly suggests that the stickleback data are most consistent with synchronous expansion after the Last Glacial Maximum (posterior probability = 0.99). The aSFS will have general application for multilevel statistical frameworks to test models involving assemblages and/or communities, and as large-scale SNP data from nonmodel species become routine, the aSFS expands the potential for powerful next-generation comparative population genomic inference.

Combining the least cost path method with population genetic data and species distribution models to identify landscape connectivity during the late Quaternary in Himalayan hemlock

Himalayan hemlock (Tsuga dumosa) experienced a recolonization event during the Quaternary period; however, the specific dispersal routes are remain unknown. Recently, the least cost path (LCP) calculation coupled with population genetic data and species distribution models has been applied to reveal the landscape connectivity. In this study, we utilized the categorical LCP method, combining species distribution of three periods (the last interglacial, the last glacial maximum, and the current period) and locality with shared chloroplast, mitochondrial, and nuclear haplotypes, to identify the possible dispersal routes of T. dumosa in the late Quaternary. Then, both a coalescent estimate of migration rates among regional groups and establishment of genetic divergence pattern were conducted. After those analyses, we found that the species generally migrated along the southern slope of Himalaya across time periods and genomic makers, and higher degree of dispersal was in the present and mtDNA haplotype. Furthermore, the direction of range shifts and strong level of gene flow also imply the existence of Himalayan dispersal path, and low area of genetic divergence pattern suggests that there are not any obvious barriers against the dispersal pathway. Above all, we inferred that a dispersal route along the Himalaya Mountains could exist, which is an important supplement for the evolutionary history of T. dumosa. Finally, we believed that this integrative genetic and geospatial method would bring new implications for the evolutionary process and conservation priority of species in the Tibetan Plateau.

In and out of refugia: historical patterns of diversity and demography in the North American Caesar's mushroom species complex

Some of the effects of past climate dynamics on plant and animal diversity make-up have been relatively well studied, but to less extent in fungi. Pleistocene refugia are thought to harbour high biological diversity (i.e. phylogenetic lineages and genetic diversity), mainly as a product of increased reproductive isolation and allele conservation. In addition, high extinction rates and genetic erosion are expected in previously glaciated regions. Some of the consequences of past climate dynamics might involve changes in range and population size that can result in divergence and incipient or cryptic speciation. Many of these dynamic processes and patterns can be inferred through phylogenetic and coalescent methods. In this study, we first delimit species within a group of closely related edible ectomycorrhizal Amanita from North America (the American Caesar's mushrooms species complex) using multilocus coalescent-based approaches; and then address questions related to effects of Pleistocene climate change on the diversity and genetics of the group. Our study includes extensive geographical sampling throughout the distribution range, and DNA sequences from three nuclear protein-coding genes. Results reveal cryptic diversity and high speciation rates in refugia. Population sizes and expansions seem to be larger at midrange latitudes (Mexican highlands and SE USA). Range shifts are proportional to population size expansions, which were overall more common during the Pleistocene. This study documents responses to past climate change in fungi and also highlights the applicability of the multispecies coalescent in comparative phylogeographical analyses and diversity assessments that include ancestral species.

Large scale patterns of genetic variation and differentiation in sugar maple from tropical Central America to temperate North America

Background

Geological events in the latter Cenozoic have influenced the distribution, abundance and genetic structure of tree populations in temperate and tropical North America. The biogeographical history of temperate vegetation that spans large ranges of latitude is complex, involving multiple latitudinal shifts that might have occurred via different migration routes. We determined the regional structuring of genetic variation of sugar maple (Acer saccharum subsp. saccharum) and its only subspecies in tropical America (Acer saccharum subsp. skutchii) using nuclear and chloroplast data. The studied populations span a geographic range from Maine, USA (46°N), to El Progreso, Guatemala (15°N). We examined genetic subdivisions, explored the locations of ancestral haplotypes, analyzed genetic data to explore the presence of a single or multiple glacial refugia, and tested whether genetic lineages are temporally consistent with a Pleistocene or older divergence.

Results

Nuclear and chloroplast data indicated that populations in midwestern USA and western Mexico were highly differentiated from populations in the rest of the sites. The time of the most recent common ancestor of the western Mexico haplotype lineage was dated to the Pliocene (5.9 Ma, 95 % HPD: 4.3–7.3 Ma). Splits during the Pleistocene separated the rest of the phylogroups. The most frequent and widespread haplotype occurred in half of the sites (Guatemala, eastern Mexico, southeastern USA, and Ohio). Our data also suggested that multiple Pleistocene refugia (tropics-southeastern USA, midwestern, and northeastern USA), but not western Mexico (Jalisco), contributed to post-glacial northward expansion of ranges. Current southern Mexican and Guatemalan populations have reduced population sizes, genetic bottlenecks and tend toward homozygosity, as indicated using nuclear and chloroplast markers.

Conclusions

The divergence of western Mexican populations from the rest of the sugar maples likely resulted from orographic and volcanic barriers to gene flow. Past connectivity among populations in the southeastern USA and eastern Mexico and Guatemala possible occurred through gene flow during the Pleistocene. The time to the most common ancestor values revealed that populations from the Midwest and Northeast USA represented different haplotype lineages, indicating major divergence of haplotypes lineages before the Last Glacial Maximum and suggesting the existence of multiple glacial refugia.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0518-7) contains supplementary material, which is available to authorized users.

A new Antaeotricha species from Florida sandhills and scrub (Lepidoptera, Depressariidae, Stenomatinae)

Antaeotricha floridella sp. n. is described and diagnosed from the closely similar Antaeotricha albulella (Walker). The species is distributed in xeric sandhill and scrub habitats in peninsular Florida, USA, and larvae feed on Quercus species. Keys are given for pale-winged Stenomatinae and similar Gelechioidea based on external characters and genitalia.