Genetic homogeneity of the invasive lionfish across the Northwestern Atlantic and the Gulf of Mexico based on Single Nucleotide Polymorphisms

Population genetic structure of Pomacea canaliculata in China based on the COI and ITS1 genes

Comprehending the phylogeography of invasive organisms enhances our insight into their distribution dynamics, which is instrumental for the development of effective prevention and management strategies. In China, Pomacea canaliculata and Pomacea maculata are the two most widespread and damaging species of the non-native Pomacea spp.. Given this species' rapid spread throughout country, it is urgent to investigate the genetic diversity and structure of its different geographic populations, a task undertaken in the current study using the COI and ITS1 mitochondrial and ribosomal DNA genes, respectively. The result of this study, based on a nationwide systematic survey, a collection of Pomacea spp., and the identification of cryptic species, showed that there is a degree of genetic diversity and differentiation in P. canaliculata, and that all of its variations are mainly due to differences between individuals within different geographical populations. Indeed, this species contains multiple haplotypes, but none of them form a systematic geographical population structure. Furthermore, the COI gene exhibits higher genetic diversity than the ITS1 gene. Our study further clarifies the invasive pathways and dispersal patterns of P. canaliculata in China to provide a theoretical basis.

The reconstruction of invasion histories with genomic data in light of differing levels of anthropogenic transport

Unravelling the history of range shifts is key for understanding past, current and future species distributions. Anthropogenic transport of species alters natural dispersal patterns and directly affects population connectivity. Studies have suggested that high levels of anthropogenic transport homogenize patterns of genetic differentiation and blur colonization pathways. However, empirical evidence of these effects remains elusive. We compared two range-shifting species (Microcosmus squamiger and Ciona robusta) to examine how anthropogenic transport affects our ability to reconstruct colonization pathways using genomic data. We first investigated shipping networks from the 18th century onwards, cross-referencing these with regions where the species have records to infer how each species has potentially been affected by different levels of anthropogenic transport. We then genotyped thousands of single-nucleotide polymorphisms from 280 M. squamiger and 190 C. robusta individuals collected across their extensive species' ranges and reconstructed colonization pathways. Differing levels of anthropogenic transport did not preclude the elucidation of population structure, though specific inferences of colonization pathways were difficult to discern in some of the considered scenario sets. We conclude that genomic data in combination with information of underlying introduction drivers provide key insights into the historic spread of range-shifting species.

This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

Genetic connectivity of lionfish (Pterois volitans) in marine protected areas of the Gulf of Mexico and Caribbean Sea

Lionfish (Pterois volitans) have rapidly invaded the tropical Atlantic and spread across the wider Caribbean in a relatively short period of time. Because of its high invasion capacity, we used it as a model to identify the connectivity among nine marine protected areas (MPAs) situated in four countries in the Gulf of Mexico and the Caribbean Sea. This study provides evidence of local genetic differentiation of P. volitans in the Gulf of Mexico and the Caribbean Sea. A total of 475 lionfish samples were characterized with 12 microsatellites, with 6-20 alleles per locus. Departures from Hardy-Weinberg equilibrium (HWE) were found in 10 of the 12 loci, all caused by heterozygous excess. Moderate genetic differentiation was observed between Chiriviche, Venezuela and Xcalak, México localities (F ST = 0.012), and between the Los Roques and the Veracruz (F ST = 0.074) sites. STRUCTURE analysis found that four genetic entities best fit our data. A unique genetic group in the Gulf of Mexico may imply that the lionfish invasion unfolded both in a counterclockwise manner in the Gulf of Mexico. In spite of the notable dispersion of P. volitans, our results show some genetic structure, as do other noninvasive Caribbean fish species, suggesting that the connectivity in some MPAs analyzed in the Caribbean is limited and caused by only a few source individuals with subsequent genetic drift leading to local genetic differentiation. This indicates that P. volitans dispersion could be caused by mesoscale phenomena, which produce stochastic connectivity pulses. Due to the isolation of some MPAs from others, these findings may hold a promise for local short-term control of by means of intensive fishing, even in MPAs, and may have regional long-term effects.

Precipitous Declines in Northern Gulf of Mexico Invasive Lionfish Populations Following the Emergence of an Ulcerative Skin Disease

Invasive Indo-Pacific lionfish Pterois volitans/miles have become well-established in many western Atlantic marine habitats and regions. However, high densities and low genetic diversity could make their populations susceptible to disease. We examined changes in northern Gulf of Mexico (nGOM) lionfish populations following the emergence of an ulcerative skin disease in August 2017, when estimated disease prevalence was as high as 40%. Ulcerated female lionfish had 9% lower relative condition compared to non-ulcerated females. Changes in lionfish size composition indicated a potential recruitment failure in early summer 2018, when the proportion of new recruits declined by >80%. Remotely operated vehicle surveys during 2016–2018 indicated lionfish population density declined in 2018 by 75% on natural reefs. The strongest declines (77–79%) in lionfish density were on high-density (>25 lionfish per 100 m²) artificial reefs, which declined to similar levels as low-density (<15 lionfish per 100 m²) artificial reefs that had prior lionfish removals. Fisheries-dependent sampling indicated lionfish commercial spearfishing landings, commercial catch per unit effort (CPUE), and lionfish tournament CPUE also declined approximately 50% in 2018. Collectively, these results provide correlative evidence for density-dependent epizootic population control, have implications for managing lionfish and impacted native species, and improve our understanding of biological invasions.

Single nucleotide polymorphisms reveal a genetic cline across the north-east Atlantic and enable powerful population assignment in the European lobster

Resolving stock structure is crucial for fisheries conservation to ensure that the spatial implementation of management is commensurate with that of biological population units. To address this in the economically important European lobster (Homarus gammarus), genetic structure was explored across the species' range using a small panel of single nucleotide polymorphisms (SNPs) previously isolated from restriction-site-associated DNA sequencing; these SNPs were selected to maximize differentiation at a range of both broad and fine scales. After quality control and filtering, 1,278 lobsters from 38 sampling sites were genotyped at 79 SNPs. The results revealed a pronounced phylogeographic break between the Atlantic and Mediterranean basins, while structure within the Mediterranean was also apparent, partitioned between lobsters from the central Mediterranean and the Aegean Sea. In addition, a genetic cline across the north-east Atlantic was revealed using both putatively neutral and outlier SNPs, but the precise driver(s) of this clinal pattern-isolation by distance, secondary contact, selection across an environmental gradient, or a combination of these factors-remains undetermined. Putatively neutral markers differentiated lobsters from Oosterschelde, an estuary on the Dutch coast, a finding likely explained by past bottlenecks and limited gene flow with adjacent North Sea populations. Building on the findings of our spatial genetic analysis, we were able to test the accuracy of assigning lobsters at various spatial scales, including to basin of origin (Atlantic or Mediterranean), region of origin and sampling location. The predictive model assembled using 79 SNPs correctly assigned 99.7% of lobsters not used to build the model to their basin of origin, but accuracy decreased to region of origin and again to sampling location. These results are of direct relevance to managers of lobster fisheries and hatcheries, and provide the basis for a genetic tool for tracing the origin of European lobsters in the food supply chain.

Genetic structure of Mexican lionfish populations in the southwest Gulf of Mexico and the Caribbean Sea

The recent expansion of the invasive lionfish throughout the Western Hemisphere is one of the most extensively studied aquatic invasions. Molecular studies have improved our understanding of larval dispersal, connectivity, and biogeographical barriers among lionfish populations, but none have included Mexican localities, an important area for the larval dispersal of Pterois volitans through the Western Caribbean and the Gulf of Mexico. Here, we present a genetic analysis of lionfishes collected along Mexican coasts, examining their connectivity with other Caribbean localities (Belize, Cuba, Puerto Rico) and the role of ocean currents on population structure. We collected 213 lionfish samples from seven locations comprising four countries. To evaluate genetic structure, mitochondrial control region and nuclear inter-simple sequence repeat markers were used. We found that lionfish collected along Mexican coasts show a similar haplotype composition (H02 followed by H01 and H04) to other Caribbean locations, and the H03 rare haplotype was not found. Haplotype composition in the southwest Gulf of Mexico suggests a discontinuity between the southern and northern areas of the Gulf of Mexico. The southern area clustered more strongly to the Caribbean region, and this is supported by the complexity of water circulation in the semi-enclosed region of the Gulf of Mexico. Mitochondrial genetic diversity parameters show small values, whereas nuclear markers produce medium to high values. Only nuclear markers highlighted significant genetic differentiation between the southwest Gulf of Mexico and Caribbean region, confirming a phylogeographic break between both regions. Separate analysis of Caribbean locations indicates restricted larval exchange between southern and northern regions of the Mesoamerican Barrier Reef System, potentially in response to regional oceanographic circulation.

Population genomics of rapidly invading lionfish in the Caribbean reveals signals of range expansion in the absence of spatial population structure

Range expansions driven by global change and species invasions may have significant genomic, evolutionary, and ecological implications. During range expansions, strong genetic drift characterized by repeated founder events can result in decreased genetic diversity with increased distance from the center of the historic range, or the point of invasion. The invasion of the Indo-Pacific lionfish, Pterois volitans, into waters off the US East Coast, Gulf of Mexico, and Caribbean Sea provides a natural system to study rapid range expansion in an invasive marine fish with high dispersal capabilities. We report results from 12,759 single nucleotide polymorphism loci sequenced by restriction enzyme-associated DNA sequencing for nine P. volitans sampling areas in the invaded range, including Florida and other sites throughout the Caribbean, as well as mitochondrial control region D-loop data. Analyses revealed low to no spatially explicit metapopulation genetic structure, which is partly consistent with previous finding of little structure within ocean basins, but partly divergent from initial reports of between-basin structure. Genetic diversity, however, was not homogeneous across all sampled sites. Patterns of genetic diversity correlate with invasion pathway. Observed heterozygosity, averaged across all loci within a population, decreases with distance from Florida while expected heterozygosity is mostly constant in sampled populations, indicating population genetic disequilibrium correlated with distance from the point of invasion. Using an F ST outlier analysis and a Bayesian environmental correlation analysis, we identified 256 and 616 loci, respectively, that could be experiencing selection or genetic drift. Of these, 24 loci were shared between the two methods.