Comparative phylogeography of unglaciated eastern North America

Temporal and spatial comparisons of angiosperm diversity between eastern Asia and North America

Eastern Asia (EA) and North America north of Mexico (NA) have comparable latitude, land area, and climate, but the overall plant diversity is much higher in EA than in NA. Despite intensive studies on disjunct taxa of the two regions, the temporal and spatial diversity patterns between the two floras remain unclear. Here we explore the floristic differences between EA and NA using the well-studied floras of China and the United States of America (USA) as exemplars, while also employing a newly generated dated phylogeny covering ∼90% of the angiosperm genera of the two countries and comprehensive spatial distribution data. We find that China possesses both higher richness and phylogenetic diversity (PD) for angiosperm genera than the USA. Notably, most lineages contribute to the PD anomaly between the two floras, with 46 of 58 lineages having higher PD in China. Temporally, China has a higher proportion of genera that originated before the Miocene than found in the USA (29.9% vs 23.2%). The eastern USA has more genera that originated during the Paleogene than does the western USA, but the reverse pattern is observed after the middle Miocene, with more genera originating in the west. Spatially, China shows a more distinct east-west deviation in diversity than the USA with eastern China possessing much higher generic richness and PD and more ancient lineages than western China. However, the eastern USA possesses lower generic richness, but higher PD and more ancient lineages than the western USA. Both the floras in China and the USA share a signature of an older east and a younger west, and this pattern may be largely driven by regional orogenic activities and climatic changes in the west of the two regions. Finally, our study indicates that more efforts are needed to enhance biodiversity conservation in southern China and the eastern USA by identifying and protecting phylodiversity hotspots.

Ecological Divergence and the History of Gene Flow in the Nearctic Milksnakes (Lampropeltis triangulum Complex)

Many phylogeographic studies on species with large ranges have found genetic-geographic structure associated with changes in habitat and physical barriers preventing or reducing gene flow. These interactions with geographic space, contemporary and historical climate, and biogeographic barriers have complex effects on contemporary population genetic structure and processes of speciation. While allopatric speciation at biogeographic barriers is considered the primary mechanism for generating species, more recently it has been shown that parapatric modes of divergence may be equally or even more common. With genomic data and better modeling capabilities, we can more clearly define causes of speciation in relation to biogeography and migration between lineages, the location of hybrid zones with respect to the ecology of parental lineages, and differential introgression of genes between taxa. Here, we examine the origins of three Nearctic milksnakes (Lampropeltis elapsoides, Lampropeltis triangulum and Lampropeltis gentilis) using genome-scale data to better understand species diversification. Results from artificial neural networks show that a mix of a strong biogeographic barrier, environmental changes, and physical space has affected genetic structure in these taxa. These results underscore conspicuous environmental changes that occur as the sister taxa L. triangulum and L. gentilis diverged near the Great Plains into the forested regions of the Eastern Nearctic. This area has been recognized as a region for turnover for many vertebrate species, but as we show here the contemporary boundary does not isolate these sister species. These two species likely formed in the mid-Pleistocene and have remained partially reproductively isolated over much of this time, showing differential introgression of loci. We also demonstrate that when L. triangulum and L. gentilis are each in contact with the much older L. elapsoides, some limited gene flow has occurred. Given the strong agreement between nuclear and mtDNA genomes, along with estimates of ecological niche, we suggest that all three lineages should continue to be recognized as unique species. Furthermore, this work emphasizes the importance of considering complex modes of divergence and differential allelic introgression over a complex landscape when testing mechanisms of speciation. [Cline; delimitation; Eastern Nearctic; Great Plains; hybrids; introgression; speciation.].

Genetic diversity in North American Cercis Canadensis reveals an ancient population bottleneck that originated after the last glacial maximum

Understanding of the present-day genetic diversity, population structure, and evolutionary history of tree species can inform resource management and conservation activities, including response to pressures presented by a changing climate. Cercis canadensis (Eastern Redbud) is an economically valuable understory tree species native to the United States (U.S.) that is also important for forest ecosystem and wildlife health. Here, we document and explain the population genetics and evolutionary history of this deciduous tree species across its distributed range. In this study, we used twelve microsatellite markers to investigate 691 wild-type trees sampled at 74 collection sites from 23 Eastern U.S. states. High genetic diversity and limited gene flow were revealed in wild, natural stands of C. canadensis with populations that are explained by two major genetic clusters. These findings indicate that an ancient population bottleneck occurred coinciding with the last glacial maximum (LGM) in North America. The structure in current populations likely originated from an ancient population in the eastern U.S. that survived LGM and then later diverged into two contemporary clusters. Data suggests that populations have expanded since the last glaciation event from one into several post-glacial refugia that now occupy this species’ current geographic range. Our enhanced understanding benchmarks the genetic variation preserved within this species and can direct future efforts in conservation, and resource utilization of adaptively resilient populations that present the greatest genetic and structural diversity.

Weak spatial-genetic structure in a native invasive, the southern pine beetle (Dendroctonus frontalis), across the eastern United States

The southern pine beetle, Dendroctonus frontalis, is a native pest of pine trees that has recently expanded its range into the northeastern United States. Understanding its colonization, dispersal, and connectivity will be critical for mitigating negative economic and ecological impacts in the newly invaded areas. Characterization of spatial-genetic structure can contribute to this; however, previous studies have reached different conclusions about regional population genetic structure, with one study reporting a weak east-west pattern, and the most recent reporting an absence of structure. Here we systematically assessed several explanations for the absence of spatial-genetic structure. To do this, we developed nine new microsatellite markers and combined them with an existing 24-locus data matrix for the same individuals. We then reanalyzed this full dataset alongside datasets in which certain loci were omitted with the goal of creating more favorable signal to noise ratios. We also partitioned the data based on the sex of D. frontalis individuals, and then employed a broad suite of genotypic clustering and isolation-by-distance (IBD) analyses. We found that neither inadequate information content in the molecular marker set, nor unfavorable signal-to-noise ratio, nor insensitivity of the analytical approaches could explain the absence of structure. Regardless of dataset composition, there was little evidence for clusters (i.e., distinct geo-genetic groups) or clines (i.e., gradients of increasing allele frequency differences over larger geographic distances), with one exception: significant IBD was repeatedly detected using an individual-based measure of relatedness whenever datasets included males (but not for female-only datasets). This is strongly indicative of broad-scale female-biased dispersal, which has not previously been reported for D. frontalis, in part owing to logistical limitations of direct approaches (e.g., capture-mark-recapture). Weak spatial-genetic structure suggests long-distance connectivity and that gene flow is high, but additional research is needed to understand range expansion and outbreak dynamics in this species using alternate approaches.

Population dynamics of caribou shaped by glacial cycles before the last glacial maximum

Pleistocene glacial cycles influenced the diversification of high‐latitude wildlife species through recurrent periods of range contraction, isolation, divergence, and expansion from refugia and subsequent admixture of refugial populations. We investigate population size changes and the introgressive history of caribou (Rangifer tarandus) in western Canada using 33 whole genome sequences coupled with larger‐scale mitochondrial data. We found that a major population expansion of caribou occurred starting around 110,000 years ago (kya), the start of the last glacial period. Additionally, we found effective population sizes of some caribou reaching ~700,000 to 1,000,000 individuals, one of the highest recorded historical effective population sizes for any mammal species thus far. Mitochondrial analyses dated introgression events prior to the LGM dating to 20–30 kya and even more ancient at 60 kya, coinciding with colder periods with extensive ice coverage, further demonstrating the importance of glacial cycles and events prior to the LGM in shaping demographic history. Reconstructing the origins and differential introgressive history has implications for predictions on species responses under climate change. Our results have implications for other whole genome analyses using pairwise sequentially Markovian coalescent (PSMC) analyses, as well as highlighting the need to investigate pre‐LGM demographic patterns to fully reconstruct the origin of species diversity, especially for high‐latitude species.

Historical biogeography of the genus Rhadinaea (Squamata: Dipsadinae)

Multiple geological and climatic events have created geographical or ecological barriers associated with speciation events, playing a role in biological diversification in North and Central America. Here, we evaluate the influence of the Neogene and Quaternary geological events, as well as the climatic changes in the diversification of the colubrid snake genus Rhadinaea using molecular dating and ancestral area reconstruction. A multilocus sequence dataset was generated for 37 individuals of Rhadinaea from most of the biogeographical provinces where the genus is distributed, representing 19 of the 21 currently recognized species, and two undescribed species. Our analyses show that the majority of the Rhadinaea species nest in two main clades, herein identified as "Eastern" and "Southern". These clades probably diverged from each other in the early Miocene, and their divergence was followed by 11 divergences during the middle to late Miocene, three divergences during the Pliocene, and six divergences in the Pleistocene. The ancestral distribution of Rhadinaea was reconstructed across the Sierra Madre del Sur. Our phylogenetic analyses do not support the monophyly of Rhadinaea. The Miocene and Pliocene geomorphology, perhaps in conjunction with climate change, appears to have triggered the diversification of the genus, while the climatic changes during the Miocene probably induced the diversification of Rhadinaea in the Sierra Madre del Sur. Our analysis suggests that the uplifting of the Trans-Mexican Volcanic Belt and Chiapan-Guatemalan highlands in this same period resulted in northward and southward colonization events. This was followed by more recent, independent colonization events in the Pliocene and Pleistocene involving the Balsas Basin, Chihuahuan Desert, Pacific Coast, Sierra Madre Occidental, Sierra Madre Oriental, Sierra Madre del Sur, Trans-Mexican Volcanic Belt, and Veracruz provinces, probably driven by the climatic fluctuations of the time.

Evolution of Castanea in North America: restriction-site-associated DNA sequencing and ecological modeling reveal a history of radiation, range shifts, and disease

PREMISE

Although chestnuts and chinquapins are some of the best known and most widely loved of any plants in North America, relatively little genomic sequencing has been done, and much is still unknown about their evolution.

METHODS

We used double-digest restriction-site-associated DNA (ddRAD) sequencing data to infer the species-level phylogeny for Castanea and assess the phylogeography of the North American species using samples collected from populations that span the full extent of the species' ranges. We also constructed species distribution models using digitized herbarium specimens and observational data from field surveys.

RESULTS

We identified strong population structure within Castanea dentata (American chestnut) that reflects a stepwise northern migration since the last glacial maximum. Our species distribution models further confirmed this scenario and matched closely with the Castanea fossil pollen record. We also found significant structure within the Castanea pumila lineage, most notably a genetic cluster that corresponds to the frequently recognized Castanea pumila var. ozarkensis.

CONCLUSIONS

The two North American Castanea species have contrasting patterns of population structure, but each is typical of plant phylogeography in North America. Within the C. pumila complex, we found novel genetic structure that provides new insights about C. pumila taxonomy. Our results also identified a series of distinctive populations that will be valuable in ongoing efforts to conserve and restore chestnuts and chinquapins in North America.

Genomic signatures of inbreeding and mutation load in a threatened rattlesnake

Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

Towards a unified framework to study causality in Earth-life systems

There is considerable interest in better understanding how earth processes shape the generation and distribution of life on Earth. This question, at its heart, is one of causation. In this article I propose that at a regional level, earth processes can be thought of as behaving somewhat deterministically and may have an organized effect on the diversification and distribution of species. However, the study of how landscape features shape biology is challenged by pseudocongruent or collinear variables. I demonstrate that causal structures can be used to depict the cause–effect relationships between earth processes and biological patterns using recent examples from the literature about speciation and species richness in montane settings. This application shows that causal diagrams can be used to better decipher the details of causal relationships by motivating new hypotheses. Additionally, the abstraction of this knowledge into structural equation metamodels can be used to formulate theory about relationships within Earth–life systems more broadly. Causal structures are a natural point of collaboration between biologists and Earth scientists, and their use can mitigate against the risk of misassigning causality within studies. My goal is that by applying causal theory through application of causal structures, we can build a systems‐level understanding of what landscape features or earth processes most shape the distribution and diversification of species, what types of organisms are most affected, and why.

Hybridization, Behavioral Patterns, and Pre- and Postzygotic Isolation Between Two Recently Derived Species of Walnut-Infesting Rhagoletis Fruit Flies in the Highlands of Mexico

Episodes of isolation and secondary contact among populations of insects of Nearctic origin during Pleistocene glacial/postglacial climatic cycles had a strong evolutionary influence on the diversity of flies in the genus Rhagoletis in mountainous areas of Mexico. As a series of experiments undertaken to gather support for phylogenetic hypotheses on the origin of three walnut-infesting species in the suavis group, we examined pre- and postzygotic isolation between Rhagoletis completa Cresson, 1929 and R. ramosae Hernández-Ortiz, 1985. Despite morphological, biological, and behavioral differences, these two species were found to be capable of hybridization. Mating experiments in large enclosures revealed asymmetric sexual isolation. There were notable differences in male sexual behavior. While R. ramosae males mated exclusively on host fruit, R. completa males used fruit and alternative mating locations. During fruit-guarding and male-male contests, R. completa and R. ramosae males adopted markedly different wing postures. R. completa females were more reluctant to copulate with heterospecific males than R. ramosae females. During no choice crosses in small enclosures, there was a reduction of egg hatch for the hybrid cross of R. completa males × R. ramosae females. Our results and previous studies on reproductive isolation between other species pairs in the suavis group support a clade in which R. ramosae, R. zoqui Bush, 1966, and R. completa are close relatives all still capable of hybridizing.

Shennongjia-Wushan Mountains-One cryptic glacial refugium introduced by the phylogeographical study of the Geometridae moth Ourapteryx szechuana Wehrli

The origin and evolution of biodiversity in the Shennongjia and Wushan Mountains, located in central China, are little known. In this study, we used Ourapteryx szechuana, which is widely distributed in China and northern Nepal, to explore whether these mountains acted as glacial refugia during climate oscillations of the Quaternary. In total, 192 samples of O. szechuana were collected throughout much of the distribution range. Phylogenetic analysis, molecular dating, demographic history reconstructions, and MAXENT were used to investigate the evolutionary history and differentiation mechanisms and predict the potential species distributions during four different periods. The phylogenetic tree and the star-like median-joining network strongly supported two reciprocally monophyletic and allopatric lineages. Lineage I was restricted to the Shennongjia and Wushan Mountains. The divergence time of O. szechuana from its sister species O. thibetaria was approximately 1.94 Ma. The differentiation processes of the two intraspecific lineages occurred at approximately 0.47 Ma. The demographic history reconstruction and the ecological niche model suggested that Lineage II experienced an expansion after the LGM (Last Glacial Maximum), whereas Lineage I did not experience any expansion. Our results suggested the Naynayxungla glaciation promoted the divergence of the two lineages by restricting them to different refugia. The valleys of the Shennongjia-Wushan Mountains may have kept stable and warm (thus ice-free) environments during Quaternary glaciations, allowing this region to act as a glacial refugia. Our studies show that the Shennongjia and Wushan Mountains are likely to be important but little studied glacial refugia for the insect and thus worthy of more attention.

Local causation

The counterfactual and regularity theories are universal accounts of causation. I argue that these should be generalized to produce local accounts of causation. A hallmark of universal accounts of causation is the assumption that apparent variation in causation between locations must be explained by differences in background causal conditions, by features of the causal-nexus or causing-complex. The local account of causation presented here rejects this assumption, allowing for genuine variation in causation to be explained by differences in location. I argue that local accounts of causation are plausible, and have pragmatic, empirical and theoretical advantages over universal accounts. I then report on the use of presheaves as models of local causation. The use of presheaves as models of local variation has precedents in algebraic geometry, category theory and physics; they are here used as models of local causal variation. The paper presents this idea as stemming from an approach using presheaves as models of local truth. Finally, I argue that a proper balance between universal and local causation can be assuaged by moving from presheaves to fully-fledged sheaf models.

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.

A new subterranean species of Anillinus Casey (Carabidae, Trechinae, Anillini) from Florida

A new species of blind subterranean ground beetle in the genus Anillinus Casey is described from Florida. Anillinus albrittonorumsp. nov. (type locality: 6 miles NW High Springs, Columbia County, Florida) has a unique structure of female genitalia and occupies an isolated position within the genus. This new species is illustrated with images of the habitus, body parts, and male and female genitalia. Relationships of A. albrittonorum to other members of the genus are discussed.

The phylogeographic history of a range disjunction in eastern North America: the role of post-glacial expansion into newly suitable habitat

PREMISE

A disjunct distribution, where a species' geographic range is discontinuous, can occur through vicariance or long-distance dispersal. Approximately 75 North American plant species exhibit a ~650 km disjunction between the Ozark and Appalachian regions. This disjunction is attributed to biogeographic forces including: (1) Eocene-Oligocene vicariance by the formation of the Mississippi embayment; (2) Pleistocene vicariance from interglacial flooding; (3) post-Pleistocene northward colonization from separate glacial refugia; (4) Hypsithermal vicariance due to climate fluctuations; and (5) recent long-distance dispersal. We investigated which of these pathways most likely gave rise to the Appalachian-Ozark disjunction in Delphinium exaltatum.

METHODS

We genotyped populations of D. exaltatum from five Ozark and seven Appalachian localities, analyzed genetic structure, tested the order and timing of divergences using DIYABC, and conducted niche reconstructions up to 21,000 years before present (YBP).

RESULTS

Populations fell into five main genetic clusters, i.e., a group in the central Appalachians, and four "lowland" groups. DIYABC analyses showed the central Appalachian and lowland lineages diverging 11,300 to17,000 YBP, and the lowland groups diverging 6800 to 10,900 YBP. Niche reconstructions showed that suitable climate for the central Appalachian lineage experienced large spatial discontinuity starting 14,000 YBP, such that divergence and persistence before this period is less plausible than divergence thereafter.

CONCLUSIONS

Our results did not fully support any of the original hypotheses. Rather, the oldest divergence likely occurred after 13,500 YBP through expansion into newly opened habitat in the Appalachians. The Appalachian-Ozark disjunction likely resulted from northward dispersal of the lowland lineage as climate warmed during the Holocene.

Genomic and phenotypic divergence informs translocation strategies for an endangered freshwater fish

Translocation, the movement of organisms for conservation purposes, can result in unintended introgression if genetic material flows between populations in new ways. The Bluemask Darter Etheostoma akatulo is a federally endangered species of freshwater fish inhabiting the Caney Fork River system and three of its tributaries (Collins River, Rocky River, and Cane Creek) in Tennessee. The current conservation strategy for Bluemask Darters involves translocating the progeny of broodstock from the Collins River (in the west) to the Calfkiller River (in the east) where the species had been extirpated. In this study, we use ddRAD sequence data from across the extant range to assess this translocation strategy in light of population structure, phylogeny, and demography. We also include museum specimen data to assess morphological variation among extant and extirpated populations. Our analyses reveal substantial genetic and phenotypic disparities between a western population in the Collins River and an eastern population encompassing the Rocky River, Cane Creek, and upper Caney Fork, the two of which shared common ancestry more than 100,000 years ago. Furthermore, morphological analyses classify 12 of 13 Calfkiller River specimens with phenotypes consistent with the eastern population. These results suggest that current translocations perturb the evolutionary boundaries between two delimited populations. Instead, we suggest that repopulating the Calfkiller River using juveniles from the Rocky River could balance conflicting signatures of demography, diversity, and divergence. Beyond conservation, the microgeographic structure of Bluemask Darter populations adds another puzzle to the phylogeography of the hyperdiverse freshwater fishes in eastern North America.

Pliocene-Early Pleistocene Geological Events Structure Pacific Martens (Martes caurina)

The complex topography, climate, and geological history of Western North America have shaped contemporary patterns of biodiversity and species distributions in the region. Pacific martens (Martes caurina) are distributed along the northern Pacific Coast of North America with disjunct populations found throughout the Northwestern Forested Mountains and Marine West Coast Forest ecoregions of the West Coast. Martes in this region have been classified into subspecies; however, the subspecific designation has been extensively debated. In this study, we use genomic data to delineate conservation units of Pacific marten in the Sierra-Cascade-Coastal montane belt in the western United States. We analyzed the mitochondrial genome for 94 individuals to evaluate the spatial distribution and divergence times of major lineages. We further genotyped 401 individuals at 13 microsatellite loci to investigate major patterns of population structure. Both nuclear and mitochondrial DNA suggest substantial genetic substructure concordant with historical subspecies designations. Our results revealed that the region contains 2 distinct mitochondrial lineages: a Cascades/Sierra lineage that diverged from the Cascades/coastal lineage 2.23 (1.48-3.14 mya), consistent with orogeny of the Cascade Mountain chain. Interestingly, Pacific Martes share phylogeographic patterns similar with other sympatric taxa, suggesting that the complex geological history has shaped the biota of this region. The information is critical for conservation and management efforts, and further investigation of adaptive diversity is warranted following appropriate revision of conservation management designations.

Integrated analysis of the ringneck snake Diadophis punctatus complex (Colubridae: Dipsadidae) in a biodiversity hotspot provides the foundation for conservation reassessment

Species classification may not reflect the underlying/cryptic genetic diversity and focusing on groups that do not represent historically independent units can misdirect conservation efforts. The identification of evolutionarily significant units (ESUs) allows cryptic genetic diversity to be accounted for when designating conservation priorities. We used multi-locus coalescent-based species delimitation methods and multivariate analyses of morphological data to examine whether the subspecies merit conservation recognition and infer the ESUs in ringneck snakes (Diadophis punctatus) throughout the California Floristic Province. Species delimitation methods failed to recover groups consistent with designated subspecies and instead inferred three well supported, mostly geographically isolated lineages. Divergence time estimates suggest that the divergences were driven by historical isolation associated with Pleistocene climate shifts. We found a correlation between increased morphological differentiation and time since divergence, and greater niche similarity between the more recently diverged eastern California and western California groups. Based on these results, we propose that the morphological similarities are due to a combination of morphological conservatism and evolutionary stasis. Our study provides the foundation necessary to re-assess the biodiversity and conservation status of ringneck snakes and offers an important step in unveiling the diversity within the western portion of the genus’ range.

Genetic diversity and population structure in Vitis species illustrate phylogeographic patterns in eastern North America

Geographical distribution and diversity of current plant species have been strongly shaped by climatic oscillations during the Quaternary. Analysing the resulting divergence among species and differentiation within species is crucial to understand the evolution of taxa like the Vitis genus, which provides very useful genetic resources for grapevine improvement and might reveal original recolonization patterns due to growth habit and dispersal mode. Here, we studied the genetic structure in natural populations of three species from eastern North America: Vitis aestivalis, V. cinerea and V. riparia using different marker types. Vitis aestivalis and V. cinerea showed higher diversity than V. riparia. The two former species are less differentiated, confirming an earlier divergence of V. riparia. V. aestivalis and V. riparia exhibited different genetic groups on both sides of the Appalachian Mountains that could mirror different recolonization routes from southern refugia. Genetic structure was stronger in V. cinerea, for which two varieties (var. berlandieri and var. cinerea) are morphologically recognized. Our results confirm this distinction and suggest the existence of three other lineages within var. cinerea. These discontinuities appear linked to adaptation of var. berlandieri to dry and limy areas of Texas and partially to the Mississippi River Valley. Rapid range expansions from refugia upon climate warming are also suggested by the low linkage disequilibrium values observed. Furthermore, large variation for downy mildew resistance was observed in the three species. Our findings appeared consistent with the vegetation history of eastern North America.

Signatures of north-eastern expansion and multiple refugia: genomic phylogeography of the Pine Barrens tree frog, Hyla andersonii (Anura: Hylidae)

Range fragmentation poses challenges for species persistence over time and can be caused by both historical and contemporary processes. We combined genomic data, phylogeographical model testing and palaeoclimatic niche modelling to infer the evolutionary history of the Pine Barrens tree frog (Hyla andersonii), a seepage bog specialist, in eastern North America to gain a better understanding of the historical context of its fragmented distribution. We sampled H. andersonii populations across the three disjunct regions of the species range: Alabama/Florida (AF), the Carolinas (CL) and New Jersey (NJ). Phylogenetic relationships within H. andersonii were consistent between the nuclear species tree and mitochondrial analyses, indicating divergence between AF and CL/NJ (Atlantic clade) ~0.9 Mya and divergence of the NJ clade ~0.15 Mya. Several predictions of north-eastern expansion along the Atlantic coast were supported by phylogeographical analyses. Model testing using genome-wide single nucleotide polymorphism data and species distribution models both provided evidence for multiple disjunct refugia. This comprehensive phylogeographical study of H. andersonii demonstrates a long history of range fragmentation within an endemic coastal plain species and highlights the influence of historical climate change on the current distribution of species and their genetic diversity.

Range and niche expansion through multiple interspecific hybridization: a genotyping by sequencing analysis of Cherleria (Caryophyllaceae)

BACKGROUND

Cherleria (Caryophyllaceae) is a circumboreal genus that also occurs in the high mountains of the northern hemisphere. In this study, we focus on a clade that diversified in the European High Mountains, which was identified using nuclear ribosomal (nrDNA) sequence data in a previous study. With the nrDNA data, all but one species was monophyletic, with little sequence variation within most species. Here, we use genotyping by sequencing (GBS) data to determine whether the nrDNA data showed the full picture of the evolution in the genomes of these species.

RESULTS

The overall relationships found with the GBS data were congruent with those from the nrDNA study. Most of the species were still monophyletic and many of the same subclades were recovered, including a clade of three narrow endemic species from Greece and a clade of largely calcifuge species. The GBS data provided additional resolution within the two species with the best sampling, C. langii and C. laricifolia, with structure that was congruent with geography. In addition, the GBS data showed significant hybridization between several species, including species whose ranges did not currently overlap.

CONCLUSIONS

The hybridization led us to hypothesize that lineages came in contact on the Balkan Peninsula after they diverged, even when those lineages are no longer present on the Balkan Peninsula. Hybridization may also have helped lineages expand their niches to colonize new substrates and different areas. Not only do genome-wide data provide increased phylogenetic resolution of difficult nodes, they also give evidence for a more complex evolutionary history than what can be depicted by a simple, branching phylogeny.

Phylogeography and demographic expansion in the widely distributed horned passalus beetle, Odontotaenius disjunctus (coleoptera: Passalidae)

Dynamic climatic oscillations during the Pleistocene had profound effects on the distribution of species across North America. Although the role of historical climate change on speciation remains controversial, the impact on genetic variation within species has been well documented. Analyses of mtDNA sequences from the cytochrome oxidase I gene (911 bp) for 115 individuals of Odontotaenius disjunctus was combined with ecological niche modelling (ENM) to infer the demographic and population differentiation scenarios under present and past conditions. We inferred three lineages that diverged during the Pleistocene and replace each other geographically across the eastern United States. One of these lineages traverses previously identified genetic barriers for terrestrial animals including the Mississippi and Apalachicola Rivers and the Appalachian Mountains. We observed overlapping ranges between two haplotype groups as well as a region of secondary contact associated with ecological transition zone in northern Florida. The two continental lineages depict a genetic signature of a recent population increase associated with expanding niche envelope, whereas the clade restricted to peninsular Florida shows stable populations in a shrinking niche envelope. Given the lack of ecological separation, overlapping distribution of haplogroups and the presence of secondary contact zones, the taxonomic status of these lineages must await robust testing using multilocus DNA data to assess species boundaries.

The choices we make and the impacts they have: Machine learning and species delimitation in North American box turtles (Terrapene spp.)

Model-based approaches that attempt to delimit species are hampered by computational limitations as well as the unfortunate tendency by users to disregard algorithmic assumptions. Alternatives are clearly needed, and machine-learning (M-L) is attractive in this regard as it functions without the need to explicitly define a species concept. Unfortunately, its performance will vary according to which (of several) bioinformatic parameters are invoked. Herein, we gauge the effectiveness of M-L-based species-delimitation algorithms by parsing 64 variably-filtered versions of a ddRAD-derived SNP data set collected from North American box turtles (Terrapene spp.). Our filtering strategies included: (i) minor allele frequencies (MAF) of 5%, 3%, 1%, and 0% (= none), and (ii) maximum missing data per-individual/per-population at 25%, 50%, 75%, and 100% (= no filtering). We found that species-delimitation via unsupervised M-L impacted the signal-to-noise ratio in our data, as well as the discordance among resolved clades. The latter may also reflect biogeographic history, gene flow, incomplete lineage sorting, or combinations thereof (as corroborated from previously observed patterns of differential introgression). Our results substantiate M-L as a viable species-delimitation method, but also demonstrate how commonly observed patterns of phylogenetic discordance can seriously impact M-L-classification.

Canada lynx (Lynx canadensis) gene flow across a mountain transition zone in western North America

Mountain ecotones have the potential to cause multiple patterns in divergence, from simple barrier effects to more fundamental ecological divergence. Most work in mountain ecotones in North America has focused on reinforcement between refugial populations, making prediction of how mountains impact species that are not restricted to separate glacial refugia remains difficult. This study focused on the Canada lynx (Lynx canadensis Kerr, 1792), a highly mobile felid considered to be a habitat and dietary specialist. Specifically, we used 14 microsatellite loci and landscape genetic tools to investigate if the Rocky Mountains and associated climatic transitions influence lynx genetic differentiation in western North America. Although lynx exhibited high gene flow across the region, analyses detected structuring of neutral genetic variation across our study area. Gene flow for lynx most strongly related to temperature and elevation compared with other landscape variables (terrain roughness, percent forest cover, and habitat suitability index) and geographic distance alone. Overall, genetic structure in lynx is most consistent with barrier effects created by the Rocky Mountains rather than ecological divergence. Furthermore, warmer temperatures had a measurable impact on gene flow, which suggests connectivity may further decrease in peripheral or fragmented populations as climate warms.

Pleistocene aridification underlies the evolutionary history of the Caribbean endemic, insular, giant Consolea (Opuntioideae)

PREMISE

The Caribbean islands are in the top five biodiversity hotspots on the planet; however, the biogeographic history of the seasonally dry tropical forest (SDTF) there is poorly studied. Consolea consists of nine species of dioecious, hummingbird-pollinated tree cacti endemic to the West Indies, which form a conspicuous element of the SDTF. Several species are threatened by anthropogenic disturbance, disease, sea-level rise, and invasive species and are of conservation concern. However, no comprehensive phylogeny yet exists for the clade.

METHODS

We reconstructed the phylogeny of Consolea, sampling all species using plastomic data to determine relationships, understand the evolution of key morphological characters, and test their biogeographic history. We estimated divergence times to determine the role climate change may have played in shaping the current diversity of the clade.

RESULTS

Consolea appears to have evolved very recently during the latter part of the Pleistocene on Cuba/Hispaniola likely from a South American ancestor and, from there, moved into the Bahamas, Jamaica, Puerto Rico, Florida, and the Lesser Antilles. The tree growth form is a synapomorphy of Consolea and likely aided in the establishment and diversification of the clade.

CONCLUSIONS

Pleistocene aridification associated with glaciation likely played a role in shaping the current diversity of Consolea, and insular gigantism may have been a key innovation leading to the success of these species to invade the often-dense SDTF. This in-situ Caribbean radiation provides a window into the generation of species diversity and the complexity of the SDTF community within the Antilles.

Resolving spatial complexities of hybridization in the context of the gray zone of speciation in North American ratsnakes (Pantherophis obsoletus complex)

Inferring the history of divergence between species in a framework that permits the presence of gene flow has been crucial for characterizing the "gray zone" of speciation, which is the period of time where lineages have diverged but have not yet achieved strict reproductive isolation. However, estimates of both divergence times and rates of gene flow often ignore spatial information, for example when considering the location and width of hybrid zones with respect to changes in the environment between lineages. Using population genomic data from the North American ratsnake complex (Pantherophis obsoletus), we connected phylogeographic estimates of lineage structure, migration, historical demography, and timing of divergence with hybrid zone dynamics. We examined the spatial context of diversification by linking migration and timing of divergence to the location and widths of hybrid zones. Artificial neural network approaches were applied to understand how landscape features and past climate have influenced population genetic structure among these lineages. We found that rates of migration between lineages were associated with the overall width of hybrid zones. Timing of divergence was not related to migration rate or hybrid zone width across species pairs but may be related to the number of alleles weakly introgressing through hybrid zones. This research underscores how incomplete reproductive isolation can be better understood by considering differential allelic introgression and the effects of historical and contemporary landscape features on the formation of lineages as well as overall genomic estimates of migration rates through time.

As old as the hills: Pliocene palaeogeographical processes influence patterns of genetic structure in the widespread, common shrub Banksia sessilis

The impact of Quaternary glaciation on the development of phylogeographic structure in plant species is well documented. In unglaciated landscapes, phylogeographic patterns tend to reflect processes relating to persistence and stochasticity, yet other factors, associated with the palaeogeographical history of the landscape, including geomorphological events, can also have a significant influence. The unglaciated landscape of south-western Western Australia is an ideal location to observe these ancient drivers of lineage diversification, with tectonic activity associated with the Darling Fault in the late Pliocene attributed to patterns of deep phylogeographic divergence in a widespread tree from this region. Interestingly, other species within this region have not shown this pattern and this palaeogeographical boundary therefore presents an opportunity to examine age and historical distribution of plant species endemic to this region. In this study, we assess patterns of genetic diversity and structure across 28 populations of the widespread shrub Banksia sessilis using three cpDNA markers and nine nuclear microsatellite markers. Sixteen cpDNA haplotypes were identified, comprising two major chloroplast DNA lineages that are estimated to have diverged in the Pliocene, approximately 3.3 million years ago. This timing coincides with major geomorphological processes in the landscape, including the separation of the Darling Plateau from the adjacent Swan Coastal Plain, as well as eustatic changes on the Swan Coastal Plain that are likely to have resulted in the physical isolation of historical plant lineages. Chloroplast lineages were broadly aligned with populations associated with older lateritic soils of the Darling Plateau and Geraldton sandplains or the younger sandy soils associated with the Swan Coastal Plain and Southern Coastline. This structural pattern of lateritic versus non-lateritic division was not observed in the nuclear microsatellite data that identified three genetic clades that roughly corresponded to populations in the North, South, and Central portions of the distributions.

Predicting amphibian intraspecific diversity with machine learning: Challenges and prospects for integrating traits, geography, and genetic data

The growing availability of genetic data sets, in combination with machine learning frameworks, offers great potential to answer long-standing questions in ecology and evolution. One such question has intrigued population geneticists, biogeographers, and conservation biologists: What factors determine intraspecific genetic diversity? This question is challenging to answer because many factors may influence genetic variation, including life history traits, historical influences, and geography, and the relative importance of these factors varies across taxonomic and geographic scales. Furthermore, interpreting the influence of numerous, potentially correlated variables is difficult with traditional statistical approaches. To address these challenges, we analysed repurposed data using machine learning and investigated predictors of genetic diversity, focusing on Nearctic amphibians as a case study. We aggregated species traits, range characteristics, and >42,000 genetic sequences for 299 species using open-access scripts and various databases. After identifying important predictors of nucleotide diversity with random forest regression, we conducted follow-up analyses to examine the roles of phylogenetic history, geography, and demographic processes on intraspecific diversity. Although life history traits were not important predictors for this data set, we found significant phylogenetic signal in genetic diversity within amphibians. We also found that salamander species at northern latitudes contained low genetic diversity. Data repurposing and machine learning provide valuable tools for detecting patterns with relevance for conservation, but concerted efforts are needed to compile meaningful data sets with greater utility for understanding global biodiversity.

Individualistic evolutionary responses of Central African rain forest plants to Pleistocene climatic fluctuations

Understanding the evolutionary dynamics of genetic diversity is fundamental for species conservation in the face of climate change, particularly in hyper-diverse biomes. Species in a region may respond similarly to climate change, leading to comparable evolutionary dynamics, or individualistically, resulting in dissimilar patterns. The second-largest expanse of continuous tropical rain forest (TRF) in the world is found in Central Africa. Here, present-day patterns of genetic structure are thought to be dictated by repeated expansion and contraction of TRFs into and out of refugia during Pleistocene climatic fluctuations. This refugia model implies a common response to past climate change. However, given the unrivalled diversity of TRFs, species could respond differently because of distinct environmental requirements or ecological characteristics. To test this, we generated genome-wide sequence data for >700 individuals of seven codistributed plants from Lower Guinea in Central Africa. We inferred species' evolutionary and demographic histories within a comparative phylogeographic framework. Levels of genetic structure varied among species and emerged primarily during the Pleistocene, but divergence events were rarely concordant. Demographic trends ranged from repeated contraction and expansion to continuous growth. Furthermore, patterns in genetic variation were linked to disparate environmental factors, including climate, soil, and habitat stability. Using a strict refugia model to explain past TRF dynamics is too simplistic. Instead, individualistic evolutionary responses to Pleistocene climatic fluctuations have shaped patterns in genetic diversity. Predicting the future dynamics of TRFs under climate change will be challenging, and more emphasis is needed on species ecology to better conserve TRFs worldwide.

Demographic Processes Linked to Genetic Diversity and Positive Selection across a Species' Range

Demography determines the strength of genetic drift, which generally reduces genetic variation and the efficacy of selection. Here, we disentangled the importance of demographic processes at a local scale (census size and mating system) and at a species-range scale (old split between population clusters, recolonization after the last glaciation cycle, and admixture) in determining within-population genomic diversity and genomic signatures of positive selection. Analyses were based on re-sequence data from 52 populations of North American Arabidopsis lyrata collected across its entire distribution. The mating system and range dynamics since the last glaciation cycle explained around 60% of the variation in genomic diversity among populations and 52% of the variation in the signature of positive selection. Diversity was lowest in selfing compared with outcrossing populations and in areas further away from glacial refugia. In parallel, reduced positive selection was found in selfing populations and in populations with a longer route of postglacial range expansion. The signature of positive selection was also reduced in populations without admixture. We conclude that recent range expansion can have a profound influence on diversity in coding and non-coding DNA, similar in magnitude to the shift toward selfing. Distribution limits may in fact be caused by reduced effective population size and compromised positive selection in recently colonized parts of the range.

Landscape and Climate Influence the Patterns of Genetic Diversity and Inbreeding in Cerrado Plant Species

The anthropization of the landscape of the Cerrado biome that has occurred over the past few decades has fragmented its natural environments, impacting the connectivity of the plant populations and altering their gene flow. Plant species may also reduce population size in response to sub-optimal climatic and environmental conditions, and observed distribution patterns may align with theoretical schemes, such as the center–periphery model, that is, it is possible that populations on the edge have lower genetic diversity than center populations, theoretically submitted to environmental conditions closer to the optimum. In this context, we evaluate whether the genetic diversity and inbreeding coefficients of Cerrado plant species are affected by landscape features and climate characteristics, and in particular, if the distribution of the genetic diversity of these plants is consistent with the center–periphery model. To do this, we conducted a literature search for genetic studies of Cerrado plant populations using Scopus, Web of Science, and Scielo databases and the species found were used as a proxy to explore patterns throughout the biome. The data were analyzed using generalized linear mixed models (GLMM) and multiple matrix regressions (MMRRs) to evaluate the effects of landscape features and climatic variables on the observed (HO) and expected heterozygosity (HE), allelic richness (AR) and inbreeding (Fis) patterns of the local populations. The landscape was evaluated in terms of the percentage land cover of agriculture (AG), forestry (FO), remnant vegetation (RV), urban areas (UA), pasture (PA), and water (WA) within buffers of 1 km, 3 km, and 5 km around the study populations. We analyzed 121 populations of 31 plant species. The GLMMs showed that HO was affected by FO regardless of buffer size, while HE was also affected by FO, but also by WA and UA. AR was affected by WA and UA in all three buffer zones while the Fis was affected by FO and AU. The MMRRs showed that WA may affect HO, HE, and Fis within the 1 km buffer, while FO affects HO and UA affects AR within the 5 km buffer. In the case of the 1 km and 3 km buffers, however, the geographic distance between populations was identified as a factor determining the genetic diversity and inbreeding indices, indicating that isolation by distance may be an important factor defining the breeding patterns of the Cerrado plant populations. The GLMMs and MMRRs also showed that the mean annual temperature (MAT) and, to a lesser extent, isothermality (ISO) can explain the variation in genetic diversity observed in the Cerrado plant populations. We also found that the center–periphery model fits the distribution pattern observed in most of the species evaluated, including Annona crassiflora,Annona coriacea, Copaifera langsdorffii, and Eugenia dysenterica. Our results indicate that changes in the climate and the landscape of Brazilian Cerrado must be considered carefully to guarantee minimizing the impacts of these processes on the genetic diversity of Cerrado plant species and ensuring the long-term conservation of these species in this biome.

Westward range expansion from middle latitudes explains the Mississippi River discontinuity in a forest herb of eastern North America

It is often expected that temperate plants have expanded their geographical ranges northward from primarily southern refugia. Evidence for this hypothesis is mixed in eastern North American species, and there is increasing support for colonization from middle latitudes. We studied genome-wide patterns of variation in RADseq loci to test hypotheses concerning range expansion in a North American forest herb (Campanula americana). First, spatial patterns of genetic differentiation were determined. Then phylogenetic relationships and divergence times were estimated. Spatial signatures of genetic drift were also studied to identify the directionality of recent range expansion and its geographical origins. Finally, spatially explicit scenarios for the spread of plants across the landscape were compared, using variation in the population mutation parameter and Tajima's D. We found strong longitudinal subdivision, with populations clustering into groups west and east of the Mississippi River. While the southeastern region was probably part of a diverse Pleistocene refugium, there is little evidence that range expansion involved founders from these southern locales. Instead, declines in genetic diversity and the loss of rare alleles support a westward colonization wave from a middle latitude refugium near the southern Appalachian Mountains, with subsequent expansion from a Pleistocene staging ground in the Mississippi River Valley (0.51-1.27 million years ago). These analyses implicate stepping stone colonization from middle latitudes as an important mechanism of species range expansion in eastern North America. This study further demonstrates the utility of population genetics as a tool to infer the routes travelled by organisms during geographical range expansion.

Do Centres of Endemism provide a spatial context for predicting and preserving plant phylogeographic patterns in the Cape Floristic Region, South Africa?

Aim

The evolutionary forces that gave rise to the exceptional plant species richness of the Cape Floristic Region (CFR) have also likely played a role at the intraspecific level (i.e. plant populations)-and thereby generating shared phylogeographic patterns among taxa. Here we test whether plant populations in the CFR exhibit phylogeographic breaks across the boundaries between Centres of Endemism (CoEs). The boundaries between CoEs (derived from the distribution ranges of endemic taxa and currently mapped at a coarse, Quarter Degree Square scale) represent a spatial proxy for the evolutionary diversifying drivers acting on plant taxa in the CFR.

Location

The CFR, located along the southern Cape of South Africa.

Methods

Published phylogeographic literature were compiled and spatial patterns of genetic divergence re-analysed to assess the frequency at which CFR plant taxa exhibit phylogeographic breaks either (1) across or (2) within CoE boundaries. Population pairs from each study were compared across and within CoEs and scored as either exhibiting a phylogeographic break or not.

Results

Phylogeographic breaks in Cape plants were found to occur across the boundaries of CoEs more often than not. Significantly more population pairs exhibited phylogeographic breaks across CoE boundaries (506 of the 540, χ² = 886, p < 0.001) and fewer breaks within CoEs (94 of 619, χ² = 300, p < 0.001) than would be expected if there was equal probability of a genetic break occurring across CoE boundaries.

Main conclusions

The evolutionary forces that have produced and maintained the exceptional plant diversity in the CFR appear to have operated at the population level, producing similar patterns of phylogeographic structuring of plant lineages regardless of life history or taxonomy. This tendency for Cape plants to exhibit shared patterns of spatially structured genetic diversity that match the distribution of endemic taxa may assist CFR phylogeographers to streamline sampling efforts and test novel hypotheses pertaining to the distribution of genetic diversity among CFR plant taxa. Additionally, the resolution at which CoEs are mapped should be refined, which may provide a valuable tool for future conservation planning and the development of precautionary guidelines for the translocation of genetic material during species reintroductions and commercial cultivation of Cape endemic crops. Thus, to answer the question 'Do Centres of Endemism provide a spatial context for predicting and preserving plant phylogeographic patterns in the Cape Floristic Region, South Africa?'-yes, CoEs do appear to be an important tool for Cape phylogeographers. However, the data is limited and more plant phylogeography work is needed in the CFR.

Genomic analysis of demographic history and ecological niche modeling in the endangered Chinese Grouse Tetrastes sewerzowi

Background

The Quaternary had worldwide consequences in forming the contemporary diversity of many populations, species and communities, which is characterized by marked climatic oscillations between glacial and interglacial periods. The origin and evolution of biodiversity in mountainous areas are highly dependent on historical orogenesis and associated climatic changes. The Chinese grouse Tetrastes sewerzowi is a forest-dwelling species endemic to the mountains to the east of the Qinghai–Tibet Plateau, which has been listed as Near Threatened with a decreasing trend by the IUCN because of ongoing deforestation and fragmentation of coniferous forests. It is important to place current population status into a broader ecological and evolutionary context to understand their demographic history.

Results

Analyses of the Chinese Grouse genome revealed fluctuations throughout the Pleistocene in effective population size. Populations decreased during early to middle Pleistocene but showed an expansion during late Pleistocene which was then followed by a sharp decline during the last glacial maximum (LGM). Ecological niche modeling indicated that a suitable habitat shift between high altitude regions to low altitude regions was due to a changing climate. This result parallels patterns of population size change in Chinese Grouse estimated from PSMC modelling, which suggested an expansion in population size from the last interglacial period (LIG) and then a peak and a bottleneck occurring at the last glacial maximum (LGM). Furthermore, the present-day distribution of Chinese Grouse is greatly reduced and fragmented. It will likely become even more fragmented in the future since coniferous forest cover is threatened in the region of their distribution and the availability of such habitat restricts their ecological niche.

Conclusions

The Chinese Grouse have experienced substantial population size changes from the beginning to the LIG and reached a peak before the LGM. A sharp decrease and bottleneck occurred during the LGM, when the coniferous forests were subjected to extensive loss. The results inferred from the whole genome sequencing and species distribution models both support historical population fluctuations. The distribution of the Chinese Grouse is strongly dependent on the coniferous forest cover. To protect the fragmented coniferous forests is an essential action to protect the Chinese Grouse.

Population Genetic Assessment of Anadromous and Resident Striped Bass (Morone saxatilis) in the Roanoke River System, Eastern United States

Striped bass is the subject of important commercial and sport fisheries in North America. The Roanoke River drainage—especially Smith Mountain Lake, Leesville Lake, and Kerr Reservoir—has popular recreational striped bass fisheries. After construction of five hydroelectric dams, populations became landlocked, declined, and have been supplemented by stocking. A key basis for responsibly augmenting populations is to characterize genetic variation and incorporate the findings into responsible hatchery and stocking practices. Genetic variation at 12 microsatellite DNA loci was evaluated among 837 striped bass representing 16 collections across the native range; populations from rivers in South Carolina, North Carolina, Chesapeake Bay, and Hudson River were screened to provide context for assessing genetic structure within the Roanoke system. Analysis of population genetic differentiation showed landlocked Roanoke River striped bass to be distinctive. Subject to genetic isolation, high M ratios, and relatively low Ne estimates suggest loss of genetic variation, and relatedness analysis showed heightened frequencies of related individuals. These insights into population genetics, demographics, and existing guidelines for broodstock acquisition and mating designs can inform genetically cognizant hatchery management and stocking for striped bass in the Roanoke River drainage. In particular, we recommend the use of larger numbers of breeders and factorial mating designs to increase the genetic diversity of propagated striped bass stocked within the Roanoke River drainage.

Spatiotemporal analyses suggest the role of glacial history and the ice-free corridor in shaping American badger population genetic variation

Recurring glacial cycles through the Quaternary period drastically altered the size and distribution of natural populations of North American flora and fauna. The "southerly refugia model" has been the longstanding framework for testing the effects of glaciation on contemporary genetic patterns; however, insights from ancient DNA have contributed to the reconstruction of more complex histories for some species. The American badger, Taxidea taxus, provides an interesting species for exploring the genetic legacy of glacial history, having been hypothesized to have postglacially emerged from a single, southerly refugium to recolonize northern latitudes. However, previous studies have lacked genetic sampling from areas where distinct glacial refugia have been hypothesized, including the Pacific Northwest and American Far North (Yukon, Alaska). In order to further investigate the phylogeographic history of American badgers, we collected mitochondrial DNA sequence data from ancient subfossil material collected within the historical range (Alaska, Yukon) and combined them with new and previously published data from across the species' contemporary distribution (n = 1,207). We reconstructed a mostly unresolved phylogenetic tree and star-like haplotype network indicative of emergence from a largely panmictic glacial refugium and recent population expansion, the latter further punctuated by significantly negative Tajima's D and Fu's Fs values. Although directionality of migration cannot be unequivocally inferred, the moderate to high levels of genetic variation exhibited by American badgers, alongside the low frequency of haplotypes with indels in the Midwest, suggest a potential recolonization into central North America after the hypothesized ice-free corridor reopened ~13,000 years ago. Overall, the expanded reconstruction of phylogeographic history of American badgers offers a broader understanding of contemporary range-wide patterns and identifies unique genetic units that can likely be used to inform conservation of at-risk populations at the northern periphery.

Extensive Genetic Connectivity and Historical Persistence Are Features of Two Widespread Tree Species in the Ancient Pilbara Region of Western Australia

Phylogeographic studies can be used as a tool to understand the evolutionary history of a landscape, including the major drivers of species distributions and diversity. Extensive research has been conducted on phylogeographic patterns of species found in northern hemisphere landscapes that were affected by glaciations, yet the body of literature for older, unaffected landscapes is still underrepresented. The Pilbara region of north-western Australia is an ancient and vast landscape that is topographically complex, consisting of plateaus, gorges, valleys, and ranges, and experiences extreme meteorological phenomena including seasonal cyclonic activity. These features are expected to influence patterns of genetic structuring throughout the landscape either by promoting or restricting the movement of pollen and seed. Whilst a growing body of literature exists for the fauna endemic to this region, less is known about the forces shaping the evolution of plant taxa. In this study we investigate the phylogeography of two iconic Pilbara tree species, the Hamersley Bloodwood (Corymbia hamersleyana) and Western Gidgee (Acacia pruinocarpa), by assessing patterns of variation and structure in several chloroplast DNA regions and nuclear microsatellite loci developed for each species. Gene flow was found to be extensive in both taxa and there was evidence of long-distance seed dispersal across the region (pollen to seed ratios of 6.67 and 2.96 for C. hamersleyana and A. pruinocarpa, respectively), which may result from flooding and strong wind gusts associated with extreme cyclonic activity. Both species possessed high levels of cpDNA genetic diversity in comparison to those from formerly glaciated landscapes (C. hamersleyana = 14 haplotypes, A. pruinocarpa = 37 haplotypes) and showed evidence of deep lineage diversification occurring from the late Miocene, a time of intensifying aridity in this landscape that appears to be a critical driver of evolution in Pilbara taxa. In contrast to another study, we did not find evidence for topographic features acting as refugia for the widely sampled C. hamersleyana.

Genotyping-by-sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest-ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome-wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal-Cascades (cluster 1), east-slope Cascades-Sierra Nevadas-Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the "stable-edge hypothesis" that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited "trailing-edge" dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the "inland dispersal hypothesis" predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific-coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable-edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific-coastal genetic lineage of quaking aspen.

Phylogeography of the Cran’s bully Gobiomorphus basalis (Gobiiformes: Eleotridae) and an analysis of species boundaries within the New Zealand radiation of Gobiomorphus

New Zealand has a complex recent history of climatic and tectonic change that has left variable signatures in the geographic distribution and genetic structure of the region’s flora and fauna. To identify concordant patterns, a broad range of taxa must be examined and compared. In New Zealand’s North Island, a consensus is forming as to the dominant biogeographic barriers in the region although obligate freshwater taxa have not been considered in this framework. We use single-nucleotide polymorphisms to investigate phylogeography in the widespread obligate freshwater fish Gobiomorphus basalis on the North Island. Phylogeographic patterns within G. basalis reveal biogeographic disjunctions that are in some ways consistent and in other ways at odds with established patterns, providing insight into the processes that have shaped the islands’ biogeography. We also use phylogeography to delineate species boundaries within the entire New Zealand radiation of Gobiomorphus and find that it contains several morphologically cryptic species. We resolve two clades within G. basalis that correspond to areas north and south of the Taupo Volcanic Zone. We confirm the distinctiveness of Gobiomorphus alpinus relative to Gobiomorphus cotidianus, as well as the presence of two lineages within Gobiomorphus breviceps that were previously identified based on mitochondrial data.

JC, McFarland K, Goffinet B. Population Genomics and Phylogeography of a Clonal Bryophyte With Spatially Separated Sexes and Extreme Sex Ratios

The southern Appalachian (SA) is one of the most biodiversity-rich areas in North America and has been considered a refugium for many disjunct plant species, from the last glacial period to the present. Our study focuses on the SA clonal hornwort, Nothoceros aenigmaticus J. C. Villarreal & K. D. McFarland. This hornwort was described from North Carolina and is widespread in the SA, growing on rocks near or submerged in streams in six and one watersheds of the Tennessee (TR) and Alabama (AR) Rivers, respectively. Males and female populations occur in different watersheds, except in the Little Tennessee (TN) River where an isolated male population exists ca. 48 km upstream from the female populations. The sex ratio of 1:0 seems extreme in each population. In this study, we use nuclear and organellar microsatellites from 250 individuals from six watersheds (seven populations) in the SA region and two populations from Mexico (23 individuals). We, then, selected 86 individuals from seven populations and used genotyping by sequencing to sample over 600 bi-allelic markers. Our results suggest that the SA N. aenigmaticus and Mexican plants are a nested within a clade of sexual tropical populations. In the US populations, we confirm an extreme sex ratio and only contiguous US watersheds share genotypes. The phylogenetic analysis of SNP data resolves four clusters: Mexican populations, male plants (Little Pigeon and Pigeon river watersheds) and two clusters of female plants; one from the Little Tennessee and Hiwassee Rivers (TR) and the other from the Ocoee (TR) and Coosa (AR) Rivers. All clusters are highly differentiated (Fst values over 0.9). In addition, our individual assignment analyses and PCAs reflect the phylogenetic results grouping the SA samples in three clades and recovering males and female plants with high genetic differentiation (Fst values between 0.5 and 0.9 using microsatellites and bi-allelic markers). Our results point to Pleistocene events shaping the biogeographical pattern seen in US populations. The extreme sex ratio reflects isolation and highlights the high vulnerability of the populations in the SA.

Different processes shape the patterns of divergence in the nuclear and chloroplast genomes of a relict tree species in East Asia

Isolation by spatial distance (IBD), environment (IBE), and historical climatic instability (IBI) are three common processes assessed in phylogeographic and/or landscape genetic studies. However, the relative contributions of these three processes with respect to spatial genetic patterns have seldom been compared. Moreover, whether the relative contribution differs in different regions or when assessed using different genetic markers has rarely been reported. Lindera obtusiloba has been found to have two sister genetic clades of chloroplast (cpDNA) and nuclear microsatellite (nSSR), both of which show discontinuous distribution in northern and southern East Asia. In this study, we used the Mantel test and multiple matrix regression with randomization (MMRR) to determine the relative contributions of IBD, IBE, and IBI with respect to L. obtusiloba populations. Independent Mantel tests and MMRR calculations were conducted for two genetic data sets (cpDNA and nSSR) and for different regions (the overall species range, and northern and southern subregions of the range). We found a significant IBI pattern in nSSR divergence for all assessed regions, whereas no clear IBI pattern was detected with respect to cpDNA. In contrast, significant (or marginal) divergent IBD patterns were detected for cpDNA in all regions, whereas although a significant IBE was apparent with respect to the overall range, the effect was not detected in the two subregions. The differences identified in nSSR and cpDNA population divergence may be related to differences in the heredity and ploidy of the markers. Compared with the southern region, the northern region showed less significant correlation patterns, which may be related to the shorter population history and restricted population range. The findings of this study serve to illustrate that comparing between markers or regions can contribute to gaining a better understanding the population histories of different genomes or within different regions of a species' range.

Oaks: an evolutionary success story

The genus Quercus is among the most widespread and species-rich tree genera in the northern hemisphere. The extraordinary species diversity in America and Asia together with the continuous continental distribution of a limited number of European species raise questions about how macro- and microevolutionary processes made the genus Quercus an evolutionary success. Synthesizing conclusions reached during the past three decades by complementary approaches in phylogenetics, phylogeography, genomics, ecology, paleobotany, population biology and quantitative genetics, this review aims to illuminate evolutionary processes leading to the radiation and expansion of oaks. From opposing scales of time and geography, we converge on four overarching explanations of evolutionary success in oaks: accumulation of large reservoirs of diversity within populations and species; ability for rapid migration contributing to ecological priority effects on lineage diversification; high rates of evolutionary divergence within clades combined with convergent solutions to ecological problems across clades; and propensity for hybridization, contributing to adaptive introgression and facilitating migration. Finally, we explore potential future research avenues, emphasizing the integration of microevolutionary and macroevolutionary perspectives.