Pattern‐oriented modelling of population genetic structure

Although several statistical approaches can be used to describe patterns of genetic variation and infer stochastic differentiation, selective responses, or interruptions of gene flow due to physical or environmental barriers, it is worthwhile to note that similar processes, controlled by several parameters in theoretical models, frequently give rise to similar patterns. Here, we develop a Pattern‐Oriented Modelling (POM) approach that allows us to determine how a complex set of parameters potentially driving empirical genetic differentiation among populations generate alternative scenarios that can be fitted to observed data. We generated 10 000 random combinations of parameters related to population size, gene flow and response to gradients (both driven by dispersal and selection) in a spatially explicit model, and analysed simulated patterns with F_ST statistics and mean correlograms using Moran's I spatial autocorrelation coefficients. These statistics were compared with observed patterns for a tree species endemic to the Brazilian Cerrado. For a best match with observed F_ST (equal to 0.182), the important parameters driving simulated scenario are mainly related to population structure, including low population size with closed populations (low N_m), strong distance decay of gene flow, in addition to a strong effect of the initial variance of allele frequencies. These scenarios present a low autocorrelation of allele frequencies. Best matching of correlograms, on the other hand, appears in simulations with a large population size, high N_m and low population differentiation and F_ST (as well as more gene flow). Thus, targeting the two statistics (correlograms and F_ST) shows that best matches with empirical data with two distinct sets of parameters in the simulations, because observed patterns involve both a relatively high F_ST and significant autocorrelation. This conflict can be resolved by assuming that initial variance in allele frequencies can be interpreted as reflecting deep‐time historical variation and evolutionary dynamics of allele frequencies, creating a relatively high level of population differentiation, whereas current patterns in gene flow creates spatial autocorrelation. This make sense in terms of the previous knowledge on population differentiation in D. alata, especially if patterns are explained by a combination of isolation‐by‐distance and allelic surfing due to range expansion after the last glacial maximum. This reveals the potential for more complex applications of POM in population genetics. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 1152–1161.     

Intercontinental genetic structure in the amphi‐Pacific peatmoss Sphagnum miyabeanum (Bryophyta: Sphagnaceae)

Unlike seed plants where global biogeographical patterns typically involve interspecific phylogenetic history, spore‐producing bryophyte species often have intercontinental distributions that are best understood from a population genetic perspective. We sought to understand how reproductive processes, especially dispersal, have contributed to the intercontinental ‘Pacific Rim’ distribution of Sphagnum miyabeanum. In total, 295 gametophyte plants from western North America (California, Oregon, British Columbia, Alaska), Russia, Japan, and China were genotyped at 12 microsatellite loci. Nucleotide sequences were obtained for seven anonymous nuclear loci plus two plastid regions from 21 plants of S. miyabeanum and two outgroup species. We detected weak but significant genetic differentiation among plants from China, Japan, Alaska, British Columbia, and the western USA. Alaskan plants are genetically most similar to Asian plants, and British Columbian plants are most similar to those in the western USA. There is detectable migration between regions, with especially high levels between Alaska and Asia (China and Japan). Migration appears to be recent and/or ongoing, and more or less equivalent in both directions. There is weak (but significant) isolation‐by‐distance within geographical regions, and the slope of the regression of genetic on geographical distance differs for Asian versus North American plants. A distinctive Vancouver Island morphotype is very weakly differentiated, and does not appear to be reproductively isolated from plants of the normal morphotype. The intercontinental geographical range of S. miyabeanum reflects recent and probably ongoing migration, facilitated by the production of tiny spores capable of effective long distance dispersal. The results of the present study are consistent with Pleistocene survival of S. miyabeanum in unglaciated Beringia, although we cannot eliminate the possibility that the species recolonized Alaska from Asia more recently. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 17–37.     

Life histories and ecotype conservation in an adaptive vertebrate: Genetic constitution of piscivorous brown trout covaries with habitat stability

Ecotype variation in species exhibiting different life history strategies may reflect heritable adaptations to optimize reproductive success, and potential for speciation. Traditionally, ecotypes have, however, been defined by morphometrics and life history characteristics, which may be confounded with individual plasticity. Here, we use the widely distributed and polytypic freshwater fish species brown trout (Salmo trutta) as a model to study piscivorous life history and its genetic characteristics in environmentally contrasting habitats; a large lake ecosystem with one major large and stable tributary, and several small tributaries. Data from 550 fish and 13 polymorphic microsatellites (H_e = 0.67) indicated ecotype‐specific genetic differentiation (θ = 0.0170, p < .0001) among Bayesian assigned small riverine resident and large, lake migrating brown trout (>35 cm), but only in the large tributary. In contrast, large trout did not constitute a distinct genetic group in small tributaries, or across riverine sites. Whereas life history data suggest a small, river resident and a large migratory piscivorous ecotype in all studied tributaries, genetic data indicated that a genetically distinct piscivorous ecotype is more likely to evolve in the large and relatively more stable river habitat. In the smaller tributaries, ecotypes apparently resulted from individual plasticity. Whether different life histories and ecotypes result from individual plasticity or define different genetic types, have important consequence for conservation strategies.     

Fine‐scale patterns of genetic divergence within and between morphologically variable subspecies of the sea urchin Heliocidaris erythrogramma (Echinometridae)

The spatial scale over which genetic divergences occur between populations and the extent that they are paralleled by morphological differences can vary greatly among marine species. In the present study, we use a hierarchical spatial design to investigate genetic structure in Heliocidaris erythrogramma occurring on near shore limestone reefs in Western Australia. These reefs are inhabited by two distinct subspecies: the thick‐spined Heliocidaris erythrogramma armigera and the thin‐spined Heliocidaris erythrogramma erythrogramma, each of which also have distinct colour patterns. In addition to pronounced morphological variation, H. erythrogramma exhibits a relatively short (3–4 days) planktonic phase before settlement and metamorphosis, which limits their capacity for dispersal. We used microsatellite markers to determine whether patterns of genetic structure were influenced more by morphological or life history limitations to dispersal. Both individual and population‐level analyses found significant genetic differentiation between subspecies, which was independent of geographical distance. Genetic diversity was considerably lower within H. e. erythrogramma than within H. e. armigera and genetic divergence was four‐fold greater between subspecies than among populations within subspecies. This pattern was consistent even at fine spatial scales (< 5 km). We did detect some evidence of gene flow between the subspecies; however, it appears to be highly restricted. Within subspecies, genetic structure was more clearly driven by dispersal capacity, although weak patterns of isolation‐by‐distance suggest that there may be other factors limiting gene exchange between populations. Our results show that spatial patterns of genetic structure in Western Australian H. erythrogramma is influenced by a range of factors but is primarily correlated with the distribution of morphologically distinct subspecies. This suggests the presence of reproductive barriers to gene exchange between them and demonstrates that morphological variation can be a good predictor of genetic divergence. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 578–592.     

Genetic population structure of spotted seatrout Cynoscion nebulosus along the south‐eastern U.S.A.

Analyses of the genetic population structure of spotted seatrout Cynoscion nebulosus along the south‐eastern U.S. coast using 13 microsatellites suggest significant population differentiation between fish in North Carolina (NC) compared with South Carolina (SC) and Georgia (GA), with New River, NC, serving as an area of integration between northern and southern C. nebulosus. Although there is a significant break in gene flow between these areas, the overall pattern throughout the sampling range represents a gradient in genetic diversification with the degree of geographic separation. Latitudinal distance and estuarine density appear to be main drivers in the genetic differentiation of C. nebulosus along the south‐eastern U.S. coast. The isolation‐by‐distance gene‐flow pattern creates fine‐scale differences in the genetic composition of proximal estuaries and dictates that stocking must be confined to within 100 km of the location of broodstock collection in order to maintain the natural gradient of genetic variation along the south‐eastern U.S. coast.     

Phylogeography of the humbug damselfish,Dascyllus aruanus (Linnaeus, 1758): evidence of Indo‐Pacific vicariance and genetic differentiation of peripheral populations

The phylogeographical structure of coral‐associated reef fishes may have been severely affected, more than species from deeper habitats, by habitat loss during periods of low sea level. The humbug damselfish, Dascyllus aruanus, is widely distributed across the Indo‐West Pacific, and exclusively inhabits branching corals. We used mitochondrial cytochrome b sequence and seven microsatellite loci on D. aruanus samples (260 individuals) from 13 locations across the Indo‐West Pacific to investigate its phylogeographical structure distribution‐wide. A major genetic partition was found between the Indian and Pacific Ocean populations, which we interpret as the result of geographical isolation on either side of the Indo‐Pacific barrier during glacial periods. The peripheral populations of the Red Sea and the Society Islands exhibited lower genetic diversity, and substantial genetic differences with the other populations, suggesting relative isolation. Thus, vicariance on either side of the Indo‐Pacific barrier and peripheral differentiation are considered to be the main drivers that have shaped the phylogeographical patterns presently observed in D. aruanus. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 931–942.     

Population genomic variation reveals roles of history,adaptation and ploidy in switchgrass

Geographic patterns of genetic variation are shaped by multiple evolutionary processes, including genetic drift, migration and natural selection. Switchgrass (Panicum virgatum L.) has strong genetic and adaptive differentiation despite life history characteristics that promote high levels of gene flow and can homogenize intraspecific differences, such as wind‐pollination and self‐incompatibility. To better understand how historical and contemporary factors shape variation in switchgrass, we use genotyping‐by‐sequencing to characterize switchgrass from across its range at 98 042 SNPs. Population structuring reflects biogeographic and ploidy differences within and between switchgrass ecotypes and indicates that biogeographic history, ploidy incompatibilities and differential adaptation each have important roles in shaping ecotypic differentiation in switchgrass. At one extreme, we determine that two Panicum taxa are not separate species but are actually conspecific, ecologically divergent types of switchgrass adapted to the extreme conditions of coastal sand dune habitats. Conversely, we identify natural hybrids among lowland and upland ecotypes and visualize their genome‐wide patterns of admixture. Furthermore, we determine that genetic differentiation between primarily tetraploid and octoploid lineages is not caused solely by ploidy differences. Rather, genetic diversity in primarily octoploid lineages is consistent with a history of admixture. This suggests that polyploidy in switchgrass is promoted by admixture of diverged lineages, which may be important for maintaining genetic differentiation between switchgrass ecotypes where they are sympatric. These results provide new insights into the mechanisms shaping variation in widespread species and provide a foundation for dissecting the genetic basis of adaptation in switchgrass.     

Selection on hybrids of ecologically divergent ecotypes of a marine snail: the relative importance of exogenous and endogenous barriers

Unravelling the form of selection acting on hybrids of ecotypes undergoing ecological speciation is essential to understand the mechanisms behind the evolution of reproductive isolation in the face of gene flow. Shell phenotype is known to be affected by natural selection and is involved in the fitness of the marine snail Littorina saxatilis. Here, we studied the association between shell traits and fitness in hybrids in order to determine the relative role of exogenous and endogenous selection in this hybrid zone of L. saxatilis. We show that directional selection is the predominant mode of selection among hybrids. We also show its heterogeneity, affecting different shell traits, within populations at the level of the microhabitat. Therefore, endogenous selection mechanisms are most probably lacking in this hybrid zone and exogenous barriers (pre‐ and post‐zygotic) are possibly one of the main forces behind the evolution of barriers to gene flow between these ecologically divergent ecotypes. This study shows how this barrier might represent an important type of reproductive isolation within ecological speciation, and this should be taken into account in future studies of speciation in hybrid zones. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 391–400.     

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater‐resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size‐assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.     

Comparative phylogeography of Australo‐Papuan mangrove‐restricted and mangrove‐associated avifaunas

The world's richest mangrove‐restricted avifauna is in Australia and New Guinea. The history of differentiation of the species involved and their patterns of intraspecific genetic variation remain poorly known. Here, we use sequence data derived from two mitochondrial protein‐coding genes to study the evolutionary history of eight co‐distributed mangrove‐restricted and mangrove‐associated birds from the Australian part of this region. Utilizing a comparative phylogeographical framework, we observed that the study species present concordantly located phylogeographical breaks across their shared geographical distribution, a plausible signature of common mechanisms of vicariance underlying this pattern. Barriers such as the Canning Gap, Bonaparte Gap, and the Carpentarian Gaps all had important but varying degrees of impact on the studied species. The Burdekin Gap along Australia's eastern seaboard probably had only a minor influence as a barrier to gene flow in mangrove birds. Statistical phylogeographical simulations were able to discriminate among alternative scenarios involving six different geographical and temporal population separations. Species exhibiting recent colonizations into mangroves include Rhipidura phasiana, Myiagra ruficollis, and Myzomela erythrocephala. By contrast, Peneoenanthe pulverulenta, Pachycephala melanura, Pachycephala lanioides, Zosterops luteus, and Colluricincla megarhyncha all had deeper histories, reflected as more marked phylogeographical divisions separating populations on the eastern seaboard/Cape York Peninsula from more western regions such as the Arnhem Land, the Pilbara, and the Kimberley. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 574–598.     

Sun skink landscape genomics: assessing the roles of micro‐evolutionary processes in shaping genetic and phenotypic diversity across a heterogeneous and fragmented landscape

Incorporating genomic data sets into landscape genetic analyses allows for powerful insights into population genetics, explicitly geographical correlates of selection, and morphological diversification of organisms across the geographical template. Here, we utilize an integrative approach to examine gene flow and detect selection, and we relate these processes to genetic and phenotypic population differentiation across South‐East Asia in the common sun skink, Eutropis multifasciata. We quantify the relative effects of geographic and ecological isolation in this system and find elevated genetic differentiation between populations from island archipelagos compared to those on the adjacent South‐East Asian continent, which is consistent with expectations concerning landscape fragmentation in island archipelagos. We also identify a pattern of isolation by distance, but find no substantial effect of ecological/environmental variables on genetic differentiation. To assess whether morphological conservatism in skinks may result from stabilizing selection on morphological traits, we perform F_ST–P_ST comparisons, but observe that results are highly dependent on the method of comparison. Taken together, this work provides novel insights into the manner by which micro‐evolutionary processes may impact macro‐evolutionary scale biodiversity patterns across diverse landscapes, and provide genomewide confirmation of classic predictions from biogeographical and landscape ecological theory.     

A cryptic genetic boundary in remnant populations of a long‐lived,bird‐pollinated shrub Banksia sphaerocarpa var. caesia (Proteaceae)

Plant species often exhibit genetic structure at multiple spatial scales. Detection of this structure may depend on the sampling strategy used. We intensively sampled a common, naturally patchy Banksia species within a 200 km² region, in order to assess patterns of genetic diversity and structure at multiple spatial scales. In total, 1321 adult shrubs from 37 geographical populations were genotyped using eight highly polymorphic microsatellite markers developed for the species. Genetic structure was detected at three spatial scales. First, we identified a stark and unexpected division of the landscape into two genetic subregions, one to the north‐east and one to the south‐west of the sampling grid. This differentiation was based on sudden, highly structured differences in common allele frequencies, the cause of which is unknown but that may relate to physical and reproductive barriers to gene flow, localised selection, and/or historical processes. Second, we observed genetic differentiation of populations within these subregions, reflecting previously described patterns of restricted pollen flow in this species. Finally, fine‐scale genetic structure, although weak, was observed within some of the populations (mean S_P = 0.01837). When feasible, intensive sampling may uncover cryptic patterns of genetic structure that would otherwise be overlooked when sampling at broader spatial scales. Further studies using a similar sampling strategy may reveal this type of discontinuity to be a feature of other south‐western Australian taxa and has implications for our understanding of evolution in this landscape as well as conservation into the future. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 115 , 241–255.     

Evolution of life‐history traits in geographically isolated populations of Vaejovis scorpions (Scorpiones: Vaejovidae)

Geographical isolation can over time accumulate life‐history variation which can eventually lead to speciation. We used five species of Vaejovis scorpions that have been isolated from one another since the Pleistocene glaciation to identify if biogeographical patterns have allowed for the accumulation of life‐history variation among species. Gravid females were captured and brought back to the lab until giving birth. Once offspring had begun to disperse, measurements of female size, reproductive investment, offspring size, offspring number, and variation in offspring size were recorded. Differences in how each species allocated energy to these variables were analysed utilizing path analysis and structural equation modelling. Female and offspring size, litter size, and total litter mass differed among species, but relative energetic investment did not. Most significant differences among species were not present after removing the effect of female size, indicating that female size is a major source of life‐history variation. Path analyses indicated that there was no size–number trade‐off within any species and that each species allocates energy toward total litter mass differently. Additionally, as offspring size increased, the variation in offspring mass decreased. These results show that each species allocates the same relative amount of energy in different ways. The variation seen could be a response to environmental variability or uncertainty, a product of maternal effects, or caused by the sufficient accumulation of genetic differences due to geographical isolation. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 715–727.     

THE RELEVANCE OF GENE FLOW IN METAPOPULATION DYNAMICS OF AN OCEANIC ISLAND ENDEMIC,OLEA EUROPAEA SUBSP. GUANCHICA

Theoretical and empirical studies suggest that geographical isolation and extinction–recolonization dynamics are two factors causing strong genetic structure in metapopulations, but their consequences in species with high dispersal abilities have not been tested at large scales. Here, we investigated the effect of population age structure and isolation by distance in the patterns of genetic diversity in a wind‐pollinated, zoochorous tree (Olea europaea subsp. guanchica) sporadically affected by volcanic events across the Canarian archipelago. Genetic variation was assessed at six nuclear microsatellites (nDNA) and six chloroplast fragments (cpDNA) in nine subpopulations sampled on four oceanic islands. Subpopulations occurring on more recent substrates were more differentiated than those on older substrates, but within‐subpopulation genetic diversity was not significantly different between age groups for any type of marker. Isolation‐by‐distance differentiation was observed for nDNA but not for cpDNA, in agreement with other metapopulation studies. Contrary to the general trend for island systems, between‐island differentiation was extremely low, and lower than differentiation between subpopulations on the same island. The pollen‐to‐seed ratio was close to one, two orders of magnitude lower than the average estimated for other wind‐pollinated, animal‐dispersed plants. Our results showed that population turnover and geographical isolation increased genetic differentiation relative to an island model at equilibrium, but overall genetic structure was unexpectedly weak for a species distributed among islands. This empirical study shows that extensive gene flow, particularly mediated by seeds, can ameliorate population subdivision resulting from extinction–recolonization dynamics and isolation by distance.     

Weak large‐scale population genetic structure in a philopatric seabird,the European Shag Phalacrocorax aristotelis

Quantifying population genetic structure is fundamental to testing hypotheses regarding gene flow, population divergence and dynamics across large spatial scales. In species with highly mobile life‐history stages, where it is unclear whether such movements translate into effective dispersal among discrete philopatric breeding populations, this approach can be particularly effective. We used seven nuclear microsatellite loci and mitochondrial DNA (ND2) markers to quantify population genetic structure and variation across 20 populations (447 individuals) of one such species, the European Shag, spanning a large geographical range. Despite high breeding philopatry, rare cross‐sea movements and recognized subspecies, population genetic structure was weak across both microsatellites and mitochondrial markers. Furthermore, although isolation‐by‐distance was detected, microsatellite variation provided no evidence that open sea formed a complete barrier to effective dispersal. These data suggest that occasional long‐distance, cross‐sea movements translate into gene flow across a large spatial scale. Historical factors may also have shaped contemporary genetic structure: cluster analyses of microsatellite data identified three groups, comprising colonies at southern, mid‐ and northern latitudes, and similar structure was observed at mitochondrial loci. Only one private mitochondrial haplotype was found among subspecies, suggesting that this current taxonomic subdivision may not be mirrored by genetic isolation.     

Genetic differentiation,speciation, and phylogeography of cactus flies (Diptera: Neriidae: Odontoloxozus) from Mexico and south‐western USA

Nucleotide sequences from the mitochondrial cytochrome c oxidase subunit I (COI) gene, comprising the standard barcode segment, were used to examine genetic differentiation, systematics, and population structure of cactus flies (Diptera: Neriidae: Odontoloxozus) from Mexico and south‐western USA. Phylogenetic analyses revealed that samples of Odontoloxozus partitioned into two distinct clusters: one comprising the widely distributed Odontoloxozus longicornis (Coquillett) and the other comprising Odontoloxozus pachycericola Mangan & Baldwin, a recently described species from the Cape Region of the Baja California peninsula, which we show is distributed northward to southern California, USA. A mean Kimura two‐parameter genetic distance of 2.8% between O. longicornis and O. pachycericola, and eight diagnostic nucleotide substitutions in the COI gene segment, are consistent with a species‐level separation, thus providing the first independent molecular support for recognizing O. pachycericola as a distinct species. We also show that the only external morphological character considered to separate adults of the two species (number of anepisternal bristles) varies with body size and is therefore uninformative for making species assignments. Analysis of molecular variance indicated significant structure among populations of O. longicornis from three main geographical areas, (1) Arizona, USA and Sonora, Mexico; (2) Santa Catalina Island, California, USA; and (3) central Mexico (Querétaro and Guanajuato), although widely‐separated populations from Arizona and Sonora showed no evidence of structure. A TCS haplotype network showed no shared haplotypes of O. longicornis among the three main regions. The potential roles of vicariance and isolation‐by‐distance in restricting gene flow and promoting genetic differentiation and speciation in Odontoloxozus are discussed. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 245–256.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

Small fish,large river: Surprisingly minimal genetic structure in a dispersal‐limited,habitat specialist fish

Genetic connectivity is expected to be lower in species with limited dispersal ability and a high degree of habitat specialization (intrinsic factors). Also, gene flow is predicted to be limited by habitat conditions such as physical barriers and geographic distance (extrinsic factors). We investigated the effects of distance, intervening pools, and rapids on gene flow in a species, the Tuxedo Darter (Etheostoma lemniscatum), a habitat specialist that is presumed to be dispersal‐limited. We predicted that the interplay between these intrinsic and extrinsic factors would limit dispersal and lead to genetic structure even at the small spatial scale of the species range (a 38.6 km river reach). The simple linear distribution of E. lemniscatum allowed for an ideal test of how these factors acted on gene flow and allowed us to test expectations (e.g., isolation‐by‐distance) of linearly distributed species. Using 20 microsatellites from 163 individuals collected from 18 habitat patches, we observed low levels of genetic structure that were related to geographic distance and rapids, though these factors were not barriers to gene flow. Pools separating habitat patches did not contribute to any observed genetic structure. Overall, E. lemniscatum maintains gene flow across its range and is comprised of a single population. Due to the linear distribution of the species, a stepping‐stone model of dispersal best explains the maintenance of gene flow across its small range. In general, our observation of higher‐than‐expected connectivity likely stems from an adaptation to disperse due to temporally unstable and patchy habitat.     

A multilocus sequencing approach reveals the cryptic phylogeographical history of Phyllodoce nipponica Makino (Ericaceae)

Discordant phylogeographical patterns among species with similar distributions may not only denote specific biogeographical histories of different species, but also could represent stochastic variance of genealogies in applied genetic markers. A multilocus investigation representing different genomes can be used to address the latter concern, allowing robust inference to biogeographical history. In the present study, we conducted a multilocus phylogeographical analysis to re‐examine the genetic structuring of Phyllodoce nipponica, in which chloroplast (cp)DNA markers exhibited a discordant pattern compared to those of other alpine plants. The geographical structure of sequence variation at five nuclear loci was not consistent with that of cpDNA and showed differentiation between the northern and southern parts of the range of this species. Its demographic history inferred from the isolation‐with‐migration model suggests that the north–south divergence originated from Pleistocene vicariance. In addition, the demographic parameters showed a lack of chloroplast‐specific gene flow, suggesting that stochastic variance in genealogy resulted in the discordant geographical structure. Thus, P. nipponica probably experienced Pleistocene vicariance between its southern and northern range parts in concordance with other alpine plants in the Japanese archipelago. The findings of the present study demonstrates the importance of using a multilocus approach for inferring population dynamics, as well as for reconciling discordant phylogeographical patterns among species. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 214–226.     

Population differentiation and speciation in the genus Characidium (Characiformes: Crenuchidae): effects of reproductive and chromosomal barriers

Both time and low gene flow are the key factors by which different biological species arise. The divergence process among lineages and the development of pre‐ or postzygotic isolation occur when gene flow events are lacking. The separation among species of the genus Characidium was analysed in relation to the geomorphological mechanisms in river courses, events of captured adjacent upland drainages in south‐eastern Brazil, and sex chromosome differences. The ZZ/ZW sex chromosomes of Characidium vary in size, morphology, degree of heterochromatinization, and presence/absence of ribosomal DNA. The goal of this study was to understand the mechanism of sex chromosome differentiation, its close association with the geological history of cladogenetic events among drainages, and reproductive isolation leading to Characidium speciation. The W‐specific probe from Characidium gomesi generated a highlighted signal on the entire W chromosome of C. gomesi, Characidium heirmostigmata, Characidium pterostictum, and Characidium sp., instead of karyotypes of three Characidium aff. zebra populations, which showed scattered signals. An evolutionary and biogeographic landscape arose by analysis of ribosomal DNA site location and differentiation of the sex chromosomes, which established mechanisms of reproductive isolation leading to meiotic barriers, keeping the biological unit distinct even if the contact among species was restored. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 541–553.