Colonization from divergent ancestors: glaciation signatures on contemporary patterns of genomic variation in Collared Pikas (Ochotona collaris)

Identifying the genetic structure of a species and the factors that drive it is an important first step in modern population management, in part because populations evolving from separate ancestral sources may possess potentially different characteristics. This is especially true for climate‐sensitive species such as pikas, where the delimitation of distinct genetic units and the characterization of population responses to contemporary and historical environmental pressures are of particular interest. We combined a restriction site‐associated DNA sequencing (RADSeq) data set containing 4156 single nucleotide polymorphisms with ecological niche models (ENMs) of present and past habitat suitability to characterize population composition and evaluate the effects of historical range shifts, contemporary climates and landscape factors on gene flow in Collared Pikas, which are found in Alaska and adjacent regions of northwestern Canada and are the lesser‐studied of North America's two pika species. The results suggest that contemporary environmental factors contribute little to current population connectivity. Instead, genetic diversity is strongly shaped by the presence of three ancestral lineages isolated during the Pleistocene (~148 and 52 kya). Based on ENMs and genetic data, populations originating from a northern refugium experienced longer‐term stability, whereas both southern lineages underwent population expansion – contradicting the southern stability and northern expansion patterns seen in many other taxa. Current populations are comparable with respect to generally low diversity within populations and little‐to‐no recent admixture. The predominance of divergent histories structuring populations implies that if we are to understand and manage pika populations, we must specifically assess and accurately account for the forces underlying genetic similarity.     

Barriers,rather than refugia,underlie the origin of diversity in toads endemic to the Brazilian Atlantic Forest

In this study, we investigated the relative contribution of geographic barriers and Pleistocene refuges in the diversification of the Rhinella crucifer species complex, a group of endemic toads with a widespread distribution in the Brazilian Atlantic Forest (AF). We used intensive sampling and multilocus DNA sequence data to compare nucleotide diversity between refuge and nonrefuge areas, investigate regional demographic patterns, estimate demographic parameters related to genetic breaks and test refuge versus barrier scenarios of diversification using approximate Bayesian computation. We did not find higher levels of genetic diversity in putative refuge areas, either at regional or biome scale. Rather, the demographic history of the species complex supports regional differences with moderate population growth in the north and central regions and stability in southern AF. Genetic breaks were dated to the Plio–Pleistocene; however, our analyses rejected the role of refuges in creating a northern and central divergence, supporting a recent colonization scenario at a smaller scale within the central AF. Overall, our data rule out massive climatically driven fragmentation and large‐scale recolonization events for populations across the biome. We confirmed the importance of geographic barriers in creating main divergences and underscored the importance of searching for cryptic discontinuities in the landscape. Comparison of our results with those of other AF taxa indicates organismal specific responses to moderate shifts in habitat and that multiple refuges may constitute a more realistic model for diversification of Atlantic Forest biota.     

East Asian allopatry and north Eurasian sympatry in Long‐tailed Tit lineages despite similar population dynamics during the late Pleistocene

Eurasia is a large continent characterized by heterogeneous environments. Glacial cycles during the late Pleistocene have had variable impacts on the avifauna across Eurasia. Bird populations from South‐East Asia show stability through the Last Glacial Maximum (LGM), while populations from Europe exhibit evidence of post‐LGM expansion. We investigated the phylogeography of the Long‐tailed Tit (Aegithalos caudatus), which spans the longitudinal breadth of Eurasia to test how climatic history and regional topographical complexity affected populations and diversification within the species complex. Our results show that two lineages from central and southern China (lineages C and D) segregate geographically, while lineages across northern Eurasia (lineage A and B) show substantial sympatry. Bayesian estimates for the timing of diversification suggest that the four lineages diverged during the middle Pleistocene, splitting in parallel and undergoing concurrent demographic histories since divergence. A. caudatus lineages experienced similar and synchronous population size dynamics during glacial cycles before the LGM. We conclude that the difference in geo‐topologic complexity may be an important factor that led to the variation in secondary admixture between northern Eurasian and eastern Asian lineages.     

Mountains too high and valleys too deep drive population structuring and demographics in a Qinghai–Tibetan Plateau frog Nanorana pleskei (Dicroglossidae)

Pleistocene glacial–interglacial climatic oscillations greatly shaped the current genetic structure of many species. However, geographic features may influence the impact of climatic cycling. Distinct geographic and environmental characters between northern and southern parts of the eastern Qinghai–Tibetan Plateau (EQTP) facilitate explorations into the impacts of geographic features on species. The northern parts of EQTP contain large areas of marsh, and the environment is rather homogeneous. In contrast, the southern EQTP harbors complex alpine valleys and a much more heterogeneous setting. We evaluate DNA sequence variation from both the mitochondrial and nuclear genomes in Nanorana pleskei, a species endemic to the EQTP. Hypothesis testing on the evolutionary history of N. pleskei indicates that northern populations can disperse freely, but alpine valleys isolate southern populations. Demographic histories between northern and southern populations also differ. Northern populations appear to have experienced population expansions, while southern frogs exhibit a far more stable demographic history. By combining climatic analyses and species' distribution models, our study suggests that geographic and environmental features drive the differences between the northern and southern EQTP.     

Late Pleistocene divergence and postglacial expansion in the Brazilian Atlantic Forest: multilocus phylogeography of Rhopias gularis (Aves: Passeriformes)

In the last decade, phylogeographic studies have revealed a complex evolutionary history of the Brazilian Atlantic Forest (AF) biota. Here, we investigated the evolutionary history of Rhopias gularis, an endemic bird of the AF, based on sequences of two mitochondrial genes and three nuclear introns from 64 specimens from 15 localities. We addressed three main questions: (1) Does the genetic diversity of R. gularis exhibit a distribution pattern congruent with the refuge hypothesis postulated for the AF? (2) Is the population genetic structure of R. gularis congruent with those observed in other AF species? (3) What were the possible historical events responsible for the population structure of this species? Our mtDNA data revealed two phylogroups: (1) phylogroup central‐south, with samples from the central and southern parts of the range; (2) and phylogroup north, which included individuals from southern Bahia. Nevertheless, nuclear loci did not reveal any evidence of population structure. Bottleneck tests indicated that the central‐south lineage experienced demographic expansion, starting around 20 kya, which coincides with the end of the last glacial maximum. However, there was no evidence of population growth in phylogroup north. Isolation with migration analysis indicated that these phylogroups split c.a. 304 kya, with limited gene flow among them. Palaeodistribution models indicated that R. gularis had a reduced distribution in the south and central AF during the last glacial maximum. Our results support a diversification scenario that is in accordance with proposed Pleistocene refugia. The phylogeographic results from our study exhibited spatial and temporal concordances and discordances with previous studies of organisms from the AF. Differences in habitat requirements of these species could be behind this complex scenario. Future studies correlating variables of the niche of these species with the observed phylogeographic patterns may help understand why there are congruent and incongruent results.     

Radiation of the polymorphic Little Devil poison frog (Oophaga sylvatica) in Ecuador

Some South American poison frogs (Dendrobatidae) are chemically defended and use bright aposematic colors to warn potential predators of their unpalatability. Aposematic signals are often frequency‐dependent where individuals deviating from a local model are at a higher risk of predation. However, extreme diversity in the aposematic signal has been documented in poison frogs, especially in Oophaga. Here, we explore the phylogeographic pattern among color‐divergent populations of the Little Devil poison frog Oophaga sylvatica by analyzing population structure and genetic differentiation to evaluate which processes could account for color diversity within and among populations. With a combination of PCR amplicons (three mitochondrial and three nuclear markers) and genome‐wide markers from a double‐digested RAD (ddRAD) approach, we characterized the phylogenetic and genetic structure of 199 individuals from 13 populations (12 monomorphic and 1 polymorphic) across the O. sylvatica distribution. Individuals segregated into two main lineages by their northern or southern latitudinal distribution. A high level of genetic and phenotypic polymorphism within the northern lineage suggests ongoing gene flow. In contrast, low levels of genetic differentiation were detected among the southern lineage populations and support recent range expansions from populations in the northern lineage. We propose that a combination of climatic gradients and structured landscapes might be promoting gene flow and phylogenetic diversification. Alternatively, we cannot rule out that the observed phenotypic and genomic variations are the result of genetic drift on near or neutral alleles in a small number of genes.     

Geographical variation in genetic structure of an Atlantic Coastal Forest frog reveals regional differences in habitat stability

Climatic oscillations throughout the Pleistocene combined with geological and topographic complexity resulted in extreme habitat heterogeneity along the Atlantic coast of Brazil. Inferring how these historic landscape patterns have structured the current diversity of the region's biota is important both for our understanding of the factors promoting diversification, as well as the conservation of this biodiversity hotspot. Here we evaluate potential historical scenarios of diversification in the Atlantic Coastal Forest of Brazil by investigating the population genetic structure of a frog endemic to the region. Using mitochondrial and nuclear sequences, we generated a Bayesian population-level phylogeny of the Thoropa miliaris species complex. We found deep genetic divergences among five geographically distinct clades. Southern clades were monophyletic and nested within paraphyletic northern clades. Analyses of historical demographic patterns suggest an overall north to south population expansion, likely associated with regional differences in habitat stability during the Pliocene and early Pleistocene. However, genetic structure among southern populations is less pronounced and likely represents more recent vicariant events resulting from Holocenic sea-level oscillations. Our analyses corroborate that the Atlantic Coastal Forest has been a biogeographically dynamic landscape and suggest that the high diversity of its fauna and flora resulted from a combination of climatic and geologic events from the Pliocene to the present.     

Evolutionary history of the Maltese wall lizard Podarcis filfolensis: insights on the ‘Expansion–Contraction’ model of Pleistocene biogeography

The expansion–contraction (EC) model predicts demographic and range contraction of temperate species during Pleistocene glaciations as a consequence of climate‐related habitat changes, and provides a paradigm for explaining the high intraspecific diversity found in refugia in terms of long‐term demographic stability. However, recent evidence has revealed a weak predictive power of this model for terrestrial species in insular and coastal settings. We investigated the Pleistocene EC dynamics and their evolutionary consequences on temperate species using the Maltese archipelago and its endemic lizard Podarcis filfolensis as a model system. The evolutionary and demographic history of P. filfolensis as inferred from mitochondrial and nuclear sequences data does not conform to the EC model predictions, supporting (i) demographic and spatial stability or expansion, rather than contraction, of the northern and southern lineages during the last glacial period; and (ii) a major role for allopatric differentiation primed by sea‐level dynamics, rather than prolonged demographic stability, in the formation of the observed genetic diversity. When combined with evidence from other Mediterranean refugia, this study shows how the incorporation of Pleistocene sea‐level variations in the EC model accounts for a reverse demographic and range response of insular and coastal temperate biotas relative to continental ones. Furthermore, this cross‐archipelago pattern in which allopatric diversity is formed and shaped by EC cycles resembles that seen between isolated populations within mainland refugia and suggests that the EC model, originally developed to explain population fluctuations into and out‐of refugia, may be appropriate for describing the demographic and evolutionary dynamics driving the high genetic diversity observed in these areas.     

Reconstructing the demographic history of orang‐utans using Approximate Bayesian Computation

Investigating how different evolutionary forces have shaped patterns of DNA variation within and among species requires detailed knowledge of their demographic history. Orang‐utans, whose distribution is currently restricted to the South‐East Asian islands of Borneo (Pongo pygmaeus) and Sumatra (Pongo abelii), have likely experienced a complex demographic history, influenced by recurrent changes in climate and sea levels, volcanic activities and anthropogenic pressures. Using the most extensive sample set of wild orang‐utans to date, we employed an Approximate Bayesian Computation (ABC) approach to test the fit of 12 different demographic scenarios to the observed patterns of variation in autosomal, X‐chromosomal, mitochondrial and Y‐chromosomal markers. In the best‐fitting model, Sumatran orang‐utans exhibit a deep split of populations north and south of Lake Toba, probably caused by multiple eruptions of the Toba volcano. In addition, we found signals for a strong decline in all Sumatran populations ~24 ka, probably associated with hunting by human colonizers. In contrast, Bornean orang‐utans experienced a severe bottleneck ~135 ka, followed by a population expansion and substructuring starting ~82 ka, which we link to an expansion from a glacial refugium. We showed that orang‐utans went through drastic changes in population size and connectedness, caused by recurrent contraction and expansion of rainforest habitat during Pleistocene glaciations and probably hunting by early humans. Our findings emphasize the fact that important aspects of the evolutionary past of species with complex demographic histories might remain obscured when applying overly simplified models.     

Population genetic structure of the endangered crayfish Austropotamobius pallipes in France based on microsatellite variation: biogeographical inferences and conservation implications

1. One important goal in conservation biology is to characterise evolutionary lineages within endangered species before management decisions are taken. Here, we assess population differentiation in the freshwater crayfish Austropotamobius pallipes, an endangered species endemic to western Europe and provide valuable information for the conservation of French populations. 2. Analysis of five microsatellite loci in 44 populations revealed very different within population levels of genetic diversity (0.000 < H₀ < 0.564). Two groups, corresponding to northern and southern French populations, showed a high degree of genetic differentiation in both allele frequencies and allele sizes. Comparison of these results with previous studies of A. pallipes strongly suggests that the divergence between northern and southern populations could have occurred during the last glaciation period of the Pleistocene from one Atlantic and one Mediterranean refuge. 3. Evidence for genetic admixture between these two lineages was revealed by correspondence analyses in southern populations, probably as the result of artificial translocations. 4. French populations appeared significantly differentiated among the different river drainages and were highly structured within rivers. The impact of population size, population bottlenecks and founder events on the population genetic differentiation are discussed. 5. Based on these results, we propose the designation of two evolutionarily significant units for A. pallipes in France. Our data also support the maintenance of separate demographic management strategies for crayfish inhabiting different river systems. However, genetic analyses will have to be combined with demographic and ecological data for sustainable conservation programmes.     

Phylogeography of the Trans‐Volcanic bunchgrass lizard (Sceloporus bicanthalis) across the highlands of south‐eastern Mexico

We quantify the population divergence processes that shaped population genetic structure in the Trans‐Volcanic bunchgrass lizard (Sceloporus bicanthalis) across the highlands of south‐eastern Mexico. Multilocus genetic data from nine nuclear loci and mitochondrial (mt)DNA were used to estimate the population divergence history for 47 samples of S. bicanthalis. Bayesian clustering methods partitioned S. bicanthalis into three populations: (1) a southern population in Oaxaca and southern Puebla; (2) a population in western Puebla; and (3) a northern population with a broad distribution across Hidalgo, Puebla, and Veracruz. The multilocus nuclear data and mtDNA both supported a Late Pleistocene increase in effective population size, and the nuclear data revealed low levels of unidirectional gene flow from the widespread northern population into the southern and western populations. Populations of S. bicanthalis experienced different demographic histories during the Pleistocene, and phylogeographical patterns were similar to those observed in many co‐distributed highland taxa. Although we recommend continuing to recognize S. bicanthalis as a single species, future research on the evolution of viviparity could gain novel insights by contrasting physiological and genomic patterns among the different populations located across the highlands of south‐eastern Mexico. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 852–865.     

Landscape genetics reveal broad and fine‐scale population structure due to landscape features and climate history in the northern leopard frog (Rana pipiens) in North Dakota

Prehistoric climate and landscape features play large roles structuring wildlife populations. The amphibians of the northern Great Plains of North America present an opportunity to investigate how these factors affect colonization, migration, and current population genetic structure. This study used 11 microsatellite loci to genotype 1,230 northern leopard frogs (Rana pipiens) from 41 wetlands (30 samples/wetland) across North Dakota. Genetic structure of the sampled frogs was evaluated using Bayesian and multivariate clustering methods. All analyses produced concordant results, identifying a major east–west split between two R. pipiens population clusters separated by the Missouri River. Substructuring within the two major identified population clusters was also found. Spatial principal component analysis (sPCA) and variance partitioning analysis identified distance, river basins, and the Missouri River as the most important landscape factors differentiating R. pipiens populations across the state. Bayesian reconstruction of coalescence times suggested the major east–west split occurred ~13–18 kya during a period of glacial retreat in the northern Great Plains and substructuring largely occurred ~5–11 kya during a period of extreme drought cycles. A range‐wide species distribution model (SDM) for R. pipiens was developed and applied to prehistoric climate conditions during the Last Glacial Maximum (21 kya) and the mid‐Holocene (6 kya) from the CCSM4 climate model to identify potential refugia. The SDM indicated potential refugia existed in South Dakota or further south in Nebraska. The ancestral populations of R. pipiens in North Dakota may have inhabited these refugia, but more sampling outside the state is needed to reconstruct the route of colonization. Using microsatellite genotype data, this study determined that colonization from glacial refugia, drought dynamics in the northern Great Plains, and major rivers acting as barriers to gene flow were the defining forces shaping the regional population structure of R. pipiens in North Dakota.     

Small body size increases the regional differentiation of populations of tropical mantellid frogs (Anura: Mantellidae)

The processes affecting species diversification may also exert an influence on patterns of genetic variability within species. We evaluated the contributions of five variables potentially influencing clade diversification (body size, reproductive mode, range size, microhabitat and skin texture) on mtDNA divergence and polymorphism among populations of 40 species of frogs (Mantellidae) from two rainforest communities in Madagascar. We report an inverse association between body size and nucleotide divergence between populations but find no influence of other variables on genetic variation. Body size explained ca. 11% of the variation in nucleotide divergence between populations and was coupled with high F_ST levels and an absence of haplotype sharing in small‐bodied and medium‐sized frogs. Low dispersal ability is likely the proximate mechanism producing higher population differentiation in small mantellids. The lack of genetic cohesion among populations establishes regional genetic fragmentation which in turn has the potential to accelerate rates of allopatric speciation in small frogs relative to large species. However, there is little evidence of increased speciation rates in these or other small‐bodied organisms. We reconcile these contradictory observations by suggesting that lower dispersal ability also curbs colonization of new areas, decelerating diversification in weak dispersers. Our results imply that the intermediate dispersal model also applies to amphibians and may explain inconsistent previous results on the correlation of body size and speciation rate.     

Signatures of Demography and Recombination at Coding Genes in Naturally-Distributed Populations of Arabidopsis Lyrata Subsp. Petraea

Demography impacts the observed standing level of genetic diversity present in populations. Distinguishing the relative impacts of demography from selection requires a baseline of expressed gene variation in naturally occurring populations. Six nuclear genes were sequenced to estimate the patterns and levels of genetic diversity in natural Arabidopsis lyrata subsp. petraea populations that differ in demographic histories since the Pleistocene. As expected, northern European populations have genetic signatures of a strong population bottleneck likely due to glaciation during the Pleistocene. Levels of diversity in the northern populations are about half of that in central European populations. Bayesian estimates of historical population size changes indicate that central European populations also have signatures of population size change since the last glacial maxima, suggesting that these populations are not as stable as previously thought. Time since divergence amongst northern European populations is higher than amongst central European populations, suggesting that the northern European populations were established before the Pleistocene and survived glaciation in small separated refugia. Estimates of demography based on expressed genes are complementary to estimates based on microsatellites and transposable elements, elucidating temporal shifts in population dynamics and confirming the importance of marker selection for tests of demography.     

Late Pleistocene human population bottlenecks, volcanic winter, and differentiation of modern humans

The “Weak Garden of Eden” model for the origin and dispersal of modern humans (Harpendinget al., 1993) posits that modern humans spread into separate regions from a restricted source, around 100 ka (thousand years ago), then passed through population bottlenecks. Around 50 ka, dramatic growth occurred within dispersed populations that were genetically isolated from each other. Population growth began earliest in Africa and later in Eurasia and is hypothesized to have been caused by the invention and spread of a more efficient Later Stone Age/Upper Paleolithic technology, which developed in equatorial Africa.Climatic and geological evidence suggest an alternative hypothesis for Late Pleistocene population bottlenecks and releases. The last glacial period was preceded by one thousand years of the coldest temperatures of the Later Pleistocene (∼71–70 ka), apparently caused by the eruption of Toba, Sumatra. Toba was the largest known explosive eruption of the Quaternary. Toba's volcanic winter could have decimated most modern human populations, especially outside of isolated tropical refugia. Release from the bottleneck could have occurred either at the end of this hypercold phase, or 10,000 years later, at the transition from cold oxygen isotope stage 4 to warmer stage 3. The largest populations surviving through the bottleneck should have been found in the largest tropical refugia, and thus in equatorial Africa. High genetic diversity in modern Africans may thus reflect a less severe bottleneck rather than earlier population growth.Volcanic winter may have reduced populations to levels low enough for founder effects, genetic drift and local adaptations to produce rapid population differentiation. If Toba caused the bottlenecks, then modern human races may have differentiated abruptly, only 70 thousand years ago.     

Strong intraspecific variation in genetic diversity and genetic differentiation in Daphnia magna: the effects of population turnover and population size

Theory predicts that genetic diversity and genetic differentiation may strongly vary among populations of the same species depending on population turnover and local population sizes. Yet, despite the importance of these predictions for evolutionary and conservation issues, empirical studies comparing high‐turnover and low‐turnover populations of the same species are scarce. In this study, we used Daphnia magna, a freshwater crustacean, as a model organism for such a comparison. In the southern/central part of its range, D. magna inhabits medium‐sized, stable ponds, whereas in the north, it occurs in small rock pools with strong population turnover. We found that these northern populations have a significantly lower genetic diversity and higher genetic differentiation compared to the southern/central populations. Total genetic diversity across populations was only about half and average within‐population diversity only about a third of that in southern/central populations. Moreover, an average southern population contains more genetic diversity than the whole metapopulation system in the north. We based our analyses both on silent sites and microsatellites. The similarity of our results despite the contrasting mutation rates of these markers suggests that the differences are caused by contemporary rather than by historical processes. Our findings show that variation in population turnover and population size may have a major impact on the genetic diversity and differentiation of populations, and hence may lead to differences in evolutionary processes like local adaptation, hybrid vigour and breeding system evolution in different parts of a species range.     

Evolutionary history of Scinax treefrogs on land‐bridge islands in south‐eastern Brazil

Aim We investigated how Pleistocene refugia and recent (c. 12,000 years ago) sea level incursions shaped genetic differentiation in mainland and island populations of the Scinax perpusillus treefrog group. Location Brazilian Atlantic Forest, São Paulo state, south‐eastern Brazil. Methods Using mitochondrial and microsatellite loci, we examined population structure and genetic diversity in three species from the S. perpusillus group, sampled from three land‐bridge islands and five mainland populations, in order to understand the roles of Pleistocene forest fragmentation and sea level incursions on genetic differentiation. We calculated metrics of relatedness and genetic diversity to assess whether island populations exhibit signatures of genetic drift and isolation. Two of the three island populations in this study have previously been described as new species based on a combination of distinct morphological and behavioural characters, thus we used the molecular datasets to determine whether phenotypic change is consistent with genetic differentiation. Results Our analyses recovered three distinct lineages or demes composed of northern mainland São Paulo populations, southern mainland São Paulo populations, and one divergent island population. The two remaining island populations clustered with samples from adjacent mainland populations. Estimates of allelic richness were significantly lower, and estimates of relatedness were significantly higher, in island populations relative to their mainland counterparts. Main conclusions Fine‐scale genetic structure across mainland populations indicates the possible existence of local refugia within São Paulo state, underscoring the small geographic scale at which populations diverge in this species‐rich region of the Atlantic Coastal Forest. Variation in genetic signatures across the three islands indicates that the populations experienced different demographic processes after marine incursions fragmented the distribution of the S. perpusillus group. Genetic signatures of inbreeding and drift in some island populations indicate that small population sizes, coupled with strong ecological selection, may be important evolutionary forces driving speciation on land‐bridge islands.     

Phylogeography and genetic differentiation along the distributional range of the orchid Epidendrum fulgens: a Neotropical coastal species not restricted to glacial refugia

Aim Phylogeographical studies in the Brazilian Atlantic Forest (BAF) have mostly included species associated with forest habitats, whereas taxa associated with grassland and sand‐dune plant communities have so far been largely overlooked. This study examines the phylogeography of the orchid Epidendrum fulgens, which occurs on coastal sand dunes and granitic outcrops, in order to identify major genetic divergences or disjunctions across the range of the species and to investigate the genetic signatures of past range contractions and expansions. Location Southern and south‐eastern seashore vegetation along the BAF biome, and granitic and arenitic outcrops that occur in the subtropical grassland plant communities located south of the BAF. Methods Nine nuclear and four plastid microsatellite loci were used to genotype 424 individuals from 16 populations across the distributional range of E. fulgens. For both sets of markers, we estimated genetic diversity and population differentiation, testing for a north–south gradient of genetic diversity. The plastid haplotype network and a Bayesian assignment analysis of nuclear markers were used to infer population structure. Past demographic changes were investigated using a coalescence approach. Results A deep disjunction was found between northern populations within the BAF and southern populations outside the BAF that occur on granitic and arenitic outcrops. Recent demographic reductions were detected in northern populations on coastal sands. Such demographic changes were not expected for those populations, as previous studies with forest species had found evidence of population expansion in the same areas. Higher genetic diversity was found in southern populations on granite, in contrast to patterns observed in previous studies of forest species. Main conclusions The results are consistent with the long‐term persistence of E. fulgens. Bottlenecks were detected in populations from areas where population expansion events have been detected in other plant (and animal) species, suggesting that forest expansion after the Last Glacial Maximum played a role in the population fragmentation and decrease in genetic diversity in E. fulgens. A substantial genetic division in E. fulgens corresponds to the ‘Portal de Torres’, a region that demarcates the northern limits of subtropical grassland plant communities and the southern limits of the BAF.     

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. F_ST outlier analysis revealed that promoter, 5′‐UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within‐population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.     

Genetic Structure and Demographic History Reveal Migration of the Diamondback Moth Plutella xylostella (Lepidoptera: Plutellidae) from the Southern to Northern Regions of China

The diamondback moth Plutella xylostella (Linnaeus) (Lepidoptera: Plutellidae) is one of the most destructive insect pests of cruciferous plants worldwide. Biological, ecological and genetic studies have indicated that this moth is migratory in many regions around the world. Although outbreaks of this pest occur annually in China and cause heavy damage, little is known concerning its migration. To better understand its migration pattern, we investigated the population genetic structure and demographic history of the diamondback moth by analyzing 27 geographical populations across China using four mitochondrial genes and nine microsatellite loci. The results showed that high haplotype diversity and low nucleotide diversity occurred in the diamondback moth populations, a finding that is typical for migratory species. No genetic differentiation among all populations and no correlation between genetic and geographical distance were found. However, pairwise analysis of the mitochondrial genes has indicated that populations from the southern region were more differentiated than those from the northern region. Gene flow analysis revealed that the effective number of migrants per generation into populations of the northern region is very high, whereas that into populations of the southern region is quite low. Neutrality testing, mismatch distribution and Bayesian Skyline Plot analyses based on mitochondrial genes all revealed that deviation from Hardy-Weinberg equilibrium and sudden expansion of the effective population size were present in populations from the northern region but not in those from the southern region. In conclusion, all our analyses strongly demonstrated that the diamondback moth migrates within China from the southern to northern regions with rare effective migration in the reverse direction. Our research provides a successful example of using population genetic approaches to resolve the seasonal migration of insects.