Tracking an invasion front with environmental DNA

Data from environmental DNA (eDNA) may revolutionize environmental monitoring and management, providing increased detection sensitivity at reduced cost and survey effort. However, eDNA data are rarely used in decision-making contexts, mainly due to uncertainty around (1) data interpretation and (2) whether and how molecular tools dovetail with existing management efforts. We address these challenges by jointly modeling eDNA detection via qPCR and traditional trap data to estimate the density of invasive European green crab (Carcinus maenas), a species for which, historically, baited traps have been used for both detection and control. Our analytical framework simultaneously quantifies uncertainty in both detection methods and provides a robust way of integrating different data streams into management processes. Moreover, the joint model makes clear the marginal information benefit of adding eDNA (or any other) additional data type to an existing monitoring program, offering a path to optimizing sampling efforts for species of management interest. Here, we document green crab eDNA beyond the previously known invasion front and find that the value of eDNA data dramatically increases with low population densities and low traditional sampling effort, as is often the case at leading-edge locations. We also highlight the detection limits of the molecular assay used in this study, as well as scenarios under which eDNA sampling is unlikely to improve existing management efforts.

Functional responses of a cosmopolitan invader demonstrate intraspecific variability in consumer-resource dynamics

BACKGROUND

Variability in the ecological impacts of invasive species across their geographical ranges may decrease the accuracy of risk assessments. Comparative functional response analysis can be used to estimate invasive consumer-resource dynamics, explain impact variability, and thus potentially inform impact predictions. The European green crab (Carcinus maenas) has been introduced on multiple continents beyond its native range, although its ecological impacts appear to vary among populations and regions. Our aim was to test whether consumer-resource dynamics under standardized conditions are similarly variable across the current geographic distribution of green crab, and to identify correlated morphological features.

METHODS

Crabs were collected from multiple populations within both native (Northern Ireland) and invasive regions (South Africa and Canada). Their functional responses to local mussels (Mytilus spp.) were tested. Attack rates and handling times were compared among green crab populations within each region, and among regions (Pacific Canada, Atlantic Canada, South Africa, and Northern Ireland). The effect of predator and prey morphology on prey consumption was investigated.

RESULTS

Across regions, green crabs consumed prey according to a Type II (hyperbolic) functional response curve. Attack rates (i.e., the rate at which a predator finds and attacks prey), handling times and maximum feeding rates differed among regions. There was a trend toward higher attack rates in invasive than in native populations. Green crabs from Canada had lower handling times and thus higher maximum feeding rates than those from South Africa and Northern Ireland. Canadian and Northern Ireland crabs had significantly larger claws than South African crabs. Claw size was a more important predictor of the proportion of mussels killed than prey shell strength.

DISCUSSION

The differences in functional response between regions reflect observed impacts of green crabs in the wild. This suggests that an understanding of consumer-resource dynamics (e.g., the per capita measure of predation), derived from simple, standardized experiments, might yield useful predictions of invader impacts across geographical ranges.

Genomewide evidence of environmentally mediated secondary contact of European green crab (Carcinus maenas) lineages in eastern North America

Genetic-environment associations are increasingly revealed through population genomic data and can occur through a number of processes, including secondary contact, divergent natural selection, or isolation by distance. Here, we investigate the influence of the environment, including seasonal temperature and salinity, on the population structure of the invasive European green crab (Carcinus maenas) in eastern North America. Green crab populations in eastern North America are associated with two independent invasions, previously shown to consist of distinct northern and southern ecotypes, with a contact zone in southern Nova Scotia, Canada. Using a RAD-seq panel of 9,137 genomewide SNPs, we detected 41 SNPs (0.49%) whose allele frequencies were highly correlated with environmental data. A principal components analysis of 25 environmental variables differentiated populations into northern, southern, and admixed sites in concordance with the observed genomic spatial structure. Furthermore, a spatial principal components analysis conducted on genomic and geographic data revealed a high degree of global structure (p < .0001) partitioning a northern and southern ecotype. Redundancy and partial redundancy analyses revealed that among the environmental variables tested, winter sea surface temperature had the strongest association with spatial structuring, suggesting that it is an important factor defining range and expansion limits of each ecotype. Understanding environmental thresholds associated with intraspecific diversity will facilitate the ability to manage current and predict future distributions of this aquatic invasive species.

RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction‐site‐associated DNA sequencing‐derived SNPs to explore fine‐scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise F_ST ~0.001 to 0.15) and spread across many loci, with a mean F_ST ~0.052 and 52% of loci examined characterized by F_ST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.

Modelling the biological invasion of Carcinus maenas (the European green crab)

This paper proposes a system of integro-difference equations to model the spread of Carcinus maenas, commonly called the European green crab, that causes severe damage to coastal ecosystems. A model with juvenile and adult classes is first studied. Here, standard theory of monotone operators for integro-difference equations can be applied and yields explicit formulas for the asymptotic spreading speeds of the juvenile and adult crabs. A second model including an infected class is considered by introducing a castrating parasite Sacculina carcini as a biological control agent. The dynamics are complicated and simulations reveal the occurrence of periodic solutions and stacked fronts. In this case, only conjectures can be made for the asymptotic spreading speeds because of the lack of mathematical theory for non-monotone operators. This paper also emphasizes the need for mathematical studies of non-monotone operators in heterogeneous environments and the existence of stacked front solutions in biological invasion models.