Montane refugia predict population genetic structure in the Large‐blotched Ensatina salamander

Understanding the biotic consequences of Pleistocene range shifts and fragmentation remains a fundamental goal in historical biogeography and evolutionary biology. Here, we combine species distribution models (SDM) from the present and two late Quaternary time periods with multilocus genetic data (mitochondrial DNA and microsatellites) to evaluate the effect of climate‐induced habitat shifts on population genetic structure in the Large‐blotched Ensatina (Ensatina eschscholtzii klauberi), a plethodontid salamander endemic to middle and high‐elevation conifer forest in the Transverse and Peninsular Ranges of southern California and northern Baja California. A composite SDM representing the range through time predicts two disjunct refugia, one in southern California encompassing the core of the species range and the other in the Sierra San Pedro Mártir of northern Baja California at the southern limit of the species range. Based on our spatial model, we would expect a pattern of high connectivity among populations within the northern refugium and, conversely, a pattern of isolation due to long‐term persistence of the Sierra San Pedro Mártir population. Our genetic results are consistent with these predictions based on the hypothetical refugia in that (i) historical measures of population connectivity among stable areas are correlated with gene flow estimates; and (ii) there is strong geographical structure between separate refugia. These results provide evidence for the role of recent climatic change in shaping patterns of population persistence and connectivity within the Transverse and Peninsular Ranges, an evolutionary hotspot.     

Genetic divergence of populations of the white char <Emphasis Type="Italic">Salvelinus Albus</Emphasis>, northern and southern forms of malma <Emphasis Type="Italic">S. malma</Emphasis> (Salmonidae), by microsatellite loci of DNA

A very small but statistically significant genetic differentiation by allele frequencies (in the Kamchatka basin F _ST = 0.016 and in Kronotskoe Lake F _ST = 0.063) is revealed between sympatric populations of the northern malma Salvelinus malma malma and of the white char S. albus by analysis of nine microsatellite loci. Factor analysis and clusterization of individual multilocal genotypes in the sum sample of sympatric S. malma malma and S. albus also demonstrate their extremely weak genetic differentiation. Divergence between geographically isolated (allopatric) populations of each of these chars exceeds by two-three times the interspecies differentiation in the zones of sympatry. Samples from allopatric populations clusterize by the geographic principle but not by the taxonomic principle. Such a situation may be explained by the absence of reproductive isolation and presence of hybridization between sympatric S. malma malma and S. albus with the simultaneous existence of some limitations of genetic exchange related to the ecological specialization of these chars. Analysis of microsatellite loci revealed highly significant genetic differences between the northern and southern forms of malma S. malma krascheninnikovi by many times exceeding genetic differences between S. malma malma and S. albus. The obtained data do not confirm the independent species status of S. albus but support the opinion of the species level of differences between the northern and southern forms of malma.     

EVOLUTIONARY AND HISTORICAL ANALYSIS OF PROTEIN VARIATION IN THE BLOTCHED FORMS OF SALAMANDERS OF THE ENSATINA COMPLEX (AMPHIBIA: PLETHODONTIDAE)

Geographic variation in 23 to 29 protein-encoding genetic loci was examined in 48 populations of the Ensatina complex, a “ring species” distributed around the Central Valley of California. The samples span two critical links in the chain of morphologically distinct units: the transition from the unblotched to blotched color pattern types in the vicinity of Lassen Peak, northeastern California, and a geographic gap in the range of the complex in the San Gabriel Mountains, southern California. A general pattern of isolation by distance with a regular buildup of genetic distance correlated with increases in geographic distance characterizes the populations studied, with the exception of a little-differentiated group of populations in the northern Sierra Nevada; this region is postulated to be a zone of genetic reticulation characterized by relatively high gene flow. An adaptively significant color pattern is thought to have spread into the northern Sierra Nevada from the south, but protein variants have been introduced both from the north and the south. Genetic distances across the San Gabriel Mountain gap match expectations from the pattern of buildup of genetic distance as a function of geographic distance elsewhere in the complex. A phylogenetic analysis of the protein data supports the reticulation hypothesis; whereas the southernmost populations currently do constitute a monophyletic assemblage, an “extinction experiment” demonstrates that the distinction could be the result of the recent extinction of populations in a present gap in our sampling. The Ensatina complex appears to be a dynamic entity representing several stages in the evolution of species. It is a ring species, and whereas various taxonomic arrangements are possible, no taxonomic changes are proposed.     

Habitat fragmentation in coastal southern California disrupts genetic connectivity in the cactus wren (Campylorhynchus brunneicapillus)

Achieving long‐term persistence of species in urbanized landscapes requires characterizing population genetic structure to understand and manage the effects of anthropogenic disturbance on connectivity. Urbanization over the past century in coastal southern California has caused both precipitous loss of coastal sage scrub habitat and declines in populations of the cactus wren (Campylorhynchus brunneicapillus). Using 22 microsatellite loci, we found that remnant cactus wren aggregations in coastal southern California comprised 20 populations based on strict exact tests for population differentiation, and 12 genetic clusters with hierarchical Bayesian clustering analyses. Genetic structure patterns largely mirrored underlying habitat availability, with cluster and population boundaries coinciding with fragmentation caused primarily by urbanization. Using a habitat model we developed, we detected stronger associations between habitat‐based distances and genetic distances than Euclidean geographic distance. Within populations, we detected a positive association between available local habitat and allelic richness and a negative association with relatedness. Isolation‐by‐distance patterns varied over the study area, which we attribute to temporal differences in anthropogenic landscape development. We also found that genetic bottleneck signals were associated with wildfire frequency. These results indicate that habitat fragmentation and alterations have reduced genetic connectivity and diversity of cactus wren populations in coastal southern California. Management efforts focused on improving connectivity among remaining populations may help to ensure population persistence.     

A biogeographic genetic approach for testing the role of reinforcement: the case of Drosophila pseudoobscura and D. persimilis

Abstract.— The role of reinforcement in speciation can be explained by two distinct models. In model I, two diverged populations hybridize and produce fertile hybrids that successfully backcross (hybridization with gene flow). In model II, two populations hybridize but succeeding backcrosses are unproductive (hybridization without gene flow). Using Drosophila persimilis and D. pseudoobscura , we have tested model I by comparing the extent of heterospecific introgression in sympatric versus allopatric populations. We show that certain expectations of this particular model of reinforcement, which is based on hybridization and gene flow between divergent populations after secondary contact, are not realized in these two species. The evidence consists of the similarity of genetic distances as well as proportions of unique/rare alleles between sympatric and allopatric heterospecific populations and a negative correlation between genetic distance and geographical distance between heterospecific populations, which suggests ecological differentiation. This approach in quantifying differential gene flow has important consequences to studies that compare sympatric and allopatric isolation using genetic distance. Following model I, one would expect a pattern of higher prezygotic isolation in sympatric species compared to allopatric species of the same genetic distance simply as a result of an underestimation of genetic distance due to introgression between sympatric populations. We suggest more parsimonious explanations such as reinforcement without genetic exchange (model II) and ecological differentiation, which require high levels of preexisting reproductive isolation between populations.     

A tree of tree frogs around the Black Sea

Speciation, the process by which one species evolves into two or more, is a major focus of ongoing debate, particularly regarding the geographic context in which it occurs. Geographic models of speciation tend to fall into discrete categories, typically referred to as allopatric, parapatric and sympatric speciation, according to whether two groups evolve reproductive isolation while geographically isolated, differentiated but connected by gene flow, or completely co‐occurring. Yet molecular studies indicate that full development of reproductive isolation can take very long compared with the timescale at which climatic oscillations occur, such that the geographic context of differentiating forms might change often during the long process to full species. Studies of genetic relationships across the ranges of organisms with low‐dispersal distances have the potential to reveal these complex histories. In a particularly elegant example in this issue, Dufresnes et al. ( 2016 ) use genetic variation and ecological niche modelling to show that a ring of populations of the eastern tree frog (Hyla orientalis) surrounding the Black Sea had a complex history of geographic differentiation. Alternating phases of geographic fragmentation and phases of gene flow between neighbouring populations have produced a pattern of gradual genetic change connecting the western, southern and eastern sides of the ring, with the northwestern and northeastern forms being most differentiated. In the north, a population in Crimea appears to have been produced through mixture of the two extreme forms. The overall genetic relationships are reminiscent of those found in ring species, which have been used as prime demonstrations of the process of speciation. The difference, however, is that the terminal forms appear to have mixed rather than be reproductively isolated, although more research is needed to infer whether there might be some reproductive isolation on the northern side of the ring.     

Differentiation of chromosome banding sequence pools and genomic dna in Holarctic natural populations of <Emphasis Type="Italic">Chironomus entis</Emphasis> Shobanov (Diptera,Chironomidae)

Polymorphism and differentiation of the chromosome banding sequence pools and genomic DNA were studied in three natural populations of Chironomus entis from Europe and North America. These populations showed a moderate level of chromosomal polymorphism and high RAPD polymorphism of genomic DNA. The Palearctic and Nearctic populations of this species did not differ significantly in the levels of chromosome and genomic DNA polymorphism. Estimation of the cytogenetic (GD_cg) and genetic (GD_DNA) distances between these C. entis populations showed that their chromosome banding sequence pools and cytogenetic structures are differentiated to a greater extent than genomic DNA. The values of cytogenetic and genetic distances between the Palearctic and Nearctic populations of C. entis are higher than the values of the corresponding distances between the Nearctic populations, but they do not reach the level of divergence between species.     

Genetic differentiation in sympatric wood ants, <Emphasis Type="Italic">Formica rufa</Emphasis> and <Emphasis Type="Italic">F. polyctena</Emphasis>

Summary Direct observations have suggested that the closely related wood ants Formica polyctena and F. rufa represent different social organizations, with high queen number in F. polyctena and a high frequency of monogynous nests in F. rufa. We examined social organization and genetic population structure in a setup where populations of the two species are sympatric and gene flow between the species is possible. Our aim was to compare social organization in the species, and study evolutionary relationships between them. The observed relatedness among colony workers suggested that the difference in the level of polygyny is quantitative rather than qualitative, with a higher queen number in F. polyctena. The observed difference in polygyny was not accompanied by a difference in spatial genetic differentiation which was weak in both species. The genetic distance between the species is consistent with limited interspecific gene flow. Identification of a few possible F. rufa migrants in F. polyctena populations suggests potential interspecific gene flow. Thus, reproductive isolation of the species may not be complete when they are sympatric.Received 14 March 2003; revised 10 October 2003; accepted 20 October 2003.     

Flies across the water: genetic differentiation and reproductive isolation in allopatric desert Drosophila

Between sister species of Drosophila, both pre- and postzygotic reproductive isolation commonly appear by the time a Nei's genetic distance of 0.5 is observed. The degree of genetic differentiation present when allopatric populations of the same Drosophila species exhibit incipient reproductive isolation has not been systematically investigated. Here we compare the relationship between genetic differentiation and pre- and postzygotic isolation among allopatric populations of three cactophilic desert Drosophila: D. mettleri, D. nigrospiracula, and D. mojavensis. The range of all three is interrupted by the Gulf of California, while two species, D. mettleri and D. mojavensis, have additional allopatric populations residing on distant Santa Catalina Island, off the coast of southern California. Significant population structure exists within all three species, but only for allopatric populations of D. mojavensis is significant isolation at the prezygotic level observed. The genetic distances for the relevant populations of D. mojavensis were in the range of 0.12, similar to that for D. mettleri whose greatest D = 0.11 was unassociated with any form of isolation. These observations suggest further investigations of Drosophila populations with genetic distances in this range be undertaken to identify any potential patterns in the relationship between degree of genetic differentiation and the appearance of pre- and/or postzygotic isolation.     

RESOURCE-ASSOCIATED POPULATION SUBDIVISION IN A SYMBIOTIC CORAL-REEF SHRIMP

The importance of sympatric speciation remains controversial. An empirical observation frequently offered in its support is the occurrence of sister taxa living in sympatry but using different resources. To examine the possibility of sympatric differentiation in producing such cases, I measured genetic, behavioral, and demographic differentiation between populations of the tropical sponge-dwelling shrimp Synalpheus brooksi occupying two alternate host species on three reefs in Caribbean Panama. This species belongs to an apparently monophyletic group of ≥ 30 species of mostly obligate, host-specific sponge-dwellers, many of which occur in sympatry. Demographic data demonstrated the potential for disruptive selection imposed by the two host species: shrimp demes from the sponge Agelas clathrodes were consistently denser, poorer in mature females, more heavily parasitized by branchial bopyrid isopods, and less parasitized by thoracic isopods, than conspecific shrimp from the sponge Spheciospongia vesparium. Laboratory assays demonstrated divergence in host preference: shrimp on all three reefs tended to choose their native sponge species more often than did conspecific shrimp from the other host. Because S. brooksi mates within the host, this habitat selection should foster assortative mating by host species. A hierarchical survey of protein-electrophoretic variation also supported host-mediated divergence, revealing the following: (1) shrimp from the two hosts are conspecific, as evidenced by absence of fixed allelic differences at any of nine allozyme loci scored; (2) strong genetic subdivision among populations of this philopatric shrimp on reefs separated by 1–3 km; and (3) significant host-associated genetic differentiation within two of the three reefs. Finally, intersexual aggression (a proxy for mating incompatibility) between shrimp from different host species was significantly elevated on the one reef where host-associated genetic differences were strongest, demonstrating concordance between genetic and behavioral estimates of divergence. Adjacent reefs appear to be semi-independent sites of host-associated differentiation, as evidenced by differences in the degree of host-associated behavioral and genetic differentiation, and in the specific loci involved, on different reefs. In philopatric organisms with highly subdivided populations, such as S. brooksi, resource-associated differentiation can occur independently in different populations, thus providing multiple “experiments” in differentiation and resulting in a mosaic pattern of polymorphism as reflected by neutral genetic markers. Several freshwater fishes, an amphipod, and a snail similarly show independent but remarkably convergent patterns of resource-associated divergence in different conspecific populations, often in the absence of obvious spatial barriers. In each case, substantial differentiation has occurred in the face of continuing gene flow.     

GENETIC POPULATION STRUCTURE AND LEVELS OF GENE FLOW IN THE STREAM DWELLING WATERSTRIDER,AQUARIUS (=GERRIS) REMIGIS (HEMIPTERA: GERRIDAE)

Gene flow, in combination with selection and drift, determines levels of differentiation among local populations. In this study we estimate gene flow in a stream dwelling, flightless waterstrider, Aquarius remigis. Twenty-eight Aquarius remigis populations from Quebec, Ontario, New Brunswick, Iowa, North Carolina, and California were genetically characterized at 15 loci using starch gel electrophoresis. Sampling over two years was designed for a hierarchical analysis of population structure incorporating variation among sites within streams, streams within watersheds, watersheds within regions, and regions within North America. Hierarchical F statistics indicated that only sites within streams maintained enough gene flow to prevent differentiation through drift (Nm = 27.5). Above the level of sites within streams gene flow is highly restricted (Nm ≤ 0.5) and no correlation is found between genetic and geographic distances. This agrees well with direct estimates of gene flow based on mark and recapture data, yielding an N_e of approximately 170 individuals. Previous assignment of subspecific status to Californian A. remigis is not supported by genetic distances between those populations and other populations in North America. Previous suggestion of specific status for south-eastern A. remigis is supported by genetic distances between North Carolina populations and other populations in North America, and a high proportion of region specific alleles in the North Carolina populations. However, because of the high degree of morphological and genetic variability throughout the range of this species, the assignment of specific or subspecific status to parts of the range may be premature.     

Interspecies mating in sympatric species of Syngnathus pipefish

Hybridization is thought to be an important source of novel genetic variation, and interspecific hybridization may increase the adaptive potential of wild populations. While hybridization has not been previously reported in syngnathid fishes (seahorses and pipefish), the sympatric occurrence of closely related species at high densities increases the probability of interspecies mating in this group. Southern California is home to five species of Syngnathus pipefish, and these species frequently co-occur in near-shore eelgrass beds along the California coast. Recent work has identified exceptionally high levels of genetic diversity in southern populations of Syngnathus leptorhynchus, a widespread species which ranges from Mexico to Alaska. Microsatellite genotyping and mitochondrial sequence data are used here to study the population genetics of S. leptorhynchus and S. auliscus at a site in San Diego Bay where they are found to co-occur at high densities. While no adult hybrids were detected in the study population, analysis of male broods indicates that interspecies mating is occurring between the two species. The lack of premating isolating mechanisms between these two relatives suggests that hybridization may be common in sympatric species of Syngnathus.     

The roles of geography and founder effects in promoting host‐associated differentiation in the generalist bogus yucca moth Prodoxus decipiens

There is ample evidence that host shifts in plant‐feeding insects have been instrumental in generating the enormous diversity of insects. Changes in host use can cause host‐associated differentiation (HAD) among populations that may lead to reproductive isolation and eventual speciation. The importance of geography in facilitating this process remains controversial. We examined the geographic context of HAD in the wide‐ranging generalist yucca moth Prodoxus decipiens. Previous work demonstrated HAD among sympatric moth populations feeding on two different Yucca species occurring on the barrier islands of North Carolina, USA. We assessed the genetic structure of P. decipiens across its entire geographic and host range to determine whether HAD is widespread in this generalist herbivore. Population genetic analyses of microsatellite and mtDNA sequence data across the entire range showed genetic structuring with respect to host use and geography. In particular, genetic differentiation was relatively strong between mainland populations and those on the barrier islands of North Carolina. Finer scale analyses, however, among sympatric populations using different host plant species only showed significant clustering based on host use for populations on the barrier islands. Mainland populations did not form population clusters based on host plant use. Reduced genetic diversity in the barrier island populations, especially on the derived host, suggests that founder effects may have been instrumental in facilitating HAD. In general, results suggest that the interplay of local adaptation, geography and demography can determine the tempo of HAD. We argue that future studies should include comprehensive surveys across a wide range of environmental and geographic conditions to elucidate the contribution of various processes to HAD.     

Sunshine versus gold: The effect of population age on genetic structure of an invasive mosquito

The genetic diversity and structure of invasive species are affected by the time since invasion, but it is not well understood how. We compare likely the oldest populations of Aedes aegypti in continental North America with some of the newest to illuminate the range of genetic diversity and structure that can be found within the invasive range of this important disease vector. Aedes aegypti populations in Florida have probably persisted since the 1600‐1700s, while populations in southern California derive from new invasions that occurred in the last 10 years. For this comparison, we genotyped 1,193 individuals from 28 sites at 12 highly variable microsatellites and a subset of these individuals at 23,961 single nucleotide polymorphisms (SNPs). This is the largest sample analyzed for genetic structure for either region, and it doubles the number of southern California populations previously analyzed. As predicted, the older populations (Florida) showed fewer indicators of recent founder effect and bottlenecks; in particular, these populations have dramatically higher genetic diversity and lower genetic structure. Geographic distance and driving distance were not good predictors of genetic distance in either region, especially southern California. Additionally, southern California had higher levels of genetic differentiation than any comparably sized documented region throughout the worldwide distribution of the species. Although population age and demographic history are likely driving these differences, differences in climate and transportation practices could also play a role.     

PATTERNS OF GENETIC VARIATION IN SALIX SECTION LONGIFOLIAE (SALICACEAE)

Enzyme electrophoresis was used to examine genetic relationships within a taxonomically complex group of diploid willows (Salix section Longifoliae). Forty-eight populations representing seven morphogeographic taxa were analyzed for 15 putative enzyme loci. Levels of genetic variation (A, P, H_e, and H_T) differed greatly among populations, but were generally lower than those reported for other species with similar life history traits. In populations and species from southern latitudes, where populations tend to be smaller and more isolated, genetic variation decreased and population differentiation increased. Genetic identity values suggest that the section is divisible into four major elements: 1) Salix interior, 2) S. taxifolia, 3) the S. exigua group, and 4) S. melanopsis. The first three occupy eastern, southern, and western portions of the continent, respectively. Salix melanopsis is sympatric with the S. exigua group in western North America, but is ecologically distinct. The S. exigua group is differentiated into four major geographic entities: S. hindsiana from west of the Sierras in California and southwest Oregon, S. sessilifolia from west of the Cascades in central Oregon to British Columbia, and a northern and southern race of S. exigua in the Intermountain West. The relationships suggested by allozyme data do not agree with traditional taxonomic concepts. Evidence of hybridization was obtained at five points where taxa come into contact, and Salix fluviatilis may be of hybrid origin. Taxonomic confusion appears to be largely the result of frequent hybridization events and overemphasis of a few highly variable morphological traits. The section appears to be a syngameon undergoing geographic speciation, with one taxon gaining reproductive isolation by a shift in habitat preference.     

Gene flow across species boundaries in sympatric, sexually deceptive Ophrys (Orchidaceae) species

Abstract Orchids of the genus Ophrys (Orchidaceae) are pollinated by male bees and wasps through sexual deception. The Ophrys sphegodes group encompasses several closely related species that differ slightly in floral morphology and are pollinated by different solitary bee species. Populations representing different species of the O. sphegodes group often flower simultaneously in sympatry. To test whether gene flow across the species boundaries occurs in these sympatric populations, or whether they are reproductively isolated, we examined the distribution of genetic variation within and among populations and species of this group. We collected at each of five different localities in southern France and Italy two sympatric, co-flowering Ophrys populations, representing six Ophrys species in total. The six microsatellite loci surveyed were highly variable. Genetic differentiation among geographically distant populations of the same species was lower than differentiation among sympatric populations of different species. However, the strength of genetic differentiation among species was among the lowest reported for orchids. Genotype assignment tests and marker-based estimates of gene flow revealed that gene flow across species boundaries occurred and may account for the low observed differentiation among species. These results suggest that sexual deceit pollination in Ophrys may be less specific than thought, or that rare mistakes occur.     

Molecular and morphological divergence in the butterfly genus Lycaeides (Lepidoptera: Lycaenidae) in North America: evidence of recent speciation

Male genital morphology, allozyme allele frequencies and mtDNA sequence variation were surveyed in the butterfly species Lycaeides idas and L. melissa from across much of their range in North America. Despite clear differences in male genital morphology, wing colour patterns and habitat characteristics, genetic variation was not taxonomically or geographically structured and the species were not identifiable by either genetic data set. Genetic distances (Nei's D=0.002–0.078, calculated from allozyme data) between all populations of both species were within the range commonly observed for conspecific populations of other butterflies. The most frequent mtDNA haplotype was present in individuals of both species in populations from southern California to Wisconsin. We conclude that speciation has probably happened recently and the lack of genetic differentiation between the species is the product of either (1) recent or ongoing gene flow at neutral loci, and/or (2) an insufficiency of time for lineage sorting. The evolution of male genital morphology, wing colour patterns and ecological characteristics has proceeded more rapidly than allozyme or mtDNA evolution.     

The Straits of Gibraltar: barrier or bridge to Ibero‐Moroccan bat diversity?

Genetic divergence in bat communities was assessed on both sides of the Straits of Gibraltar and cryptic diversity was examined. Screening was carried out using partial sequences of the mitochondrial (mt)DNA cytochrome b gene on 399 individual bats belonging to the 18 species found on both sides of the Straits of Gibraltar. For those bats that showed important genetic discontinuities, molecular markers (ND1 and nuclear RAG2 genes) were added to expand the sampling process. Phylogenetic reconstructions were obtained using maximum parsinomy, genetic distances, maximum likelihood, and Bayesian criteria. As an estimate of bats' flight performance, we measured for each species the wing aspect ratio and wing loading indexes, and correlated them with the maximum pairwise genetic distances obtained between southern Iberian and northern Moroccan populations. Genetic mtDNA distances between populations on both continents exceed 5% in seven out of 18 bat species analysed and unknown lineages were uncovered within the species complexes Myotis nattereri and Myotis mystacinus . We did not find a general pattern in the degree of permeability of the Straits of Gibraltar for bats. Genetic distances were not correlated with the ability to cross the Straits. Our study shows that the cryptic diversity uncovered among bats continues to increase as more regions are studied.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 434–450.