Interannual variation in the at‐sea behavior of California sea lions (Zalophus californianus)

To be successful, marine predators must alter their foraging behavior in response to changes in their environment. To understand the impact and severity of environmental change on a population it is necessary to first describe typical foraging patterns and identify the underlying variability that exists in foraging behavior. Therefore, we characterized the at‐sea behavior of adult female California sea lions (n = 32) over three years (2003, 2004, and 2005) using satellite transmitters and time‐depth recorders and examined how foraging behavior varied among years. In all years, sea lions traveled on average 84.7 ± 11.1 km from the rookery during foraging trips that were 3.2 ± 0.3 d. Sea lions spent 42.7% ± 1.9% of their time at sea diving and displayed short (2.2 ± 0.2 min), shallow dives (58.5 ± 8.5 m). Among individuals, there was significant variation in both dive behavior and movement patterns, which was found in all years. Among years, differences were found in trip durations, distances traveled, and some dive variables (e.g., dive duration and bottom time) as sea lions faced moderate variability in their foraging habitat (increased sea‐surface temperatures, decreased upwelling, and potential decreased prey abundance). The flexibility we found in the foraging behavior of California sea lions may be a mechanism to cope with environmental variability among years and could be linked to the continuing growth of sea lion populations.     

Population structure and genetic variability of pearl oyster Pinctada mazatlanica along Pacific coasts from Mexico to Panama

Populations of the Calafia pearl oyster Pinctada mazatlanica ofthe American Pacific coasts have been considered endangeredbecause of overfishing and/or alteration to coastal areas. Weassessed genetic variability and the pattern of populationstructure among 9 samples collected from Mexico to Panama, usingmtDNA RFLP analysis of two genes: 12S rRNA and subunit one ofCytochrome oxydase (COI). Haplotype diversity varied from 0.000to 0.856. The Panama population appeared to be monomorphic, whilethe other samples exhibited a level of haplotypic variabilitysimilar to those reported in the literature for the same kind ofanalysis on other bivalves species. A test for the impact ofdemographic history on genetic diversity was applied on thesequence data, and the results were congruent with a recentdecline of population sizes. Genetic differentiation was shown tofollow a scheme of isolation by distance, with low levels ofdifferentiation at the scales of ten to one hundred kilometres,whereas stronger and significant genetic structure was detectedat a larger scale. Three significantly distinct groups could thenbe defined, which correspond to Northern Mexico, Southern Mexico,and Panama.     

Environmental gradients predict the genetic population structure of a coral reef fish in the Red Sea

The relatively recent fields of terrestrial landscape and marine seascape genetics seek to identify the influence of biophysical habitat features on the spatial genetic structure of populations or individuals. Over the last few years, there has been accumulating evidence for the effect of environmental heterogeneity on patterns of gene flow and connectivity in marine systems. Here, we investigate the population genetic patterns of an anemonefish, Amphiprion bicinctus, along the Saudi Arabian coast of the Red Sea. We collected nearly one thousand samples from 19 locations, spanning approximately 1500 km, and genotyped them at 38 microsatellite loci. Patterns of gene flow appeared to follow a stepping‐stone model along the northern and central Red Sea, which was disrupted by a distinct genetic break at a latitude of approximately 19°N. The Red Sea is characterized by pronounced environmental gradients along its axis, roughly separating the northern and central from the southern basin. Using mean chlorophyll‐a concentrations as a proxy for this gradient, we ran tests of isolation by distance (IBD, R² = 0.52) and isolation by environment (IBE, R² = 0.64), as well as combined models using partial Mantel tests and multiple matrix regression with randomization (MMRR). We found that genetic structure across our sampling sites may be best explained by a combined model of IBD and IBE (Mantel: R² = 0.71, MMRR: R² = 0.86). Our results highlight the potential key role of environmental patchiness in shaping patterns of gene flow in species with pelagic larval dispersal. We support growing calls for the integration of biophysical habitat characteristics into future studies of population genetic structure.     

Isolation by oceanographic distance explains genetic structure for Macrocystis pyrifera in the Santa Barbara Channel

Ocean currents are expected to be the predominant environmental factor influencing the dispersal of planktonic larvae or spores; yet, their characterization as predictors of marine connectivity has been hindered by a lack of understanding of how best to use oceanographic data. We used a high-resolution oceanographic model output and Lagrangian particle simulations to derive oceanographic distances (hereafter called transport times) between sites studied for Macrocystis pyrifera genetic differentiation. We build upon the classical isolation-by-distance regression model by asking how much additional variability in genetic differentiation is explained when adding transport time as predictor. We explored the extent to which gene flow is dependent upon seasonal changes in ocean circulation. Because oceanographic transport between two sites is inherently asymmetric, we also compare the explanatory power of models using the minimum or the mean transport times. Finally, we compare the direction of connectivity as estimated by the oceanographic model and genetic assignment tests. We show that the minimum transport time had higher explanatory power than the mean transport time, revealing the importance of considering asymmetry in ocean currents when modelling gene flow. Genetic assignment tests were much less effective in determining asymmetry in gene flow. Summer-derived transport times, in particular for the month of June, which had the strongest current speed, greatest asymmetry and highest spore production, resulted in the best-fit model explaining twice the variability in genetic differentiation relative to models that use geographic distance or habitat continuity. The best overall model also included habitat continuity and explained 65% of the variation in genetic differentiation among sites.     

Crossing the uncrossable: novel trans‐valley biogeographic patterns revealed in the genetic history of low‐dispersal mygalomorph spiders (Antrodiaetidae,Antrodiaetus) from California

Antrodiaetus riversi is a dispersal‐limited, habitat‐specialized mygalomorph spider species endemic to mesic woodlands of northern and central California. Here, we build upon prior phylogeographic research using a much larger geographic sample and include additional nuclear genes, providing more detailed biogeographic insights throughout the range of this complex. Of particular interest is the uncovering of unexpected and replicated trans‐valley biogeographic patterns, where in two separate genetic clades western haplotypes in the California south Coast Ranges are phylogenetically closely related to eastern haplotypes from central and northern Sierran foothills. In both instances, these trans‐valley phylogenetic patterns are strongly supported by multiple genes. These western and eastern populations are currently separated by the Central Valley, a well‐recognized modern‐day and historical biogeographic barrier in California. For one clade, the directionality is clearly northeast to southwest, and all available evidence is consistent with a jump dispersal event estimated at 1.2–1.3 Ma. During this time period, paleogeographic data indicate that northern Sierran rivers emptied to the ocean in the south Coast Ranges, rather than at the San Francisco Bay. For the other trans‐valley clade genetic evidence is less conclusive regarding the mechanism and directionality of biogeographic exchange, although the estimated timeframe is similar (approximately 1.8 Ma). Despite the large number of biogeographic studies previously conducted in central California, to the best of our knowledge no prior studies have discussed or revealed a northern Sierran to south Coast Range biogeographic connection. This uniqueness may reflect the low‐dispersal biology of mygalomorph spiders, where ‘post‐event’ gene exchange rarely erases historical biogeographic signal.     

Drivers for genetic structure at different geographic scales for Pacific red snapper (Lutjanus peru) and yellow snapper (Lutjanus argentiventris) in the tropical eastern Pacific

The tropical eastern Pacific (TEP) is a highly dynamic region and a model system to study how habitat discontinuities affect the distribution of shorefishes, particularly for species that display ontogenetic habitat shifts, including snappers (Lutjanidae). To evaluate the genetic structure of the Pacific red snapper (Lutjanus peru) and the yellow snapper (Lutjanus argentiventris) throughout their distribution range along the TEP, 13 and 11 microsatellite loci were analysed, respectively. The genetic diversity of L. peru (N = 446) and L. argentiventris (N = 170) was evaluated in 10 and 5 localities, respectively, showing slightly higher but non-significant values in the Gulf of California for both species. The genetic structure analysis identified the presence of significant genetic structure in both species, but the locations of the identified barriers for the gene flow differed between species. The principal driver for the genetic structure at large scales >2500 km was isolation by distance. At smaller scales (<250 km), the habitat discontinuity for juveniles and adults and the environmental differences throughout the distribution range represented potential barriers to gene flow between populations for both species.     

A synthesis of genetic connectivity in deep‐sea fauna and implications for marine reserve design

With anthropogenic impacts rapidly advancing into deeper waters, there is growing interest in establishing deep‐sea marine protected areas (MPAs) or reserves. Reserve design depends on estimates of connectivity and scales of dispersal for the taxa of interest. Deep‐sea taxa are hypothesized to disperse greater distances than shallow‐water taxa, which implies that reserves would need to be larger in size and networks could be more widely spaced; however, this paradigm has not been tested. We compiled population genetic studies of deep‐sea fauna and estimated dispersal distances for 51 studies using a method based on isolation‐by‐distance slopes. Estimates of dispersal distance ranged from 0.24 km to 2028 km with a geometric mean of 33.2 km and differed in relation to taxonomic and life‐history factors as well as several study parameters. Dispersal distances were generally greater for fishes than invertebrates with the Mollusca being the least dispersive sampled phylum. Species that are pelagic as adults were more dispersive than those with sessile or sedentary lifestyles. Benthic species from soft‐substrate habitats were generally less dispersive than species from hard substrate, demersal or pelagic habitats. As expected, species with pelagic and/or feeding (planktotrophic) larvae were more dispersive than other larval types. Many of these comparisons were confounded by taxonomic or other life‐history differences (e.g. fishes being more dispersive than invertebrates) making any simple interpretation difficult. Our results provide the first rough estimate of the range of dispersal distances in the deep sea and allow comparisons to shallow‐water assemblages. Overall, dispersal distances were greater for deeper taxa, although the differences were not large (0.3–0.6 orders of magnitude between means), and imbalanced sampling of shallow and deep taxa complicates any simple interpretation. Our analyses suggest the scales of dispersal and connectivity for reserve design in the deep sea might be comparable to or slightly larger than those in shallow water. Deep‐sea reserve design will need to consider the enormous variety of taxa, life histories, hydrodynamics, spatial configuration of habitats and patterns of species distributions. The many caveats of our analyses provide a strong impetus for substantial future efforts to assess connectivity of deep‐sea species from a variety of habitats, taxonomic groups and depth zones.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

Inverse relationship between F and microsatellite polymorphism in the marine fish, walleye pollock (Theragra chalcogramma): implications for resolving weak population structure

Microsatellites have proved to be useful for the detection of weak population structure in marine fishes and other species characterized by large populations and high gene flow. None the less, uncertainty remains about the net effects of the particular mutational properties of these markers, and the wide range of locus polymorphism they exhibit, on estimates of differentiation. We examined the effect of varying microsatellite polymorphism on the magnitude of observed differentiation in a population survey of walleye pollock, Theragra chalcogramma. Genetic differentiation at 14 microsatellite loci among six putative populations from across the North Pacific Ocean and Bering Sea was weak but significant on large geographical scales and conformed to an isolation-by-distance pattern. A negative relationship was found between locus variability and the magnitude of estimated population subdivision. Estimates of F(ST) declined with locus polymorphism, resulting in diminished power to discriminate among samples, and we attribute this loss to the effects of size homoplasy. This empirical result suggests that mutation rates of some microsatellite loci are sufficiently high to limit resolution of weak genetic structure typical of many marine fishes.     

Hierarchical components of genetic variation at a species boundary: population structure in two sympatric varieties of Lupinus microcarpus (Leguminosae)

Lupinus microcarpus is a self-compatible annual plant that forms a species complex of morphologically variable but indeterminate varieties. In order to examine the hypothesis that varieties of L. microcarpus comprise genetically differentiated and reproductively isolated species, populations of L. microcarpus var. horizontalis and var. densiflorus were sampled from an area of sympatry in central California and genotyped using six microsatellite loci. Bayesian clustering divided the total sample into two groups corresponding to the named varieties with extremely low levels of inferred coancestry. Similarly, maximum likelihood and distance methods for genetic assignment placed individuals in two nonoverlapping groups. Evidence for isolation by distance (IBD) within each variety was found at shorter distance classes, but varieties remained differentiated in sympatry. Furthermore, coalescent estimates of divergence time indicate separation within the past 950-5050 generations, with minimal gene flow after divergence. A four-level hierarchical analysis of molecular variance (amova) found significant levels of genetic differentiation among varieties (theta(P) = 0.292), populations within varieties (theta(S) = 0.449), subpopulations within populations (theta(SS) = 0.623), and individuals within subpopulations (f = 0.421); but the greatest degree of differentiation was at the subpopulation level. Although it is sometimes assumed that the magnitude of genetic differences (e.g. F(ST)) should be greater between species than among populations or subpopulations of the same species, shared ancestral polymorphism may lead to relatively low levels of differentiation at the species level, even as the stochastic effects of genetic drift generate higher levels of differentiation at lower hierarchical levels. These results suggest that L. microcarpus var. horizontalis and var. densiflorus are recently diverged yet reproductively isolated species, with high levels of inbreeding resulting from the combined effects of limited gene flow, demographic bottlenecks, and partial selfing in finite, geographically structured populations.     

History,geography and host use shape genomewide patterns of genetic variation in the redheaded pine sawfly (Neodiprion lecontei)

Divergent host use has long been suspected to drive population differentiation and speciation in plant‐feeding insects. Evaluating the contribution of divergent host use to genetic differentiation can be difficult, however, as dispersal limitation and population structure may also influence patterns of genetic variation. In this study, we use double‐digest restriction‐associated DNA (ddRAD) sequencing to test the hypothesis that divergent host use contributes to genetic differentiation among populations of the redheaded pine sawfly (Neodiprion lecontei), a widespread pest that uses multiple Pinus hosts throughout its range in eastern North America. Because this species has a broad range and specializes on host plants known to have migrated extensively during the Pleistocene, we first assess overall genetic structure using model‐based and model‐free clustering methods and identify three geographically distinct genetic clusters. Next, using a composite‐likelihood approach based on the site frequency spectrum and a novel strategy for maximizing the utility of linked RAD markers, we infer the population topology and date divergence to the Pleistocene. Based on existing knowledge of Pinus refugia, estimated demographic parameters and patterns of diversity among sawfly populations, we propose a Pleistocene divergence scenario for N. lecontei. Finally, using Mantel and partial Mantel tests, we identify a significant relationship between genetic distance and geography in all clusters, and between genetic distance and host use in two of three clusters. Overall, our results indicate that Pleistocene isolation, dispersal limitation and ecological divergence all contribute to genomewide differentiation in this species and support the hypothesis that host use is a common driver of population divergence in host‐specialized insects.     

Population genetic structure of the cyclic snowshoe hare (Lepus americanus) in southwestern Yukon,Canada

Spatial population structure has important ecological and evolutionary consequences. Little is known about the population structure of snowshoe hares (Lepus americanus), despite their ecological importance in North American boreal forests. We used seven variable microsatellite DNA loci to determine the spatial genetic structure of snowshoe hares near Kluane Lake, Yukon during a cyclic population peak. We sampled 317 hares at 12 sites separated by distances ranging from 3 to 140 km, and used 46 additional samples from Alaska and Montana. The level of genetic variation was high (13.4 alleles/locus, 0.67 expected heterozygosity) and the distribution of alleles and genotypes was not homogeneous across the sites. The degree of differentiation was low among Yukon sites (FST = 0.015) and between Yukon and Alaska (FST = 0.012), but the Montana site was highly differentiated (FST = 0.20). A weak pattern of isolation by distance was found over the Yukon study area, with an indication that local genetic drift may be important in shaping the regional genetic structure. Landscape barriers expected to influence gene flow did not consistently affect genetic structure, although there was evidence for a partial barrier effect of Kluane Lake. The high level of inferred gene flow confirms that snowshoe hare dispersal is widespread, successful and equal between the sexes. A stepping-stone model of gene flow, potentially influenced by the synchronous density cycle, appears to best explain the observed genetic structure. Our results suggest that despite their dramatic fluctuations in density, snowshoe hares in the northern boreal forest have a large evolutionary effective population size and are not strongly subdivided by either physical or social barriers to gene flow.     

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.     

Living on the edge: the role of geography and environment in structuring genetic variation in the southernmost populations of a tropical oak

Understanding the factors determining genetic diversity and structure in peripheral populations is a long‐standing goal of evolutionary biogeography, yet little empirical information is available for tropical species. In this study, we combine information from nuclear microsatellite markers and niche modelling to analyse the factors structuring genetic variation across the southernmost populations of the tropical oak Quercus segoviensis. First, we tested the hypothesis that genetic variability decreases with population isolation and increases with local habitat suitability and stability since the Last Glacial Maximum (LGM). Second, we employed a recently developed multiple matrix regression with randomisation (MMRR) approach to study the factors associated with genetic divergence among the studied populations and test the relative contribution of environmental and geographic isolation to contemporary patterns of genetic differentiation. We found that genetic diversity was negatively correlated with average genetic differentiation with other populations, indicating that isolation and limited gene flow have contributed to erode genetic variability in some populations. Considering the relatively small size of the study area (<120 km), analyses of genetic structure indicate a remarkable inter‐population genetic differentiation. Environmental dissimilarity and differences in current and past climate niche suitability and their additive effects were not associated with genetic differentiation after controlling for geographic distance, indicating that local climate does not contribute to explain spatial patterns of genetic structure. Overall, our data indicate that geographic isolation, but not current or past climate, is the main factor determining contemporary patterns of genetic diversity and structure within the southernmost peripheral populations of this tropical oak.     

Isolation by distance and vicariance drive genetic structure of a coral reef fish in the Pacific Ocean

We studied the genetic diversity of a coral reef fish species to investigate the origin of the differentiation. A total of 727 Acanthurus triostegus collected from 15 locations throughout the Pacific were analyzed for 20 polymorphic loci. The genetic structure showed limited internal disequilibrium within each population; 3.7% of the loci showed significant Hardy-Weinberg disequilibrium, mostly associated with Adh*, and we subsequently removed this locus from further analysis of geographic pattern. The genetic structure of A. triostegus throughout the tropical Pacific Ocean revealed a strong geographic pattern. Overall, there was significant population differentiation (multilocus F(ST) = 0.199), which was geographically structured according to bootstraps of neighbor-joining analysis on Nei's unbiased genetic distances and AMOVA analysis. The genetic structure revealed five geographic groups in the Pacific Ocean: western Pacific (Guam, Philippines, Palau, and Great Barrier Reef); central Pacific (Solomons, New Caledonia, and Fiji); and three groups made up of the eastern populations, namely Hawaiian Archipelago (north), Marquesas (equatorial), and southern French Polynesia (south) that incorporates Clipperton Island located in the northeastern Pacific. In addition, heterozygosity values were found to be geographically structured with higher values grouped within Polynesian and Clipperton populations, which exhibited lower population size. Finally, the genetic differentiation (F(ST)) was significantly correlated with geographic distance when populations from the Hawaiian and Marquesas archipelagos were separated from all the other locations. These results show that patterns of differentiation vary within the same species according to the spatial scale, with one group probably issued from vicariance, whereas the other followed a pattern of isolation by distance. The geographic pattern for A. triostegus emphasizes the diversity of the evolutionary processes that lead to the present genetic structure with some being more influential in certain areas or according to a particular spatial scale.     

Population genetic structure of the lemon shark (Negaprion brevirostris) in the western Atlantic: DNA microsatellite variation

DNA microsatellite markers were used to characterize the population genetic structure of the lemon shark, Negaprion brevirostris, in the western Atlantic. This study demonstrates for the first time the usefulness of microsatellites to study population genetic structure and mating systems in the Chondricthyes. Lemon sharks (mostly juveniles) were sampled non-destructively from four locations, Gullivan Bay and Marquesas Key in Florida, Bimini, Bahamas, and Atol das Rocas, Brazil. At least 545 individuals were genotyped at each of four dinucleotide loci. The number of alleles per locus ranged from 19 to 43, and expected heterozygosities ranged from 0.69 to 0.90. Relatively little genetic structure was found in the western Atlantic, with small but significant values for estimators of F(ST) and R(ST) among populations, theta (0.016) and rho (0.026), respectively. No sharp discontinuities were found between the Caribbean sites and Brazil, and most alleles were found at all four sites, indicating that gene flow occurs throughout the western Atlantic with no evidence for distinct stocks.