Phylogeography of the mottled spinefoot Siganus fuscescens: Pleistocene divergence and limited genetic connectivity across the Philippine archipelago

Historical isolation during Pleistocene low sea level periods is thought to have contributed to divergence among marine basin populations across the Coral Triangle. In the Philippine archipelago, populations in the South China Sea, Sulu Sea–inland seas, and Philippine Sea‐Celebes Sea basins might have been partially isolated. Meanwhile, present‐day broadscale oceanographic circulation patterns suggest connectivity between these basins. To evaluate hypotheses regarding the influence of historical and contemporary factors on genetic structure, phylogeographic patterns based on mitochondrial control region sequences for a reef‐associated fish, Siganus fuscescens, were analysed. Three distinct lineages were recovered. One lineage was identified as the morphologically similar species Siganus canaliculatus, while two lineages are monophyletic with S. fuscescens. Clade divergence and demographic expansion in S. fuscescens occurred during the Pleistocene. A strong signal of latitudinal structure was detected (Φ_CT = 0.188), driven by marked differences in clade distribution: one clade is widely distributed (clade A), while a second clade (clade B) has a restricted northern distribution. Regional structure of clade A is consistent with the basin isolation hypothesis (Φ_CT = 0.040) and suggests isolation of the South China Sea (Φ_CT = 0.091). Fine‐scale structure was observed in the South China Sea and south Philippine Sea, while Sulu Sea and inland seas were unstructured. Genetic structure across multiple spatial scales (archipelagic, regional, and fine‐scale within basins) suggests the influence of vicariant barriers and contemporary limits to gene flow in S. fuscescens that may be influenced by oceanographic circulation, geographical distance between available habitats, and latitudinal temperature differences.     

Phylogeographic mitogenomics of Atlantic cod Gadus morhua: Variation in and among trans‐Atlantic,trans‐Laurentian,Northern cod,and landlocked fjord populations

The historical phylogeography, biogeography, and ecology of Atlantic cod (Gadus morhua) have been impacted by cyclic Pleistocene glaciations, where drops in sea temperatures led to sequestering of water in ice sheets, emergence of continental shelves, and changes to ocean currents. High‐resolution, whole‐genome mitogenomic phylogeography can help to elucidate this history. We identified eight major haplogroups among 153 fish from 14 populations by Bayesian, parsimony, and distance methods, including one that extends the species coalescent back to ca. 330 kya. Fish from the Barents and Baltic Seas tend to occur in basal haplogroups versus more recent distribution of fish in the Northwest Atlantic. There was significant differentiation in the majority of trans‐Atlantic comparisons (Φ_ST = .029–.180), but little or none in pairwise comparisons within the Northwest Atlantic of individual populations (Φ_ST = .000–.060) or defined management stocks (Φ_ST = .000–.023). Monte Carlo randomization tests of population phylogeography showed significantly nonrandom trans‐Atlantic phylogeography versus absence of such structure within various partitions of trans‐Laurentian, Northern cod (NAFO 2J3KL) and other management stocks, and Flemish Cap populations. A landlocked meromictic fjord on Baffin Island comprised multiple identical or near‐identical mitogenomes in two major polyphyletic clades, and was significantly differentiated from all other populations (Φ_ST = .153–.340). The phylogeography supports a hypothesis of an eastern origin of genetic diversity ca. 200–250 kya, rapid expansion of a western superhaplogroup comprising four haplogroups ca. 150 kya, and recent postglacial founder populations.     

Limited connectivity and a phylogeographic break characterize populations of the pink anemonefish,Amphiprion perideraion,in the Indo‐Malay Archipelago: inferences from a mitochondrial and microsatellite loci

To enhance the understanding of larval dispersal in marine organisms, species with a sedentary adult stage and a pelagic larval phase of known duration constitute ideal candidates, because inferences can be made about the role of larval dispersal in population connectivity. Members of the immensely diverse marine fauna of the Indo‐Malay Archipelago are of particular importance in this respect, as biodiversity conservation is becoming a large concern in this region. In this study, the genetic population structure of the pink anemonefish, Amphiprion perideraion, is analyzed by applying 10 microsatellite loci as well as sequences of the mitochondrial control region to also allow for a direct comparison of marker‐derived results. Both marker systems detected a strong overall genetic structure (Φ_ST = 0.096, P < 0.0001; mean D_est = 0.17; F_ST = 0.015, P < 0.0001) and best supported regional groupings (Φ_CT = 0.199 P < 0.0001; F_CT = 0.018, P < 0.001) that suggested a differentiation of the Java Sea population from the rest of the archipelago. Differentiation of a New Guinea group was confirmed by both markers, but disagreed over the affinity of populations from west New Guinea. Mitochondrial data suggest higher connectivity among populations with fewer signals of regional substructure than microsatellite data. Considering the homogenizing effect of only a few migrants per generation on genetic differentiation between populations, marker‐specific results have important implications for conservation efforts concerning this and similar species.     

The presence of a cryptic barrier in the West Pacific Ocean suggests the effect of glacial climate changes on a widespread sea‐dispersed plant,Vigna marina (Fabaceae)

Ocean currents are an important driver of evolution for sea‐dispersed plants, enabling them to maintain reciprocal gene flow via sea‐dispersed diaspores and obtain wide distribution ranges. Although geographic barriers are known to be the primary factors shaping present genetic structure of sea‐dispersed plants, cryptic barriers which form clear genetic structure within oceanic regions are poorly understood. To test the presence of a cryptic barrier, we conducted a phylogeographic study together with past demographic inference for a widespread sea‐dispersed plant, Vigna marina, using 308 individuals collected from the entire Indo‐West Pacific (IWP) region. Chloroplast DNA variation showed strong genetic structure that separated populations into three groups: North Pacific (NP), South Pacific (SP) and Indian Ocean (IN) (F′_CT among groups = 0.954–1.000). According to the Approximate Bayesian computation inference, splitting time between NP and SP was approximately 20,200 years (95%HPD, 4,530–95,400) before present. Moreover, a signal of recent population expansion was detected in the NP group. This study clearly showed the presence of a cryptic barrier in the West Pacific region of the distributional range of V. marina. The locations of the cryptic barrier observed in V. marina corresponded to the genetic breaks found in other plants, suggesting the presence of a common cryptic barrier for sea‐dispersed plants. Demographic inference suggested that genetic structure related to this cryptic barrier has been present since the last glacial maximum and may reflect patterns of past population expansion from refugia.     

Multilocus phylogeography of Australian teals (Anas spp.): a case study of the relationship between vagility and genetic structure

Biogeographic barriers potentially restrict gene flow but variation in dispersal or vagility can influence the effectiveness of these barriers among different species and produce characteristic patterns of population genetic structure. The objective of this study was to investigate interspecific and intraspecific genetic structure in two closely related species that differ in several life‐history characteristics. The grey teal Anas gracilis is geographically widespread throughout Australia with a distribution that crosses several recognized biogeographic barriers. This species has high vagility as its extensive movements track broad‐scale patterns in rainfall. In contrast, the closely related chestnut teal A. castanea is endemic to the mesic southeastern and southwestern regions of Australia and is more sedentary. We hypothesized that these differences in life‐history characteristics would result in more pronounced population structuring in the chestnut teal. We sequenced five nuclear loci (nuDNA) for 49 grey teal and 23 chestnut teal and compared results to published mitochondrial DNA (mtDNA) sequences. We used analysis of molecular variance to examine population structure, and applied coalescent based approaches to estimate demographic parameters. As predicted, chestnut teal were more strongly structured at both mtDNA and nuDNA (Φ_ST= 0.163 and 0.054, respectively) than were grey teal (Φ_ST < 0.0001 for both sets of loci). Surprisingly, a greater proportion of the total genetic variation was partitioned among populations within species (Φ_SC= 0.014 and 0.047 for nuDNA and mtDNA, respectively) than between the two species (Φ_CT < 0.0001 for both loci). The ‘Isolation with Migration’ coalescent model suggested a late Pleistocene divergence between the taxa, but remarkably, a deeper divergence between the southeastern and southwestern populations of chestnut teal. We conclude that dispersal potential played a prominent role in the structuring of populations within these species and that divergent selection associated with ecology and life history traits likely contributed to rapid and recent speciation in this pair.     

Exploring the role of Micronesian islands in the maintenance of coral genetic diversity in the Pacific Ocean

Understanding how genetic diversity is maintained across patchy marine environments remains a fundamental problem in marine biology. The Coral Triangle, located in the Indo‐West Pacific, is the centre of marine biodiversity and has been proposed as an important source of genetic diversity for remote Pacific reefs. Several studies highlight Micronesia, a scattering of hundreds of small islands situated within the North Equatorial Counter Current, as a potentially important migration corridor. To test this hypothesis, we characterized the population genetic structure of two ecologically important congeneric species of reef‐building corals across greater Micronesia, from Palau to the Marshall Islands. Genetic divergences between islands followed an isolation‐by‐distance pattern, with Acropora hyacinthus exhibiting greater genetic divergences than A. digitifera, suggesting different migration capabilities or different effective population sizes for these closely related species. We inferred dispersal distance using a biophysical larval transport model, which explained an additional 15–21% of the observed genetic variation compared to between‐island geographical distance alone. For both species, genetic divergence accumulates and genetic diversity diminishes with distance from the Coral Triangle, supporting the hypothesis that Micronesian islands act as important stepping stones connecting the central Pacific with the species‐rich Coral Triangle. However, for A. hyacinthus, the species with lower genetic connectivity, immigration from the subequatorial Pacific begins to play a larger role in shaping diversity than input from the Coral Triangle. This work highlights the enormous dispersal potential of broadcast‐spawning corals and identifies the biological and physical drivers that influence coral genetic diversity on a regional scale.     

Genetic structure of Octopus vulgaris (Cephalopoda,Octopodidae) in the central Mediterranean Sea inferred from the mitochondrial COIII gene

The polymorphism of the mitochondrial gene cytochrome oxidase III was studied in the Mediterranean octopus, Octopus vulgaris Cuvier, 1797. A total of 202 specimens from seven sampling sites were analysed with the aim of elucidating patterns of genetic structure in the central Mediterranean Sea and to give an insight into the phylogeny of the Octopus genus. Phylogenetic analyses showed that individuals from the central Mediterranean belong to the O. vulgaris species whose limits should nevertheless be clarified. Concerning genetic structure, two high-frequency haplotypes were present in all locations. The overall genetic divergence (Φ_ST = 0.05, P < 0.05) indicated a significant genetic structuring in the study area and an AMOVA highlighted a significant break between western and eastern Mediterranean basins (Φ_CT = 0.094, P < 0.05). Possible explanations for the observed patterns of genetic structuring are discussed with reference to their relevance for fisheries management.     

Seascape continuity plays an important role in determining patterns of spatial genetic structure in a coral reef fish

Detecting patterns of spatial genetic structure (SGS) can help identify intrinsic and extrinsic barriers to gene flow within metapopulations. For marine organisms such as coral reef fishes, identifying these barriers is critical to predicting evolutionary dynamics and demarcating evolutionarily significant units for conservation. In this study, we adopted an alternative hypothesis‐testing framework to identify the patterns and predictors of SGS in the Caribbean reef fish Elacatinus lori. First, genetic structure was estimated using nuclear microsatellites and mitochondrial cytochrome b sequences. Next, clustering and network analyses were applied to visualize patterns of SGS. Finally, logistic regressions and linear mixed models were used to identify the predictors of SGS. Both sets of markers revealed low global structure: mitochondrial Φ_ST = 0.12, microsatellite F_ST = 0.0056. However, there was high variability among pairwise estimates, ranging from no differentiation between sites on contiguous reef (Φ_ST = 0) to strong differentiation between sites separated by ocean expanses ≥ 20 km (maximum Φ_ST = 0.65). Genetic clustering and statistical analyses provided additional support for the hypothesis that seascape discontinuity, represented by oceanic breaks between patches of reef habitat, is a key predictor of SGS in E. lori. Notably, the estimated patterns and predictors of SGS were consistent between both sets of markers. Combined with previous studies of dispersal in E. lori, these results suggest that the interaction between seascape continuity and the dispersal kernel plays an important role in determining genetic connectivity within metapopulations.     

Population structure of the dusky pipefish (Syngnathus floridae) from the Atlantic and Gulf of Mexico,as revealed by mitochondrial DNA and microsatellite analyses

Aim To elucidate the historical phylogeography of the dusky pipefish (Syngnathus floridae) in the North American Atlantic and Gulf of Mexico ocean basins. Location Southern Atlantic Ocean and northern Gulf of Mexico within the continental United States. Methods A 394‐bp fragment of the mitochondrial cytochrome b gene and a 235‐bp fragment of the mitochondrial control region were analysed from individuals from 10 locations. Phylogenetic reconstruction, haplotype network, mismatch distributions and analysis of molecular variance were used to infer population structure between ocean basins and time from population expansion within ocean basins. Six microsatellite loci were also analysed to estimate population structure and gene flow among five populations using genetic distance methods (F_ST, Nei’s genetic distance), isolation by distance (Mantel’s test), coalescent‐based estimates of genetic diversity and migration patterns, Bayesian cluster analysis and bottleneck simulations. Results Mitochondrial analyses revealed significant structuring between ocean basins in both cytochrome b (Φ_ST = 0.361, P < 0.0001; Φ_CT = 0.312, P < 0.02) and control region (Φ_ST = 0.166, P < 0.0001; Φ_CT = 0.128, P < 0.03) sequences. However, phylogenetic reconstructions failed to show reciprocal monophyly in populations between ocean basins. Microsatellite analyses revealed significant population substructuring between all locations sampled except for the two locations that were in closest proximity to each other (global F_ST value = 0.026). Bayesian analysis of microsatellite data also revealed significant population structuring between ocean basins. Coalescent‐based analyses of microsatellite data revealed low migration rates among all sites. Mismatch distribution analysis of mitochondrial loci supports a sudden population expansion in both ocean basins in the late Pleistocene, with the expansion of Atlantic populations occurring more recently. Main conclusions Present‐day populations of S. floridae do not bear the mitochondrial DNA signature of the strong phylogenetic discontinuity between the Atlantic and Gulf coasts of North America commonly observed in other species. Rather, our results suggest that Atlantic and Gulf of Mexico populations of S. floridae are closely related but nevertheless exhibit local and regional population structure. We conclude that the present‐day phylogeographic pattern is the result of a recent population expansion into the Atlantic in the late Pleistocene, and that life‐history traits and ecology may play a pivotal role in shaping the realized geographical distribution pattern of this species.     

New Zealand triplefin fishes (family Tripterygiidae): contrasting population structure and mtDNA diversity within a marine species flock

Triplefin fishes (Family Tripterygiidae) dominate the New Zealand temperate coastal fish fauna in diversity (26 endemic species, 14 genera). Most species appear to have evolved as a local radiation and mostly occupy sympatric distributions throughout New Zealand. To investigate the forces driving current gene‐flow patterns and past evolutionary histories, we searched for common patterns of population genetic subdivision within eight species sampled throughout their distributions [mitchochondrial DNA (mtDNA) control region, n = 1086]. We hypothesised that common phylogeographical and population differentiation patterns would reveal past or ongoing physical processes, with differences reflecting stochastic or species‐specific processes. Striking differences between species were apparent, with strong phylogeographical structure detected in Grahamina capito and the estuarine species G. nigripenne. G. capito fell into three distinct geographically restricted lineages. G. nigripenne largely separated into northern and southern lineages (Φ_ST 0.834). Strong population structuring and isolation by distance was evident in Bellapiscis medius, B. lesleyae and Forsterygion lapillum (Φ_ST 0.686, 0.436 and 0.115, respectively). High gene flow was apparent in G. gymnota and Ruanoho whero, and F. varium. However, for the latter species, isolation was apparent with samples collected from the offshore Three Kings Islands. Overall, a strong relationship was found between habitat depth and population structure among species, and species inhabiting shallower water habitats showed lower genetic diversity with higher levels of population subdivision. High‐latitude populations generally showed low haplotype and nucleotide diversity. These data suggest that processes resulting from intraspecific differences in habitat preference, climatic histories and/or larval ecologies have subdivided populations of shallow water triplefin species.     

Long‐term panmixia in a cosmopolitan Indo‐Pacific coral reef fish and a nebulous genetic boundary with its broadly sympatric sister species

Phylogeographical studies have shown that some shallow‐water marine organisms, such as certain coral reef fishes, lack spatial population structure at oceanic scales, despite vast distances of pelagic habitat between reefs and other dispersal barriers. However, whether these dispersive widespread taxa constitute long‐term panmictic populations across their species ranges remains unknown. Conventional phylogeographical inferences frequently fail to distinguish between long‐term panmixia and metapopulations connected by gene flow. Moreover, marine organisms have notoriously large effective population sizes that confound population structure detection. Therefore, at what spatial scale marine populations experience independent evolutionary trajectories and ultimately species divergence is still unclear. Here, we present a phylogeographical study of a cosmopolitan Indo‐Pacific coral reef fish Naso hexacanthus and its sister species Naso caesius, using two mtDNA and two nDNA markers. The purpose of this study was two‐fold: first, to test for broad‐scale panmixia in N. hexacanthus by fitting the data to various phylogeographical models within a Bayesian statistical framework, and second, to explore patterns of genetic divergence between the two broadly sympatric species. We report that N. hexacanthus shows little population structure across the Indo‐Pacific and a range‐wide, long‐term panmictic population model best fit the data. Hence, this species presently comprises a single evolutionary unit across much of the tropical Indian and Pacific Oceans. Naso hexacanthus and N. caesius were not reciprocally monophyletic in the mtDNA markers but showed varying degrees of population level divergence in the two nuclear introns. Overall, patterns are consistent with secondary introgression following a period of isolation, which may be attributed to oceanographic conditions of the mid to late Pleistocene, when these two species appear to have diverged.     

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of F_ST and to evaluate whether highly variable microsatellite loci resulted in deflated F_ST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (F_CT = 0.010) and within (F_ST = 0.014–0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (P_CT > F_CT). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (P_ST > F_ST), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated F_ST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for P_ST–F_ST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.     

Whole‐genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole‐genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N_e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long‐term population dynamics of Branchiostoma. The inferred ancestral N_e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10^?9) or third glaciation (mutation rate = 9 × 10^?9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long‐term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.     

Cranial morphometrics and mitochondrial DNA sequences distinguish cryptic species of the longface emperor (Lethrinus olivaceus), an emblematic fish of Indo-West Pacific coral reefs

Range-wide morphometric variability (cranial measurements) and genetic variability (nucleotide sequences of the cytochrome b gene) were investigated in the longface emperor, Lethrinus olivaceus (Lethrinidae), an emblematic large predatory fish of Indo-West Pacific coral reefs. Two cranial morphotypes were observed, one present from the Indian Ocean to the Coral Triangle and the other one, from the Coral Triangle to the western Central Pacific. The two morphotypes are concordant with reciprocally monophyletic mitochondrial lineages separated by 9.5% net nucleotide distance. These results suggest an old evolutionary history for L. olivaceus, which consists of two distinct species (Lethrinus sp. A in the Indian Ocean and Coral Triangle, Lethrinus sp. B in the western Pacific Ocean), whose distribution ranges meet or overlap in the eastern part of the Coral Triangle, in Taiwan and in West Papua. Lethrinus sp. A comprises two distinct mitochondrial lineages separated by 1.7% net nucleotide distance, one exclusive to the populations from the Indian Ocean, the other exclusive to Coral Triangle populations. The latter observation might be explained by vicariance, whereby the two lineages have been isolated from one another on either side of the Sunda Shelf because of low sea level in the Pleistocene. To clarify the nomenclature of this species complex, we recommend sequencing a fragment of the cytochrome b gene of the holotypes of L. olivaceus and of its first junior synonyms L. rostratus and L. waigiensis.     

How sea level change mediates genetic divergence in coastal species across regions with varying tectonic and sediment processes

Plate tectonics and sediment processes control regional continental shelf topography. We examine the genetic consequences of how glacial‐associated sea level change interacted with variable nearshore topography since the last glaciation. We reconstructed the size and distribution of areas suitable for tidal estuary formation from the last glacial maximum, ~20 thousand years ago, to present from San Francisco, California, USA (~38°N) to Reforma, Sinaloa, Mexico (~25°N). We assessed range‐wide genetic structure and diversity of three codistributed tidal estuarine fishes (California Killifish, Shadow Goby, Longjaw Mudsucker) along ~4,600 km using mitochondrial control region and cytB sequence, and 16–20 microsatellite loci from a total of 524 individuals. Results show that glacial‐associated sea level change limited estuarine habitat to few, widely separated refugia at glacial lowstand, and present‐day genetic clades were sourced from specific refugia. Habitat increased during postglacial sea level rise and refugial populations admixed in newly formed habitats. Continental shelves with active tectonics and/or low sediment supply were steep and hosted fewer, smaller refugia with more genetically differentiated populations than on broader shelves. Approximate Bayesian computation favoured the refuge–recolonization scenarios from habitat models over isolation by distance and seaway alternatives, indicating isolation at lowstand is a major diversification mechanism among these estuarine (and perhaps other) coastal species. Because sea level change is a global phenomenon, we suggest this top‐down physical control of extirpation–isolation–recolonization may be an important driver of genetic diversification in coastal taxa inhabiting other topographically complex coasts globally during the Mid‐ to Late Pleistocene and deeper timescales.     

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (Φ_ST = 0.088) were 5.2 times more differentiated than autosomal DNA (Φ_ST = 0.017) but comparable to mtDNA (Φ_ST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (Φ_ST) and absolute (d_XY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (Φ_ST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.     

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.     

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep‐sea shark (Centroscymnus coelolepis)

Knowledge of the mechanisms limiting connectivity and gene flow in deep‐sea ecosystems is scarce, especially for deep‐sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep‐sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497‐bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ_ST = ?0.0091, F_ST = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ_ST = 0.0744, P < 0.01; F_ST = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ_ST = 0.3808, F_ST = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.     

Phylogeography of Ophiorrhiza japonica (Rubiaceae) in continental islands,the Ryukyu Archipelago,Japan

Aim Phylogeographical patterns in the Ryukyu Archipelago have been explained primarily by landbridge formation and the opening of two straits in the Pliocene, namely the Tokara and Kerama gaps. These old straits have been considered to be the barriers most likely to determine genetic boundaries. To test this, we conducted a molecular analysis of the herb Ophiorrhiza japonica. We discuss the causes of and processes involved in its phylogeographical structure and explore aspects of island separation other than the duration of the straits to explain genetic boundaries at the gaps. Location Ryukyu Archipelago, Japan. Methods Plants were collected from 40 localities in the archipelago and vicinity. Non‐coding regions of chloroplast DNA were sequenced. The genealogical relationships among haplotypes were estimated using a statistical parsimony network. To examine the phylogeographical structure, we compared two parameters of population differentiation, namely G_ST and N_ST, and conducted correlation analysis of genetic and geographical distances. Genetic boundaries were identified using Monmonier’s maximum difference algorithm. To test vicariance–dispersal hypotheses, that is, vicariance after migration via the Pliocene landbridge or over‐sea dispersal in the Pleistocene, molecular dating analysis was conducted. Results A statistical parsimony network revealed that the haplotypes from the Ryukyu Archipelago and northwards coalesce to one ancestral haplotype in Taiwan. A clear phylogeographical structure was observed: plants within the same population and populations in geographical proximity were phylogenetically close. A genetic boundary was recognized across the Kerama Gap, but not across the Tokara Gap. Dating analysis suggested that population divergence across the Kerama Gap occurred in the early to late Pleistocene. Main conclusions The statistical parsimony network suggests migration from Taiwan and northward range expansion in the archipelago. Based on the divergence time, over‐sea dispersal in the Pleistocene is likely, although migration via a Pliocene landbridge is not totally rejected. Negligible genetic differentiation across the Tokara Gap suggests recent over‐sea dispersal, possibly facilitated by the small geographical width of the gap. Conversely, the large genetic differentiation across the Kerama Gap is probably explained by the large geographical distance across it. The past splitting of a landbridge would have had a significant influence on population differentiation after a certain geographical distance was reached.