Patterns of invasion and colonization of the sea lamprey (Petromyzon marinus) in North America as revealed by microsatellite genotypes

Upstream migrant sea lamprey (Petromyzon marinus) show signs of increasing oxidative stress but maintain aerobic capacity with age

Following the parasitic juvenile phase of their life cycle, sea lamprey (Petromyzon marinus) mature into a reproductive but rapidly aging and deteriorating adult, and typically die shortly after spawning in May or June. However, pre-spawning upstream migrant sea lamprey can be maintained for several months beyond their natural lifespan when held in cold water (∼4-8 °C) under laboratory conditions. We exploited this feature to investigate the interactions between senescence, oxidative stress, and metabolic function in this phylogenetically ancient fish. We investigated how life history traits and mitochondria condition, as indicated by markers of oxidative stress (catalase activity, lipid peroxidation) and aerobic capacity (citrate synthase activity), changed in adult sea lamprey from June to December after capture during their upstream spawning migration. Body mass but not liver mass declined with age, resulting in an increase in hepatosomatic index. Both effects were most pronounced in males, which also tended to have larger livers than females. Lamprey experienced greater oxidative stress with age, as reflected by increasing activity of the antioxidant enzyme catalase and increasing levels of lipid peroxidation in liver mitochondrial isolates over time. Surprisingly, the activity of citrate synthase also increased with age in both sexes. These observations implicate mitochondrial dysfunction and oxidative stress in the senescence of sea lamprey. Due to their unique evolutionary position and the technical advantage of easily delaying the onset of senescence in lampreys using cold water, these animals could represent an evolutionary unique and tractable model to investigate senescence in vertebrates.

Genetic diversity and structure of a recent fish invasion: Tench (Tinca tinca) in eastern North America

Introduced and geographically expanding populations experience similar eco-evolutionary challenges, including founder events, genetic bottlenecks, and novel environments. Theory predicts that reduced genetic diversity resulting from such phenomena limits the success of introduced populations. Using 1900 SNPs obtained from restriction-site-associated DNA sequencing, we evaluated hypotheses related to the invasion history and connectivity of an invasive population of Tench (Tinca tinca), a Eurasian freshwater fish that has been expanding geographically in eastern North America for three decades. Consistent with the reported history of a single introduction event, our findings suggest that multiple introductions from distinct genetic sources are unlikely as Tench had a small effective population size (~114 [95% CI = 106-123] individuals), no strong population subdivision across time and space, and evidence of a recent genetic bottleneck. The large genetic neighbourhood size (220 km) and weak within-population genetic substructure suggested high connectivity across the invaded range, despite the relatively large area occupied. There was some evidence for a small decay in genetic diversity as the species expanded northward, but not southward, into new habitats. As eradicating the species within a ~112 km radius would be necessary to prevent recolonization, eradicating Tench is likely not feasible at watershed-and possibly local-scales. Management should instead focus on reducing abundance in priority conservation areas to mitigate adverse impacts. Our study indicates that introduced populations can thrive and exhibit relatively high levels of genetic diversity despite severe bottlenecks (<1.5% of the ancestral effective population size) and suggests that landscape heterogeneity and population demographics can generate variability in spatial patterns of genetic diversity within a single range expansion.

Population genomics of New Zealand pouched lamprey (kanakana; piharau; Geotria australis)

Pouched lamprey (Geotria australis) or kanakana/piharau is a culturally and ecologically significant jawless fish that is distributed throughout Aotearoa New Zealand. Despite its importance, much remains unknown about historical relationships and gene flow between populations of this enigmatic species within New Zealand. To help inform management, we assembled a draft G. australis genome and completed the first comprehensive population genomics analysis of pouched lamprey within New Zealand using targeted gene sequencing (Cyt-b and COI) and restriction site-associated DNA sequencing (RADSeq) methods. Employing 16 000 genome-wide single nucleotide polymorphisms (SNPs) derived from RADSeq (n = 186) and sequence data from Cyt-b (766 bp, n = 94) and COI (589 bp, n = 20), we reveal low levels of structure across 10 sampling locations spanning the species range within New Zealand. F-statistics, outlier analyses, and STRUCTURE suggest a single panmictic population, and Mantel and EEMS tests reveal no significant isolation by distance. This implies either ongoing gene flow among populations or recent shared ancestry among New Zealand pouched lamprey. We can now use the information gained from these genetic tools to assist managers with monitoring effective population size, managing potential diseases, and conservation measures such as artificial propagation programs. We further demonstrate the general utility of these genetic tools for acquiring information about elusive species.

Pedigree analysis and estimates of effective breeding size characterize sea lamprey reproductive biology

The sea lamprey (Petromyzon marinus) is an invasive species in the Great Lakes and the focus of a large control and assessment program. Current assessment methods provide information on the census size of spawning adult sea lamprey in a small number of streams, but information characterizing reproductive success of spawning adults is rarely available. We used RAD-capture sequencing to genotype single nucleotide polymorphism (SNP) loci for ~1600 sea lamprey larvae collected from three streams in northern Michigan (Black Mallard, Pigeon, and Ocqueoc Rivers). Larval genotypes were used to reconstruct family pedigrees, which were combined with Gaussian mixture analyses to identify larval age classes for estimation of spawning population size. Two complementary estimates of effective breeding size (N b), as well as the extrapolated minimum number of spawners (N s), were also generated for each cohort. Reconstructed pedigrees highlighted inaccuracies of cohort assignments from traditionally used mixture analyses. However, combining genotype-based pedigree information with length-at-age assignment of cohort membership greatly improved cohort identification accuracy. Population estimates across all three streams sampled in this study indicate a small number of successfully spawning adults when barriers were in operation, implying that barriers limited adult spawning numbers but were not completely effective at blocking access to spawning habitats. Thus, the large numbers of larvae present in sampled systems were a poor indicator of spawning adult abundance. Overall, pedigree-based N b and N s estimates provide a promising and rapid assessment tool for sea lamprey and other species.

An assessment of terminology for intraspecific diversity in fishes, with a focus on "ecotypes" and "life histories"

Understanding and preserving intraspecific diversity (ISD) is important for species conservation. However, ISD units do not have taxonomic standards and are not universally recognized. The terminology used to describe ISD is varied and often used ambiguously. We compared definitions of terms used to describe ISD with use in recent studies of three fish taxa: sticklebacks (Gasterosteidae), Pacific salmon and trout (Oncorhynchus spp., "PST"), and lampreys (Petromyzontiformes). Life history describes the phenotypic responses of organisms to environments and includes biological parameters that affect population growth or decline. Life-history pathway(s) are the result of different organismal routes of development that can result in different life histories. These terms can be used to describe recognizable life-history traits. Life history is generally used in organismal- and ecology-based journals. The terms paired species/species pairs have been used to describe two different phenotypes, whereas in some species and situations a continuum of phenotypes may be expressed. Our review revealed overlapping definitions for race and subspecies, and subspecies and ecotypes. Ecotypes are genotypic adaptations to particular environments, and this term is often used in genetic- and evolution-based journals. "Satellite species" is used for situations in which a parasitic lamprey yields two or more derived, nonparasitic lamprey species. Designatable Units, Evolutionary Significant Units (ESUs), and Distinct Population Segments (DPS) are used by some governments to classify ISD of vertebrate species within distinct and evolutionary significant criteria. In situations where the genetic or life-history components of ISD are not well understood, a conservative approach would be to call them phenotypes.

Rapid genetic adaptation to recently colonized environments is driven by genes underlying life history traits

Background

Uncovering the mechanisms underlying rapid genetic adaptation can provide insight into adaptive evolution and shed light on conservation, invasive species control, and natural resource management. However, it can be difficult to experimentally explore rapid adaptation due to the challenges associated with propagating and maintaining species in captive environments for long periods of time. By contrast, many introduced species have experienced strong selection when colonizing environments that differ substantially from their native range and thus provide a “natural experiment” for studying rapid genetic adaptation. One such example occurred when sea lamprey (Petromyzon marinus), native to the northern Atlantic, naturally migrated into Lake Champlain and expanded their range into the Great Lakes via man-made shipping canals.

Results

Utilizing 368,886 genome-wide single nucleotide polymorphisms (SNPs), we calculated genome-wide levels of genetic diversity (i.e., heterozygosity and π) for sea lamprey collected from native (Connecticut River), native but recently colonized (Lake Champlain), and invasive (Lake Michigan) populations, assessed genetic differentiation between all populations, and identified candidate genes that responded to selection imposed by the novel environments. We observed a 14 and 24% reduction in genetic diversity in Lake Michigan and Lake Champlain populations, respectively, compared to individuals from the Connecticut River, suggesting that sea lamprey populations underwent a genetic bottleneck during colonization. Additionally, we identified 121 and 43 outlier genes in comparisons between Lake Michigan and Connecticut River and between Lake Champlain and Connecticut River, respectively. Six outlier genes that contained synonymous SNPs in their coding regions and two genes that contained nonsynonymous SNPs may underlie the rapid evolution of growth (i.e., GHR), reproduction (i.e., PGR, TTC25, STARD10), and bioenergetics (i.e., OXCT1, PYGL, DIN4, SLC25A15).

Conclusions

By identifying the genomic basis of rapid adaptation to novel environments, we demonstrate that populations of invasive species can be a useful study system for understanding adaptive evolution. Furthermore, the reduction in genome-wide levels of genetic diversity associated with colonization coupled with the identification of outlier genes underlying key life history traits known to have changed in invasive sea lamprey populations (e.g., growth, reproduction) illustrate the utility in applying genomic approaches for the successful management of introduced species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07553-x.

Incipient resistance to an effective pesticide results from genetic adaptation and the canalization of gene expression

The resistance of pest species to chemical controls has vast ecological, economic, and societal costs. In most cases, resistance is only detected after spreading throughout an entire population. Detecting resistance in its incipient stages, by comparison, provides time to implement preventative strategies. Incipient resistance can be detected by coupling standard toxicology assays with large-scale gene expression experiments. We apply this approach to a system where an invasive parasite, sea lamprey (Petromyzon marinus), has been treated with the highly effective pesticide 3-trifluoromethyl-4-nitrophenol (TFM) for 60 years. Toxicological experiments revealed that lamprey from treated populations did not have higher survival to TFM exposure than lamprey from untreated populations, demonstrating that full-fledged resistance has not yet evolved. In contrast, we find hundreds of genes differentially expressed in response to TFM in the population with the longest history of exposure, many of which relate to TFM's primary mode of action, the uncoupling of oxidative phosphorylation, and subsequent depletion of ATP. Three genes critical to oxidative phosphorylation, ATP5PB, PLCB1, and NDUFA9, were nearly fixed for alternative alleles in comparisons of SNPs between treated and untreated populations (FST  > 5 SD from the mean). ATP5PB encodes subunit b of ATP synthase and an additional subunit, ATP5F1B, was canalized for high expression in treated populations, but remained plastic in response to TFM treatment in individuals from the untreated population. These combined genomic and transcriptomic results demonstrate that an adaptive, genetic response to TFM is likely driving incipient resistance in a damaging pest species.

Achieving Sea Lamprey Control in Lake Champlain

The control of parasitic sea lamprey in Lake Champlain has been a necessary component of its fishery restoration and recovery goals for 30 years. While adopting the approach of the larger and established sea lamprey control program of the Laurentian Great Lakes, local differences emerged that shifted management focus and effort as the program evolved. Increased investment in lamprey assessment and monitoring revealed under-estimations of population density and distribution in the basin, where insufficient control efforts went unnoticed. As control efforts improved in response to a better understanding of the population, the effects of lamprey on the fishery lessened. A long-term evaluation of fishery responses when lamprey control was started, interrupted, delayed, and enhanced provided evidence of a recurring relationship between the level of control effort applied and the measured suppression of the parasitic sea lamprey population. Changes in levels of control efforts over time showed repeatedly that measurable suppression of the parasitic population required effective control of 80% of the known larval population. Understanding the importance of assessment and monitoring and the relationship between control effort and population suppression has led to recognition that a comprehensive, not incremental, approach is needed to achieve effective control of sea lamprey in Lake Champlain.

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site-associated DNA sequencing (RAD-Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00-0.018) were concordant with previous microsatellite-based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age-1 and age-2 families of full and half-siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non-native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.

Rapid resistance to pesticide control is predicted to evolve in an invasive fish

Xenobiotic resistance is commonly found in species with short generation times such as bacteria, annual plants, and insects. Nevertheless, the fundamental evolutionary principles that govern the spread of resistance alleles hold true for species with longer generation times. One such example could occur with sea lamprey (Petromyzon marinus), a parasitic invasive species in the Laurentian Great Lakes that decimated native fish populations prior to its control with the pesticide 3-trifluoromethyl-4-nitrophenol (TFM). Since the 1950s, tributaries have been treated annually with TFM, where treatments effectively remove most, but not all, larval sea lamprey. We developed an eco-genetic model of sea lamprey to examine factors affecting the evolution of resistance and found that resistance alleles rapidly rise to fixation after 40-80 years of treatment, despite the species' relatively long generation time (4-7 years). The absence of natal homing allows resistant individuals to spread quickly throughout the entire system, but also makes the early detection of resistance challenging. High costs of resistance and density independent reproduction can delay, but not prevent, the onset of resistance. These results illustrate that sea lamprey have the potential to evolve resistance to their primary control agent in the near future, highlighting the urgent need for alternative controls.

Control of invasive sea lampreys using the piscicides TFM and niclosamide: Toxicology, successes & future prospects

The invasion of the Laurentian Great Lakes of North America by sea lampreys (Petromyzon marinus) in the early 20th century contributed to the depletion of commercial, recreational and culturally important fish populations, devastating the economies of communities that relied on the fishery. Sea lamprey populations were subsequently controlled using an aggressive integrated pest-management program which employed barriers and traps to prevent sea lamprey from migrating to their spawning grounds and the use of the piscicides (lampricides) 3-trifluoromethyl-4-nitrophenol (TFM) and niclosamide to eliminate larval sea lampreys from their nursery streams. Although sea lampreys have not been eradicated from the Great Lakes, populations have been suppressed to less than 10% of their peak numbers in the mid-1900s. The ongoing use of lampricides provides the foundation for sea lamprey control in the Great Lakes, one of the most successful invasive species control programs in the world. Yet, significant gaps remain in our understanding of how lampricides are taken-up and handled by sea lampreys, how lampricides exert their toxic effects, and how they adversely affect non-target invertebrate and vertebrates species. In this review we examine what has been learned about the uptake, handling and elimination, and the mode of TFM and niclosamide toxicity in lampreys and in non-target animals, particularly in the last 10 years. It is now clear that the mode of TFM toxicity is the same in non-target fishes and lampreys, in which TFM interferes with oxidative phosphorylation by the mitochondria leading to decreased ATP production. Vulnerability to TFM is related to abiotic factors such as water pH and alkalinity, which we propose changes the relative amounts of the bioavailable un-ionized form of TFM in the gill microenvironment. Niclosamide, which is also a molluscicide used to control snails in areas prone to schistosomiasis infections of humans, also likely works by uncoupling oxidative phosphorylation, but less is known about other aspects of its toxicology. The effects of TFM include reductions in energy stores, particularly glycogen and high energy phosphagens. However, non-target fishes readily recover from sub-lethal TFM exposure as demonstrated by the rapid restoration of energy stores and clearance of TFM. Although both TFM and niclosamide are non-persistent in the environment and critical for sea lamprey control, increasing public and institutional concerns about pesticides in the environment makes it imperative to explore other means of sea lamprey control. Accordingly, we also address possible "next-generation" strategies of sea lamprey control including genetic tools such as RNA interference and CRISPR-Cas9 to impair critical physiological processes (e.g. reproduction, digestion, metamorphosis) in lamprey, and the use of green chemistry to develop more environmentally benign chemical methods of sea lamprey control.

Development of microsatellite markers for the Siberian weasel Mustela sibirica

The Siberian weasel (Mustela sibirica) is widely distributed in mainland Asia, but its introduction into Japan and subsequent expansion have affected the Japanese weasel (M. itatsi). To provide a useful tool for population genetic studies and control of M. sibirica, we developed 10 polymorphic microsatellite markers. Among 40 individuals of M. sibirica collected in Hubei Province, China, the number of alleles per locus varied from 2 to 19, with the observed heterozygosity ranging from 0.050 to 1.000 and the expected heterozygosity ranging from 0.049 to 0.920. None of the loci deviated from Hardy-Weinberg equilibrium. These markers will be useful in further studies investigating the population structure and natural history of M. sibirica, and may thus provide new insights for the efficient management of this species.

Genetic variation and population structure among larval Lethenteron spp. within the Yukon River drainage, Alaska

The absence of information on genetic variation and population structure of lampreys Lethenteron spp. in the eastern part of their distribution limits our understanding of the migration ecology and spatial population genetic structure of the species. We examined genetic variation within and among three aggregations of Lethenteron spp. larvae in the Yukon River drainage, Alaska, using microsatellite genotypes. A total of 120 larval lampreys were genotyped at eight microsatellite loci. Global F_ST was 0.053 (95% CI 0.021-0.086), while pairwise F_ST values ranged from 0.048-0.057. Model-based Bayesian clustering analyses with sample locality priors (LOCPRIOR) identified three distinct, but admixed, genetic clusters that corresponded with the three aggregations. Estimates of contemporary gene flow indicate substantial reciprocal migration among sites consistent with no or low-fidelity natal homing. These results are largely in agreement with previous reports of historic and contemporary gene flow among Lethenteron spp. in other parts of their geographic distribution.

THE OPTIMAL IMPLEMENTATION OF THETROJAN Y CHROMOSOMEERADICATION STRATEGY OF AN INVASIVE SPECIES

Invasive aquatic species continue to be a persistent problem around the world. The Trojan Y Chromosome (TYC) eradication strategy has recently been developed to help fight the problem in aquatic systems by targeting only the invasive species, sparing native marine stock. It involves rearing genetically modified samples of the invasive species and introducing them into the environment to alter the sex ratio of the invasive population. The paper is devoted to finding the optimal implementation of the TYC eradication strategy of an invasive species as well as a modified, potentially more cost-effective strategy. The modified TYC strategy (MTYC) eliminates one round of exposure to sex hormones compared to the TYC strategy. After introducing both strategies, the optimal control problems for each are formulated. The two strategies are compared through numerical simulations. Our results illustrate that the MTYC strategy, with lower implementation costs, is a better strategy option when trying to minimize the overall effective cost in most scenarios.

Population genetic analysis informs the invasion history of the emerging trematode Dicrocoelium dendriticum into Canada

Parasite distributions are constantly changing due to climate change, local and global movement of animals and humans, as well as land use and habitat change. The trematode Dicrocoelium dendriticum is a relatively recent invader of Canada, being first reported in eastern Canada in the 1930s and western Canada in the 1970s. However, historical records are scarce and its emergence is poorly understood. The establishment of this parasite in Canada provides an interesting opportunity to explore the use of population genetic approaches to help elucidate the invasion history of a relatively recently established helminth parasite. In this study, we compare the genetic diversity and population structure of a number of D. dendriticum populations from western and eastern Canada, and compare these with much longer established European populations. Two independent genetic marker systems were used; a microsatellite marker panel and a cytochrome c oxidase 1 (cox1) mitochondrial (mt)DNA sequence marker. We found distinct differences in both genetic diversity and population structure of the different Canadian populations that provide insights into their invasion histories compared with the European populations. Two populations from British Columbia, Canada - Salt Spring and Vancouver Islands - are of low diversity, show evidence of a population bottleneck and are closely related to each other, suggesting a shared recent history of establishment. These west coast populations are otherwise most closely related to those from eastern Canada and western Europe, and in contrast are genetically divergent from those in Cypress Hills, Alberta, Canada. Although the Alberta parasite population is the most recently reported in Canada, being first identified there in the early 1990s, it was the most genetically diverse of those examined and showed a strong pattern of admixture of genotypes present in western and eastern Europe. Overall, our results are consistent with a model in which western Europe is likely the source of flukes on the east coast of Canada, which were then subsequently translocated to the west coast of Canada. The most recently reported D. dendriticum population in Canada appears to have a different history and likely has multiple origins.

Population structure and dispersal routes of an invasive parasite, Fascioloides magna, in North America and Europe

Background

Fascioloides magna (Trematoda: Fasciolidae) is an important liver parasite of a wide range of free-living and domestic ruminants; it represents a remarkable species due to its large spatial distribution, invasive character, and potential to colonize new territories. The present study provides patterns of population genetic structure and admixture in F. magna across all enzootic regions in North America and natural foci in Europe, and infers migratory routes of the parasite on both continents.

Methods

In total, 432 individuals from five North American enzootic regions and three European foci were analysed by 11 microsatellite loci. Genetic data were evaluated by several statistical approaches: (i) the population genetic structure of F. magna was inferred using program STRUCTURE; (ii) the genetic interrelationships between populations were analysed by PRINCIPAL COORDINATES ANALYSIS; and (iii) historical dispersal routes in North America and recent invasion routes in Europe were explored using MIGRATE.

Results

The analysis of dispersal routes of the parasite in North America revealed west-east and south-north lineages that partially overlapped in the central part of the continent, where different host populations historically met. The exact origin of European populations of F. magna and their potential translocation routes were determined. Flukes from the first European focus, Italy, were related to F. magna from northern Pacific coast, while parasites from the Czech focus originated from south-eastern USA, particularly South Carolina. The Danube floodplain forests (third and still expanding focus) did not display relationship with any North American population; instead the Czech origin of the Danube population was indicated. A serial dilution of genetic diversity along the dispersion route across central and eastern Europe was observed. The results of microsatellite analyses were compared to previously acquired outputs from mitochondrial haplotype data and correlated with past human-directed translocations and natural migration of the final cervid hosts of F. magna.

Conclusions

The present study revealed a complex picture of the population genetic structure and interrelationships of North American and European populations, global distribution and migratory routes of F. magna and an origin of European foci.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1811-z) contains supplementary material, which is available to authorized users.

Stochastic models for the Trojan Y-Chromosome eradication strategy of an invasive species

The Trojan Y-Chromosome (TYC) strategy, an autocidal genetic biocontrol method, has been proposed to eliminate invasive alien species. In this work, we develop a Markov jump process model for this strategy, and we verify that there is a positive probability for wild-type females going extinct within a finite time. Moreover, when sex-reversed Trojan females are introduced at a constant population size, we formulate a stochastic differential equation (SDE) model as an approximation to the proposed Markov jump process model. Using the SDE model, we investigate the probability distribution and expectation of the extinction time of wild-type females by solving Kolmogorov equations associated with these statistics. The results indicate how the probability distribution and expectation of the extinction time are shaped by the initial conditions and the model parameters.

Genetic models reveal historical patterns of sea lamprey population fluctuations within Lake Champlain

The origin of sea lamprey (Petromyzon marinus) in Lake Champlain has been heavily debated over the past decade. Given the lack of historical documentation, two competing hypotheses have emerged in the literature. First, it has been argued that the relatively recent population size increase and concomitant rise in wounding rates on prey populations are indicative of an invasive population that entered the lake through the Champlain Canal. Second, recent genetic evidence suggests a post-glacial colonization at the end of the Pleistocene, approximately 11,000 years ago. One limitation to resolving the origin of sea lamprey in Lake Champlain is a lack of historical and current measures of population size. In this study, the issue of population size was explicitly addressed using nuclear (nDNA) and mitochondrial DNA (mtDNA) markers to estimate historical demography with genetic models. Haplotype network analysis, mismatch analysis, and summary statistics based on mtDNA noncoding sequences for NCI (479 bp) and NCII (173 bp) all indicate a recent population expansion. Coalescent models based on mtDNA and nDNA identified two potential demographic events: a population decline followed by a very recent population expansion. The decline in effective population size may correlate with land-use and fishing pressure changes post-European settlement, while the recent expansion may be associated with the implementation of the salmonid stocking program in the 1970s. These results are most consistent with the hypothesis that sea lamprey are native to Lake Champlain; however, the credibility intervals around parameter estimates demonstrate that there is uncertainty regarding the magnitude and timing of past demographic events.

Reconnaissance of contaminants in larval Pacific lamprey (Entosphenus tridentatus) tissues and habitats in the Columbia River Basin, Oregon and Washington, USA

Pacific lampreys (Entosphenus tridentatus) have resided in the Columbia River Basin for millennia and have great ecological and cultural importance. The role of habitat contamination in the recent decline of the species has rarely been studied and was the main objective of this effort. A wide range of contaminants (115 analytes) was measured in sediments and tissues at 27 sites across a large geographic area of diverse land use. This is the largest dataset of contaminants in habitats and tissues of Pacific lamprey in North America and the first study to compare contaminant bioburden during the larval life stage and the anadromous, adult portion of the life cycle. Bioaccumulation of pesticides, flame retardants, and mercury was observed at many sites. Based on available data, contaminants are accumulating in larval Pacific lamprey at levels that are likely detrimental to organism health and may be contributing to the decline of the species.

Of paleo-genes and Perch: what if an "alien" is actually a native?

Documenting whether a biotic taxon is native or alien to an ecosystem has theoretical value for ecological and evolutionary studies, and has practical value because it can potentially identify a taxon as a desirable component of an ecosystem or target it for removal. In some cases, however, such background information is inadequate or unavailable. Here we use paleo-DNA to re-evaluate the historical status of yellow perch in the 6 million acre Adirondack State Park of northern New York. Yellow perch DNA in a 2200-year sediment record reveals a long-term native status for these supposedly alien fish and challenges assumptions that they necessarily exclude native trout from upland lakes. Similar approaches could be applied to other species with uncertain historical distributions and could help to identify unrecognized pockets of biodiversity.

Contrasting population genetic structure among freshwater-resident and anadromous lampreys: the role of demographic history, differential dispersal and anthropogenic barriers to movement

The tendency of many species to abandon migration remains a poorly understood aspect of evolutionary biology that may play an important role in promoting species radiation by both allopatric and sympatric mechanisms. Anadromy inherently offers an opportunity for the colonization of freshwater environments, and the shift from an anadromous to a wholly freshwater life history has occurred in many families of fishes. Freshwater-resident forms have arisen repeatedly among lampreys (within the Petromyzontidae and Mordaciidae), and there has been much debate as to whether anadromous lampreys, and their derived freshwater-resident analogues, constitute distinct species or are divergent ecotypes of polymorphic species. Samples of 543 European river lamprey Lampetra fluviatilis (mostly from anadromous populations) and freshwater European brook lamprey Lampetra planeri from across 18 sites, primarily in the British Isles, were investigated for 13 polymorphic microsatellite DNA loci, and 108 samples from six of these sites were sequenced for 829 bp of mitochondrial DNA (mtDNA). We found contrasting patterns of population structure for mtDNA and microsatellite DNA markers, such that low diversity and little structure were seen for all populations for mtDNA (consistent with a recent founder expansion event), while fine-scale structuring was evident for nuclear markers. Strong differentiation for microsatellite DNA loci was seen among freshwater-resident L. planeri populations and between L. fluviatilis and L. planeri in most cases, but little structure was evident among anadromous L. fluviatilis populations. We conclude that postglacial colonization founded multiple freshwater-resident populations with strong habitat fidelity and limited dispersal tendencies that became highly differentiated, a pattern that was likely intensified by anthropogenic barriers.

Integrative neuro-endocrine pathways in the control of reproduction in lamprey: a brief review

The gonadotropin-releasing hormone (GnRH) system is well known as the main regulator of reproductive physiology in vertebrates. It is also part of a network of brain structures and pathways that integrate information from the internal and external milieu and coordinate the adaptive behavioral and physiological responses to social and reproductive survival needs. In this paper we review the state of knowledge of the GnRH system in relation to the behavior, external, and internal factors that control reproduction in one of the oldest lineage of vertebrates, the lampreys.

Analysis of the Trojan Y-Chromosome eradication strategy for an invasive species

The Trojan Y-Chromosome (TYC) strategy, an autocidal genetic biocontrol method, has been proposed to eliminate invasive alien species. In this work, we analyze the dynamical system model of the TYC strategy, with the aim of studying the viability of the TYC eradication and control strategy of an invasive species. In particular, because the constant introduction of sex-reversed trojan females for all time is not possible in practice, there arises the question: What happens if this injection is stopped after some time? Can the invasive species recover? To answer that question, we perform a rigorous bifurcation analysis and study the basin of attraction of the recovery state and the extinction state in both the full model and a certain reduced model. In particular, we find a theoretical condition for the eradication strategy to work. Additionally, the consideration of an Allee effect and the possibility of a Turing instability are also studied in this work. Our results show that: (1) with the inclusion of an Allee effect, the number of the invasive females is not required to be very low when the introduction of the sex-reversed trojan females is stopped, and the remaining Trojan Y-Chromosome population is sufficient to induce extinction of the invasive females; (2) incorporating diffusive spatial spread does not produce a Turing instability, which would have suggested that the TYC eradication strategy might be only partially effective, leaving a patchy distribution of the invasive species.

Functional and evolutionary implications of the cellular composition of the gill epithelium of feeding adults of a freshwater parasitic species of lamprey, Ichthyomyzon unicuspis

This paper provides the first description of the cellular composition of the gill epithelium of feeding adults of Ichthyomyzon unicuspis Hubbs and Trautman, 1937 (silver lamprey), a parasitic species of lamprey that is confined to fresh water. The surface layer of this epithelium consists solely of pavement cells and intercalated mitochondria-rich cells, which are the only cell types found in all freshwater stages of lampreys and thus considered responsible for the uptake of Na+ and Cl– in hypotonic environments. This epithelium does not contain, however, the chloride cells present during the marine parasitic phase of anadromous lamprey species, such as Petromyzon marinus L., 1758 (sea lamprey), and which are responsible for secreting excess Na+ and Cl–. The absence of this cell type in parasitic adults of I. unicuspis also differs from its presence in parasitic adults of landlocked P. marinus and metamorphosing individuals of the exclusively freshwater nonparasitic species Lethenteron appendix (DeKay, 1842) (American brook lamprey), and which thus reflects the retention of a cell type that was crucial for osmoregulation during the marine phase of their respective anadromous parasitic ancestors. The absence of chloride cells in I. unicuspis is consistent with the hypothesis that Ichthyomyzon, which is at or close to the base of the phylogenetic tree for Northern Hemisphere lampreys (Petromyzontidae), evolved in fresh water or has been confined to fresh water for a very long period.

Neither philopatric nor panmictic: microsatellite and mtDNA evidence suggests lack of natal homing but limits to dispersal in Pacific lamprey

Most species with lengthy migrations display some degree of natal homing; some (e.g. migratory birds and anadromous salmonids) show spectacular feats of homing. However, studies of the sea lamprey (Petromyzon marinus) indicate that this anadromous species locates spawning habitat based on pheromonal cues from larvae rather than through philopatry. Previous genetic studies in the anadromous Pacific lamprey (Entosphenus tridentatus) have both supported and rejected the hypothesis of natal homing. To resolve this, we used nine microsatellite loci to examine the population structure in 965 Pacific lamprey from 20 locations from central British Columbia to southern California and supplemented this analysis with mitochondrial DNA restriction fragment length polymorphism analysis on a subset of 530 lamprey. Microsatellite analysis revealed (i) relatively low but often statistically significant genetic differentiation among locations (97% pairwise F(ST) values were <0.04 but 73.7% were significant); and (ii) weak but significant isolation by distance (r(2) = 0.0565, P = 0.0450) but no geographic clustering of samples. The few moderate F(ST) values involved comparisons with sites that were geographically distant or far upstream. The mtDNA analysis--although providing less resolution among sites (only 4.7%F(ST) values were significant)--was broadly consistent with the microsatellite results: (i) the southernmost site and some sites tributary to the Salish Sea were genetically distinct; and (ii) southern sites showed higher haplotype and private haplotype richness. These results are inconsistent with philopatry, suggesting that anadromous lampreys are unusual among species with long migrations, but suggest that limited dispersal at sea precludes panmixia in this species.

Formation of a fluvial non-parasitic population of Lethenteron camtschaticum as the first step in petromyzontid speciation

To elucidate the petromyzontid speciation process, the genetic independence of the fluvial non-parasitic populations within the anadromous parasitic Lethenteron camtschaticum was estimated by using polymorphic microsatellite loci. Abundant gene flow was revealed in multitemporal scales between potentially sympatric populations, suggesting ongoing gene flow resulting from imperfect size-assortative mating between them and plastic determination of life histories. On the contrary, landlocked fluvial non-parasitic populations in the upper region of dams were genetically divergent from anadromous parasitic populations. The temporal heterogeneity of gene flow, i.e. contemporary little gene flow but significant gene flow over the long-term between the landlocked fluvial non-parasitic and anadromous parasitic populations was elucidated. In addition, the divergence time of isolation of the landlocked populations from the ancestral anadromous parasitic population was estimated to have occurred 17.9-428.2 years ago, which includes the construction times of an initial dam c. 90 years ago. These instances indicate that the landlocked populations should have very recently been established, and subsequent accumulation of divergence and development of reproductive isolation are predicted. The present landlocked fluvial non-parasitic populations should be analogous to the founder populations in terms of petromyzontid speciation. The data also strongly support the hypothesis of multitemporal and multispatial speciation in the petromyzontid stem-satellite species complex.

Processing and translation initiation of non-long terminal repeat retrotransposons by hepatitis delta virus (HDV)-like self-cleaving ribozymes

Many non-long terminal repeat (non-LTR) retrotransposons lack internal promoters and are co-transcribed with their host genes. These transcripts need to be liberated before inserting into new loci. Using structure-based bioinformatics, we show that several classes of retrotransposons in phyla-spanning arthropods, nematodes, and chordates utilize self-cleaving ribozymes of the hepatitis delta virus (HDV) family for processing their 5' termini. Ribozyme-terminated retrotransposons include rDNA-specific R2, R4, and R6, telomere-specific SART, and Baggins and RTE. The self-scission of the R2 ribozyme is strongly modulated by the insertion site sequence in the rDNA, with the most common insertion sequences promoting faster processing. The ribozymes also promote translation initiation of downstream open reading frames in vitro and in vivo. In some organisms HDV-like and hammerhead ribozymes appear to be dedicated to processing long and short interspersed elements, respectively. HDV-like ribozymes serve several distinct functions in non-LTR retrotransposition, including 5' processing, translation initiation, and potentially trans-templating.

Changes in the cellular composition of the gill epithelium during the life cycle of a nonparasitic lamprey: functional and evolutionary implications

Several lamprey species form pairs, comprising an anadromous parasitic species and a derivative nonparasitic species that neither leaves fresh water nor feeds as an adult. This paper provides the first description of the radical changes undergone by the cellular composition of the gill epithelium during the major phases in the life cycle of a nonparasitic lamprey (American brook lamprey, Lethenteron appendix (DeKay, 1842) (=  Lampetra appendix (DeKay, 1842)) and discusses their potential functional and evolutionary significance. The gill epithelium of the larva of L. appendix contains ammocoete mitochondrion-rich cells (AMRCs), intercalated mitochondrion-rich cells, and pavement cells, as does that of the larva of anadromous parasitic species which likewise lives in fresh water. By the completion of metamorphosis, the AMRCs have disappeared and well-developed chloride cells have been produced, the latter cell type being essential for osmoregulation by its closely related anadromous species in hypertonic environments. By the attainment of sexual maturity, the chloride cells have been lost. Such changes in the timing of chloride cell representation could help account for the ability of some metamorphosing, but not mature individuals of another nonparasitic species ( Lampetra planeri (Bloch, 1784)), to osmoregulate in up to 70% of seawater. The well-developed chloride cells in the nonparasitic L. appendix represent the retention of an ancestral character.

Population genetics and management units of invasive common carp Cyprinus carpio in the Murray-Darling Basin, Australia

Common carp Cyprinus carpio were introduced into Australia on several occasions and are now the dominant fish in the Murray-Darling Basin (MDB), the continent's largest river system. In this study, variability at 14 microsatellite loci was examined in C. carpio (n = 1037) from 34 sites throughout the major rivers in the MDB, from 3 cultured populations, from Prospect Reservoir in the Sydney Basin and from Lake Sorrell in Tasmania. Consistent with previous studies, assignment testing indicated that the Boolara, Yanco and koi strains of C. carpio are present in the MDB. Unique to this study, however, the Prospect strain was widely distributed throughout the MDB. Significant genetic structuring of populations (Fisher's exact test, AMOVA and distribution of the different strains) amongst the MDB sub-drainages was detected, and was strongly associated with contemporary barriers to dispersal and population history. The distributions of the strains were used to infer the history of introduction and spread of C. carpio in the MDB. Fifteen management units are proposed for control programmes that have high levels of genetic diversity, contain multiple interbreeding strains and show no evidence of founder effects or recent population bottlenecks.

Sea lamprey Petromyzon marinus: an exception to the rule of homing in anadromous fishes

Anadromous fishes are believed to make regular circuits of migration in the sea before homing to their natal rivers. Sea lamprey Petromyzon marinus is an anadromous fish that is an exception to this life-history pattern. It also differs from other anadromous fishes in that its adult phase is parasitic, a feeding strategy that should make homing problematic for lamprey cohorts that become widely dispersed through transport by the diverse hosts they parasitize. We sequenced a portion of the mitochondrial DNA control region from sea lampreys collected from 11 North American east coast rivers to test for genetic evidence of homing. There were no significant differences (chi2=235.1, p=0.401) in haplotype frequencies among them, with almost 99 per cent of haplotypic diversity occurring within populations. These findings, together with concordant genetic results from other geographical regions and ancillary information on pheromonal communication, suggest that sea lamprey does not home but rather exhibits regional panmixia while using a novel 'suitable river' strategy to complete its life cycle.

Invasion genetics of the Eurasian round goby in North America: tracing sources and spread patterns

The Eurasian round goby Neogobius melanostomus (Apollonia melanostoma) invaded the North American Great Lakes in 1990 through ballast water, spread rapidly, and now is widely distributed and moving through adjacent tributaries. We analyse its genetic diversity and divergence patterns among 25 North American (N = 744) and 22 Eurasian (N = 414) locations using mitochondrial DNA cytochrome b gene sequences and seven nuclear microsatellite loci in order to: (i) identify the invasion's founding source(s), (ii) test for founder effects, (iii) evaluate whether the invasive range is genetically heterogeneous, and (iv) determine whether fringe and central areas differ in genetic diversity. Tests include F(ST) analogues, neighbour-joining trees, haplotype networks, Bayesian assignment, Monmonier barrier analysis, and three-dimensional factorial correspondence analysis. We recovered 13 cytochrome b haplotypes and 232 microsatellite alleles in North America and compared these to variation we previously described across Eurasia. Results show: (i) the southern Dnieper River population was the primary Eurasian donor source for the round goby's invasion of North America, likely supplemented by some alleles from the Dniester and Southern Bug rivers, (ii) the overall invasion has high genetic diversity and experienced no founder effect, (iii) there is significant genetic structuring across North America, and (iv) some expansion areas show reduced numbers of alleles, whereas others appear to reflect secondary colonization. Sampling sites in Lake Huron's Saginaw Bay and Lake Ontario significantly differ from all others, having unique alleles that apparently originated from separate introductions. Substantial genetic variation, multiple founding sources, large number of propagules, and population structure thus likely aided the goby's ecological success.