Tri-oceanic connectivity of the supposedly cosmopolitan polychaete,Harmothoe imbricata (Annelida: Polynoidae): insights from the COI marker

ABSTRACT

In this study, we investigated the global population genetic structure of the polychaete, Harmothoe imbricata, to determine connectivity patterns within the species. We sequenced the mtDNA marker, cytochrome c oxidase I (COI) from 29 specimens sampled across three sites from the New England coast of the United States. These were supplemented with 145 archived sequences from GenBank and the Barcode of Life Database, representing 16 global populations, which encompassed three broadly defined marine biogeographic regions: the Atlantic, Pacific and Arctic. The resulting haplotype network and pairwise AMOVA results showed marked structure across all the major biogeographic regions and also provides evidence for cryptic diversity in the species. Haplotypes from Arctic populations were more closely related to each other than those from the northwestern and northeastern Atlantic. Two evolutionary divergent lineages were recovered from Los Angeles, California and Manitoba, Canada. The highest genetic diversity was observed in the Arctic populations, providing evidence for an Arctic origin for H. imbricata. While human-mediated introductions may have likely contributed to some of the genetic patterns observed in this study, future work should incorporate a nuclear DNA component which could shed more light on contemporary movement of this species across large spatial scales.     

Eelgrass meadows in the California Channel Islands and adjacent coast reveal a mosaic of two species, evidence for introgression and variable clonality

BACKGROUND AND AIMS: Seagrasses are important facilitator species in shallow, soft-bottom marine environments worldwide and, in many places, are threatened by coastal development and eutrophication. One narrow-leaved species (Zostera marina) and one wide-leaved species, variously designated as Z. marina, Z. pacifica or Z. asiatica, are found off the California Channel Islands and adjacent California-Mexico coast. The aim of the present study was to confirm species identification genetically and to link patterns of genetic diversity, connectivity and hybridization among and within the populations with historical sea levels (Ice Age) or the contemporary environment. METHODS: Samples (n = 11-100) were collected from 28 sites off five California Channel Islands and six sites off the adjacent coast of southern California and Baja California, Mexico. DNA polymorphisms of the rDNA-ITS (internal transcribed spacer) cistron (nuclear), the matK intron (chloroplast) and nine microsatellite loci (nuclear) were examined in a population genetic and phylogeographic context. KEY RESULTS: All wide-leaved individuals were Z. pacifica, whereas narrow-leaved forms were Z. marina. Microsatellite genotypes were consistent with hybridization between the two species in three populations. The present distribution of Z. pacifica follows a glacial age land mass rather than present oceanographic regimes, but no link was observed between the present distribution of Z. marina and past or present environments. Island populations of Z. marina often were clonal and characterized by low genotypic diversity compared with populations along the Baja California coast. The high level of clonal connectivity around Santa Catalina Island indicated the importance of dispersal and subsequent re-establishment of vegetative fragments. CONCLUSIONS: The pristine environmental conditions of offshore islands do not guarantee maximum genetic diversity. Future restoration and transplantation efforts of seagrasses must recognize cryptic species and consider the degree of both genetic and genotypic variation in candidate donor populations.     

Mitochondrial DNA variation in Northeast Atlantic and Mediterranean populations of Norway lobster, Nephrops norvegicus

Analysis of the genetic structure of the Norway lobster (Nephrops norvegicus), a marine crustacean with high commercial value, was undertaken to gain information regarding the differentiation of Atlantic from Mediterranean populations of marine invertebrates. Restriction fragment length polymorphism analysis of two mitochondrial DNA segments, 3.6 kilobases in total, was performed. Twelve populations from the North Sea, Irish Sea, Portuguese coast and Aegean Sea were analysed. Low levels of differentiation were found among them (F(ST) = 0.018, P < 0.001) and there were no signs of an Atlantic-Mediterranean divide or of an isolation-by-distance scheme of differentiation. Possible reasons for these low levels of differentiation can be found in the recent expansion of N. norvegicus populations. This is supported by the mismatch distribution of pairwise haplotype differences, as well as by the high mean haplotype diversity (h = 0.93) combined with medium nucleotide diversity (pi = 0.0057) (in comparison to values for marine crustaceans or teleosts) found in this study. This combination of high levels of haplotype diversity with moderate to low levels of nucleotide diversity has also been frequently attributed to a recent time of divergence for various marine species. No evidence was found for a Mediterranean refugium during glaciation periods, separate from the Atlantic, as has been reported for some marine species. The Irish Sea population was the most differentiated as a result of reduced levels of diversity. Results are also discussed in the light of future management of N. norvegicus stocks.     

Seascape genomics reveals fine‐scale patterns of dispersal for a reef fish along the ecologically divergent coast of Northwestern Australia

Understanding the drivers of dispersal among populations is a central topic in marine ecology and fundamental for spatially explicit management of marine resources. The extensive coast of Northwestern Australia provides an emerging frontier for implementing new genomic tools to comparatively identify patterns of dispersal across diverse and extreme environmental conditions. Here, we focused on the stripey snapper (Lutjanus carponotatus), which is important to recreational, charter‐based and customary fishers throughout the Indo‐West Pacific. We collected 1,016 L. carponotatus samples at 51 locations in the coastal waters of Northwestern Australia ranging from the Northern Territory to Shark Bay and adopted a genotype‐by‐sequencing approach to test whether realized connectivity (via larval dispersal) was related to extreme gradients in coastal hydrodynamics. Hydrodynamic simulations using CONNIE and a more detailed treatment in the Kimberley Bioregion provided null models for comparison. Based on 4,402 polymorphic single nucleotide polymorphism loci shared across all individuals, we demonstrated significant genetic subdivision between the Shark Bay Bioregion in the south and all locations within the remaining, more northern bioregions. More importantly, we identified a zone of admixture spanning a distance of 180 km at the border of the Kimberley and Canning bioregions, including the Buccaneer Archipelago and adjacent waters, which collectively experiences the largest tropical tidal range and some of the fastest tidal currents in the world. Further testing of the generality of this admixture zone in other shallow water species across broader geographic ranges will be critical for our understanding of the population dynamics and genetic structure of marine taxa in our tropical oceans.     

Genetic structure of the shrub Daphne laureola across the Baetic Ranges, a Mediterranean glacial refugium and biodiversity hotspot

Distribution margins constitute areas particularly prone to random and/or adaptive intraspecific differentiation in plants. This trend may be particularly marked in species discontinuously distributed across mountain ranges, where sharp geographic isolation gradients and habitat boundaries will enhance genetic isolation among populations. In this study, we analysed the level of neutral genetic differentiation among populations of the long-lived shrub Daphne laureola (Thymelaeaceae) across the Baetic Ranges, a glacial refugium and biodiversity hotspot in the western Mediterranean Basin. Within this area, core and marginal populations of D. laureola were compared with regard to their spatial isolation, size, genetic diversity and differentiation. A spatially explicit analysis conducted on the vast majority of the species' known populations in the study area (N = 111) showed that marginal populations (western and eastern) present larger spatial isolation than core populations, but are not smaller. We compared genetic diversity and differentiation between core and marginal populations using a subsample of 15 populations and 225 amplified fragment length polymorphism (AFLP) markers. Core and marginal populations did not differ in genetic diversity, probably because of the occurrence of large populations on the local margins. Western populations were strongly differentiated from the other populations. In addition, spatial and genetic differentiation among populations was larger on the western margin. Eastern populations constituted a genetically homogeneous group closely related to core populations, despite their greater spatial isolation. Results suggest that studies on phenotypic differentiation between core and marginal populations of D. laureola, and presumably other species having discontinuous distributions across the Baetic ranges, should take into account geographical differences in levels of genetic differentiation between the different distribution borders.     

Ecological and historical determinants of population genetic structure and diversity in the Mediterranean shrub Rosmarinus officinalis (Lamiaceae)

Population genetic studies of widespread Mediterranean shrubs are scarce compared with those of trees and narrow endemics or studies from phylogeographical perspectives, despite the key role these species may play in Mediterranean ecosystems. Knowledge on the effect of ecological factors in shaping their genetic patterns is also limited. In this study we investigate genetic diversity and population structure across 18 populations of Rosmarinus officinalis, a Mediterranean shrubland plant. Populations were sampled along two elevational gradients, one each on calcareous and siliceous soils in a mountain system in the eastern Iberian Peninsula, to decipher the effect of ecological factors on the genetic diversity and structure based on 11 microsatellite loci. We found overall high levels of genetic diversity and weak population structure. Genetic diversity increased with elevation, whereas population differentiation was stronger among populations growing on siliceous soils. The nested analysis of elevational gradients within soil types revealed that these general patterns were mostly driven by siliceous populations, whereas calcareous populations were more homogeneous along elevational belts. Bayesian analysis of population structure revealed genetic membership of lowland and high‐elevation populations to different genetic clusters and a higher admixture of intermediate‐elevation populations to both clusters. High‐elevation populations were less differentiated from a hypothetical ancestral cluster, suggesting the persistence of their gene pool during the Pleistocene glaciations. In contrast, lowland populations resulted from more recent divergence. We propose that life‐history and reproductive traits mostly contribute to explain the high levels of genetic diversity and weak population structure, whereas ecological and historical factors mostly contribute to the stronger differentiation of siliceous populations and a rapid expansion of R. officinalis on calcareous soils possibly mediated by human landscape transformations, © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2016, 180 , 50–63.     

Population genetic structure of annual and perennial populations of Zostera marina L. along the Pacific coast of Baja California and the Gulf of California

The Baja California peninsula represents a biogeographical boundary contributing to regional differentiation among populations of marine animals. We investigated the genetic characteristics of perennial and annual populations of the marine angiosperm, Zostera marina, along the Pacific coast of Baja California and in the Gulf of California, respectively. Populations of Z. marina from five coastal lagoons along the Pacific coast and four sites in the Gulf of California were studied using nine microsatellite loci. Analyses of variance revealed significant interregional differentiation, but no subregional differentiation. Significant spatial differentiation, assessed using θ_ST values, was observed among all populations within the two regions. Z. marina populations along the Pacific coast are separated by more than 220 km and had the greatest θ_ST (0.13–0.28) values, suggesting restricted gene flow. In contrast, lower but still significant genetic differentiation was observed among populations within the Gulf of California (θ_ST = 0.04–0.18), even though populations are separated by more than 250 km. This suggests higher levels of gene flow among Gulf of California populations relative to Pacific coast populations. Direction of gene flow was predominantly southward among Pacific coast populations, whereas no dominant polarity in the Gulf of California populations was observed. The test for isolation by distance (IBD) showed a significant correlation between genetic and geographical distances in Gulf of California populations, but not in Pacific coast populations, perhaps because of shifts in currents during El Niño Southern Oscillation (ENSO) events along the Pacific coast.     

The evolution of Hessian fly from the Old World to the New World: Evidence from molecular markers

Eighteen polymorphic microsatellite loci and 11 single‐nucleotide polymorphisms were genotyped in 1 095 individual Hessian fly specimens representing 23 populations from North America, southern Europe, and southwest Asia. The genotypes were used to assess genetic diversity and interrelationship of Hessian fly populations. While phylogenetic analysis indicates that the American populations most similar to Eurasian populations come from the east coast of the United States, genetic distance is least between (Alabama and California) and (Kazakhstan and Spain). Allelic diversity and frequency vary across North America, but they are not correlated with distance from the historically documented point of introduction in New York City or with temperature or precipitation. Instead, the greatest allelic diversity mostly occurs in areas with Mediterranean climates. The microsatellite data indicate a general deficiency for heterozygotes in Hessian fly. The North American population structure is consistent with multiple introductions, isolation by distance, and human‐abetted dispersal by bulk transport of puparia in infested straw or on harvesting equipment.     

Maritime climate influence on chaparral composition and diversity in the coast range of central California

We investigated the hypothesis that maritime climatic factors associated with summer fog and low cloud stratus (summer marine layer) help explain the compositional diversity of chaparral in the coast range of central California. We randomly sampled chaparral species composition in 0.1‐hectare plots along a coast‐to‐interior gradient. For each plot, climatic variables were estimated and soil samples were analyzed. We used Cluster Analysis and Principle Components Analysis to objectively categorize plots into climate zone groups. Climate variables, vegetation composition and various diversity measures were compared across climate zone groups using ANOVA and nonmetric multidimensional scaling. Differences in climatic variables that relate to summer moisture availability and winter freeze events explained the majority of variance in measured conditions and coincided with three chaparral assemblages: maritime (lowland coast where the summer marine layer was strongest), transition (upland coast with mild summer marine layer influence and greater winter precipitation), and interior sites that generally lacked late summer water availability from either source. Species turnover (β‐diversity) was higher among maritime and transition sites than interior sites. Coastal chaparral differs from interior chaparral in having a higher obligate seeder to facultative seeder (resprouter) ratio and by being dominated by various Arctostaphylos species as opposed to the interior dominant, Adenostoma fasciculatum. The maritime climate influence along the California central coast is associated with patterns of woody plant composition and β‐diversity among sites. Summer fog in coastal lowlands and higher winter precipitation in coastal uplands combine to lower late dry season water deficit in coastal chaparral and contribute to longer fire return intervals that are associated with obligate seeders and more local endemism. Soil nutrients are comparatively less important in explaining plant community composition, but heterogeneous azonal soils contribute to local endemism and promote isolated chaparral patches within the dominant forest vegetation along the coast.     

CONCORDANCE BETWEEN GENETIC AND SPECIES DIVERSITY IN CORAL REEF FISHES ACROSS THE PACIFIC OCEAN BIODIVERSITY GRADIENT

The relationship between genetic diversity and species diversity provides insights into biogeography and historic patterns of evolution and is critical for developing contemporary strategies for biodiversity conservation. Although concordant large‐scale clines in genetic and species diversity have been described for terrestrial organisms, whether these parameters co‐vary in marine species remains largely unknown. We examined patterns of genetic diversity for 11 coral reef fish species sampled at three locations across the Pacific Ocean species diversity gradient (Australia: ～1600 species; New Caledonia: ～1400 species; French Polynesia: ～800 species). Combined genetic diversity for all 11 species paralleled the decline in species diversity from West to East, with French Polynesia exhibiting lowest total haplotype and nucleotide diversities. Haplotype diversity consistently declined toward French Polynesia in all and nucleotide diversity in the majority of species. The French Polynesian population of most species also exhibited significant genetic differentiation from populations in the West Pacific. A number of factors may have contributed to the general positive correlation between genetic and species diversity, including location and time of species origin, vicariance events, reduced gene flow with increasing isolation, and decreasing habitat area from West to East. However, isolation and habitat area, resulting in reduced population size, are likely to be the most influential.     

Strong genetic clines and geographical variation in gene flow in the rocky intertidal barnacle Balanus glandula

A long-standing issue in marine biology is identifying spatial scales at which populations of sessile adults are connected by planktonic offspring. We examined the genetic continuity of the acorn barnacle Balanus glandula, an abundant member of rocky intertidal communities of the northeastern Pacific Ocean, and compared these genetic patterns to the nearshore oceanography described by trajectories of surface drifters. Consistent with its broad dispersal potential, barnacle populations are genetically similar at both mitochondrial (cytochrome oxidase I) and nuclear (elongation factor 1-alpha) loci across broad swaths of the species' range. In central California, however, there is a striking genetic cline across 475 km of coastline between northern and southern populations. These patterns indicate that gene flow within central California is far more restricted spatially than among other populations. Possible reasons for the steep cline include the slow secondary introgression of historically separated populations, a balance between diversifying selection and dispersal, or some mix of both. Geographic trajectories of oceanic drifters closely parallel geographical patterns of gene flow. Drifters placed to the north (Oregon; approximately 44 degrees N) and south (Santa Barbara, California; approximately 34 degrees N) of the cline disperse hundreds of kilometers within 40 days, yet over the long-term their trajectories never overlapped. The lack of communication between waters originating in Oregon and southern California probably helps to maintain strong genetic differentiation between these regions. More broadly, the geographical variation in gene flow implies that focusing on species-level averages of gene flow can mask biologically important variance within species which reflects local environmental conditions and historical events.     

Genetic Structure Among 50 Species of the Northeastern Pacific Rocky Intertidal Community

Comparing many species'' population genetic patterns across the same seascape can identify species with different levels of structure, and suggest hypotheses about the processes that cause such variation for species in the same ecosystem. This comparative approach helps focus on geographic barriers and selective or demographic processes that define genetic connectivity on an ecosystem scale, the understanding of which is particularly important for large-scale management efforts. Moreover, a multispecies dataset has great statistical advantages over single-species studies, lending explanatory power in an effort to uncover the mechanisms driving population structure. Here, we analyze a 50-species dataset of Pacific nearshore invertebrates with the aim of discovering the most influential structuring factors along the Pacific coast of North America. We collected cytochrome c oxidase I (COI) mtDNA data from populations of 34 species of marine invertebrates sampled coarsely at four coastal locations in California, Oregon, and Alaska, and added published data from 16 additional species. All nine species with non-pelagic development have strong genetic structure. For the 41 species with pelagic development, 13 show significant genetic differentiation, nine of which show striking F_ST levels of 0.1–0.6. Finer scale geographic investigations show unexpected regional patterns of genetic change near Cape Mendocino in northern California for five of the six species tested. The region between Oregon and Alaska is a second focus of intraspecific genetic change, showing differentiation in half the species tested. Across regions, strong genetic subdivision occurs more often than expected in mid-to-high intertidal species, a result that may reflect reduced gene flow due to natural selection along coastal environmental gradients. Finally, the results highlight the importance of making primary research accessible to policymakers, as unexpected barriers to marine dispersal break the coast into separate demographic zones that may require their own management plans.     

Looking at both sides of the invasion: patterns of colonization in the violet tunicate Botrylloides violaceus

Understanding the ecological and evolutionary forces that shape the genetic structure of invasive populations and facilitate their expansion across a large spectrum of environments is critical for the prediction of spread and management of ongoing invasions. Here, we study the dynamics of postestablishment colonization in the colonial ascidian Botrylloides violaceus, a notorious marine invader. After its initial introduction from the Northwest Pacific, B. violaceus spread rapidly along the Pacific and Atlantic coasts of North America, impacting both aquaculture facilities and natural ecosystems. We compare genetic diversity and patterns of gene flow among 25 populations (N=679) from the West and East coasts, and evaluate the contribution of sexual vs. asexual reproduction to this species' invasion success using data from the mitochondrial cytochrome c oxidase subunit I (COI) gene and 13 nuclear polymorphic microsatellite loci. Our results reveal contrasting patterns of spread in the coastal waters of North America. While the West coast was colonized by noncontiguous (long-distance) dispersal, the East coast invasion appears to have occurred through contiguous (stepping-stone) spread. Molecular data further indicate that although dispersal in colonial ascidians is predominantly achieved through sexually produced propagules, aquaculture practices such as high-pressure washing can facilitate fragmentation and potentially exacerbate infestations and spread via asexual propagules. The results presented here suggest that caution should be used against the general assumption that all invasions, even within a single species, exhibit similar patterns of colonization, as highly contrasting dynamics may transpire in different invaded ranges.     

Genetic diversity in native and introduced populations of the amethyst gem clam Gemma gemma (Totten, 1834) from the U.S. east and west coasts

Reduced genetic diversity due to founder effects often is expected for invasive populations. The present study examined two nuclear gene regions and one mitochondrial gene to evaluate the origins and genetic diversity of Gemma gemma, a ‘stow-away’ that was introduced to California more than 100 years ago with the importation of the Eastern oyster, Crassostrea virginica, from the United States’ Atlantic coast. A previous investigation involving mitochondrial DNA cytochrome-c-oxidase subunit I sequences reported no significant difference in haplotype diversity between the native and introduced populations; however, estimates of allelic (or haplotypic) variability are insensitive to losses of rare alleles that may accompany founder events and population bottlenecks. Estimates of allele richness and the distribution of rare alleles provide more sensitive indicators of such events. The present investigation of introduced and potential source populations identified lower allele richness and number of singleton alleles in California samples. Atlantic coast Gemma exhibit a sharp phylogeographic transition between northeastern (New York through New England) and mid-Atlantic (southern New Jersey through Virginia) subpopulations, which appear latitudinally inverted for the California Gemma populations. These genetic results, and information from the transportation history of the Eastern oyster, help to clarify processes involved in the introduction of this invasive species.     

Nearshore subtidal community structure compared between inner coast and outer coast sites in Southeast Alaska

Processes that structure subarctic marine communities, particularly in glaciated regions, are not well understood. This understanding is needed as a baseline and to manage these communities in the face of future climate-driven changes. This study investigates two coastal regions of Southeast Alaska with the goals to (a) identify and compare patterns of subtidal community structure for macroalgal, fish, macroinvertebrate (>5?cm), and small epibenthic invertebrate (<5?cm) communities between inner coast and outer coast sites and (b) link patterns of community structure to habitat and environmental parameters. Species assemblage and benthic habitat data were used to compare species diversity and community composition at 6?m and 12?m depths at nine inner coast and nine outer coast sites. Multivariate analysis was applied to reduce environmental variables to major gradients, to resolve community structure, and to relate community structure to habitat and environmental variables. Increased salinity and decreased temperature at outer coast sites compared with inner coast sites were associated with community structure, with greater species diversity at outer coast sites at 6?m depth. Invertebrate community composition was associated with benthic habitat, including crust and coralline algae for macroinvertebrates, and algal cover and substrate for small epibenthic invertebrates. This research suggests that marine communities in glaciated regions are strongly influenced by freshwater input and that future climate-driven changes in freshwater input will likely result in marine community composition changes.     

Latitude‐wide genetic patterns reveal historical effects and contrasting patterns of turnover and nestedness at the range peripheries of a tropical marine fish

Few studies have examined core–periphery genetic patterns in tropical marine taxa. The core–periphery hypothesis (CPH) predicts that core populations will have higher genetic diversity and lower genetic differentiation than peripheral populations as a consequence of greater population sizes and population connectivity in the core. However, the applicability of the CPH to many tropical marine taxa may be confounded by their complex population histories and/or high (asymmetric) population connectivity. In this study we investigated genetic patterns (based on mtDNA) across the latitudinal range of the neon damselfish Pomacentrus coelestis (36°N, Japan – 37°S, east Australia). We suggest a novel hypothetical framework for core–periphery genetic patterns and extend typical analyses to include genealogical analyses, partitioned β‐diversity measures (total β_SOR, turnover β_SIM, and nestedness‐resultant β_SNE), and analyses of nestedness. We found that the existence of two divergent lineages of the neon damselfish led levels of genetic diversity to deviate from CPH expectations. When focusing on the widespread lineage (Pacific clade) nucleotide diversity was higher in the core, supporting the CPH. However, genetic patterns differed toward the northern and southern peripheries of the Pacific clade. The turnover of haplotypes (pairwise‐β_sim) increased over distance in the north, indicative of historical colonization with little contemporary migration. In contrast, although turnover was still dominant in the south (β_SIM), there was no relationship to distance (pairwise‐β_sim), suggesting the influence of more contemporary processes. Moreover, the haplotype compositions of populations in the south were nested according to latitude, indicating immigration from lower latitudes toward the southern periphery. By extending the typical characterizations of core–periphery genetic patterns we were able to identify the effects of lineage sympatry on measures of genetic diversity and contrasting demographic histories toward the latitudinal peripheries of the neon damselfish's range.     

Among population differentiation at nuclear genes in native Scots pine (Pinus sylvestris L.) in Scotland

In the Scottish Highlands, Scots pine is at the north-western extreme of its wide natural distribution. Here, the remaining native populations are patchily distributed in highly variable environments, from the more continental, drier eastern Highlands to the milder, wetter Atlantic Ocean coast. As these pinewoods are the remnants of a naturally established forest, they form a valuable system for analysis of genetic and adaptive variation in heterogeneous environments. Using samples from across the Scottish population, we analysed data from nuclear and mitochondrial genes to assess patterns of within and between population genetic variation. Within population diversity levels were high, and significant genetic differentiation among pairs of Scottish populations at relatively small spatial scales was present at several nuclear loci. At these loci, no differentiation had been found among continental populations, even those separated by large geographic distances. Overall, no clear clustering of Scottish samples was found in population structure analysis suggesting that geographically distant populations with high intra-population nucleotide diversity are not strongly isolated or diverged from each other. Scottish populations lacked a mitotype that is widespread in eastern and north-eastern Europe, indicating that pines from that area may not have participated in the most recent colonisation of the British Isles.     

Geographic proximity not a prerequisite for invasion: Hawaii not the source of California invasion by light brown apple moth (Epiphyas postvittana)

Background

The light brown apple moth (LBAM), Epiphyas postvittana (Walker), is native to Australia but invaded England, New Zealand, and Hawaii more than 100 years ago. In temperate climates, LBAM can be a major agricultural pest. In 2006 LBAM was discovered in California, instigating eradication efforts and quarantine against Hawaiian agriculture, the assumption being that Hawaii was the source of the California infestation. Genetic relationships among populations in Hawaii, California, and New Zealand are crucial to understanding LBAM invasion dynamics across the Pacific.

Methodology/Principal Findings

We sequenced mitochondrial DNA (mtDNA) from 1293 LBAM individuals from California (695), Hawaii (448), New Zealand (147), and Australia (3) to examine haplotype diversity and structure among introduced populations, and evaluate the null hypothesis that invasive populations are from a single panmictic source. However, invasive populations in California and New Zealand harbor deep genetic diversity, whereas Hawaii shows low level, shallow diversity.

Conclusions/Significance

LBAM recently has established itself in California, but was in Hawaii and New Zealand for hundreds of generations, yet California and New Zealand show similar levels of genetic diversity relative to Hawaii. Thus, there is no clear relationship between duration of invasion and genetic structure. Demographic statistics suggest rapid expansion occurring in California and past expansions in New Zealand; multiple introductions of diverse, genetically fragmented lineages could contribute to these patterns. Hawaii and California share no haplotypes, therefore, Hawaii is not the source of the California introduction. Paradoxically, Hawaii and California share multiple haplotypes with New Zealand. New Zealand may be the source for the California and Hawaii infestations, but the introductions were independent, and Hawaii was invaded only once. This has significant implications for quarantine, and suggests that probability of invasion is not directly related to geographic distance. Surprisingly, Hawaiian LBAM populations have much lower genetic diversity than California, despite being older.