Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

Gene flow from multiple sources maintains high genetic diversity and stable population history of Common Moorhen Gallinula chloropus in China

In this study, we used mitochondrial control sequences and microsatellite data from 231 Common Moorhen Gallinula chloropus individuals sampled from 19 sites in China to analyse their genetic structure and evolutionary history. High genetic diversity was found for all populations, although microsatellite analysis showed that the genetic diversity in non‐migratory populations was significantly higher than in migratory populations. High gene flow occurred between neighbouring populations, although long‐distance gene flow also occurred. The Huazhong population was the single greatest genetic source for other populations. High gene flow probably led to the shallow genetic structure that we observed. Demographic expansion was found in migratory populations, non‐migratory populations and with all individuals combined. The expansion time for all populations combined was estimated to be 221 000 years ago. The Common Moorhen population grew rapidly during the interglacial before the last glacial maximum (LGM), then remained generally stable from the LGM to the present.     

Postglacial range shift and demographic expansion of the marine intertidal snail Batillaria attramentaria

To address the impacts of past climate changes, particularly since the last glacial period, on the history of the distribution and demography of marine species, we investigated the evolutionary and demographic responses of the intertidal batillariid gastropod, Batillaria attramentaria, to these changes, using the snail as a model species in the northwest Pacific. We applied phylogeographic and divergence population genetic approaches to mitochondrial COI sequences from B. attramentaria. To cover much of its distributional range, 197 individuals collected throughout Korea and 507 publically available sequences (mostly from Japan) were used. Finally, a Bayesian skyline plot (BSP) method was applied to reconstruct the demographic history of this species. We found four differentiated geographic groups around Korea, confirming the presence of two distinct, geographically subdivided haplogroups on the Japanese coastlines along the bifurcated routes of the warm Tsushima and Kuroshio Currents. These two haplogroups were estimated to have begun to split approximately 400,000 years ago. Population divergence analysis supported the hypothesis that the Yellow Sea was populated by a northward range expansion of a small fraction of founders that split from a southern ancestral population since the last glacial maximum (LGM: 26,000–19,000 years ago), when the southern area became re‐submerged. BSP analyses on six geographically and genetically defined groups in Korea and Japan consistently demonstrated that each group has exponentially increased approximately since the LGM. This study resolved the phylogeography of B. attramentaria as a series of events connected over space and time; while paleoceanographic conditions determining the connectivity of neighboring seas in East Asia are responsible for the vicariance of this species, the postglacial sea‐level rise and warming temperatures have played a crucial role in rapid range shifts and broad demographic expansions of its populations.     

Genetic patterns and changes in availability of suitable habitat support a colonisation history of a North American perennial plant

• Climatic fluctuations during the Pleistocene influenced the geographical distribution of plant species across the southern region of California. Following an integrative approach, we combined genetic data analysis with Environmental Niche Models (ENMs) to assess the historical range expansion of Yucca schidigera, a long‐lived desert perennial native of the Baja California Peninsula.
• We genotyped 240 individuals with seven nuclear microsatellite to investigate genetic diversity distribution across 13 populations. Indeed, we used Environmental Niche Models to examine the changes on the distribution of suitable climatic conditions for this species during the LIG (~120 ka), LGM (~22 ka) and Mid Holocene (~6 ka).
• We detected high genetic diversity across Y. schidigera populations (A_R = 9.94 ± 0.38 SE; H_exp = 0.791 ± 0.011 SE) with genetic variation decreasing significantly with latitude (allelic richness: R ² = 0.38, P = 0.023; expected heterocigosity: R² = 0.32, P = 0.042). We observed low, but significant genetic differentiation (F_ST = 0.0678; P < 0.001) which was consistent with the parapatric distribution of the three genetic groupings detected by the Bayesian clustering algorithm. The ENMs suggest that suitable habitat for this species increased since the LGM.
• Our results support a range expansion of Y. schidigera across northwestern Baja California during the late Quaternary. Genetic data suggest that colonization of the current distribution followed a southward directionality as suitable climatic conditions became widely available in this region. High genetic variation across our sample suggests large historic effective population sizes for this section of the geographical range.

     

Winter climate change and coastal wetland foundation species: salt marshes vs. mangrove forests in the southeastern United States

We live in an era of unprecedented ecological change in which ecologists and natural resource managers are increasingly challenged to anticipate and prepare for the ecological effects of future global change. In this study, we investigated the potential effect of winter climate change upon salt marsh and mangrove forest foundation species in the southeastern United States. Our research addresses the following three questions: (1) What is the relationship between winter climate and the presence and abundance of mangrove forests relative to salt marshes; (2) How vulnerable are salt marshes to winter climate change‐induced mangrove forest range expansion; and (3) What is the potential future distribution and relative abundance of mangrove forests under alternative winter climate change scenarios? We developed simple winter climate‐based models to predict mangrove forest distribution and relative abundance using observed winter temperature data (1970–2000) and mangrove forest and salt marsh habitat data. Our results identify winter climate thresholds for salt marsh–mangrove forest interactions and highlight coastal areas in the southeastern United States (e.g., Texas, Louisiana, and parts of Florida) where relatively small changes in the intensity and frequency of extreme winter events could cause relatively dramatic landscape‐scale ecosystem structural and functional change in the form of poleward mangrove forest migration and salt marsh displacement. The ecological implications of these marsh‐to‐mangrove forest conversions are poorly understood, but would likely include changes for associated fish and wildlife populations and for the supply of some ecosystem goods and services.     

Comparative phylogeography and demographic history of two marsupials of the Atlantic Forest in eastern Brazil

The Atlantic Forest biodiversity hotspot in eastern South America has been the focus of several phylogeographic studies concerning relationships between populations and areas and how taxa respond to environmental changes. We infer and compare the demographic and biogeographic histories of two didelphid marsupial species, Gracilinanus microtarsus and Marmosops incanus, from the Atlantic Forest of eastern Brazil to determine how these species responded to environmental changes over time, using mitochondrial and nuclear DNA sequences. We found great intraspecific genetic divergence in both species and a strong geographic structure related to similar and spatially cohesive groups within each species. These groups are consistent with the same topographical barriers, such as mountains and river valleys. Intraspecific clades are very old, dating back to a period of tectonic activities in the Neogene (5.39–8.57 Mya). Changes in the environment over the last 7 million years lead to fairly concordant demographic changes in both marsupial species, including population expansion during the last glacial maximum (ca. 21,000 years ago) or last interglacial (ca. 120,000 years ago) or both. These results do not fit the Pleistocene refuge hypothesis as an explanation of the historical biogeography and diversification of both species in the Atlantic Forest, but are compatible with the Atlantis Forest hypothesis.     

Microsatellite‐based analysis of genetic structure and gene flow of Mythimna separata (Walker) (Lepidoptera: Noctuidae) in China

The oriental armyworm, Mythimna separata, is a serious agricultural pest in China. Seasonal and roundtrip migration has recently led to sudden, localized outbreaks and crop losses. To evaluate genetic differentiation between populations in eastern and western China and elucidate gene flow, the genetic structure of 20 natural populations from nine provinces was examined using seven microsatellite markers. The results indicated high genetic diversity. However, little to moderate (0 < F_ST < 0.15) genetic differentiation was detected, and there was no correlation between genetic distance and geographical distance. Bayesian clustering analysis identified three groups whereas discriminant analysis of principal components identified ten clusters that were considered as two clear‐cut clusters and one admixed group. Gene flow occurred frequently in most population pairs, and an asymmetrical migration rate was detected in several pairwise population comparisons. The bottleneck test showed that few populations had experienced recent bottlenecks. Correspondingly, large‐scale and long‐distance migration of M. separata has caused low genetic differentiation and frequent gene exchange. Our findings are important for studying genetic evolution and help to improve predictions of M. separata outbreaks in China.     

Unravelling the evolutionary history and future prospects of endemic species restricted to former glacial refugia

The contemporary distribution and genetic composition of biodiversity bear a signature of species’ evolutionary histories and the effects of past climatic oscillations. For many European species, the Mediterranean peninsulas of Iberia, Italy and the Balkans acted as glacial refugia and the source of range recolonization, and as a result, they contain disproportionately high levels of diversity. As these areas are particularly threatened by future climate change, it is important to understand how past climatic changes affected their biodiversity. We use an integrated approach, combining markers with different evolutionary rates and combining phylogenetic analysis with approximate Bayesian computation and species distribution modelling across temporal scales. We relate phylogeographic processes to patterns of genetic variation in Myotis escalerai, a bat species endemic to the Iberian Peninsula. We found a distinct population structure at the mitochondrial level with a strong geographic signature, indicating lineage divergence into separate glacial refugia within the Iberian refugium. However, microsatellite markers suggest higher levels of gene flow resulting in more limited structure at recent time frames. The evolutionary history of M. escalerai was shaped by the effects of climatic oscillations and changes in forest cover and composition, while its future is threatened by climatically induced range contractions and the role of ecological barriers due to competition interactions in restricting its distribution. This study warns that Mediterranean peninsulas, which provided refuge for European biodiversity during past glaciation events, may become a trap for limited dispersal and ecologically limited endemic species under future climate change, resulting in loss of entire lineages.     

Differential threshold effects of habitat fragmentation on gene flow in two widespread species of bush crickets

Effects of habitat fragmentation on genetic diversity vary among species. This may be attributed to the interacting effects of species traits and landscape structure. While widely distributed and abundant species are often considered less susceptible to fragmentation, this may be different if they are small sized and show limited dispersal. Under intensive land use, habitat fragmentation may reach thresholds at which gene flow among populations of small-sized and dispersal-limited species becomes disrupted. Here, we studied the genetic diversity of two abundant and widespread bush crickets along a gradient of habitat fragmentation in an agricultural landscape. We applied traditional (G(ST), θ) and recently developed (G'ST', D) estimators of genetic differentiation on microsatellite data from each of twelve populations of the grassland species Metrioptera roeselii and the forest-edge species Pholidoptera griseoaptera to identify thresholds of habitat fragmentation below which genetic population structure is affected. Whereas the grassland species exhibited a uniform genetic structuring (G(ST) = 0.020-0.033; D = 0.085-0.149) along the whole fragmentation gradient, the forest-edge species' genetic differentiation increased significantly from D < 0.063 (G(ST) < 0.018) to D = 0.166 (G(ST) = 0.074), once the amount of suitable habitat dropped below a threshold of 20% and its proximity decreased substantially at the landscape scale. The influence of fragmentation on genetic differentiation was qualitatively unaffected by the choice of estimators of genetic differentiation but quantitatively underestimated by the traditional estimators. These results indicate that even for widespread species in modern agricultural landscapes fragmentation thresholds exist at which gene flow among suitable habitat patches becomes restricted.     

Phylogeography and demographic history of Gyrodactylus konovalovi (Monogenoidea: Gyrodactylidae), an ectoparasite on the East Asia Amur minnow (Cyprinidae) in Central China

Gyrodactylus konovalovi is an ectoparasite on the Amur minnow (Rhynchocypris lagowskii) that is widely distributed in the cold fresh waters of East Asia. In the present study, the phylogeography and demographic history of G. konovalovi and the distribution of its host in the Qinling Mountains are examined. A total of 79 individual parasites was sequenced for a 528 bp region of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene, and 25 haplotypes were obtained. The substitution rate (dN/dS) was 0.068 and indicated purifying selection. Haplotype diversity (h) and nucleotide diversity (π) varied widely in the Qinling Mountains. Phylogenetic trees based on Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods and network analysis revealed that all haplotypes were consistently well‐supported in three different lineages, indicating a significant geographic distribution pattern. There was a significant positive correlation between genetic differentiation (F_st) and geographic distance. The results of mismatch distribution, neutrality test and Bayesian skyline plot analyses showed that whole populations underwent population contraction during the Pleistocene. Based on the molecular clock calibration, the most common ancestor was estimated to have emerged in the middle Pleistocene. Our study suggests for the first time that a clearly phylogeography of G. konovalovi was shaped by geological events and climate fluctuations, such as orogenesis, drainage capture changes, and vicariance, during the Pleistocene in the Qinling Mountains.     

Population structure and gene flow in the global pest,Helicoverpa armigera

Helicoverpa armigera is a major agricultural pest that is distributed across Europe, Asia, Africa and Australasia. This species is hypothesized to have spread to the Americas 1.5 million years ago, founding a population that is at present, a distinct species, Helicoverpa zea. In 2013, H. armigera was confirmed to have re‐entered South America via Brazil and subsequently spread. The source of the recent incursion is unknown and population structure in H. armigera is poorly resolved, but a basic understanding would highlight potential biosecurity failures and determine the recent evolutionary history of region‐specific lineages. Here, we integrate several end points derived from high‐throughput sequencing to assess gene flow in H. armigera and H. zea from populations across six continents. We first assemble mitochondrial genomes to demonstrate the phylogenetic relationship of H. armigera with other Heliothine species and the lack of distinction between populations. We subsequently use de novo genotyping‐by‐sequencing and whole‐genome sequences aligned to bacterial artificial chromosomes, to assess levels of admixture. Primarily, we find that Brazilian H. armigera are derived from diverse source populations, with strong signals of gene flow from European populations, as well as prevalent signals of Asian and African ancestry. We also demonstrate a potential field‐caught hybrid between H. armigera and H. zea, and are able to provide genomic support for the presence of the H. armigera conferta subspecies in Australasia. While structure among the bulk of populations remains unresolved, we present distinctions that are pertinent to future investigations as well as to the biosecurity threat posed by H. armigera.     

Roads to isolation: Similar genomic history patterns in two species of freshwater crabs with contrasting environmental tolerances and range sizes

Freshwater species often show high levels of endemism and risk of extinction owing to their limited dispersal abilities. This is exemplified by the stenotopic freshwater crab, Johora singaporensis which is one of the world's 100 most threatened species, and currently inhabits less than 0.01 km² of five low order hill streams within the highly urbanized island city‐state of Singapore. We compared populations of J. singaporensis with that of the non‐threatened, widespread, abundant, and eurytopic freshwater crab, Parathelphusa maculata, and found surprisingly high congruence between their population genomic histories. Based on 2,617 and 2,470 genome‐wide SNPs mined via the double‐digest restriction‐associated DNA sequencing method for ~90 individuals of J. singaporensis and P. maculata, respectively, the populations are strongly isolated (F_ST = 0.146–0.371), have low genetic diversity for both species (also for COI), and show signatures of recent genetic bottlenecks. The most genetically isolated populations for both species are separated from other populations by one of the oldest roads in Singapore. These results suggest that anthropogenic developments may have impacted stream‐dependent species in a uniform manner, regardless of ubiquity, habitat preference, or dispersal modes of the species. While signs of inbreeding were not detected for the critically endangered species, the genetic distinctiveness and low diversity of the populations call for genetic rescue and connecting corridors between the remaining fragments of the natural habitat.     

No evidence of divergence at neutral genetic markers between the two morphologically different subspecies of the most numerous Arctic seabird

Identifying natural populations that might be considered separate units using morphology, genotype or both is important in understanding the process of speciation and for conservation. We examined the relationships between the only two subspecies of the most numerous Arctic seabird, the Little Auk Alle alle, using both morphological (wing and head‐bill lengths) and genetic data (482 base pairs of the mitochondrial control region and seven nuclear microsatellite loci). We found significant morphological differences between the subspecies, A. a. polaris being significantly larger than the nominate A. a. alle. However, we did not find the subspecies to be differentiated at either mitochondrial DNA or at microsatellite loci. Consequently, one evolutionary significant unit is proposed. The similarity of the two subspecies at neutral genetic markers may be due to contemporary gene flow between populations, as well as large population sizes both in the present and in the past, combined with recent post‐glacial colonization of the Artic.     

Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter‐ and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.