High genetic diversity and stable Pleistocene distributional ranges in the widespread Mexican red oak Quercus castanea Née (1801) (Fagaceae)

The Mexican highlands are areas of high biological complexity where taxa of Nearctic and Neotropical origin and different population histories are found. To gain a more detailed view of the evolution of the biota in these regions, it is necessary to evaluate the effects of historical tectonic and climate events on species. Here, we analyzed the phylogeographic structure, historical demographic processes, and the contemporary period, Last Glacial Maximum (LGM) and Last Interglacial (LIG) ecological niche models of Quercus castanea, to infer the historical population dynamics of this oak distributed in the Mexican highlands. A total of 36 populations of Q. castanea were genotyped with seven chloroplast microsatellite loci in four recognized biogeographic provinces of Mexico: the Sierra Madre Occidental (western mountain range), the Central Plateau, the Trans‐Mexican Volcanic Belt (TMVB, mountain range crossing central Mexico from west to east) and the Sierra Madre del Sur (SMS, southern mountain range). We obtained standard statistics of genetic diversity and structure and tested for signals of historical demographic expansions. A total of 90 haplotypes were identified, and 29 of these haplotypes were restricted to single populations. The within‐population genetic diversity was high (mean h_S = 0.72), and among‐population genetic differentiation showed a strong phylogeographic structure (N_ST = 0.630 > G_ST = 0.266; p < .001). Signals of demographic expansion were identified in the TMVB and the SMS. The ecological niche models suggested a considerable percentage of stable distribution area for the species during the LGM and connectivity between the TMVB and the SMS. High genetic diversity, strong phylogeographic structure, and ecological niche models suggest in situ permanence of Q. castanea populations with large effective population sizes. The complex geological and climatic histories of the TMVB help to explain the origin and maintenance of a large proportion of the genetic diversity in this oak species.     

Species distribution modelling and cpSSR reveal population history of the Neotropical annual herb Tithonia rotundifolia (Asteraceae)

• The impacts of the historical geologic and climatic events on the diversity and genetic structure of Neotropical taxa have recently become a subject of study. However, annual plants associated with tropical dry forests remain under‐studied. The exploration of additional taxa in contrasting environments will improve the current understanding of responses of the Neotropical biota to these events. Here, we explore the species distribution and geographic structure of the annual herb Tithonia rotundifolia.
• We sampled 175 individuals from 19 populations of T. rotundifolia. Species distribution modelling and six microsatellite chloroplast loci were used to infer its population history. We identified areas of historical climate suitability and then tested if there is genetic structuring among these areas.
• Haplotypes showed strong phylogeographic structure. Historical climatic suitability areas were found along the Pacific coast; however, a gap was found at the Isthmus of Tehuantepec (IT). Although Bayesian analysis showed population structuring, amova revealed that the IT is not its main driver. Instead, a subdivision into a higher number of regions had higher F_CT values. Also, populations to the east of the IT showed evidence of recent population expansion and migration in a south–north direction.
• Pleistocene climate fluctuations partially explain the geographic structure of T. rotundifolia. However, life‐history characteristics such as limited seed dispersal and the patchy distribution of suitable habitats explain the high haplotype diversity and population sub‐structuring and diversity. Lastly, the absence of geographic structure of some haplotypes may indicate long‐distance dispersal, or hybridisation with the closely related T. tubaeformis.

     

The Quaternary evolutionary history,potential distribution dynamics,and conservation implications for a Qinghai–Tibet Plateau endemic herbaceous perennial,Anisodus tanguticus (Solanaceae)

Various hypotheses have been proposed about the Quaternary evolutionary history of plant species on the Qinghai–Tibet Plateau (QTP), yet only a handful of studies have considered both population genetics and ecological niche context. In this study, we proposed and compared climate refugia hypotheses based on the phylogeographic pattern of Anisodus tanguticus (three plastid DNA fragments and nuclear internal transcribed spacer regions from 32 populations) and present and past species distribution models (SDMs). We detected six plastid haplotypes in two well‐differentiated lineages. Although all haplotypes could be found in its western (sampling) area, only haplotypes from one lineage occurred in its eastern area. Meanwhile, most genetic variations existed between populations (F_ST = 0.822). The SDMs during the last glacial maximum and last interglacial periods showed range fragmentation in the western area and significant range contraction in the eastern area, respectively, in comparison with current potential distribution. This species may have undergone intraspecific divergence during the early Quaternary, which may have been caused by survival in different refugia during the earliest known glacial in the QTP, rather than geological isolation due to orogenesis events. Subsequently, climate oscillations during the Quaternary resulted in a dynamic distribution range for this species as well as the distribution pattern of its plastid haplotypes and nuclear genotypes. The interglacial periods may have had a greater effect on A. tanguticus than the glacial periods. Most importantly, neither genetic data nor SDM alone can fully reveal the climate refugia history of this species. We also discuss the conservation implications for this important Tibetan folk medicine plant in light of these findings and SDMs under future climate models. Together, our results underline the necessity to combine phylogeographic and SDM approaches in future investigations of the Quaternary evolutionary history of species in topographically complex areas, such as the QTP.     

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid‐Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least‐cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (Φ_ST(n) = 0.463 and Φ_ST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations’ connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long‐term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.     

Genetic structure and demographic history of Lymantria dispar (Linnaeus, 1758) (Lepidoptera: Erebidae) in its area of origin and adjacent areas

We analyzed the population genetic structure and demographic history of 20 Lymantria dispar populations from Far East Asia using microsatellite loci and mitochondrial genes. In the microsatellite analysis, the genetic distances based on pairwise F_ST values ranged from 0.0087 to 0.1171. A NeighborNet network based on pairwise F_ST genetic distances showed that the 20 regional populations were divided into five groups. Bayesian clustering analysis (K = 3) demonstrated the same groupings. The populations in the Korean Peninsula and adjacent regions, in particular, showed a mixed genetic pattern. In the mitochondrial genetic analysis based on 98 haplotypes, the median‐joining network exhibited a star shape that was focused on three high‐frequency haplotypes (Haplotype 1: central Korea and adjacent regions, Group 1; Haplotype 37: southern Korea, Group 2; and Haplotype 90: Hokkaido area, Group 3) connected by low‐frequency haplotypes. The mismatch distribution dividing the three groups was unimodal. In the neutral test, Tajima's D and Fu's FS tests were negative. We can thus infer that the Far East Asian populations of L. dispar underwent a sudden population expansion. Based on the age expansion parameter, the expansion time was inferred to be approximately 53,652 years before present (ybp) for Group 1, approximately 65,043 ybp for Group 2, and approximately 76,086 ybp for Group 3. We propose that the mixed genetic pattern of the inland populations of Far East Asia is due to these expansions and that the inland populations of the region should be treated as valid subspecies that are distinguishable from other subspecies by genetic traits.     

Genetic differentiation among European samples of the arctic‐alpine leaf beetle,Chrysomela lapponica

Genetic differentiation within and among isolated populations of the arctic‐alpine leaf beetle, Chrysomela lapponica L. (Coleoptera: Chrysomelidae), specialized on either Salix L. (Salicaceae) or Betula L. (Betulaceae) species, was assessed by F‐statistics analysis at seven allozyme loci. Beetles were collected on Salix spp. in lowland Finland (four samples), at mid elevation in the Black Forest in Germany (450 m) and the Massif Central in France (two samples, 930–1 300 m), and at high elevation in the French Alps (2 300 m). Beetles sampled in the Czech Republic (650 m) fed on Betula pubescens Ehrh. Larvae feeding on Salix spp. secreted host‐derived salicylaldehyde as major toxin; those feeding on B. pubescens secreted isobutyrates and 2‐methylbutyrates of mixed plant–insect origin. In all samples, a heterozygote deficit was observed (0.120<F_IS<0.568), in particular in populations collected at mid or high elevation (F_IS>0.4). The estimated mean level of genetic differentiation among all populations was high (F_ST = 0.276). Differentiation was highly variable between pairs of populations (F_ST = 0.093–0.455, all significant) without any correlation with geographic distance. The sample collected in the Czech Republic, from the only population on B. pubescens, was not the most divergent.     

Phylogeography of Chinese cherry (Prunus pseudocerasus Lindl.) inferred from chloroplast and nuclear DNA: insights into evolutionary patterns and demographic history

Chinese cherry (Prunus pseudocerasus Lindl.) is a commercially valuable fruit crop in China. In order to obtain new insights into its evolutionary history and provide valuable recommendations for resource conservation, phylogeographic patterns of 26 natural populations (305 total individuals) from six geographic regions were analyzed using chloroplast and nuclear DNA fragments. Low levels of haplotype and nucleotide diversity were found in these populations, especially in landrace populations. It is likely that a combined effect of botanical characteristics impact the effective population size, such as inbreeding mating system, long life span, as well as vegetative reproduction. In addition, strong bottleneck effect caused by domestication, together with founder effect after dispersal and subsequent demographic expansion, might also accelerate the reduction of the genetic variation in landrace populations. Interestingly, populations from Longmen Mountain (LMM) and Daliangshan Mountain (DLSM) exhibited relatively higher levels of genetic diversity, inferring the two historical genetic diversity centers of the species. Moreover, moderate population subdivision was also detected by both chloroplast DNA (G_ST = 0.215; N_ST = 0.256) and nuclear DNA (G_ST = 0.146; N_ST = 0.342), respectively. We inferred that the episodes of efficient gene flow through seed dispersal, together with features of long generation cycle and inbreeding mating system, were likely the main contributors causing the observed phylogeographic patterns. Finally, factors that led to the present demographic patterns of populations from these regions and taxonomic varieties were also discussed.     

Phylogeography of Podocarpus matudae (Podocarpaceae): pre‐Quaternary relicts in northern Mesoamerican cloud forests

Aim Cloud forests of northern Mesoamerica represent the northern and southern limit of the contact zone between species otherwise characteristic of North or South America. Several phylogeographic studies featuring temperate conifer species have improved our understanding of species responses to environmental changes. In contrast, conifer species that presumably colonized northern Mesoamerica from South America are far less studied. A phylogeographic study of Podocarpus matudae (Podocarpaceae) was conducted to identify any major evolutionary divergences or disjunctions across its range and to determine if its current distribution is associated with pre‐Quaternary climatic and/or long‐distance dispersal events. Location Northern Mesoamerica (Mexico and Guatemala). Methods Sixteen populations (157 individuals) of P. matudae were screened for variation at two plastid DNA markers. The intra‐specific phylogenetic relationships among haplotypes were reconstructed using Bayesian inference. Population genetic analyses were undertaken to gain insight into the evolutionary history of these populations. To test whether genetic divergence among populations occurred at different time‐scales plastid DNA sequence data and fossil‐ and coalescent‐based calibrations were integrated. Results The combination of plastid markers yielded 11 haplotypes. Differentiation among populations based on DNA variation (G_ST) (0.707, SE 0.0807) indicated a clear population structure in P. matudae. Differentiation for ordered alleles (N_ST) (0.811, SE 0.0732) was higher than that for G_ST, indicating phylogeographical structure in P. matudae. Most of the total variation (81.3%, P < 0.0001) was explained by differences among populations. The estimated divergence time between the unique haplotypes from a Guatemalan population and the two most common haplotypes from the Sierra Madre Oriental in Mexico was between 10 and 20 Ma, and further haplotype divergence in the poorly resolved clade of the Sierra Madre Oriental occurred between 3 and 0.5 Ma. Main conclusions Divergence estimations support the hypothesis that extant Podocarpus matudae populations are pre‐Quaternary relicts. This finding is consistent with fossil and pollen data that support a Miocene age for temperate floristic elements in Mesoamerican cloud forests, whereas further haplotype divergence within the Sierra Madre Oriental, Chiapas and Guatemala occurred more recently, coinciding with Pleistocene cloud forest refugia.     

Elevation as a barrier: genetic structure for an Atlantic rain forest tree (Bathysa australis) in the Serra do Mar mountain range,SE Brazil

Distance and discrete geographic barriers play a role in isolating populations, as seed and pollen dispersal become limited. Nearby populations without any geographic barrier between them may also suffer from ecological isolation driven by habitat heterogeneity, which may promote divergence by local adaptation and drift. Likewise, elevation gradients may influence the genetic structure and diversity of populations, particularly those marginally distributed. Bathysa australis (Rubiaceae) is a widespread tree along the elevation gradient of the Serra do Mar, SE Brazil. This self‐compatible species is pollinated by bees and wasps and has autochoric seeds, suggesting restricted gene dispersal. We investigated the distribution of genetic diversity in six B. australis populations at two extreme sites along an elevation gradient: a lowland site (80–216 m) and an upland site (1010–1100 m.a.s.l.). Nine microsatellite loci were used to test for genetic structure and to verify differences in genetic diversity between sites. We found a marked genetic structure on a scale as small as 6 km (F_ST = 0.21), and two distinct clusters were identified, each corresponding to a site. Although B. australis is continuously distributed along the elevation gradient, we have not observed a gene flow between the extreme populations. This might be related to B. australis biological features and creates a potential scenario for adaptation to the different conditions imposed by the elevation gradient. We failed to find an isolation‐by‐distance pattern; although on the fine scale, all populations showed spatial autocorrelation until ~10‐20 m. Elevation difference was a relevant factor though, but we need further sampling effort to check its correlation with genetic distance. The lowland populations had a higher allelic richness and showed higher rare allele counts than the upland ones. The upland site may be more selective, eliminating rare alleles, as we did not find any evidence for bottleneck.     

Detecting a hierarchical genetic population structure: the case study of the Fire Salamander (Salamandra salamandra) in Northern Italy

The multistep method here applied in studying the genetic structure of a low dispersal and philopatric species, such as the Fire Salamander Salamandra salamandra, was proved to be effective in identifying the hierarchical structure of populations living in broad‐leaved forest ecosystems in Northern Italy. In this study, 477 salamander larvae, collected in 28 sampling populations (SPs) in the Prealpine and in the foothill areas of Northern Italy, were genotyped at 16 specie‐specific microsatellites. SPs showed a significant overall genetic variation (Global F_ST = 0.032, P < 0.001). The genetic population structure was assessed by using STRUCTURE 2.3.4. We found two main genetic groups, one represented by SPs inhabiting the Prealpine belt, which maintain connections with those of the Eastern foothill lowland (PEF), and a second group with the SPs of the Western foothill lowland (WF). The two groups were significantly distinct with a Global F_ST of 0.010 (P < 0.001). While the first group showed a moderate structure, with only one divergent SP (Global F_ST = 0.006, P < 0.001), the second group proved more structured being divided in four clusters (Global F_ST = 0.017, P = 0.058). This genetic population structure should be due to the large conurbations and main roads that separate the WF group from the Prealpine belt and the Eastern foothill lowland. The adopted methods allowed the analysis of the genetic population structure of Fire Salamander from wide to local scale, identifying different degrees of genetic divergence of their populations derived from forest fragmentation induced by urban and infrastructure sprawl.     

Bayesian analysis of molecular variance in pyrosequences quantifies population genetic structure across the genome of Lycaeides butterflies

The distribution of genetic variation within and among populations is commonly used to infer their demographic and evolutionary histories. This endeavour has the potential to benefit substantially from high‐throughput next‐generation sequencing technologies through a rapid increase in the amount of data available and a corresponding increase in the precision of parameter estimation. Here we report the results of a phylogeographic study of the North American butterfly genus Lycaeides using 454 sequence data. This study serves the dual purpose of demonstrating novel molecular and analytical methods for population genetic analyses with 454 sequence data and expanding our knowledge of the phylogeographic history of Lycaeides. We obtained 341 045 sequence reads from 12 populations that we were able to assemble into 15 262 contigs (most of which were variable), representing one of the largest population genetic data sets for a non‐model organism to date. We examined patterns of genetic variation using a hierarchical Bayesian analysis of molecular variance model, which provides precise estimates of genome‐level φ_ST while appropriately modelling uncertainty in locus‐specific φ_ST. We found that approximately 36% of sequence variation was partitioned among populations, suggesting historical or current isolation among the sampled populations. Estimates of pairwise genome‐level φ_ST were largely consistent with a previous phylogeographic model for Lycaeides, suggesting fragmentation into two to three refugia during Pleistocene glacial cycles followed by post‐Pleistocene range expansion and secondary contact leading to introgressive hybridization. This study demonstrates the potential of using genome‐level data to better understand the phylogeographic history of populations.     

Phylogeography of two closely related species of Allium endemic to East Asia: Population evolution in response to climate oscillations

This study investigated the effects of climate oscillations on the evolution of two closely related Allium species, A. neriniflorum and A. tubiflorum. We sequenced three cp DNA (cpDNA) fragments (rps16, rpl32‐trnL, and trnD‐trnT, together approximately 2,500 bp in length) of two closely related Allium species, with samples from 367 individuals in 47 populations distributed across the total range of these species. The interspecific and intraspecific divergence times of the two species were in the Quaternary glaciation. The population divergence was high for the cpDNA variation, suggesting a significant phylogeographic structure (N_ST = 0.844, G_ST = 0.798, p < 0.05). Remarkable ecological differentiation was also revealed by Niche models and statistical analyses. Our results suggest the speciation event of the two species was triggered by violent climatic changes during the Quaternary glaciation.     

Influence of drainage divides versus arid corridors on genetic structure and demography of a widespread freshwater turtle,Emydura macquarii krefftii,from Australia

The influence of Pleistocene climatic cycles on Southern Hemisphere biotas is not yet well understood. Australia's eastern coastal margin provides an ideal setting for examining the relative influence of landscape development, sea level fluctuation, and cyclic climatic aridity on the evolution of freshwater biodiversity. We examined the impact of climatic oscillations and physical biogeographic barriers on the evolutionary history of the wide‐ranging Krefft's river turtle (Emydura macquarii krefftii), using range‐wide sampling (649 individuals representing 18 locations across 11 drainages) and analysis of mitochondrial sequences (~1.3‐kb control region and ND4) and nuclear microsatellites (12 polymorphic loci). A range of phylogeographic (haplotype networks, molecular dating), demographic (neutrality tests, mismatch distributions), and population genetic analyses (pairwise F_ST, analysis of molecular variance, Bayesian clustering analysis) were implemented to differentiate between competing demographic (local persistence vs. range expansion) and biogeographic (arid corridor vs. drainage divide) scenarios. Genetic data reveal population genetic structure in Krefft's river turtles primarily reflects isolation across drainage divides. Striking north‐south regional divergence (2.2% ND4 p‐distance; c. 4.73 Ma, 95% higher posterior density (HPD) 2.08–8.16 Ma) was consistent with long‐term isolation across a major drainage divide, not an adjacent arid corridor. Ancient divergence among regional lineages implies persistence of northern Krefft's populations despite the recurrent phases of severe local aridity, but with very low contemporary genetic diversity. Stable demography and high levels of genetic diversity are inferred for southern populations, where aridity was less extreme. Range‐wide genetic structure in Krefft's river turtles reflects contemporary and historical drainage architecture, although regional differences in the extent of Plio–Pleistocene climatic aridity may be reflected in current levels of genetic diversity.     

Genetic evidence for the persistence of the critically endangered Sierra Nevada red fox in California

California is home to both the native state-threatened Sierra Nevada red fox (Vulpes vulpes necator), which historically inhabited high elevations of the Sierra Nevada and Cascade mountains, and to multiple low-elevation red fox populations thought to be of exotic origin. During the past few decades the lowland populations have dramatically expanded their distribution, and possibly moved into the historic range of the native high-elevation fox. To determine whether the native red fox persists in its historic range in California, we compared mitochondrial cytochrome-b haplotypes of the only currently-known high-elevation population (n = 9 individuals) to samples from 3 modern lowland populations (n = 35) and historic (1911–1941) high-elevation (n = 22) and lowland (n = 7) populations. We found no significant population differentiation among the modern and historic high-elevation populations (average pairwise F _ST = 0.06), but these populations differed substantially from all modern and historic lowland populations (average pairwise F _ST = 0.52). Among lowland populations, the historic and modern Sacramento Valley populations were not significantly differentiated from one another (F _ST = −0.06), but differed significantly from recently founded populations in the San Francisco Bay region and in southern California (average pairwise F _ST = 0.42). Analysis of molecular variance indicated that 3 population groupings (mountain, Sacramento Valley, and other lowland regions) explained 45% of molecular variance (F _CT = 0.45) whereas only 4.5% of the variance was partitioned among populations within these groupings (F _SC = 0.08). These findings provide strong evidence that the native Sierra Nevada red fox has persisted in northern California. However, all nine samples from this population had the same haplotype, suggesting that several historic haplotypes may have become lost. Unidentified barriers have apparently prevented gene flow from the Sacramento Valley population to other eastern or southern populations in California. Future studies involving nuclear markers are needed to assess the origin of the Sierra Nevada red fox and to quantify levels of nuclear gene flow.     

Bioclimatic regions influence genetic structure of four Jordanian Stipa species

Strong environmental gradients can affect the genetic structure of plant populations, but little is known as to whether closely related species respond similarly or idiosyncratically to ecogeographic variation. We analysed the extent to which gradients in temperature and rainfall shape the genetic structure of four Stipa species in four bioclimatic regions in Jordan. Genetic diversity, differentiation and structure of Stipa species were investigated using amplified fragment length polymorphism (AFLP) molecular markers. For each of the four study species, we sampled 120 individuals from ten populations situated in distinct bioclimatic regions and assessed the degree of genetic diversity and genetic differentiation within and among populations. The widespread ruderals Stipa capensis and S. parviflora had higher genetic diversity than the geographically restricted semi‐desert species S. arabica and S. lagascae. In three of the four species, genetic diversity strongly decreased with precipitation, while genetic diversity increased with temperature in S. capensis. Most genetic diversity resided among populations in the semi‐desert species (Φ_ST = 0.572/0.595 in S. arabica/lagascae) but within populations in the ruderal species (Φ_ST = 0.355/0.387 S. capensis/parviflora). Principal coordinate analysis ( PCoA) and STRUCTURE analysis showed that Stipa populations of all species clustered ecogeographically. A genome scan revealed that divergent selection at particular AFLP loci contributed to genetic differentiation. Irrespective of their different life histories, Stipa species responded similarly to the bioclimatic gradient in Jordan. We conclude that, in addition to predominant random processes, steep climatic gradients might shape the genetic structure of plant populations.     

Molecular phylogeography of Fagus engleriana (Fagaceae) in subtropical China: limited admixture among multiple refugia

Beech is one of the most important trees in the temperate and subtropical forests of the Northern Hemisphere. Despite Chinese beeches have the particularity that only grow in subtropical areas, they have received few phylogeographic research. In this study, we sampled 25 populations of the northernmost-distributed Chinese beech, Fagus engleriana, and detected six haplotypes across 350 individuals by using sequences of two chloroplast intergenic spacers. The chloroplast genetic diversity was relatively low (h _T?=?0.659), with most genetic variance residing among populations (G _ST?=?0.831, N _ST?=?0.855, G _ST??N _ST). SAMOVA analysis indicated that populations clustered into six groups with little admixture among them (most groups were characterized by a unique hapotype). Pairwise difference among haplotypes and Fu??s Fs statistic indicated that populations of F. engleriana have not experienced recent sudden expansions. Both the phylogeographic and demographic patterns found in this study suggest that F. engleriana remained fragmented in multiple refugia throughout the Pleistocene climatic changes, and experienced limited both glacial and interglacial/postglacial expansion. The results of this study imply that long-term isolation among multiple refugia, coupled with little admixture among populations of different refugia provided numerous opportunities for population divergence and allopatric speciation, which might be a driving factor for the exceptionally broad temperate species diversity in southern China.     

Congruent population genetic structures and divergence histories in anther‐smut fungi and their host plants Silene italica and the Silene nutans species complex

Study of the congruence of population genetic structure between hosts and pathogens gives important insights into their shared phylogeographical and coevolutionary histories. We studied the population genetic structure of castrating anther‐smut fungi (genus Microbotryum) and of their host plants, the Silene nutans species complex, and the morphologically and genetically closely related Silene italica, which can be found in sympatry. Phylogeographical population genetic structure related to persistence in separate glacial refugia has been recently revealed in the S. nutans plant species complex across Western Europe, identifying several distinct lineages. We genotyped 171 associated plant–pathogen pairs of anther‐smut fungi and their host plant individuals using microsatellite markers and plant chloroplastic single nucleotide polymorphisms. We found clear differentiation between fungal populations parasitizing S. nutans and S. italica plants. The population genetic structure of fungal strains parasitizing the S. nutans plant species complex mirrored the host plant genetic structure, suggesting that the pathogen was isolated in glacial refugia together with its host and/or that it has specialized on the plant genetic lineages. Using random forest approximate Bayesian computation (ABC‐RF), we found that the divergence history of the fungal lineages on S. nutans was congruent with that previously inferred for the host plant and probably occurred with ancient but no recent gene flow. Genome sequences confirmed the genetic structure and the absence of recent gene flow between fungal genetic lineages. Our analyses of individual host–pathogen pairs contribute to a better understanding of co‐evolutionary histories between hosts and pathogens in natural ecosystems, in which such studies remain scarce.     

Genetic reconstruction of a bullfrog invasion to elucidate vectors of introduction and secondary spread

Reconstructing historical colonization pathways of an invasive species is critical for uncovering factors that determine invasion success and for designing management strategies. The American bullfrog (Lithobates catesbeianus) is endemic to eastern North America, but now has a global distribution and is considered to be one of the worst invaders in the world. In Montana, several introduced populations have been reported, but little is known of their sources and vectors of introduction and secondary spread. We evaluated the genetic composition of introduced populations at local (Yellowstone River floodplain) and regional (Montana and Wyoming) scales in contrast to native range populations. Our objectives were to (1) estimate the number of introductions, (2) identify probable native sources, (3) evaluate genetic variation relative to sources, and (4) characterize properties of local‐ and regional‐scale spread. We sequenced 937 bp of the mitochondrial cytochrome b locus in 395 tadpoles collected along 100 km of the Yellowstone River, from three additional sites in MT and a proximate site in WY. Pairwise Φ_ST revealed high divergence among nonnative populations, suggesting at least four independent introductions into MT from diverse sources. Three cyt b haplotypes were identical to native haplotypes distributed across the Midwest and Great Lakes regions, and AMOVA confirmed the western native region as a likely source. While haplotype (H_d = 0.69) and nucleotide diversity (π = 0.005) were low in introduced bullfrogs, the levels of diversity did not differ significantly from source populations. In the Yellowstone, two identified haplotypes implied few introduction vectors and a significant relationship between genetic and river distance was found. Evidence for multiple invasions and lack of subsequent regional spread emphasizes the importance of enforcing legislation prohibiting bullfrog importation and the need for continuing public education to prevent transport of bullfrogs in MT. More broadly, this study demonstrates how genetic approaches can reveal key properties of a biological invasion to inform management strategies.     

European wildcat populations are subdivided into five main biogeographic groups: consequences of Pleistocene climate changes or recent anthropogenic fragmentation?

Extant populations of the European wildcat are fragmented across the continent, the likely consequence of recent extirpations due to habitat loss and over‐hunting. However, their underlying phylogeographic history has never been reconstructed. For testing the hypothesis that the European wildcat survived the Ice Age fragmented in Mediterranean refuges, we assayed the genetic variation at 31 microsatellites in 668 presumptive European wildcats sampled in 15 European countries. Moreover, to evaluate the extent of subspecies/population divergence and identify eventual wild × domestic cat hybrids, we genotyped 26 African wildcats from Sardinia and North Africa and 294 random‐bred domestic cats. Results of multivariate analyses and Bayesian clustering confirmed that the European wild and the domestic cats (plus the African wildcats) belong to two well‐differentiated clusters (average Ф_ST = 0.159, R_st  = 0.392, P > 0.001; Analysis of molecular variance [AMOVA]). We identified from c. 5% to 10% cryptic hybrids in southern and central European populations. In contrast, wild‐living cats in Hungary and Scotland showed deep signatures of genetic admixture and introgression with domestic cats. The European wildcats are subdivided into five main genetic clusters (average Ф_ST = 0.103, R_st  = 0.143, P > 0.001; AMOVA) corresponding to five biogeographic groups, respectively, distributed in the Iberian Peninsula, central Europe, central Germany, Italian Peninsula and the island of Sicily, and in north‐eastern Italy and northern Balkan regions (Dinaric Alps). Approximate Bayesian Computation simulations supported late Pleistocene–early Holocene population splittings (from c. 60 k to 10 k years ago), contemporary to the last Ice Age climatic changes. These results provide evidences for wildcat Mediterranean refuges in southwestern Europe, but the evolution history of eastern wildcat populations remains to be clarified. Historical genetic subdivisions suggest conservation strategies aimed at enhancing gene flow through the restoration of ecological corridors within each biogeographic units. Concomitantly, the risk of hybridization with free‐ranging domestic cats along corridor edges should be carefully monitored.     

Phylogeography and genetic connectivity of the marine macro‐alga Sargassum ilicifolium (Phaeophyceae,Ochrophyta) in the northwestern Pacific1

The evolutionary influences of historical and contemporary factors on the population connectivity and phylogeographic structure of a brown seaweed, Sargassum ilicifolium, were elucidated using the nuclear ITS2 and mitochondrial COI markers for the collections newly sampled within its distribution range in the northwestern Pacific (NWP). Significant genetic structure at variable levels was identified between populations (pairwise F_ST) and among populations grouped by geographical proximity (Φ_CT among regions). The adjacent groups of populations with moderate structure revealed from AMOVA appeared to have high genetic connectivity. However, a lack of genealogical concordance with the geographic distribution was uncovered for S. ilicifolium from the NWP. Such genetic homogeneity is interpreted as a result of the interaction between postglacial recolonization and dynamic oceanic current regimes in the region. Two separated glacial refugia, the South China Sea and the Okinawa Trough, in the marginal seas of east China were recognized based on the presence of endemic haplotypes and high haplotype diversity in the populations at southern China and northeast of Taiwan. Populations persisting in these refugia may have served as the source for recolonization in the NWP with the rise of sea level during the warmer interglacial periods. The role of oceanic currents in maintaining genetic connectivity of S. ilicifolium in the region was further corroborated by the coherence between the direction of oceanic currents and that of gene flow, especially along the eastern coast of Taiwan. This study underlines the interaction between historical postglacial recolonization and contemporary coastal hydrodynamics in contributing to population connectivity and distribution for this tropical seaweed in the NWP.