Genetic patterns across an invasion's history: a test of change versus stasis for the Eurasian round goby in North America

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.

The Genomic Processes of Biological Invasions: From Invasive Species to Cancer Metastases and Back Again

The concept of invasion is useful across a broad range of contexts, spanning from the fine scale landscape of cancer tumors up to the broader landscape of ecosystems. Invasion biology provides extraordinary opportunities for studying the mechanistic basis of contemporary evolution at the molecular level. Although the field of invasion genetics was established in ecology and evolution more than 50 years ago, there is still a limited understanding of how genomic level processes translate into invasive phenotypes across different taxa in response to complex environmental conditions. This is largely because the study of most invasive species is limited by information about complex genome level processes. We lack good reference genomes for most species. Rigorous studies to examine genomic processes are generally too costly. On the contrary, cancer studies are fortified with extensive resources for studying genome level dynamics and the interactions among genetic and non-genetic mechanisms. Extensive analysis of primary tumors and metastatic samples have revealed the importance of several genomic mechanisms including higher mutation rates, specific types of mutations, aneuploidy or whole genome doubling and non-genetic effects. Metastatic sites can be directly compared to primary tumor cell counterparts. At the same time, clonal dynamics shape the genomics and evolution of metastatic cancers. Clonal diversity varies by cancer type, and the tumors’ donor and recipient tissues. Still, the cancer research community has been unable to identify any common events that provide a universal predictor of “metastatic potential” which parallels findings in evolutionary ecology. Instead, invasion in cancer studies depends strongly on context, including order of events and clonal composition. The detailed studies of the behavior of a variety of human cancers promises to inform our understanding of genome level dynamics in the diversity of invasive species and provide novel insights for management.

Evolutionary trajectory of fish Piscine novirhabdovirus (=Viral Hemorrhagic Septicemia Virus) across its Laurentian Great Lakes history: Spatial and temporal diversification

Piscine novirhabdovirus = Viral Hemorrhagic Septicemia Virus (VHSV) first appeared in the Laurentian Great Lakes with large outbreaks from 2005 to 2006, as a new and novel RNA rhabdovirus subgenogroup (IVb) that killed >30 fish species. Interlude periods punctuated smaller more localized outbreaks in 2007, 2010, and 2017, although some fishes tested positive in the intervals. There have not been reports of outbreaks or positives from 2018, 2019, or 2020. Here, we employ a combined population genetics and phylogenetic approach to evaluate spatial and temporal evolutionary trajectory on its G-gene sequence variation, in comparison with whole-genome sequences (11,083 bp) from a subset of 44 individual isolates (including 40 newly sequenced ones). Our results show that IVb (N = 184 individual fish isolates) diversified into 36 G-gene haplotypes from 2003 to 2017, stemming from two originals ("a" and "b"). G-gene haplotypes "a" and "b" differed by just one synonymous single-nucleotide polymorphism (SNP) substitution, remained the most abundant until 2011, then disappeared. Group "a" descendants (14 haplotypes) remained most prevalent in the Upper and Central Great Lakes, with eight (51%) having nonsynonymous substitutions. Group "b" descendants primarily have occurred in the Lower Great Lakes, including 22 haplotypes, of which 15 (68%) contained nonsynonymous changes. Evolutionary patterns of the whole-genome sequences (which had 34 haplotypes among 44 isolates) appear congruent with those from the G-gene. Virus populations significantly diverged among the Upper, Central, and Lower Great Lakes, diversifying over time. Spatial divergence was apparent in the overall patterns of nucleotide substitutions, while amino acid changes increased temporally. VHSV-IVb thus significantly differentiated across its less than two decades in the Great Lakes, accompanied by declining outbreaks and virulence. Continuing diversification likely allowed the virus to persist at low levels in resident fish populations, and may facilitate its potential for further and future spread to new habitats and nonacclimated hosts.

Increasing confidence for discerning species and population compositions from metabarcoding assays of environmental samples: case studies of fishes in the Laurentian Great Lakes and Wabash River

Community composition data are essential for conservation management, facilitating identification of rare native and invasive species, along with abundant ones. However, traditional capture-based morphological surveys require considerable taxonomic expertise, are time consuming and expensive, can kill rare taxa and damage habitats, and often are prone to false negatives. Alternatively, metabarcoding assays can be used to assess the genetic identity and compositions of entire communities from environmental samples, comprising a more sensitive, less damaging, and relatively time- and cost-efficient approach. However, there is a trade-off between the stringency of bioinformatic filtering needed to remove false positives and the potential for false negatives. The present investigation thus evaluated use of four mitochondrial (mt) DNA metabarcoding assays and a customized bioinformatic Bioinformatic pipeline to increase confidence in species identifications by removing false positives, while achieving high detection probability. Positive controls were used to calculate sequencing error, and results that fell below those cutoff values were removed, unless found with multiple assays. The performance of this approach was tested to discern and identify North American freshwater fishes using lab experiments (mock communities and aquarium experiments) and processing of a bulk ichthyoplankton sample. The method then was applied to field environmental (e) DNA water samples taken concomitant with electrofishing surveys and morphological identifications. This protocol detected 100% of species present in concomitant electrofishing surveys in the Wabash River and an additional 21 that were absent from traditional sampling. Using single 1 L water samples collected from just four locations, the metabarcoding assays discerned 73% of the total fish species that were discerned during four months of an extensive electrofishing river survey in the Maumee River, along with an additional nine species. In both rivers, total fish species diversity was best resolved when all four metabarcoding assays were used together, which identified 35 additional species missed by electrofishing. Ecological distinction and diversity levels among the fish communities also were better resolved with the metabarcoding assays than with morphological sampling and identifications, especially using all four assays together. At the population-level, metabarcoding analyses targeting the invasive round goby Neogobius melanostomus and the silver carp Hypophthalmichthys molitrix identified all population haplotype variants found using Sanger sequencing of morphologically sampled fish, along with additional intra-specific diversity, meriting further investigation. Overall findings demonstrated that the use of multiple metabarcoding assays and custom bioinformatics that filter potential error from true positive detections improves confidence in evaluating biodiversity.

Coalescent models characterize sources and demographic history of recent round goby colonization of Great Lakes and inland waters

The establishment and spread of aquatic invasive species are ecologically and economically harmful and a source of conservation concern internationally. Processes of species invasion have traditionally been inferred from observational data of species presence/absence and relative abundance. However, genetic-based approaches can provide valuable sources of inference. Restriction site-associated DNA sequencing was used to identify and genotype single nucleotide polymorphism (SNP) loci for Round Gobies (Neogobius melanostomus) (N = 440) from 18 sampling locations in the Great Lakes and in three Michigan, USA, drainages (Flint, Au Sable, and Cheboygan River basins). Sampled rivers differed in size, accessibility, and physical characteristics including man-made dispersal barriers. Population levels of genetic diversity and interpopulation variance in SNP allele frequency were used in coalescence-based approximate Bayesian computation (ABC) to statistically compare models representing competing hypotheses regarding source population, postcolonization dispersal, and demographic history in the Great Lakes and inland waters. Results indicate different patterns of colonization across the three drainages. In the Flint River, models indicate a strong population bottleneck (<3% of contemporary effective population size) and a single founding event from Saginaw Bay led to the colonization of inland river segments. In the Au Sable River, analyses could not distinguish potential source populations, but supported models indicated multiple introductions from one source population. In the Cheboygan River, supported models indicated that colonization likely proceeded from east (Lake Huron source) to west among inland locales sampled in the system. Despite the recent occupancy of Great Lakes and inland habitats, large numbers of loci analyzed in an ABC framework enable statistically supported identification of source populations and reconstruction of the direction of inland spread and demographic history following establishment. Information from analyses can direct management actions to limit the spread of invasive species from identified sources and most probable vectors into additional inland aquatic habitats.

Invasion genetics of the silver carp Hypophthalmichthys molitrix across North America: Differentiation of fronts, introgression, and eDNA metabarcode detection

In the 1970s, the introduced silver carp Hypophthalmichthys molitrix (which is indigenous to eastern Asia) escaped from southern U.S. aquaculture to spread throughout the Mississippi River basin, and since has steadily moved northward. This large, prolific filter-feeder reduces food availability for other fishes. It now has reached the threshold of the Laurentian Great Lakes, where it likely will significantly impact food chains and fisheries. Our study evaluates population genetic variability and differentiation of the silver carp using 10 nuclear DNA microsatellite loci, and sequences of two mitochondrial genes-cytochrome b and cytochrome c oxidase subunit 1, along with the nuclear ribosomal protein S7 gene intron 1. We analyze population samples from: two primary Great Lakes' invasion fronts (at the Illinois River outside of Chicago, IL in Lake Michigan and in the Wabash River, which leads into the Maumee River and western Lake Erie), the original establishment "core" in the Lower Mississippi River, and expansion areas in the Upper Mississippi and Missouri rivers. We analyze and compare our results with bighead and other invasive carps, and cyprinid relatives. Results reveal that the silver carp invasion possesses moderate levels of genetic diversity, with more mtDNA haplotypes and unique microsatellite alleles in the "core" Lower Mississippi River population, which also diverges the most. The two invasion fronts also significantly genetically differ. About 3% of individuals (including all populations except the Illinois River) contain a unique and very divergent mtDNA haplotype, which likely stems from historic introgression in Asia with female largescale silver carp H. harmandi. The nuclear microsatellites and S7 sequences of the introgressed individuals do not differ from silver carp and are very distant from bighead carp. These sequence variation data are employed to design and evaluate a targeted high-throughput metabarcoding sequence assay that identifies and distinguishes among species of invasive carps (i.e., silver, bighead, grass, black, and common carps, along with goldfish), as well as native cyprinids, using cytochrome b. Our assay further differentiates among selected silver carp haplotypes (including between H. molitrix and H. harmandi), for use in population genetics and future analyses of spread pathways. We test and evaluate this assay on environmental (e)DNA water samples from 48 bait shops in the Great Lakes' region (along the Lake Erie, Lake St. Clair, and Wabash River watersheds), using positive and negative controls and custom bioinformatic processing. Test results discern silver carp eDNA in four of the shops-three in Lake Erie and one in the Wabash River watershed-and bighead carp from one of the same Lake Erie venues, suggesting that retailers (who often source from established southerly populations) comprise another introduction vector. Our overall findings thus provide key population genetic and phylogenetic data for understanding and tracing introductions, vectors, and spread pathways for silver carp, their variants, and their relatives.

The ghost of introduction past: Spatial and temporal variability in the genetic diversity of invasive smallmouth bass

Understanding the demographic history of introduced populations is essential for unravelling their invasive potential and adaptability to a novel environment. To this end, levels of genetic diversity within the native and invasive range of a species are often compared. Most studies, however, focus solely on contemporary samples, relying heavily on the premise that the historic population structure within the native range has been maintained over time. Here, we assess this assumption by conducting a three-way comparison of the genetic diversity of native (historic and contemporary) and invasive (contemporary) smallmouth bass (Micropterus dolomieu) populations. Analyses of a total of 572 M. dolomieu samples, representing the contemporary invasive South African range, contemporary and historical native USA range (dating back to the 1930s when these fish were first introduced into South Africa), revealed that the historical native range had higher genetic diversity levels when compared to both contemporary native and invasive ranges. These results suggest that both contemporary populations experienced a recent genetic bottleneck. Furthermore, the invasive range displayed significant population structure, whereas both historical and contemporary native US populations revealed higher levels of admixture. Comparison of contemporary and historical samples showed both a historic introduction of M. dolomieu and a more recent introduction, thereby demonstrating that undocumented introductions of this species have occurred. Although multiple introductions might have contributed to the high levels of genetic diversity in the invaded range, we discuss alternative factors that may have been responsible for the elevated levels of genetic diversity and highlight the importance of incorporating historic specimens into demographic analyses.

Human-mediated and natural dispersal of an invasive fish in the eastern Great Lakes

The globally invasive Round Goby (Neogobius melanostomus) was introduced to the Great Lakes around 1990, spreading widely and becoming the dominant benthic fish in many areas. The speed and scope of this invasion is remarkable and calls into question conventional secondary spread models and scenarios. We utilized nine microsatellites to identify large-scale genetic structure in Round Goby populations in the eastern Great Lakes, and assessed the role of colonization vs. secondary transport and dispersal in developing this structure. We identified three clusters, corresponding with Lake Huron, eastern Lake Erie, and western Lake Erie plus eastern Lake Ontario, along with three highly divergent populations. Bottleneck analysis identified founder effects in two divergent populations. Regression analyses of isolation by distance and allelic richness vs. distance from the initial invasion site were consistent with limited migration. However, some populations in eastern Lake Erie and Lake Ontario showed anomalously low genetic distance from the original site of colonization, consistent with secondary transport of large numbers of individuals via ballast water. We conclude that genetic structure of Round Goby in the Great Lakes principally resulted from long-distance secondary transport via ballast water with additional movement of individual via bait buckets and natural dispersal. The success of Round Gobies represents an interesting model for colonization characterization; however, those same attributes present significant challenges for conservation and fisheries management. Current management likely prevents many new species from arriving in the Great Lakes, but fails to address the transport of species within the lakes after they arrive; this is an issue of clear and pressing importance.

Similar yet different: co-analysis of the genetic diversity and structure of an invasive nematode parasite and its invasive mammalian host

Animal parasitic nematodes can cause serious diseases and their emergence in new areas can be an issue of major concern for biodiversity conservation and human health. Their ability to adapt to new environments and hosts is likely to be affected by their degree of genetic diversity, with gene flow between distinct populations counteracting genetic drift and increasing effective population size. The raccoon roundworm (Baylisascaris procyonis), a gastrointestinal parasite of the raccoon (Procyon lotor), has increased its global geographic range after being translocated with its host. The raccoon has been introduced multiple times to Germany, but not all its populations are infected with the parasite. While fewer introduced individuals may have led to reduced diversity in the parasite, admixture between different founder populations may have counteracted genetic drift and bottlenecks. Here, we analyse the population genetic structure of the roundworm and its raccoon host at the intersection of distinct raccoon populations infected with B. procyonis. We found evidence for two parasite clusters resulting from independent introductions. Both clusters exhibited an extremely low genetic diversity, suggesting small founding populations subjected to inbreeding and genetic drift with no, or very limited, genetic influx from population admixture. Comparison of the population genetic structures of both host and parasite suggested that the parasite spread to an uninfected raccoon founder population. On the other hand, an almost perfect match between cluster boundaries also suggested that the population genetic structure of B. procyonis has remained stable since its introduction, mirroring that of its raccoon host.

Genetic structure of brown pelicans (Pelecanus occidentalis) in the northern Gulf of Mexico in the context of human management and disturbance

Environmental disturbances, both natural and anthropogenic, have the capacity to substantially impact animal behavior and abundance, which can in turn influence patterns of genetic diversity and gene flow. However, little empirical information is available on the nature and degree of such changes due to the relative rarity of longitudinal genetic sampling of wild populations at appropriate intervals. Addressing this knowledge gap is therefore of interest to evolutionary biologists, policy makers, and managers. In the past half century, populations of the brown pelican (Pelecanus occidentalis) in the southeastern United States have been exposed to regional extirpations, translocations, colony losses, and oil spills, but potential impacts on genetic diversity and population structure remain unknown. To investigate the cumulative impacts of recent disturbances and management actions, we analyzed seven microsatellite loci using genetic samples collected from 540 nestlings across twelve pelican colonies from two time periods, corresponding to before (n = 305) and after (n = 235) the 2010 Deepwater Horizon oil spill. Pre-2010 populations in Texas were significantly differentiated from Louisiana, Alabama, and Florida populations to the east, with reintroduced populations in southeastern Louisiana having less genetic diversity than sites in Texas, consistent with a recent bottleneck. In contrast, there was no evidence of a geographic component to genetic structure among colonies sampled after the spill, consistent with increased dispersal among sites following the event. This pattern may be associated with reduced philopatry in response to colony abandonment in the areas most heavily impacted by the Deepwater Horizon event, though other factors (e.g., rehabilitation and translocation of oiled birds or colony loss due to erosion and tropical storms) were likely also involved. Future monitoring is necessary to determine if bottlenecks and loss of genetic variation are associated with the oil spill over time, and is recommended for other systems in which disturbance effects may be inferred via repeated genetic sampling.