Comparative phylogeography of unglaciated eastern North America

At the confluence of vicariance and dispersal: Phylogeography of cavernicolous springtails (Collembola: Arrhopalitidae, Tomoceridae) codistributed across a geologically complex karst landscape in Illinois and Missouri

The processes of vicariance and dispersal are central to our understanding of diversification, yet determining the factors that influence these processes remains a significant challenge in evolutionary biology. Caves offer ideal systems for examining the mechanisms underlying isolation, divergence, and speciation. Intrinsic ecological differences among cavernicolous organisms, such as the degree of cave dependence, are thought to be major factors influencing patterns of genetic isolation in caves. Using a comparative phylogeographic approach, we employed mitochondrial and nuclear markers to assess the evolutionary history of two ecologically distinct groups of terrestrial cave-dwelling springtails (Collembola) in the genera Pygmarrhopalites (Arrhopalitidae) and Pogonognathellus (Tomoceridae) that are codistributed in caves throughout the Salem Plateau-a once continuous karst region, now bisected by the Mississippi River Valley in Illinois and Missouri. Contrasting phylogeographic patterns recovered for troglobiotic Pygmarrhopalites sp. and eutroglophilic Pogonognathellus sp. suggests that obligate associations with cave habitats can restrict dispersal across major geographic barriers such as rivers and valleys, but may also facilitate subterranean dispersal between neighboring cave systems. Pygmarrhopalites sp. populations spanning the Mississippi River Valley were estimated to have diverged 2.9-4.8 Ma, which we attribute to vicariance resulting from climatic and geological processes involved in Mississippi River Valley formation beginning during the late Pliocene/early Pleistocene. Lastly, we conclude that the detection of many deeply divergent, morphologically cryptic, and microendemic lineages highlights our poor understanding of microarthropod diversity in caves and exposes potential conservation concerns.

Molecular variation across populations of a widespread North American firefly, Photinus pyralis, reveals that coding changes do not underlie flash color variation or associated visual sensitivity

Background

Genes underlying signal production and reception are expected to evolve to maximize signal detection in specific environments. Fireflies vary in their light signal color both within and between species, and thus provide an excellent system in which to study signal production and reception in the context of signaling environments. Differences in signal color have been hypothesized to be due to variation in the sequence of luciferase, the enzyme that catalyzes the light reaction. Similarly, differences in visual sensitivity, which are expected to match signal color, have been hypothesized to be due to variation in the sequence of opsins, the protein component of visual pigments. Here we investigated (1) whether sequence variation in luciferase correlates with variation in signal color and (2) whether sequence variation in opsins correlates with inferred matching visual sensitivity across populations of a widespread North American firefly species, Photinus pyralis. We further tested (3) whether selection has acted on these loci by examining their population-level differentiation relative to the distribution of differentiation derived from a genome-wide sample of loci generated by double-digest RADseq.

Results

We found virtually no coding variation in luciferase or opsins. However, there was extreme divergence in non-coding variation in luciferase across populations relative to a panel of random genomic loci.

Conclusions

The absence of protein variation at both loci challenges the paradigm that variation in signal color and visual sensitivity in fireflies is exclusively due to coding variation in luciferase and opsin genes. Instead, flash color variation within species must involve other mechanisms, such as abdominal pigmentation or regulation of light organ physiology. Evidence for selection at non-coding variation in luciferase suggests that selection is targeting luciferase regulation and may favor differ expression levels across populations.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1251-9) contains supplementary material, which is available to authorized users.

A survey of molecular diversity and population genetic structure in North American clearwing moths (Lepidoptera: Sesiidae) using cytochrome c oxidase I

The phylogeographic structure of insect species in North America is poorly understood. The moth family Sesiidae (Lepidoptera) contains many economically important pests of agriculture and forestry, as well as beneficial species used in biological control. Despite their significance, this study constitutes the first broad-ranging population genetic study on North American sesiids. It probes the population structure of eight species of sesiid moths based on sequence variation in cytochrome c oxidase I (N = 191). Haplotype diversity levels were high in seven of the eight species, while nucleotide diversity varied considerably. Patterns ranged from limited structure and a starburst pattern in the raspberry crown borer Pennisetia marginata to highly geographically structured populations in the peachtree borer Synanthedon exitiosa and the maple callus borer Synanthedon acerni. These varied patterns suggest differing evolutionary histories and dispersal abilities. By elucidating population genetic structure and barriers to dispersal we can begin to devise conservation and management plans.

Approximate Bayesian computation analysis of EST-associated microsatellites indicates that the broadleaved evergreen tree Castanopsis sieboldii survived the Last Glacial Maximum in multiple refugia in Japan

Climatic changes have played major roles in plants' evolutionary history. Glacial oscillations have been particularly important, but some of their effects on plants' populations are poorly understood, including the numbers and locations of refugia in Asian warm temperate zones. In the present study, we investigated the demographic history of the broadleaved evergreen tree species Castanopsis sieboldii (Fagaceae) during the last glacial period in Japan. We used approximate Bayesian computation (ABC) for model comparison and parameter estimation for the demographic modeling using 27 EST-associated microsatellites. We also performed the species distribution modeling (SDM). The results strongly support a demographic scenario that the Ryukyu Islands and the western parts in the main islands (Kyushu and western Shikoku) were derived from separate refugia and the eastern parts in the main islands and the Japan Sea groups were diverged from the western parts prior to the coldest stage of the Last Glacial Maximum (LGM). Our data indicate that multiple refugia survived at least one in the Ryukyu Islands, and the other three regions of the western and eastern parts and around the Japan Sea of the main islands of Japan during the LGM. The SDM analysis also suggests the potential habitats under LGM climate conditions were mainly located along the Pacific Ocean side of the coastal region. Our ABC-based study helps efforts resolve the demographic history of a dominant species in warm temperate broadleaved forests during and after the last glacial period, which provides a basic model for future phylogeographical studies using this approach.

A biogeographical profile of the sand cockroach Arenivaga floridensis and its bearing on origin hypotheses for Florida scrub biota

Florida scrub is a xeric ecosystem associated with the peninsula's sand ridges, whose intermittent Pliocene-Pleistocene isolation is considered key to scrub endemism. One scrub origin hypothesis posits endemics were sourced by the Pliocene dispersal of arid-adapted taxa from southwestern North America; a second invokes Pleistocene migration within eastern North America. Only one study to date has explicitly tested these competing hypotheses, supporting an eastern origin for certain scrub angiosperms. For further perspective, we conducted a genetic analysis of an endemic arthropod, the Florida sand cockroach (Arenivaga floridensis), with two aims: (1) to reconstruct the peninsular colonization and residence history of A. floridensis and (2) determine whether its biogeographic profile favors either origin hypothesis. We sequenced the cox2 mitochondrial gene for 237 specimens (65 populations) as well as additional loci (cox1, nuclear H3) for a subset of Florida roaches and congeners. Using Network and Bayesian inference methods, we identified three major lineages whose genetic differentiation and phylogeographical structure correspond with late Pliocene peninsula insularization, indicating Arenivaga was present and broadly distributed in Florida at that time. Stem and crown divergence estimates (6.36 Ma; 2.78 Ma) between A. floridensis and western sister taxa span a period of extensive dispersal by western biota along an arid Gulf Coast corridor. These phylogeographical and phylogenetic results yield a biogeographic profile consistent with the western origin hypothesis. Moreover, age estimates for the roach's peninsular residence complement those of several other endemics, favoring a Pliocene (or earlier) inception of the scrub ecosystem. We argue that eastern versus western hypotheses are not mutually exclusive; rather, a composite history of colonization involving disparate biotas better explains the diverse endemism of Florida scrub.

Phylogeography and Ecological Niche Modeling Reveal Reduced Genetic Diversity and Colonization Patterns of Skunk Cabbage (Symplocarpus foetidus; Araceae) From Glacial Refugia in Eastern North America

Alternating glacial and interglacial periods during the Quaternary have dramatically affected the distribution and population genetic structure of plant and animal species throughout the northern hemisphere. Surprisingly, little is known about the post-glacial recolonization history of wetland herbaceous perennials that are widely distributed in the understory of deciduous or mixed deciduous-evergreen forests in eastern North America. In this study, we investigated infraspecific variation among 32 populations of skunk cabbage, Symplocarpus foetidus, to test the hypothesis that the extant species diversity of skunk cabbage is the result of a post-glacial range expansion from southern refugia during the Quaternary Ice Age. A total of 4041 base pairs (bp) of the chloroplast intergenic spacer region (cpDNA) was sequenced from 485 individuals sampled from glaciated (18 populations, 275 individuals) and unglaciated (14 populations, 210 individuals) regions east and west of the Appalachian Mountains. Haplotype number, haplotype diversity, and nucleotide diversity were calculated, and genetic variation within and among populations was assessed by analysis of molecular variance (AMOVA). The geographic pattern of genetic differentiation was further investigated with a spatial analysis of molecular variance (SAMOVA). A total of eight haplotypes and three genetic groups (SAMOVA) were recovered and a much higher haplotype number (eight haplotypes) and haplotype diversity (0.7425) was observed in unglaciated compared to glaciated populations (five haplotypes, haplotype diversity = 0.6099). All haplotypes found in glaciated regions represented a subset of haplotypes found in unglaciated regions. Haplotypes of S. foetidus likely diverged during the Tertiary (mid-Miocene and late Pliocene), predating the last glacial maximum (LGM). Predictions based on ecological niche modeling (ENM) suggested that there was considerably less suitable habitat for skunk cabbage during the LGM, and the habitat range was further south compared to the current distribution. Reduced variation and a subset of haplotypes in glaciated regions suggest a founder effect associated with range expansion via long-distance seed dispersal. Our results do not support the "Driftless Area" scenario for the northern refugium, rather the data suggest a "Northeastern" refugium near the southernmost extent of the LGM.

Evolutionary history of Mexican domesticated and wild Meleagris gallopavo

Background

The distribution of the wild turkey (Meleagris gallopavo) extends from Mexico to southeastern Canada and to the eastern and southern regions of the USA. Six subspecies have been described based on morphological characteristics and/or geographical variations in wild and domesticated populations. In this paper, based on DNA sequence data from the mitochondrial D-loop, we investigated the genetic diversity and structure, genealogical relationships, divergence time and demographic history of M. gallopavo populations including domesticated individuals.

Results

Analyses of 612 wild and domesticated turkey mitochondrial D-loop sequences, including 187 that were collected for this study and 425 from databases, revealed 64 haplotypes with few mutations, some of which are shared between domesticated and wild turkeys. We found a high level of haplotype and nucleotide diversity, which suggests that the total population of this species is large and stable with an old evolutionary history. The results of genetic differentiation, haplotype network, and genealogical relationships analyses revealed three main genetic groups within the species: mexicana as a population relict (C1), merriami (C2), and mexicana/intermedia/silvestris/osceola (C3). Haplotypes detected in domesticated turkeys belong to group C3. Estimates of divergence times agree with range expansion and diversification events of the relict population of M. gallopavo in northwestern Mexico during the Pliocene–Pleistocene and Pleistocene–Holocene boundaries. Demographic reconstruction showed that an expansion of the population occurred 110,000 to 130,000 years ago (Kya), followed by a stable period 100 Kya and finally a decline ~ 10 Kya (Pleistocene–Holocene boundary). In Mexico, the Trans-Mexican Volcanic Belt may be responsible for the range expansion of the C3 group. Two haplotypes with different divergence times, MGMDgoB/MICH1 and MICH2, are dominant in domesticated and commercial turkeys.

Conclusions

During the Pleistocene, a large and stable population of M. gallopavo covered a wide geographic distribution from the north to the center of America (USA and Mexico). The mexicana, merriami, and mexicana/intermedia/silvestris/osceola genetic groups originated after divergence and range expansion from northwestern Mexico during the Pliocene–Pleistocene and Pleistocene–Holocene boundaries. Old and new maternal lines of the mexicana/intermedia/silvestris/osceola genetic group were distributed within the Trans-Mexican Volcanic Belt where individuals were captured for domestication. Two haplotypes are the main founder maternal lines of domesticated turkeys.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0388-8) contains supplementary material, which is available to authorized users.

Resolving relationships and phylogeographic history of the Nyssa sylvatica complex using data from RAD-seq and species distribution modeling

Nyssa sylvatica complex consists of several woody taxa occurring in eastern North America. These taxa were recognized as two or three species including three or four varieties by different authors. Due to high morphological similarities and complexity of morphological variation, classification and delineation of taxa in the group have been difficult and controversial. Here we employ data from RAD-seq to elucidate the genetic structure and phylogenetic relationships within the group. Using the genetic evidence, we evaluate previous classifications and delineate species. We also employ Species Distribution Modeling (SDM) to evaluate impacts of climatic changes on the ranges of the taxa and to gain insights into the relevant refugia in eastern North America. Results from Molecular Variance Analysis (AMOVA), STRUCTURE, phylogenetic analyses using Maximum likelihood, Bayesian Inference, and Splittree methods of RAD-seq data strongly support a two-clade pattern, largely separating samples of N. sylvatica from those of N. biflora-N. ursina mix. Divergence time analysis with BEAST suggests the two clades diverged in the mid Miocene. The ancestor of the present trees of N. sylvatica was suggested to be in the Pliocene and that of N. biflora-N. ursina mix in the end of the Miocene. Results from SDM predicted a smaller range in the southern part of the species present range of each clade during the Last Glacial Maximum (LGM). A northward expansion of the ranges during interglacial period and a northward shift of the ranges in the future under a model of global warming were also predicted. Our results support the recognition of two species in the complex, N. sylvatica and N. biflora, following the phylogenetic species concept. We found no genetic evidence supporting recognitions of intraspecific taxa. However, we propose subsp. ursina and subsp. biflora within N. biflora due to their distinction in habits, distributions, and habitats. Our results further support movements of trees in eastern North America in response to climatic changes. Finally, our study demonstrates that RAD-seq data and a combination of population genomics and SDM are valuable in resolving relationship and biogeographic history of closely related species that are taxonomically difficult.

Phylogeographic structure in three North American tent caterpillar species (Lepidoptera: Lasiocampidae): Malacosoma americana, M. californica, and M. disstria

While phylogeographic structure has been examined in many North American vertebrate species, insects have received much less attention despite their central ecological roles. The moth genus Malacosoma (Hübner, 1820), is an important group of forestry pests responsible for large-scale defoliation across much of the Nearctic and Palearctic. The present study uses sequence variation in the mitochondrial cytochrome c oxidase 1 (COI) gene to examine the population genetic structure of the three widespread Malacosoma species (M. americana, M. californica, and M. disstria). Populations of all three species showed highest diversity in the south, suggesting that modern populations derived from southern refugia with loss of variation as these lineages dispersed northwards. However, despite similar life histories and dispersal abilities, the extent of regional variation varied among the taxa. M. americana, a species restricted to eastern North America, showed much less genetic structure than the western M. californica or the widespread M. disstria. The regional differentiation in the latter reflects the likely derivation of modern lineages from several refugia, as well as taxonomic uncertainty in M. californica. In these respects, the three species of Malacosoma share phylogeographic patterns similar to those detected in vertebrates which are characterised by greater phylogeographic breaks in the western half of the continent and limited structure in the east.

The role of diversification in community assembly of the oaks (Quercus L.) across the continental U.S

PREMISE OF THE STUDY

Evolutionary and biogeographic history, including past environmental change and diversification processes, are likely to have influenced the expansion, migration, and extinction of populations, creating evolutionary legacy effects that influence regional species pools and the composition of communities. We consider the consequences of the diversification process in shaping trait evolution and assembly of oak-dominated communities throughout the continental United States (U.S.).

METHODS

Within the U.S. oaks, we tested for phylogenetic and functional trait patterns at different spatial scales, taking advantage of a dated phylogenomic analysis of American oaks and the U.S. Forest Service (USFS) Forest Inventory and Analysis (FIA).

KEY RESULTS

We find (1) phylogenetic overdispersion at small grain sizes throughout the U.S. across all spatial extents and (2) a shift from overdispersion to clustering with increasing grain sizes. Leaf traits have evolved in a convergent manner, and these traits are clustered in communities at all spatial scales, except in the far west, where species with contrasting leaf types co-occur.

CONCLUSIONS

Our results support the hypotheses that (1) interspecific interactions were important in parallel adaptive radiation of the genus into a range of habitats across the continent and (2) that the diversification process is a critical driver of community assembly. Functional convergence of complementary species from distinct clades adapted to the same local habitats is a likely mechanism that allows distantly related species to coexist. Our findings contribute to an explanation of the long-term maintenance of high oak diversity and the dominance of the oak genus in North America.

Pervasive migration across rainforest and sandy coastal plain Aechmea nudicaulis (Bromeliaceae) populations despite contrasting environmental conditions

Understanding the colonization of extreme marginal habitats and the relative roles of space and environment in maintaining peripheral populations remains challenging. Here, we leverage a system of pairs of rainforest and sandy coastal plain communities that allow us to decouple spatial and environmental effects in the population structure and migration rates of the bromeliad Aechmea nudicaulis. Structure and gene flow between populations were estimated from Bayesian clustering and coalescent-based migration models applied to chloroplast sequence and nuclear microsatellite data. Contrary to our initial expectation, the sharp environmental gradient between rainforest and sandy plains does not seem to have affected the colonization and migration dynamics in A. nudicaulis. Our analyses uncover pervasive gene flow between neighbouring habitats in both chloroplast and nuclear data despite the striking differences in environmental conditions. This result is consistent with a scenario of repeated colonization of the sandy coastal plains from forest populations through seed dispersal, as well as the maintenance of gene flow between habitats through pollination. We also recovered a broad north/south population structure that has been found in other Atlantic rainforest groups and possibly reflects older phylogeographic dynamics.

Challenges of comprehensive taxon sampling in comparative biology: Wrestling with rosids

Using phylogenetic approaches to test hypotheses on a large scale, in terms of both species sampling and associated species traits and occurrence data-and doing this with rigor despite all the attendant challenges-is critical for addressing many broad questions in evolution and ecology. However, application of such approaches to empirical systems is hampered by a lingering series of theoretical and practical bottlenecks. The community is still wrestling with the challenges of how to develop species-level, comprehensively sampled phylogenies and associated geographic and phenotypic resources that enable global-scale analyses. We illustrate difficulties and opportunities using the rosids as a case study, arguing that assembly of biodiversity data that is scale-appropriate-and therefore comprehensive and global in scope-is required to test global-scale hypotheses. Synthesizing comprehensive biodiversity data sets in clades such as the rosids will be key to understanding the origin and present-day evolutionary and ecological dynamics of the angiosperms.

Spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus revealed by museum genomics

Analyzing genetic variation through time and space is important to identify key evolutionary and ecological processes in populations. However, using contemporary genetic data to infer the dynamics of genetic diversity may be at risk of a bias, as inferences are performed from a set of extant populations, setting aside unavailable, rare, or now extinct lineages. Here, we took advantage of new developments in next-generation sequencing to analyze the spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus, a steppic Southwestern-Palearctic species. We applied a recently developed hybridization capture (hyRAD) protocol that allows retrieving orthologous sequences even from degraded DNA characteristic of museum specimens. We identified single nucleotide polymorphisms in 68 historical and 51 modern samples in order to (i) unravel the spatial genetic structure across part of the species distribution and (ii) assess the loss of genetic diversity over the past century in Swiss populations. Our results revealed (i) the presence of three potential glacial refugia spread across the European continent and converging spatially in the Alpine area. In addition, and despite a limited population sample size, our results indicate (ii) a loss of allelic richness in contemporary Swiss populations compared to historical populations, whereas levels of expected heterozygosities were not significantly different. This observation is compatible with an increase in the bottleneck magnitude experienced by central European populations of O. decorus following human-mediated land-use change impacting steppic habitats. Our results confirm that application of hyRAD to museum samples produces valuable information to study genetic processes across time and space.

Comparison of genetic variation between northern and southern populations of Lilium cernuum (Liliaceae): Implications for Pleistocene refugia

The so-called "Baekdudaegan" (BDDG), a mountain range that stretches along the Korean Peninsula, has been recently proposed as a major "southern" glacial refugium for boreal or temperate plant species based on palaeoecological and, especially, genetic data. Genetic studies comparing genetic variation between population occurring on the BDDG and more northern ones (i.e. in NE China and/or in Russian Far East) are, however, still too few to draw firm conclusions on the role of the BDDG as a refugium and a source for possible northward post-glacial recolonizations. In order to fill this gap, we selected a boreal/temperate herb, Lilium cernuum, and compared levels of allozyme-based genetic diversity of five populations from NE China with five populations from South Korea (home of its hypothesized refuge areas). As a complementary tool, we used the maximum entropy algorithm implemented in MaxEnt to infer the species' potential distribution for the present time, which was projected to different past climate scenarios for the Last Glacial Maximum (LGM). Permutation tests revealed that Korean populations harbored significantly higher levels of within-population genetic variation than those from NE China (expected heterozygosity = 0.173 vs. 0.095, respectively). Our results suggest that the lowered levels of genetic diversity in NE Chinese populations might be due to founder effects associated with post-glacial migration from southern regions. Congruent with genetic data, past distribution models showed higher probability of occurrence in southern ranges than in northern ones during the LGM. In addition, a positive correlation was detected between the expected heterozygosity and environmental LGM suitability. From a conservation perspective, our results further suggest that the southern populations in South Korea may be particularly worthy of protection.

Climatic suitability, isolation by distance and river resistance explain genetic variation in a Brazilian whiptail lizard

Spatial patterns of genetic variation can help understand how environmental factors either permit or restrict gene flow and create opportunities for regional adaptations. Organisms from harsh environments such as the Brazilian semiarid Caatinga biome may reveal how severe climate conditions may affect patterns of genetic variation. Herein we combine information from mitochondrial DNA with physical and environmental features to study the association between different aspects of the Caatinga landscape and spatial genetic variation in the whiptail lizard Ameivula ocellifera. We investigated which of the climatic, environmental, geographical and/or historical components best predict: (1) the spatial distribution of genetic diversity, and (2) the genetic differentiation among populations. We found that genetic variation in A. ocellifera has been influenced mainly by temperature variability, which modulates connectivity among populations. Past climate conditions were important for shaping current genetic diversity, suggesting a time lag in genetic responses. Population structure in A. ocellifera was best explained by both isolation by distance and isolation by resistance (main rivers). Our findings indicate that both physical and climatic features are important for explaining the observed patterns of genetic variation across the xeric Caatinga biome.

Adding loci improves phylogeographic resolution in red mangroves despite increased missing data: comparing microsatellites and RAD-Seq and investigating loci filtering

The widespread adoption of RAD-Seq data in phylogeography means genealogical relationships previously evaluated using relatively few genetic markers can now be addressed with thousands of loci. One challenge, however, is that RAD-Seq generates complete genotypes for only a small subset of loci or individuals. Simulations indicate that loci with missing data can produce biased estimates of key population genetic parameters, although the influence of such biases in empirical studies is not well understood. Here we compare microsatellite data (8 loci) and RAD-Seq data (six datasets ranging from 239 to 25,198 loci) from red mangroves (Rhizophora mangle) in Florida to evaluate how different levels of data filtering influence phylogeographic inferences. For all datasets, we calculated population genetic statistics and evaluated population structure, and for RAD-Seq datasets, we additionally examined population structure using coalescence. We found higher F_ST using microsatellites, but that RAD-Seq-based estimates approached those based on microsatellites as more loci with more missing data were included. Analyses of RAD-Seq datasets resolved the classic Gulf-Atlantic coastal phylogeographic break, which was not significant in the microsatellite analyses. Applying multiple levels of filtering to RAD-Seq datasets can provide a more complete picture of potential biases in the data and elucidate subtle phylogeographic patterns.

Crop domestication facilitated rapid geographical expansion of a specialist pollinator, the squash bee Peponapis pruinosa

Squash was first domesticated in Mexico and is now found throughout North America (NA) along with Peponapis pruinosa, a pollen specialist bee species of the squash genus Cucurbita The origin and spread of squash cultivation is well-studied archaeologically and phylogenetically; however, no study has documented how cultivation of this or any other crop has influenced species in mutualistic interactions. We used molecular markers to reconstruct the demographic range expansion and colonization routes of P. pruinosa from its native range into temperate NA. Populations east of the Rocky Mountains expanded from the wild host plant's range in Mexico and were established by a series of founder events. Eastern North America was most likely colonized from squash bee populations in the present-day continental Midwest USA and not from routes that followed the Gulf and Atlantic coasts from Mexico. Populations of P. pruinosa west of the Rockies spread north from the warm deserts much more recently, showing two genetically differentiated populations with no admixture: one in California and the other one in eastern Great Basin. These bees have repeatedly endured severe bottlenecks as they colonized NA, following human spread of their Cucurbita pollen hosts during the Holocene.

Evolutionary lessons from California plant phylogeography

Phylogeography documents the spatial distribution of genetic lineages that result from demographic processes, such as population expansion, population contraction, and gene movement, shaped by climate fluctuations and the physical landscape. Because most phylogeographic studies have used neutral markers, the role of selection may have been undervalued. In this paper, we contend that plants provide a useful evolutionary lesson about the impact of selection on spatial patterns of neutral genetic variation, when the environment affects which individuals can colonize new sites, and on adaptive genetic variation, when environmental heterogeneity creates divergence at specific loci underlying local adaptation. Specifically, we discuss five characteristics found in plants that intensify the impact of selection: sessile growth form, high reproductive output, leptokurtic dispersal, isolation by environment, and the potential to evolve longevity. Collectively, these traits exacerbate the impact of environment on movement between populations and local selection pressures-both of which influence phylogeographic structure. We illustrate how these unique traits shape these processes with case studies of the California endemic oak, Quercus lobata, and the western North American lichen, Ramalina menziesii Obviously, the lessons we learn from plant traits are not unique to plants, but they highlight the need for future animal, plant, and microbe studies to incorporate its impact. Modern tools that generate genome-wide sequence data are now allowing us to decipher how evolutionary processes affect the spatial distribution of different kinds of genes and also to better model future spatial distribution of species in response to climate change.

Comparative phylogeography clarifies the complexity and problems of continental distribution that drove A. R. Wallace to favor islands

Deciphering the geographic context of diversification and distributional dynamics in continental biotas has long been an interest of biogeographers, ecologists, and evolutionary biologists. Thirty years ago, the approach now known as comparative phylogeography was introduced in a landmark study of a continental biota. Here, I use a set of 455 studies to explore the current scope of continental comparative phylogeography, including geographic, conceptual, temporal, ecological, and genomic attributes. Geographically, studies are more frequent in the northern hemisphere, but the south is catching up. Most studies focus on a Quaternary timeframe, but the Neogene is well represented. As such, explanations for geographic structure and history include geological and climatic events in Earth history, and responses include vicariance, dispersal, and range contraction-expansion into and out of refugia. Focal taxa are biased toward terrestrial or semiterrestrial vertebrates, although plants and invertebrates are well represented in some regions. The use of various kinds of nuclear DNA markers is increasing, as are multiple locus studies, but use of organelle DNA is not decreasing. Species distribution models are not yet widely incorporated into studies. In the future, continental comparative phylogeographers will continue to contribute to erosion of the simple vicariance vs. dispersal paradigm, including exposure of the widespread nature of temporal pseudocongruence and its implications for models of diversification; provide new templates for addressing a variety of ecological and evolutionary traits; and develop closer working relationships with earth scientists and biologists in a variety of disciplines.

Effects of Mountain Uplift and Climatic Oscillations on Phylogeography and Species Divergence in Four Endangered Notopterygium Herbs

Mountain uplift and climatic fluctuations are important driving forces that have affected the geographic distribution and population dynamics history of organisms. However, it is unclear how geological and climatic events might have affected the phylogeographic history and species divergence in high-alpine herbal plants. In this study, we analyzed the population demographic history and species differentiation of four endangered Notopterygium herbs on the high-altitude Qinghai-Tibetan Plateau (QTP) and adjacent areas. We combined phylogeographic analysis with species distribution modeling to detect the genetic variations in four Notopterygium species (N. incisum, N. franchetii, N. oviforme, and N. forrestii). In total, 559 individuals from 74 populations of the four species were analyzed based on three maternally inherited chloroplast fragments (matK, rbcL, and trnS-trnG) and one nuclear DNA region (internal transcribed spacer, ITS). Fifty-five chloroplast DNA (cpDNA) and 48 ITS haplotypes were identified in the four species. All of the cpDNA and ITS haplotypes were species-specific, except N. franchetii and N. oviforme shared one cpDNA haplotype, H32. Phylogenetic analysis suggested that all four species formed a monophyletic clade with high bootstrap support, where N. franchetii and N. oviforme were sisters. In addition, each Notopterygium species generated an individual clade that corresponded to their respective species in the ITS tree. Population dynamics analyses and species distribution modeling showed that the two widely distributed herbs N. incisum and N. franchetii exhibited obvious demographic expansions during the Pleistocene ice ages. Molecular dating suggested that the divergence of the four Notopterygium species occurred approximately between 3.6 and 1.2 Mya, and it was significantly associated with recent extensive uplifts of the QTP. Our results support the hypothesis that mountain uplift and Quaternary climatic oscillations profoundly shaped the population genetic divergence and demographic dynamics of Notopterygium species. The findings of this and previous studies provide important insights into the effects of QTP uplifts and climatic changes on phylogeography and species differentiation in high altitude mountainous areas. Our results may also facilitate the conservation of endangered herbaceous medicinal plants in the genus Notopterygium.

A geographic cline in the ability to self-fertilize is unrelated to the pollination environment

The reproductive assurance (RA) hypothesis predicts that the ability to autonomously self-fertilize should be favored in environments where a lack of mates or pollinators limits outcross reproduction. Because such limits to outcrossing are predicted to be most severe at range edges, elevated autonomy in peripheral populations is often attributed to RA. We test this hypothesis in 24 populations spanning the range of Campanula americana, including sampling at the range interior and three geographic range edges. We scored autonomous fruit set in a pollinator-free environment and detected clinal variation-autonomy increased linearly from the southern to the northern edge, and from the eastern to the western edge. We then address whether the cline reflects the contemporary pollination environment. We measured population size, plant density, pollinator visitation, outcross pollen limitation and RA in natural populations over two years. Most populations were pollen limited, and those that experienced higher visitation rates by bumblebees had reduced pollen limitation. Reproductive assurance, however, was generally low across populations and was unrelated to pollen limitation or autonomy. Neither pollen limitation nor RA displayed geographic clines. Finally, autonomy was not associated with pollinator visitation rates or mate availability. Thus, the data do not support the RA hypothesis; clinal variation in autonomy is unrelated to the current pollination environment. Therefore, geographic patterns of autonomy are likely the result of historical processes rather than contemporary natural selection for RA.

Population genomics of the eastern cottonwood (Populus deltoides)

Despite its economic importance as a bioenergy crop and key role in riparian ecosystems, little is known about genetic diversity and adaptation of the eastern cottonwood (Populus deltoides). Here, we report the first population genomics study for this species, conducted on a sample of 425 unrelated individuals collected in 13 states of the southeastern United States. The trees were genotyped by targeted resequencing of 18,153 genes and 23,835 intergenic regions, followed by the identification of single nucleotide polymorphisms (SNPs). This natural P. deltoides population showed low levels of subpopulation differentiation (F_ST = 0.022–0.106), high genetic diversity (θ_W = 0.00100, π = 0.00170), a large effective population size (N_e ≈ 32,900), and low to moderate levels of linkage disequilibrium. Additionally, genomewide scans for selection (Tajima's D), subpopulation differentiation (X^TX), and environmental association analyses with eleven climate variables carried out with two different methods (LFMM and BAYENV2) identified genes putatively involved in local adaptation. Interestingly, many of these genes were also identified as adaptation candidates in another poplar species, Populus trichocarpa, indicating possible convergent evolution. This study constitutes the first assessment of genetic diversity and local adaptation in P. deltoides throughout the southern part of its range, information we expect to be of use to guide management and breeding strategies for this species in future, especially in the face of climate change.

Integrative taxonomy of the ornamental 'peppermint' shrimp public market and population genetics of Lysmata boggessi, the most heavily traded species worldwide

The ornamental trade is a worldwide industry worth >15 billion USD with a problem of rampant product misidentification. Minimizing misidentification is critical in the face of overexploitation of species in the trade. We surveyed the peppermint shrimp ornamental marketplace in the southeastern USA, the most intense market for peppermint shrimps worldwide, to characterize the composition of species in the trade, reveal the extent of misidentification, and describe the population genetics of the true target species. Shrimps were bought from aquarium shops in FL, GA, SC, and NC. We demonstrated, contrary to popular belief (information from dealers), that the most heavily traded species in the market was Lysmata boggessi, an endemic species to the eastern Gulf of Mexico, and not Lysmata wurdemanni. Importantly, only when color pattern or genetic markers in conjunction with morphological traits were employed, was it was possible to unequivocally identify L. boggessi as the only species in the trade. The intensity of the market for peppermint shrimps in the USA has led to L. boggessi being the most traded species worldwide. Misidentification in the shrimp aquarium trade is accidental and involuntary, and is explained by remarkable similarity among congeneric species. Using sequences of the 16S-mt-DNA marker, we found no indication of population genetic structure in the endemic L. boggessi across  550 km of linear coast. Therefore, this species can be considered genetically homogeneous and a single fished stock. Still, we argue in favor of additional studies using more powerful markers (e.g., SNPs) capable of revealing genetic structure at a finer spatial-scale. Our results will help advance management and conservation policies in this lucrative yet understudied fishery. Future studies of other ornamental fisheries will benefit from using an integrative taxonomic approach, as we demonstrate here.

Ice ages leave genetic diversity 'hotspots' in Europe but not in Eastern North America

After the last glacial cycle, temperate European trees migrated northward, experiencing genetic bottlenecks and founder effects, which left high haplotype endemism in southern populations and clines in genetic diversity northward. These patterns are thought to be ubiquitous across temperate forests, and are therefore used to anticipate the potential genetic consequences of future warming. We compared existing and new phylogeographic data sets (chloroplast DNA) from 14 woody taxa in Eastern North America (ENA) to data sets from 21 ecologically similar European species to test for common impacts of Quaternary climate swings on genetic diversity across diverse taxa and between continents. Unlike their European counterparts, ENA taxa do not share common southern centres of haplotype endemism and they generally maintain high genetic diversity even at their northern range limits. Differences between the genetic impacts of Quaternary climate cycles across continents suggest refined lessons for managing genetic diversity in today's warming world.

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

Phylogeographic structure, cryptic speciation and demographic history of the sharpbelly (Hemiculter leucisculus), a freshwater habitat generalist from southern China

BACKGROUND

Species with broad distributions frequently divide into multiple genetic forms and may therefore be viewed as "cryptic species". Here, we used the mitochondrial cytochrome b (Cytb) and 12 nuclear DNA loci to investigate phylogeographic structures of the sharpbelly (Hemiculter leucisculus) in rivers in southern China and explored how the geological and climatic factors have shaped the genetic diversity and evolutionary history of this species.

RESULTS

Our mitochondrial phylogenetic analysis identified three major lineages (lineages A, B, and C). Lineages B and C showed a relatively narrower geographic distribution, whereas lineage A was widely distributed in numerous drainages. Divergence dates suggested that H. leucisculus populations diverged between 1.61-2.38 Ma. Bayesian species delimitation analysis using 12 nuclear DNA loci indicated the three lineages probably represented three valid taxa. Isolation-with-migration (IM) analysis found substantial gene flow has occurred among the three lineages. Demographic analyses showed that lineages B and C have experienced rapid demographic expansion at 0.03 Ma and 0.08 Ma, respectively.

CONCLUSIONS

Hemiculter leucisculus populations in drainages in southern China comprise three mtDNA lineages, and each of which may represent a separate species. Intense uplift of the Qinghai-Tibetan Plateau, evolution of Asian monsoons, changes in paleo-drainages, and poor dispersal ability may have driven the divergence of the three putative species. However, gene flow occurs among the three lineages. Climatic fluctuations have a prominent impact on the populations from the lineages B and C, but exerted little influence on the lineage A.

Habitat Association Predicts Genetic Diversity and Population Divergence in Amazonian Birds

The ecological traits of organisms may predict their genetic diversity and population genetic structure and mediate the action of evolutionary processes important for speciation and adaptation. Making these ecological-evolutionary links is difficult because it requires comparable genetic estimates from many species with differing ecologies. In Amazonian birds, habitat association is an important component of ecological diversity. Here, we examine the link between habitat association and genetic parameters using 20 pairs of closely related Amazonian bird species in which one member of the pair occurs primarily in forest edge and floodplains and the other occurs in upland forest interior. We use standardized geographic sampling and data from 2,416 genomic markers to estimate genetic diversity, population genetic structure, and statistics reflecting demographic and evolutionary processes. We find that species of upland forest have greater genetic diversity and divergence across the landscape as well as signatures of older histories and less gene flow than floodplain species. Our results reveal that species ecology in the form of habitat association is an important predictor of genetic diversity and population divergence and suggest that differences in diversity between floodplain and upland avifaunas in the Amazon may be driven by differences in the demographic and evolutionary processes at work in the two habitats.

A time and a place for everything: phylogenetic history and geography as joint predictors of oak plastome phylogeny

Owing to high rates of introgressive hybridization, the plastid genome is poorly suited to fine-scale DNA barcoding and phylogenetic studies of the oak genus (Quercus, Fagaceae). At the tips of the oak plastome phylogeny, recent gene migration and reticulation generally cause topology to reflect geographic structure, while deeper branches reflect lineage divergence. In this study, we quantify the simple and partial effects of geographic proximity and nucleome-inferred phylogenetic history on oak plastome phylogeny at different evolutionary scales. Our study compares pairwise phylogenetic distances based on complete plastome sequences, pairwise phylogenetic distances from nuclear restriction site-associated DNA sequences (RADseq), and pairwise geographic distances for 34 individuals of the white oak clade representing 24 North American and Eurasian species. Within the North American white oak clade alone, phylogenetic history has essentially no effect on plastome variation, while geography explains 11%–21% of plastome phylogenetic variance. However, across multiple continents and clades, phylogeny predicts 30%–41% of plastome variation, geography 3%–41%. Tipwise attenuation of phylogenetic informativeness in the plastome means that in practical terms, plastome data has little use in solving phylogenetic questions, but can still be a useful barcoding or phylogenetic marker for resolving questions among major clades.

Turtle soup, Prohibition, and the population genetic structure of Diamondback Terrapins (Malaclemys terrapin)

Diamondback terrapins (Malaclemys terrapin) were a popular food item in early twentieth century America, and were consumed in soup with sherry. Intense market demand for terrapin meat resulted in population declines, notably along the Atlantic seaboard. Efforts to supply terrapins to markets resulted in translocation events, as individuals were moved about to stock terrapin farms. However, in 1920 the market for turtle soup buckled with the enactment of the eighteenth amendment to the United States’ Constitution—which initiated the prohibition of alcoholic drinks—and many terrapin fisheries dumped their stocks into local waters. We used microsatellite data to show that patterns of genetic diversity along the terrapin’s coastal range are consistent with historical accounts of translocation and cultivation activities. We identified possible instances of human-mediated dispersal by estimating gene flow over historical and contemporary timescales, Bayesian model testing, and bottleneck tests. We recovered six genotypic clusters along the Gulf and Atlantic coasts with varying degrees of admixture, including increased contemporary gene flow from Texas to South Carolina, from North Carolina to Maryland, and from North Carolina to New York. In addition, Bayesian models incorporating translocation events outperformed stepping-stone models. Finally, we were unable to detect population bottlenecks, possibly due to translocation reintroducing genetic diversity into bottlenecked populations. Our data suggest that current patterns of genetic diversity in the terrapin were altered by the demand for turtle soup followed by the enactment of alcohol prohibition. In addition, our study shows that population genetic tools can elucidate metapopulation dynamics in taxa with complex genetic histories impacted by anthropogenic activities.

Phylogenetics of Kingsnakes, Lampropeltis getula Complex (Serpentes: Colubridae), in Eastern North America

Kingsnakes of the Lampropeltis getula complex range throughout much of temperate and subtropical North America. Studies over the last century have used morphology and color pattern to describe numerous subspecies. More recently, DNA analyses have made invaluable contributions to our understanding of their evolution and taxonomy. We use genetic and ecological methods to test previous hypotheses of distinct evolutionary lineages by examining 66 total snakes and 1) analyzing phylogeographic structure using 2 mtDNA loci and 1 nuclear locus, 2) estimating divergence dates and historical demography among lineages in a Bayesian coalescent framework, and 3) applying ecological niche modeling (ENM). Our molecular data and ENMs illustrate that 3 previously recognized subspecies in the eastern United States comprise well-supported monophyletic lineages that diverged during the Pleistocene. The geographic boundaries of these 3 lineages correspond closely to known biogeographic barriers (Florida peninsula, Appalachian Mountains, and Apalachicola River) previously identified for other plants and animals, indicating shared geographic influences on evolutionary history. We conclude that genetic, ecological, and morphological data support recognition of these 3 lineages as distinct species (Lampropeltis floridana, Lampropeltis getula, and Lampropeltis meansi).