Non-native species spread in a complex network: the interaction of global transport and local population dynamics determines invasion success.

The number of released individuals, which is a component of propagule pressure, is considered to be a major driver for the establishment success of non-native species. However, propagule pressure is often assumed to result from single or few release events, which does not necessarily apply to the frequent releases of invertebrates or other taxa through global transport. For instance, the high intensity of global shipping may result in frequent releases of large numbers of individuals, and the complexity of shipping dynamics impedes predictions of invasion dynamics. Here, we present a mathematical model for the spread of planktonic organisms by global shipping, using the history of movements by 33 566 ships among 1477 ports to simulate population dynamics for the comb jelly Mnemiopsis leidyi as a case study. The degree of propagule pressure at one site resulted from the coincident arrival of individuals from other sites with native or non-native populations. Key to sequential spread in European waters was a readily available source of propagules and a suitable recipient environment. These propagules were derived from previously introduced 'bridgehead' populations supplemented with those from native sources. Invasion success is therefore determined by the complex interaction of global shipping and local population dynamics. The general findings probably hold true for the spread of species in other complex systems, such as insects or plant seeds exchanged via commercial trade or transport.

Population attenuation in zooplankton communities during transoceanic transfer in ballast water.

Successful biological invasion requires introduction of a viable population of a nonindigenous species (NIS). Rarely have ecologists assessed changes in populations while entrained in invasion pathways. Here, we investigate how zooplankton communities resident in ballast water change during transoceanic voyages. We used next-generation sequencing technology to sequence a nuclear small subunit ribosomal DNA fragment of zooplankton from ballast water during initial, middle, and final segments as a vessel transited between Canada and Brazil. Operational taxonomic unit (OTU) diversity decreased as voyage duration increased, indicating loss of community-based genetic diversity and development of bottlenecks for zooplankton taxa prior to discharge of ballast water. On average, we observed 47, 26, and 24 OTUs in initial, middle, and final samples, respectively. Moreover, a comparison of genetic diversity within taxa indicated likely attenuation of OTUs in final relative to initial samples. Abundance of the most common taxa (copepods) declined in all final relative to initial samples. Some taxa (e.g., Copepoda) were represented by a high number of OTUs throughout the voyage, and thus had a high level of intraspecific genetic variation. It is not clear whether genotypes that were most successful in surviving transit in ballast water will be the most successful upon introduction to novel environments. This study highlights that population bottlenecks may be common prior to introduction of NIS to new ecosystems.

Are genetic databases sufficiently populated to detect non-indigenous species?

Correct species identifications are of tremendous importance for invasion ecology, as mistakes could lead to misdirecting limited resources against harmless species or inaction against problematic ones. DNA barcoding is becoming a promising and reliable tool for species identifications, however the efficacy of such molecular taxonomy depends on gene region(s) that provide a unique sequence to differentiate among species and on availability of reference sequences in existing genetic databases. Here, we assembled a list of aquatic and terrestrial non-indigenous species (NIS) and checked two leading genetic databases for corresponding sequences of six genome regions used for DNA barcoding. The genetic databases were checked in 2010, 2012, and 2016. All four aquatic kingdoms (Animalia, Chromista, Plantae and Protozoa) were initially equally represented in the genetic databases, with 64, 65, 69, and 61 % of NIS included, respectively. Sequences for terrestrial NIS were present at rates of 58 and 78 % for Animalia and Plantae, respectively. Six years later, the number of sequences for aquatic NIS increased to 75, 75, 74, and 63 % respectively, while those for terrestrial NIS increased to 74 and 88 % respectively. Genetic databases are marginally better populated with sequences of terrestrial NIS of plants compared to aquatic NIS and terrestrial NIS of animals. The rate at which sequences are added to databases is not equal among taxa. Though some groups of NIS are not detectable at all based on available data-mostly aquatic ones-encouragingly, current availability of sequences of taxa with environmental and/or economic impact is relatively good and continues to increase with time.

Invasion pathway of the Ctenophore Mnemiopsis leidyi in the Mediterranean Sea.

Gelatinous zooplankton outbreaks have increased globally owing to a number of human-mediated factors, including food web alterations and species introductions. The invasive ctenophore Mnemiopsis leidyi entered the Black Sea in the early 1980s. The invasion was followed by the Azov, Caspian, Baltic and North Seas, and, most recently, the Mediterranean Sea. Previous studies identified two distinct invasion pathways of M. leidyi from its native range in the western Atlantic Ocean to Eurasia. However, the source of newly established populations in the Mediterranean Sea remains unclear. Here we build upon our previous study and investigate sequence variation in both mitochondrial (Cytochrome c Oxidase subunit I) and nuclear (Internal Transcribed Spacer) markers in M. leidyi, encompassing five native and 11 introduced populations, including four from the Mediterranean Sea. Extant genetic diversity in Mediterranean populations (n = 8, N a = 10) preclude the occurrence of a severe genetic bottleneck or founder effects in the initial colonizing population. Our mitochondrial and nuclear marker surveys revealed two possible pathways of introduction into Mediterranean Sea. In total, 17 haplotypes and 18 alleles were recovered from all surveyed populations. Haplotype and allelic diversity of Mediterranean populations were comparable to populations from which they were likely drawn. The distribution of genetic diversity and pattern of genetic differentiation suggest initial colonization of the Mediterranean from the Black-Azov Seas (pairwise F ST = 0.001-0.028). However, some haplotypes and alleles from the Mediterranean Sea were not detected from the well-sampled Black Sea, although they were found in Gulf of Mexico populations that were also genetically similar to those in the Mediterranean Sea (pairwise F ST = 0.010-0.032), raising the possibility of multiple invasion sources. Multiple introductions from a combination of Black Sea and native region sources could be facilitated by intense local and transcontinental shipping activity, respectively.

Genetic diversity in introduced golden mussel populations corresponds to vector activity.

We explored possible links between vector activity and genetic diversity in introduced populations of Limnoperna fortunei by characterizing the genetic structure in native and introduced ranges in Asia and South America. We surveyed 24 populations: ten in Asia and 14 in South America using the mitochondrial cytochrome c oxidase subunit I (COI) gene, as well as eight polymorphic microsatellite markers. We performed population genetics and phylogenetic analyses to investigate population genetic structure across native and introduced regions. Introduced populations in Asia exhibit higher genetic diversity (H(E) = 0.667-0.746) than those in South America (H(E) =  0.519-0.575), suggesting higher introduction effort for the former populations. We observed pronounced geographical structuring in introduced regions, as indicated by both mitochondrial and nuclear markers based on multiple genetic analyses including pairwise Ф(ST), F(ST), bayesian clustering method, and three-dimensional factorial correspondence analyses. Pairwise F(ST) values within both Asia (F(ST) = 0.017-0.126, P = 0.000-0.009) and South America (F(ST) =0.004-0.107, P = 0.000-0.721) were lower than those between continents (F(ST) = 0.180-0.319, P = 0.000). Fine-scale genetic structuring was also apparent among introduced populations in both Asia and South America, suggesting either multiple introductions of distinct propagules or strong post-introduction selection and demographic stochasticity. Higher genetic diversity in Asia as compared to South America is likely due to more frequent propagule transfers associated with higher shipping activities between source and donor regions within Asia. This study suggests that the intensity of human-mediated introduction vectors influences patterns of genetic diversity in non-indigenous species.