Divergence,gene flow,and speciation in eight lineages of trans‐Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans‐Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split‐migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.     

Illumination of cryptic species boundaries in long-tailed shrew tenrecs (Mammalia: Tenrecidae; Microgale), with new insights into geographic variation and distributional constraints

The increasing use of mitochondrial DNA (mtDNA) to explore and test species limits among morphologically similar species is potentially compromised by phenomena poorly reflective of organismal history and speciation, including (but not limited to) stochastic lineage sorting and gene flow. In situations where molecular data are only available from a single gene or linkage partition (e.g. mtDNA), corroboration of suspected species boundaries should be sought from independent lines of evidence, such as morphology. Recent attempts to delimit species using mtDNA and morphology have either implicitly or explicitly ignored the possibility that distinct species can occur in direct sympatry throughout much of their range, presumably because such situations are believed to be rare. We examined phylogenetic relationships within the long‐tailed shrew tenrecs (Mammalia: Tenrecidae; Microgale spp.) from Madagascar. Current taxonomy recognizes two broadly sympatric species, though as many as six have been described. Given that alpha taxonomy within shrew tenrecs has been controversial, and that patterns of morphological variation can be especially difficult to assess for this group, some authors have suggested that additional cryptic species may exist. To examine this possibility, we conducted a phylogenetic study using the mitochondrial NADH dehydrogenase subunit 2 gene and a morphometric analysis of 29 craniodental, postcranial, and external measurements from a broad geographical sample of long‐tailed shrew tenrecs. The two data sets were nearly perfectly congruent in identifying four groups that can be classified as species, thereby doubling the currently recognized number of species. We present previously unrecognized distributional evidence consistent with our conclusions and provide an empirical example of how a revised understanding of species limits alters inferences of geographic variation and species coexistence, particularly with respect to fine‐scale habitat partitioning. The results of this study suggest that certain species pairs, previously assumed to be single species occupying broad elevational ranges, are actually reproductively isolated units that are partitioning their environment along elevational lines. © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 83 , 1–22.     

Speciation with gene flow: Evidence from a complex of alpine butterflies (Coenonympha,Satyridae)

Until complete reproductive isolation is achieved, the extent of differentiation between two diverging lineages is the result of a dynamic equilibrium between genetic isolation and mixing. This is especially true for hybrid taxa, for which the degree of isolation in regard to their parental species is decisive in their capacity to rise as a new and stable entity. In this work, we explored the past and current patterns of hybridization and divergence within a complex of closely related butterflies in the genus Coenonympha in which two alpine species, C. darwiniana and C. macromma, have been shown to result from hybridization between the also alpine C. gardetta and the lowland C. arcania. By testing alternative scenarios of divergence among species, we show that gene flow has been uninterrupted throughout the speciation process, although leading to different degrees of current genetic isolation between species in contact zones depending on the pair considered. Nonetheless, at broader geographic scale, analyses reveal a clear genetic differentiation between hybrid lineages and their parental species, pointing out to an advanced stage of the hybrid speciation process. Finally, the positive correlation observed between ecological divergence and genetic isolation among these butterflies suggests a potential role for ecological drivers during their speciation processes.     

Recent divergence with gene flow in Tennessee cave salamanders (Plethodontidae: Gyrinophilus) inferred from gene genealogies

Cave organisms occupy a special place in evolutionary biology because convergent morphologies of many species demonstrate repeatability in evolution even as they obscure phylogenetic relationships. The origin of specialized cave-dwelling species also raises the issue of the relative importance of isolation vs. natural selection in speciation. Two alternative hypotheses describe the origin of subterranean species. The 'climate-relict' model proposes allopatric speciation after populations of cold-adapted species become stranded in caves due to climate change. The 'adaptive-shift' model proposes parapatric speciation driven by divergent selection between subterranean and surface habitats. Our study of the Tennessee cave salamander complex shows that the three nominal forms (Gyrinophilus palleucus palleucus, G. p. necturoides, and G. gulolineatus) arose recently and are genealogically nested within the epigean (surface-dwelling) species, G. porphyriticus. Short branch lengths and discordant gene trees were consistent with a complex history involving gene flow between diverging forms. Results of coalescent-based analysis of the distribution of haplotypes among groups reject the allopatric speciation model and support continuous or recurrent genetic exchange during divergence. These results strongly favour the hypothesis that Tennessee cave salamanders originated from spring salamanders via divergence with gene flow.     

Phylogeographic structure, gene flow and species status in blue grouse (Dendragapus obscurus)

We investigated the genetic population structure and species status of a relatively sedentary bird that is a permanent resident of western North American forests, the blue grouse (Dendragapus obscurus). Phylogenetic analysis of complete mitochondrial control region DNA sequences resulted in the identification of three basal clades of haplotypes that were largely congruent with well-known biogeographical regions. These clades corresponded to the parapatric sooty (D. o. fuliginosus) and dusky (D. o. obscurus) subspecies groups of blue grouse plus a previously unrecognized division between northern and southern dusky grouse populations; the latter does not correspond closely to any currently recognized subspecies boundary. Approximately 66% of the total genetic variation was distributed among these three regions. Maximum likelihood estimates of gene flow between the regions were low or asymmetric; gene flow has been insufficient to prevent genetic divergence between dusky and sooty grouse. Estimates of gene flow among populations within sooty grouse were large except across the Columbia River valley. Among populations of dusky grouse, estimates of gene flow were heterogeneous and asymmetrical, reflecting large-scale fragmentation of the distribution due to landscape features and associated vegetation. Genetic, morphological and behavioural evidence suggest that sooty and dusky grouse are species-level taxa; the specific status of a third clade remains ambiguous.     

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.     

Ecological divergence in the presence of gene flow in two closely related Oryza species (Oryza rufipogon and O. nivara)

Ecological divergence plays a prominent role in the process of speciation, but how divergence occurs in the face of gene flow is still less clear, and remains controversial among evolutionists. Here we investigated the nucleotide diversity, divergence and gene flow between Oryza nivara and O. rufipogon using sequences of seven chloroplast and nuclear loci. By analysing samples from 26 wild populations across the geographic ranges of the two species, we showed that both species were highly structured and O. rufipogon maintained a higher level of species‐wide diversity than O. nivara. Notably, phylogenetic, amova and F_ST analyses were unable to detect significant nucleotide differentiation between the two species. We estimated that the two species began to diverge at c. 0.16 million years ago. Our coalescent‐based simulations strongly rejected the simple isolation model of zero migration between species, but rather provided unambiguous evidence of bidirectional gene flow between species, particularly from O. rufipogon to O. nivara. Our simulations also indicated that gene flow was recurrent during the divergence process rather than arising from secondary contact after allopatric divergence. In conjunction with different morphological and life‐history traits and habitat preference in the two species, this study supports the hypothesis that these Oryza species are better treated as ecotypes that diverged quite recently and are still under the process of divergence. Importantly, we demonstrate the ecological divergence between O. rufipogon and O. nivara in the presence of significant gene flow, implying that natural selection plays a primary role in driving the divergence of the two Oryza species.     

Genomic footprints of speciation in Atlantic eels (Anguilla anguilla and A. rostrata)

The importance of speciation‐with‐geneflow scenarios is increasingly appreciated. However, the specific processes and the resulting genomic footprints of selection are subject to much discussion. We studied the genomics of speciation between the two panmictic, sympatrically spawning sister species; European (Anguilla anguilla) and American eel (A. rostrata). Divergence is assumed to have initiated more than 3 Ma, and although low gene flow still occurs, strong postzygotic barriers are present. Restriction‐site‐associated DNA (RAD) sequencing identified 328 300 SNPs for subsequent analysis. However, despite the presence of 3757 strongly differentiated SNPs (F_ST > 0.8), sliding window analyses of F_ST showed no larger genomic regions (i.e. hundreds of thousands to millions of bases) of elevated differentiation. Overall F_ST was 0.041, and linkage disequilibrium was virtually absent for SNPs separated by more than 1000 bp. We suggest this to reflect a case of genomic hitchhiking, where multiple regions are under directional selection between the species. However, low but biologically significant gene flow and high effective population sizes leading to very low genetic drift preclude accumulation of strong background differentiation. Genes containing candidate SNPs for positive selection showed significant enrichment for gene ontology (GO) terms relating to developmental processes and phosphorylation, which seems consistent with assumptions that differences in larval phase duration and migratory distances underlie speciation. Most SNPs under putative selection were found outside coding regions, lending support to emerging views that noncoding regions may be more functionally important than previously assumed. In total, the results demonstrate the necessity of interpreting genomic footprints of selection in the context of demographic parameters and life‐history features of the studied species.     

High gene flow in two thrips‐pollinated South‐East Asian pioneer trees: genetic diversity and population structure of Macaranga hypoleuca and M. beccariana (Euphorbiaceae)

As a result of intensive exploitation, disturbed forests now dominate large areas of lowland tropical rainforest in South‐East Asia. The genus Macaranga comprises some of the most important pioneer tree species of the region, among them M. beccariana and M. hypoleuca, two closely related obligate ant‐plants pollinated by thrips. We used nuclear and plastid DNA markers to address questions of genetic diversity and population structure. Twelve plastid haplotypes were detected among 281 samples, three of which were shared between the two study species. Hybrids between the two species appear to be rare. Overall, genetic diversity in both species was moderate to high, with low levels of population differentiation, consistent with other tropical pioneer trees. Genetic structure was generally more pronounced in plastid than in nuclear data, indicating that gene flow via pollen may be more efficient than via seeds. Thrips apparently also serve as efficient pollinators over long distances, perhaps through a combination of passive dispersal by wind and active search for inflorescences in the target area. Our results indicate that M. beccariana and M. hypoleuca populations from recently disturbed habitats do not yet suffer from reduced genetic diversity or increased inbreeding. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 606–621.     

Geographic patterns of population genetic structure in Mytilocypris (Ostracoda: Cyprididae): interpreting breeding systems,gene flow and history in species with differing distributions

Samples from 83 populations of salt lake Ostracods belonging to the genus Mytilocypris were collected from 74 saline lakes and ponds in the semi-arid regions of Australia. These populations were examined for variation at six polymorphic enzyme loci to diagnose breeding systems and to measure population structure, to investigate relative levels of gene flow in species with differing distributions and hence different presumed dispersal capabilities. Despite the occurrence of some populations in disjunct, peripheral, and recent ephemeral habitats, all populations of each species were found to reproduce sexually. Gene flow does occur on a local basis and appears to be facilitated by occupation in the same drainage basin for some species. There was considerable gene pool fragmentation among peripheral populations of four of the five species. Only one species, M. mytiloides, was relatively homogeneous across its range. It may be that gene flow is non-existent into peripheral populations because of poor dispersal abilities, or it may not be frequent enough to overcome local selective pressures. Regardless of these possibilities, the observed gene pool fragmentation has implications for allopatric speciation.     

Oligoryzomys flavescens (Rodentia, Muridae): gene flow among populations from central-eastern Argentina

In species acting as hosts of infectious agents, the extent of gene flow between populations is of particular interest because the expansion of different infectious diseases is usually related to the dispersal of the host. We have estimated levels of gene flow among populations of the sigmodontine rodent Oligoryzomys flavescens, in which high titers of antibodies have been detected for a Hantavirus in Argentina that produces a severe pulmonary syndrome. Enzyme polymorphism was studied by means of starch gel electrophoresis in 10 populations from the area where human cases of Hantavirus have occurred. Genetic differentiation between populations was calculated from F_ST values with the equation Nm = [(1/F_ST−1]/4. To assess the relative importance of current gene flow and historical associations between populations, the relationship of population pairwise log Nm and log geographic distance was examined. Low F_ST (mean = 0.038) and high Nm (15.27) values suggest high levels of gene flow among populations. The lack of an isolation by distance pattern would indicate that this species has recently colonized the area. The northernmost population, located on the margin of a great river, shows very high levels of gene flow with the downstream populations despite the large geographic distances. Passive transport of animals down the river by floating plants would promote unidirectional gene flow. This fact and the highest mean heterozygosity of that northernmost population suggest it is a center of dispersal within the species' range. This revised version was published online in July 2006 with corrections to the Cover Date.     

Morphological and genetic differentiation in anadromous smelt Osmerus mordax (Mitchill): disentangling the effects of geography and morphology on gene flow

Morphological analyses were combined with genetic analyses at nine microsatellite loci to examine the determinants of gene flow at 21 spawning locations of rainbow smelt Osmerus mordax along the east coast of Canada. Associations between morphology, geography and gene flow were examined using a computational geometric approach and partial Mantel tests. Significant barriers to gene flow and discontinuities in morphology were observed between Newfoundland and mainland Canada, as well as within Newfoundland samples. On regional scales, contrasting patterns were present with restricted gene flow between Newfoundland populations ( F _ST= c . 0·11) and high gene flow between mainland populations ( F _ST= c . 0·017). Within Newfoundland populations, geographic distance was significantly associated with gene flow ( r = 0·85, P < 0·001) contrasting mainland samples where gene flow was most associated with phenotypic divergence ( r = 0·33, P < 0·001). At large spatial scales, weak ( r = 0·19, P = 0·02) associations between gene flow and geographic distance were observed, and moderate associations were also observed between gene flow and morphology ( r = 0·28, P < 0·001). The presence of significant genetic isolation by distance in Newfoundland samples and the clear discontinuity associated with the Cabot Strait suggest geography may be the primary determinant of gene flow. Interestingly, the association between genetic and morphological divergence within mainland samples and overall, supports the hypothesis that gene flow may be moderated by morphological divergence at larger spatial scales even in high gene flow environments.     

Reproductive division of labour and thelytoky result in sympatric barriers to gene flow in honeybees (Apis mellifera L.)

Determining the extent and causes of barriers to gene flow is essential for understanding sympatric speciation, but the practical difficulties of quantifying reproductive isolation remain an obstacle to analysing this process. Social parasites are common in eusocial insects and tend to be close phylogenetic relatives of their hosts (= Emery's rule). Sympatric speciation caused by reproductive isolation between host and parasite is a possible evolutionary pathway. Socially parasitic workers of the Cape honeybee, Apis mellifera capensis, produce female clonal offspring parthenogenetically and invade colonies of the neighbouring subspecies A. m. scutellata. In the host colony, socially parasitic workers can become pseudoqueens, an intermediate caste with queenlike pheromone secretion. Here, we show that over an area of approximately 275.000 km2, all parasitic workers bear the genetic signature of a clone founded by a single ancestral worker genotype. Any gene flow from the host to the parasite is impossible because honeybee workers cannot mate. Gene flow from the parasite to the host is possible, as parasitic larvae can develop into queens. However, we show that despite sympatric coexistence for more than a decade, gene flow between host and social parasite (F(st) = 0.32) and hybridizations (0.71%) are rare, resulting in reproductive isolation. Our data suggest a new barrier to gene flow in sympatry, which is not based on assortative matings but on thelytoky and reproductive division of labour in eusocial insects, thereby suggesting a new potential pathway to Emery's rule.     

Gene flow,divergent selection and resistance to introgression in two species of morning glories (Ipomoea)

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.     

Dynamics of drift, gene flow, and selection during speciation in Silene

The mechanics of speciation with gene flow are still unclear. Disparity among genes in population differentiation (F(ST)) between diverging species is often interpreted as evidence for semipermeable species boundaries, with selection preventing "key" genes from introgressing despite ongoing gene flow. However, F(ST) can remain high before it reaches equilibrium between the lineage sorting of species divergence and the homogenizing effects of gene flow (via secondary contact). Thus, when interpreting F(ST), the dynamics of drift, gene flow, and selection need to be taken into account. We illustrate this view with a multigenic analyses of gene flow and selection in three closely related Silene species, S. latifolia, S. dioica, and S. diclinis. We report that although S. diclinis appears to have evolved in allopatry, isolation with (bidirectional) gene flow between S. latifolia and S. dioica is likely, perhaps as a result of parapatric speciation followed by more extensive sympatry. Interestingly, we detected the signatures of apparently independent instances of positive selection at the same locus in S. latifolia and S. dioica. Despite gene flow between the species, the adaptive alleles have not crossed the species boundary, suggesting that this gene has independently undergone species-specific (diversifying or parallel) selection.     

Species-specific usage of sleeping sites in two sympatric mouse lemur species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar

We investigated the sleeping site ecology of two sympatric mouse lemur species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar during the second half of the dry season with respect to the type, quality, and usage pattern of the sleeping sites, as well as to social sleeping habits and response to potential threats. The type and quality of the sleeping sites differed between the two species. M. murinus used protected wooden shelters (tree holes) more frequently than M. ravelobensis, and M. ravelobensis used a broader variety of less protected sites (e.g., branches, lianas, and leaves) than M. murinus. Whereas male M. murinus usually slept alone, and female M. murinus mostly slept in groups, both sexes of M. ravelobensis slept in mixed-sex sleeping groups. M. murinus relied on crypsis in their sleeping sites, whereas M. ravelobensis regularly showed a flight response to the approach of an observer. This behavioral difference could indicate an adaptation to a higher predation risk in less protected sleeping sites. Whereas female M. murinus showed a high site fidelity, male M. murinus and both sexes of M. ravelobensis frequently changed their sleeping sites, which may also be interpreted as an antipredator strategy. The results are discussed with respect to three possible ecological explanations: interspecies competition, restricted resource availability, and niche differentiation. The latter is the most likely explanation for these interspecific differences.