Genetic reconstructions of invasion history

Complex European invasion history of Anoplophora glabripennis (Motschulsky): new insights in its population genomic differentiation using genotype-by-sequencing

Anthropogenic activities like trade facilitate increasing rates of biological invasions. Asian long-horned beetle (ALB), which is naturally distributed in eastern Asia (China, Korean peninsula), was introduced via wood packing materials (WPM) used in trade to North America (1996) and Europe (2001). We used 7810 single nucleotide polymorphisms (SNPs) derived by a genotype-by-sequencing (GBS) approach to decipher the introduction patterns into Europe. This is applied for the first time on European ALB outbreaks from Germany, Switzerland, and Italy, both from still active and already eradicated infestations. The genome-wide SNPs detected signs of small and highly structured populations within Europe, showing clear founder effects. The very high population differentiation is presumably derived from multiple independent introductions to Europe, which are spatially restricted in mating. By admixture and phylogenetic analyses, some cases of secondary dispersal were observed. Furthermore, some populations suggest admixture, which might have been originated by either multiple introductions from different sources into the new sites or recurrent introductions from an admixed source population. Our results confirmed a complex invasion history of the ALB into Europe and the usability of GBS obtained SNPs in invasion science even without source populations.

The natural and human-mediated expansion of a human-commensal lizard into the fringes of Southeast Asia

BACKGROUND

Human-commensal species often display deep ancestral genetic structure within their native range and founder-effects and/or evidence of multiple introductions and admixture in newly established areas. We investigated the phylogeography of Eutropis multifasciata, an abundant human-commensal scincid lizard that occurs across Southeast Asia, to determine the extent of its native range and to assess the sources and signatures of human introduction outside of the native range. We sequenced over 350 samples of E. multifasciata for the mitochondrial ND2 gene and reanalyzed a previous RADseq population genetic dataset in a phylogenetic framework.

RESULTS

Nuclear and mitochondrial trees are concordant and show that E. multifasciata has retained high levels of genetic structure across Southeast Asia despite being frequently moved by humans. Lineage boundaries in the native range roughly correspond to several major biogeographic barriers, including Wallace's Line and the Isthmus of Kra. Islands at the outer fringe of the range show evidence of founder-effects and multiple introductions.

CONCLUSIONS

Most of enormous range of E. multifasciata across Southeast Asia is native and it only displays signs of human-introduction or recent expansion along the eastern and northern fringe of its range. There were at least three events of human-introductions to Taiwan and offshore islands, and several oceanic islands in eastern Indonesia show a similar pattern. In Myanmar and Hainan, there is a founder-effect consistent with post-warming expansion after the last glacial maxima or human introduction.

Genomic Tools in Biological Invasions: Current State and Future Frontiers

Human activities are accelerating rates of biological invasions and climate-driven range expansions globally, yet we understand little of how genomic processes facilitate the invasion process. Although most of the literature has focused on underlying phenotypic correlates of invasiveness, advances in genomic technologies are showing a strong link between genomic variation and invasion success. Here, we consider the ability of genomic tools and technologies to (i) inform mechanistic understanding of biological invasions and (ii) solve real-world issues in predicting and managing biological invasions. For both, we examine the current state of the field and discuss how genomics can be leveraged in the future. In addition, we make recommendations pertinent to broader research issues, such as data sovereignty, metadata standards, collaboration, and science communication best practices that will require concerted efforts from the global invasion genomics community.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

Stronger evidence for genetic ancestry than environmental conditions in shaping the evolution of a complex signalling trait during biological invasion

Introductions of invasive species to new environments often result in rapid rates of trait evolution. While in some cases these evolutionary transitions are adaptive and driven by natural selection, they can also result from patterns of genetic and phenotypic variation associated with the invasion history. Here, we examined the brown anole (Anolis sagrei), a widespread invasive lizard for which genetic data have helped trace the sources of non-native populations. We focused on the dewlap, a complex signalling trait known to be subject to multiple selective pressures. We measured dewlap reflectance, pattern and size in 30 non-native populations across the southeastern United States. As well, we quantified environmental variables known to influence dewlap signal effectiveness, such as canopy openness. Further, we used genome-wide data to estimate genetic ancestry, perform association mapping and test for signatures of selection. We found that among-population variation in dewlap characteristics was best explained by genetic ancestry. This result was supported by genome-wide association mapping, which identified several ancestry-specific loci associated with dewlap traits. Despite the strong imprint of this aspect of the invasion history on dewlap variation, we also detected significant relationships between dewlap traits and local environmental conditions. However, we found limited evidence that dewlap-associated genetic variants have been subject to selection. Our study emphasizes the importance of genetic ancestry and admixture in shaping phenotypes during biological invasion, while leaving the role of selection unresolved, likely due to the polygenic genetic architecture of dewlaps and selection acting on many genes of small effect.

The Invasion History of Diaphorina citri (Hemiptera: Liviidae) in Puerto Rico: Past, Present, and Future Perspectives

The Asian citrus psyllid, Diaphorina citri Kuwayama, 1908 (Hemiptera: Liviidae), is a phloem-feeding insect pest of Citrus L. ( Sapindales: Rutaceae), Murraya J. Koenig ex L. (Sapindales: Rutaceae), and other Rutaceae, native to the Indian subcontinent and transmits the citrus greening pathogens (huanglongbing). We herein evaluate the occurrence of D. citri in Puerto Rico at a fine scale and the potential invasion history of the Caribbean from the native range. We found that the genetic diversity of D. citri in Puerto Rico, using the mitochondrial genetic marker CO1, indicated that the Caribbean haplotype is unique and more related to Southwest Asia haplotypes than other groups. Results then suggest that the invasion of the Caribbean did not occur from the United States. Also, we found that isothermality was the abiotic variable that influenced D. citri species distribution model performance the most in current and future climate change scenarios. Given the data available from the Global Biodiversity Information Facility (GBIF) for present D.citri reported occurrence, future projections suggest that climate change might promote an expansion of the pest to other Caribbean islands and Latin American countries. Still, it might elicit a range contraction across the Cordillera Central in Puerto Rico, possibly promoting citriculture. Regarding pest management strategies, we must consider a holistic approach to understanding D. citri invasion and resilience in novel habitats by including broader and more robust genetic population structure studies and if its current distribution reflects projected models across the neotropics.

Genetic diversity and population structure among native, naturalized, and invasive populations of the common yellow monkeyflower, Mimulus guttatus (Phrymaceae)

An ongoing controversy in invasion biology is the prevalence of colonizing plant populations that are able to establish and spread, while maintaining limited amounts of genetic variation. Invasive populations can be established through several routes including from a single source or from multiple introductions. The aim of this study was to examine genetic diversity in populations of Mimulus guttatus in the United Kingdom, where the species is considered invasive, and compare this diversity to that in native populations on the west coast of North America. Additionally, we looked at diversity in non-native populations that have not yet become invasive (naturalized populations) in eastern North America. We investigated population structure among populations in these three regions and attempted to uncover the sources for populations that have established in the naturalized and invasive regions. We found that genetic diversity was, on average, relatively high in populations from the invasive UK region and comparable to native populations. Contrastingly, two naturalized M. guttatus populations were low in both genetic and genotypic diversity, indicating a history of asexual reproduction and self-fertilization. A third naturalized population was found to be a polyploid Mimulus hybrid of unknown origin. Our results demonstrate that M. guttatus has likely achieved colonization success outside of its native western North America distribution by a variety of establishment pathways, including those with genetic and demographic benefits resulting from multiple introductions in the UK, reproductive assurance through selfing, and asexual reproduction in eastern North America, and possible polyploidization in one Canadian population.

Revisiting the historical scenario of a disease dissemination using genetic data and Approximate Bayesian Computation methodology: The case of Pseudocercospora fijiensis invasion in Africa

The reconstruction of geographic and demographic scenarios of dissemination for invasive pathogens of crops is a key step toward improving the management of emerging infectious diseases. Nowadays, the reconstruction of biological invasions typically uses the information of both genetic and historical information to test for different hypotheses of colonization. The Approximate Bayesian Computation framework and its recent Random Forest development (ABC-RF) have been successfully used in evolutionary biology to decipher multiple histories of biological invasions. Yet, for some organisms, typically plant pathogens, historical data may not be reliable notably because of the difficulty to identify the organism and the delay between the introduction and the first mention. We investigated the history of the invasion of Africa by the fungal pathogen of banana Pseudocercospora fijiensis, by testing the historical hypothesis against other plausible hypotheses. We analyzed the genetic structure of eight populations from six eastern and western African countries, using 20 microsatellite markers and tested competing scenarios of population foundation using the ABC-RF methodology. We do find evidence for an invasion front consistent with the historical hypothesis, but also for the existence of another front never mentioned in historical records. We question the historical introduction point of the disease on the continent. Crucially, our results illustrate that even if ABC-RF inferences may sometimes fail to infer a single, well-supported scenario of invasion, they can be helpful in rejecting unlikely scenarios, which can prove much useful to shed light on disease dissemination routes.

Multiple introductions, polyploidy and mixed reproductive strategies are linked to genetic diversity and structure in the most widespread invasive plant across Southern Ocean archipelagos

Biological invasions in remote areas that experience low human activity provide unique opportunities to elucidate processes responsible for invasion success. Here we study the most widespread invasive plant species across the isolated islands of the Southern Ocean, the annual bluegrass, Poa annua. To analyse geographical variation in genome size, genetic diversity and reproductive strategies, we sampled all major sub-Antarctic archipelagos in this region and generated microsatellite data for 470 individual plants representing 31 populations. We also estimated genome sizes for a subset of individuals using flow cytometry. Occasional events of island colonization are expected to result in high genetic structure among islands, overall low genetic diversity and increased self-fertilization, but we show that this is not the case for P. annua. Microsatellite data indicated low population genetic structure and lack of isolation by distance among the sub-Antarctic archipelagos we sampled, but high population structure within each archipelago. We identified high levels of genetic diversity, low clonality and low selfing rates in sub-Antarctic P. annua populations (contrary to rates typical of continental populations). In turn, estimates of selfing declined in populations as genetic diversity increased. Additionally, we found that most P. annua individuals are probably tetraploid and that only slight variation exists in genome size across the Southern Ocean. Our findings suggest multiple independent introductions of P. annua into the sub-Antarctic, which promoted the establishment of genetically diverse populations. Despite multiple introductions, the adoption of convergent reproductive strategies (outcrossing) happened independently in each major archipelago. The combination of polyploidy and a mixed reproductive strategy probably benefited P. annua in the Southern Ocean by increasing genetic diversity and its ability to cope with the novel environmental conditions.

There and Back Again: The Unexpected Journeys of Metridium de Blainville, 1824 between the Old Oceans and throughout the Modern World

AbstractMembers of the sea anemone genus Metridium are abundant in temperate rocky habitats and fouling communities. Their biogeographic history is expected to reflect changes in currents and habitats that have influenced benthic communities, such as the climate-influenced changes that occurred during the Last Glacial Maximum. More recently, however, anthropogenic influences such as shipping transportation and the creation of artificial habitat have altered and affected the composition of modern-day marine communities. Here we use sequence-capture data to examine the genetic structure of Metridium across its shallow-water distribution to (1) evaluate species boundaries within Metridium, (2) elucidate the dispersal history of Metridium between and among oceans, and (3) assess the influence of anthropogenic movement on modern-day populations. We find strong evidence for two species within Metridium: M. farcimen and M. senile. Dispersal from the Pacific to the Atlantic included a subsequent isolation of a small population in or above the Bering Sea, which has presumably moved southward. Within the native range of M. senile, admixture is prevalent even between oceans as a result of anthropogenic activities. The nonnative populations in Chile and the Falkland Islands came from at least two distinct introduction events originating from both coasts of the United States in the North Pacific and North Atlantic Oceans. Hybridization between M. senile and M. farcimen is documented as occurring in anthropogenically influenced habitats. The heavy influence from anthropogenic activities will continue to impact our understanding of marine organisms, particularly within the native range and for those that are easily transported across long distances.

Genetic Signature of Pinctada fucata Inferred from Population Genomics: Source Tracking of the Invasion in Mischief Reef of Nansha Islands

Among the anthropogenic stresses that marine ecosystems face, biological invasions are one of the major threats. Recently, as a result of increasingly intense anthropogenic disturbance, numerous marine species have been introduced to their non-native ranges. However, many introduced species have uncertain original sources. This prevents the design and establishment of methods for controlling or preventing these introduced species. In the present study, genomic sequencing and population genetic analysis were performed to detect the geographic origin of the introduced Pinctada fucata population in the Mischief Reef of the South China Sea. The results of population genetic structure analysis showed a close relationship between the Mischief Reef introduced population and the Lingshui population, indicating that Lingshui may be the potential geographical origin. Furthermore, lower heterozygosity and nucleotide diversity were observed in the introduced population in Mischief Reef, indicating lower genetic diversity than in other native populations. We also identified some selected genomic regions and genes of the introduced population, including genes related to temperature and salinity tolerance. These genes may play important roles in the adaptation of the introduced population. Our study will improve our understanding of the invasion history of the P. fucata population. Furthermore, the results of the present study will also facilitate further control and prevention of invasion in Mischief Reef, South China Sea.

Advances in understanding and predicting the spread of invading insect populations

Understanding and predicting the spread of invading insects is a critical challenge in management programs that aim to minimize ecological and economic harm to native ecosystems. Although efforts to quantify spread rates have been well studied over the past several decades, opportunities to improve our ability to estimate rates of spread, and identify the factors, such as habitat suitability and climate, that influence spread, remain. We review emerging sources of data that can be used to delineate distributional boundaries through time and thus serve as a basis for quantifying spread rates. We then address advances in modeling methods that facilitate our understanding of factors that drive invasive insect spread. We conclude by highlighting some remaining challenges in understanding and predicting invasive insect spread, such as the role of climate change and biotic similarity between the native and introduced ranges, particularly as it applies to decision-making in management programs.

Genetics reveals shifts in reproductive behaviour of the invasive bird parasite Philornis downsi collected from Darwin’s finch nests

Due to novel or dynamic fluctuations in environmental conditions and resources, host and parasite relationships can be subject to diverse selection pressures that may lead to significant changes during and after invasion of a parasite. Genomic analyses are useful for elucidating evolutionary processes in invasive parasites following their arrival to a new area and host. Philornis downsi (Diptera: Muscidae), the avian vampire fly, was introduced to the Galápagos Islands circa 1964 and has since spread across the archipelago, feeding on the blood of developing nestlings of endemic land birds. Since its discovery, there have been significant changes to the dynamics of P. downsi and its novel hosts, such as shifting mortality rates and changing oviposition behaviour, however no temporal genetic studies have been conducted. We collected P. downsi from nests and traps from a single island population over a 14-year period, and genotyped flies at 469 single nucleotide polymorphisms (SNPs) using restriction-site associated DNA sequencing (RADSeq). Despite significant genetic differentiation (FST) between years, there was no evidence for genetic clustering within or across four sampling years between 2006 and 2020, suggesting a lack of population isolation. Sibship reconstructions from P. downsi collected from 10 Darwin’s finch nests sampled in 2020 showed evidence for shifts in reproductive behaviour compared to a similar genetic analysis conducted in 2004–2006. Compared with this previous study, females mated with fewer males, individual females oviposited fewer offspring per nest, but more unique females oviposited per nest. These findings are important to consider within reproductive control techniques, and have fitness implications for both parasite evolution and host fitness.

Population structure and pattern of geographic differentiation of Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae) in Turkey

BACKGROUND

The Colorado potato beetle (CPB) is the most harmful pest of potato in potato cultivation regions globally. Although it is an economically important agricultural pest, the population structure and colonization route of this species in Turkey are uncertain. We used microsatellite and mitochondrial DNA (mtDNA) markers to obtain information about the population source, structure and bio-invasion route of CPB populations in Turkey.

RESULTS

The common single mtDNA haplotype in European CPB populations was obtained in all Turkish CPB populations based on mtDNA data analysis. However, microsatellites revealed a low level of genetic variation in CPB populations. The results of microsatellite analysis [factorial correspondence analysis (FCA), Bayesian analysis of genetic population structure (BAPS), unweighted pair group method with arithmetic mean (UPGMA) dendrogram, F-statistics and Nei's distances] indicated three groups for invasive CPB: Thrace-Marmara and Aegean; Black Sea, Central Anatolia and Mediterranean; Northeastern Anatolia. Region-specific alleles have been identified in regions, where commercial potato cultivation and insecticide use are intensive.

CONCLUSION

The detection of a single fixed European haplotype in all Turkish populations has proved that CPB in Turkey originated from Europe as a result of a founder event occurred in European populations. Low genetic variation was due to the short time period since the spread of CPB from America to Europe. The highest number of private alleles were found in the top commercial potato cultivation region-Central Anatolia from where the CPB populations spread to other parts of Turkey. © 2021 Society of Chemical Industry.

The fast invasion of Europe by the box tree moth: an additional example coupling multiple introduction events, bridgehead effects and admixture events

Identifying the invasion routes of non-native species is crucial to understanding invasions and customizing management strategies. The box tree moth, Cydalima perspectalis, is native to Asia and was recently accidentally introduced into Europe as a result of the ornamental plant trade. Over the last 15 years, it has spread across the continent and has reached the Caucasus and Iran. It is threatening Buxus trees in both urban areas and forests. To investigate the species’ invasion routes, native and invasive box tree moth populations were sampled, and moth’s genetic diversity and structure were compared using microsatellite markers. Our approximate Bayesian computation analyses strongly suggest that invasion pathways were complex. Primary introductions originating from eastern China probably occurred independently twice in Germany and once in the Netherlands. There were also possibly bridgehead effects, where at least three invasive populations may have served as sources for other invasive populations within Europe, with indication of admixture between the two primary invasive populations. The bridgehead populations were likely those in the countries that play a major role in the ornamental plant trade in Europe, notably Germany, the Netherlands, and Italy. All these invasion processes likely facilitated its fast expansion across Europe and illustrate the role played by the ornamental plant trade not only in the moth’s introduction from China but also in the species’ spread across Europe, leading to an invasion with a complex pattern.

Genomics reveals the history of a complex plant invasion and improves the management of a biological invasion from the South African-Australian biotic exchange

Many plants exchanged in the global redistribution of species in the last 200 years, particularly between South Africa and Australia, have become threatening invasive species in their introduced range. Refining our understanding of the genetic diversity and population structure of native and alien populations, introduction pathways, propagule pressure, naturalization, and initial spread, can transform the effectiveness of management and prevention of further introductions. We used 20,221 single nucleotide polymorphisms to reconstruct the invasion of a coastal shrub, Chrysanthemoides monilifera ssp. rotundata (bitou bush) from South Africa, into eastern Australia (EAU), and Western Australia (WAU). We determined genetic diversity and population structure across the native and introduced ranges and compared hypothesized invasion scenarios using Bayesian modeling. We detected considerable genetic structure in the native range, as well as differentiation between populations in the native and introduced range. Phylogenetic analysis showed the introduced samples to be most closely related to the southern‐most native populations, although Bayesian analysis inferred introduction from a ghost population. We detected strong genetic bottlenecks during the founding of both the EAU and WAU populations. It is likely that the WAU population was introduced from EAU, possibly involving an unsampled ghost population. The number of private alleles and polymorphic SNPs successively decreased from South Africa to EAU to WAU, although heterozygosity remained high. That bitou bush remains an invasion threat in EAU, despite reduced genetic diversity, provides a cautionary biosecurity message regarding the risk of introduction of potentially invasive species via shipping routes.

Next generation biosecurity: Towards genome based identification to prevent spread of agronomic pests and pathogens using nanopore sequencing

The unintentional movement of agronomic pests and pathogens is steadily increasing due to the intensification of global trade. Being able to identify accurately and rapidly early stages of an invasion is critical for developing successful eradication or management strategies. For most invasive organisms, molecular diagnostics is today the method of choice for species identification. However, the currently implemented tools are often developed for certain taxa and need to be adapted for new species, making them ill-suited to cope with the current constant increase in new invasive species. To alleviate this impediment, we developed a fast and accurate sequencing tool allowing to modularly obtain genetic information at different taxonomical levels. Using whole genome amplification (WGA) followed by Oxford nanopore MinION sequencing, our workflow does not require any a priori knowledge on the investigated species and its classification. While mainly focusing on harmful plant pathogenic insects, we also demonstrate the suitability of our workflow for the molecular identification of bacteria (Erwinia amylovora and Escherichia coli), fungi (Cladosporium herbarum, Colletotrichum salicis, Neofabraea alba) and nematodes (Globodera rostochiensis). On average, the pairwise identity between the generated consensus sequences and best GenBank BLAST matches was 99.6 ± 0.6%. Additionally, assessing the generated insect genomic dataset, the potential power of the workflow to detect pesticide resistance genes, as well as arthropod-infecting viruses and endosymbiotic bacteria is demonstrated.

Multiple introductions and overwintering shape the progressive invasion of Aedes albopictus beyond the Alps

Aedes albopictus originates from Southeast Asia and is considered one of the most invasive species globally. This mosquito is a nuisance and a disease vector of significant public health relevance. In Europe, Ae. albopictus is firmly established and widespread south of the Alps, a mountain range that forms a formidable biogeographic barrier to many organisms. Recent reports of Ae. albopictus north of the Alps raise questions of (1) the origins of its recent invasion, and (2) if this mosquito has established overwintering populations north of the Alps. To answer these questions, we analyzed population genomic data from >4000 genome‐wide SNPs obtained through double‐digest restriction site‐associated DNA sequencing. We collected SNP data from specimens from six sites in Switzerland, north and south of the Alps, and analyzed them together with specimens from other 33 European sites, five from the Americas, and five from its Asian native range. At a global level, we detected four genetic clusters with specimens from Indonesia, Brazil, and Japan as the most differentiated, whereas specimens from Europe, Hong Kong, and USA largely overlapped. Across the Alps, we detected a weak genetic structure and high levels of genetic admixture, supporting a scenario of rapid and human‐aided dispersal along transportation routes. While the genetic pattern suggests frequent re‐introductions into Switzerland from Italian sources, the recovery of a pair of full siblings in two consecutive years in Strasbourg, France, suggests the presence of an overwintering population north of the Alps. The suggestion of overwintering populations of Ae. albopictus north of the Alps and the expansion patterns identified points to an increased risk of further northward expansion and the need for increased surveillance of mosquito populations in Northern Europe.

Genome-scale phylogeography resolves the native population structure of the Asian longhorned beetle, Anoplophora glabripennis (Motschulsky)

Human-assisted movement has allowed the Asian longhorned beetle (ALB, Anoplophora glabripennis (Motschulsky)) to spread beyond its native range and become a globally regulated invasive pest. Within its native range of China and the Korean peninsula, human-mediated dispersal has also caused cryptic translocation of insects, resulting in population structure complexity. Previous studies used genetic methods to detangle this complexity but were unable to clearly delimit native populations which is needed to develop downstream biosurveillance tools. We used genome-wide markers to define historical population structure in native ALB populations and contemporary movement between regions. We used genotyping-by-sequencing to generate 6102 single-nucleotide polymorphisms (SNPs) and amplicon sequencing to genotype 53 microsatellites. In total, we genotyped 712 individuals from ALB's native distribution. We observed six distinct population clusters among native ALB populations, with a clear delineation between northern and southern groups. Most of the individuals from South Korea were distinct from populations in China. Our results also indicate historical divergence among populations and suggest limited large-scale admixture, but we did identify a restricted number of cases of contemporary movement between regions. We identified SNPs under selection and describe a clinal allele frequency pattern in a missense variant associated with glycerol kinase, an important enzyme in the utilization of an insect cryoprotectant. We further demonstrate that small numbers of SNPs can assign individuals to geographic regions with high probability, paving the way for novel ALB biosurveillance tools.

Endemism, invasion, and overseas dispersal: the phylogeographic history of the Lesser Antillean frog, Eleutherodactylus johnstonei

Cryptogenic species are those whose native and introduced ranges are unknown. The extent and long history of human migration rendered numerous species cryptogenic. Incomplete knowledge regarding the origin and native habitat of a species poses problems for conservation management and may confound ecological and evolutionary studies. The Lesser Antilles pose a particular challenge with regard to cryptogenic species because these islands have been anthropogenically connected since before recorded history. Here, we use population genetic and phylogeographic tools in an attempt to determine the origin of Eleutherodactylus johnstonei, a frog species with a potentially widespread introduced range and whose native range within the Lesser Antilles is unknown. Based on elevated estimates of genetic diversity and within-island geographic structure not present elsewhere in the range, we identify Montserrat as the native island of E. johnstonei. We also document two major clades within E. johnstonei, only one of which is the primary source of introduced populations throughout the Americas. Our results demonstrate the utility of genetic tools for resolving cryptogenic species problems and highlight E. johnstonei as a potential system for understanding differences in invasive potential among sister lineages.

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)’.

Complex genetic patterns and distribution limits mediated by native congeners of the worldwide invasive red-eared slider turtle

Non-native (invasive) species offer a unique opportunity to study the geographic distribution and range limits of species, wherein the evolutionary change driven by interspecific interactions between native and non-native closely related species is a key component. The red-eared slider turtle, Trachemys scripta elegans (TSE), has been introduced and successfully established worldwide. It can coexist with its native congeners T. cataspila, T. venusta and T. taylori in Mexico. We performed comprehensive fieldwork, executed a battery of genetic analyses and applied a novel species distribution modeling approach to evaluate their historical lineage relationships and contemporary population genetic patterns. Our findings support the historical common ancestry between native TSE and non-native (TSE_alien ), while also highlighting the genetic differentiation of the exotic lineage. Genetic patterns are associated with their range size/endemism gradient, the microendemic T. taylori showed significant reduced genetic diversity and high differentiation, whereas TSE_alien showed the highest diversity and signals of population size expansion. Counter to our expectations, lower naturally occurring distribution overlap and little admixture patterns were found between TSE and its congeners, exhibiting reduced gene flow and clear genetic separation across neighboring species despite having zones of contact. We demonstrate that these native Trachemys species have distinct climatic niche suitability, likely preventing establishment of and displacement by the TSE_alien . Moreover, we found major niche overlap between TSE_alien and native species worldwide, supporting our prediction that sites with closer ecological optima to the invasive species have higher establishment risk than those that are closer to the niche-center of the native species.

Global high-throughput genotyping of organellar genomes reveals insights into the origin and spread of invasive starry stonewort (Nitellopsis obtusa)

Aquatic invasive species are damaging to native ecosystems. Preventing their spread and achieving comprehensive control measures requires an understanding of the genetic structure of an invasive population. Organellar genomes (plastid and mitochondrial) are useful for population level analyses of invasive plant distributions. In this study we generate complete organellar reference genomes using PacBio sequencing, then use these reference sequences for SNP calling of high-throughput, multiplexed, Illumina based organellar sequencing of fresh and historical samples from across the native and invasive range of Nitellopsis obtusa (Desv. in Loisel.) J.Groves, an invasive macroalgae. The data generated by the analytical pipeline we develop indicate introduction to North America from Western Europe. A single nucleotide transversion in the plastid genome separates a group of five samples from Michigan and Wisconsin that either resulted from introductions of two closely related genotypes or a mutation that has arisen in the invasive range. This transversion will serve as a useful tool to understand how Nitellopsis obtusa moves across the landscape. The methods and analyses described here are broadly applicable to invasive and native plant and algae species, and allow efficient genotyping of variable quality samples, including 100-year-old herbarium specimens, to determine population structure and geographic distributions.

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.

Allopatric differentiation in the Enteromius anoplus complex in South Africa, with the revalidation of Enteromius cernuus and Enteromius oraniensis, and description of a new species, Enteromius mandelai (Teleostei: Cyprinidae)

The chubbyhead barb, Enteromius anoplus, as currently described, is the most widely distributed freshwater fish in South Africa. The species occurs in almost all the major river systems across the country, with the exception of the small coastal drainages on the south coast. The use of a comprehensive data set of mitochondrial (mtDNA) cytochrome b (cyt b) sequences uncovered the presence of four distinct lineages or operational taxonomic units (OTUs) within E. anoplus: (a) the Gouritz lineage endemic to the Gouritz River system, (b) the Olifants lineage endemic to the Olifants-Doring River system, (c) the Orange lineage endemic to the Orange River system and (d) the Eastern Cape lineage which is distributed across six river systems from the Great Fish to the Buffalo rivers in the Eastern Cape Province of South Africa. The present study provides a new diagnosis for E. anoplus s.s. (the Gouritz lineage), revalidates Enteromius cernuus (the Olifants lineage) and Enteromius oraniensis (the Orange lineage) as distinct species. The study also provides a description for a new species, Enteromius mandelai sp. nov (the Eastern Cape lineage). E. cernuus and E. mandelai differ from both E. anoplus and E. oraniensis by having a complete lateral line (vs. an incomplete lateral line in the latter species). E. cernuus further differs from the other three species by having long maxillary barbels which reach or exceed the vertical through the middle of the eye and the lowest number of circumpendicular scales (10-12 vs. 12-16 for the other three species). E. oraniensis is distinctive from the other three species by having inconspicuous barbels. These findings add to the growing body of literature that shows that the freshwater fishes of the Cape Fold and adjacent freshwater ecoregions have narrow geographic ranges. This has ramifications for conservation planning and management, as well as the understanding of the evolutionary history of the stream fishes in these global endemic hotspots.

Revealing the Introduction History and Phylogenetic Relationships of Passiflora foetida sensu lato in Australia

Genomic analysis can be a valuable tool to assistmanagement of non-native invasive species, through determining source and number of introductions as well as clarifying phylogenetic relationships. Here, we used whole chloroplast sequencing to investigate the introduction history of Passiflora foetida sensu lato in Australia and clarify its relationship with other Passiflora species present. Phylogenetic analysis of chloroplast genome data identified three separate genetic lineages of P. foetida s. l. present in Australia, indicating multiple introductions. These lineages had affinities to samples from three separate areas within the native range in Central and South America that represented phylogenetically distinct lineages. These results provide a basis for a targeted search of the native range of P. foetida s. l. for candidate biological control agents that have co-evolved with this species and are thus better adapted to the lineages that are present in Australia. Results also indicated that the Passiflora species native to Australia are in a separate clade to that of P. foetida s. l. and other introduced Passiflora species cultivated in Australia. This knowledge is important to assess the likelihood of finding biological control agents for P. foetida s. l. that will be sufficiently host-specific for introduction in Australia. As P. foetida s. l. is a widespread non-native invasive species across many regions of the world, outcomes from this work highlight the importance of first evaluating the specific entities present in a country before the initiation of a biological control program.

Extensive human-mediated jump dispersal within and across the native and introduced ranges of the invasive termite Reticulitermes flavipes

As native ranges are often geographically structured, invasive species originating from a single source population only carry a fraction of the genetic diversity present in their native range. The invasion process is thus often associated with a drastic loss of genetic diversity resulting from a founder event. However, the fraction of diversity brought to the invasive range may vary under different invasion histories, increasing with the size of the propagule, the number of reintroduction events, and/or the total genetic diversity represented by the various source populations in a multiple-introduction scenario. In this study, we generated a SNP data set for the invasive termite Reticulitermes flavipes from 23 native populations in the eastern United States and six introduced populations throughout the world. Using population genetic analyses and approximate Bayesian computation random forest, we investigated its worldwide invasion history. We found a complex invasion pathway with multiple events out of the native range and bridgehead introductions from the introduced population in France. Our data suggest that extensive long-distance jump dispersal appears common in both the native and introduced ranges of this species, probably through human transportation. Overall, our results show that similar to multiple introduction events into the invasive range, admixture in the native range prior to invasion can potentially favour invasion success by increasing the genetic diversity that is later transferred to the introduced range.

Genomic data support multiple introductions and explosive demographic expansions in a highly invasive aquatic insect

The study of the genetic makeup and demographic fate of alien species is essential to understand their capacity to recover from founder effects, adapt to new environmental conditions and, ultimately, become invasive and potentially damaging. Here, we employ genomic data to gain insights into key demographic processes that might help to explain the extraordinarily successful invasion of the Western Mediterranean region by the North American boatman Trichocorixa verticalis (Hemiptera: Corixidae). Our analyses revealed the genetic distinctiveness of populations from the main areas comprising the invasive range and coalescent-based simulations supported that they originated from independent introductions events probably involving different source populations. Testing of alternative demographic models indicated that all populations experienced a strong bottleneck followed by a recent and instantaneous demographic expansion that restored a large portion (>30%) of their ancestral effective population sizes shortly after introductions took place (<60 years ago). Considerable genetic admixture of some populations suggest that hypothetical barriers to dispersal (i.e., land and sea water) are permeable to gene flow and/or that they originated from introductions involving multiple lineages. This study demonstrates the repeated arrival of propagules with different origins and short time lags between arrival and establishment, emphasizing the extraordinary capacity of the species to recover from founder effects and genetically admix in invaded areas. This can explain the demonstrated capacity of this aquatic insect to spread and outcompete native species once it colonizes new suitable regions. Future genomic analyses of native range populations could help to infer the genetic makeup of introduced populations and track invasion routes.

Human activity strongly influences genetic dynamics of the most widespread invasive plant in the sub-Antarctic

The link between the successful establishment of alien species and propagule pressure is well-documented. Less known is how humans influence the post-introduction dynamics of invasive alien populations. The latter requires studying parallel invasions by the same species in habitats that are differently impacted by humans. We analysed microsatellite and genome size variation, and then compared the genetic diversity and structure of invasive Poa annua L. on two sub-Antarctic islands: human-occupied Marion Island and unoccupied Prince Edward Island. We also carried out niche modelling to map the potential distribution of the species on both islands. We found high levels of genetic diversity and evidence for extensive admixture between genetically distinct lineages of P. annua on Marion Island. By contrast, the Prince Edward Island populations showed low genetic diversity, no apparent admixture, and had smaller genomes. On both islands, high genetic diversity was apparent at human landing sites, and on Marion Island, also around human settlements, suggesting that these areas received multiple introductions and/or acted as initial introduction sites and secondary sources (bridgeheads) for invasive populations. More than 70 years of continuous human activity associated with a meteorological station on Marion Island led to a distribution of this species around human settlements and along footpaths, which facilitates ongoing gene flow among geographically separated populations. By contrast, this was not the case for Prince Edward Island, where P. annua populations showed high genetic structure. The high levels of genetic variation and admixture in P. annua facilitated by human activity, coupled with high habitat suitability on both islands, suggest that P. annua is likely to increase its distribution and abundance in the future.

The evolutionary dynamics of biological invasions: A multi-approach perspective

Biological invasions, the establishment and spread of non-native species in new regions, can have extensive economic and environmental consequences. Increased global connectivity accelerates introduction rates, while climate and land-cover changes may decrease the barriers to invasive populations spread. A detailed knowledge of the invasion history, including assessing source populations, routes of spread, number of independent introductions, and the effects of genetic bottlenecks and admixture on the establishment success, adaptive potential, and further spread, is crucial from an applied perspective to mitigate socioeconomic impacts of invasive species, as well as for addressing fundamental questions on the evolutionary dynamics of the invasion process. Recent advances in genomics together with the development of geographic information systems provide unprecedented large genetic and environmental datasets at global and local scales to link population genomics, landscape ecology, and species distribution modeling into a common framework to study the invasion process. Although the factors underlying population invasiveness have been extensively reviewed, analytical methods currently available to optimally combine molecular and environmental data for inferring invasive population demographic parameters and predicting further spreading are still under development. In this review, we focus on the few recent insect invasion studies that combine different datasets and approaches to show how integrating genetic, observational, ecological, and environmental data pave the way to a more integrative biological invasion science. We provide guidelines to study the evolutionary dynamics of invasions at each step of the invasion process, and conclude on the benefits of including all types of information and up-to-date analytical tools from different research areas into a single framework.

Insights into invasive species from whole-genome resequencing

Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.

Effect of social structure and introduction history on genetic diversity and differentiation

Invasive species are a global threat to biodiversity, and understanding their history and biology is a major goal of invasion biology. Population-genetic approaches allow insights into these features, as population structure is shaped by factors such as invasion history (number, origin and age of introductions) and life-history traits (e.g., mating system, dispersal capability). We compared the relative importance of these factors by investigating two closely related ants, Tetramorium immigrans and Tetramorium tsushimae, that differ in their social structure and invasion history in North America. We used mitochondrial DNA sequences and microsatellite alleles to estimate the source and number of introduction events of the two species, and compared genetic structure among native and introduced populations. Genetic diversity of both species was strongly reduced in introduced populations, which also differed genetically from native populations. Genetic differentiation between ranges and the reduction in microsatellite diversity were more severe in the more recently introduced and supercolonial T. tsushimae. However, the loss of mitochondrial haplotype diversity was more pronounced in T. immigrans, which has single-queen colonies and was introduced earlier. Tetramorium immigrans was introduced at least twice from Western Europe to North America and once independently to South America. Its monogyny might have limited genetic diversity per introduction, but new mutations and successive introductions over a long time may have added to the gene pool in the introduced range. Polygyny in T. tsushimae probably facilitated the simultaneous introduction of several queens from a Japanese population to St. Louis, USA. In addition to identifying introduction pathways, our results reveal how social structure can influence the population-genetic consequences of founder events.

Founder effects shape linkage disequilibrium and genomic diversity of a partially clonal invader

The genomic variation of an invasive species may be affected by complex demographic histories and evolutionary changes during the invasion. Here, we describe the relative influence of bottlenecks, clonality, and population expansion in determining genomic variability of the widespread red macroalga Agarophyton vermiculophyllum. Its introduction from mainland Japan to the estuaries of North America and Europe coincided with shifts from predominantly sexual to partially clonal reproduction and rapid adaptive evolution. A survey of 62,285 SNPs for 351 individuals from 35 populations, aligned to 24 chromosome-length scaffolds indicate that linkage disequilibrium (LD), observed heterozygosity (H_o ), Tajima's D, and nucleotide diversity (Pi) were greater among non-native than native populations. Evolutionary simulations indicate LD and Tajima's D were consistent with a severe population bottleneck. Also, the increased rate of clonal reproduction in the non-native range could not have produced the observed patterns by itself but may have magnified the bottleneck effect on LD. Elevated marker diversity in the genetic source populations could have contributed to the increased H_o and Pi observed in the non-native range. We refined the previous invasion source region to a ~50 km section of northeastern Honshu Island. Outlier detection methods failed to reveal any consistently differentiated loci shared among invaded regions, probably because of the complex A. vermiculophyllum demographic history. Our results reinforce the importance of demographic history, specifically founder effects, in driving genomic variation of invasive populations, even when localized adaptive evolution and reproductive system shifts are observed.

Bridgehead effect and multiple introductions shape the global invasion history of a termite

Native to eastern Asia, the Formosan subterranean termite Coptotermes formosanus (Shiraki) is recognized as one of the 100 worst invasive pests in the world, with established populations in Japan, Hawaii and the southeastern United States. Despite its importance, the native source(s) of C. formosanus introductions and their invasive pathway out of Asia remain elusive. Using ~22,000 SNPs, we retraced the invasion history of this species through approximate Bayesian computation and assessed the consequences of the invasion on its genetic patterns and demography. We show a complex invasion history, where an initial introduction to Hawaii resulted from two distinct introduction events from eastern Asia and the Hong Kong region. The admixed Hawaiian population subsequently served as the source, through a bridgehead, for one introduction to the southeastern US. A separate introduction event from southcentral China subsequently occurred in Florida showing admixture with the first introduction. Overall, these findings further reinforce the pivotal role of bridgeheads in shaping species distributions in the Anthropocene and illustrate that the global distribution of C. formosanus has been shaped by multiple introductions out of China, which may have prevented and possibly reversed the loss of genetic diversity within its invasive range.

Using Genomics to Link Populations of an Invasive Species to Its Potential Sources

The introduction and subsequent range expansion of the Northern snakehead (Channa argus: Channidae, Anabantiformes) is one of a growing number of problematic biological invasions in the United States. This harmful aquatic invasive species is a predatory freshwater fish native to northeastern Asia that, following deliberate introduction, has established itself in multiple water basins in the eastern United States, as well as expanding its range into the Midwest. Previous work assessed the population structure and estimated the long-term effective population sizes of the populations present in the United States, but the source of the initial introduction(s) to the U.S. remains unidentified. Building on earlier work, we used whole genome scans (2b-RAD genomic sequencing) to analyze single nucleotide polymorphisms (SNPs) from C. argus to screen the genomes of these invasive fish from United States waters and from three sites in their native range in China. We recovered 2,822 SNP loci from genomic DNA extracted from 164 fish sampled from the eastern United States and Arkansas (Mississippi River basin), plus 30 fish sampled from three regions of the Yangtze River basin in China (n = 10 individuals per basin). Our results provide evidence supporting the Yangtze River basin in China, specifically the Bohu and/or Liangzi lakes, is a likely source of the C. argus introductions in multiple regions of the U.S., including the Lower Hudson River basin, Upper Hudson River basin and Philadelphia (Lower Delaware River basin). This information, in conjunction with additional sampling from the native range, will help to determine the source(s) of introduction for the other U.S. populations. Additionally, this work will provide valuable information for management to help prevent and manage future introductions into United States waterways, as well as aid in the development of more targeted strategies to regulate established populations and inhibit further spread.

Spatial population genomics of a recent mosquito invasion

Population genomic approaches can characterize dispersal across a single generation through to many generations in the past, bridging the gap between individual movement and intergenerational gene flow. These approaches are particularly useful when investigating dispersal in recently altered systems, where they provide a way of inferring long-distance dispersal between newly established populations and their interactions with existing populations. Human-mediated biological invasions represent such altered systems which can be investigated with appropriate study designs and analyses. Here we apply temporally restricted sampling and a range of population genomic approaches to investigate dispersal in a 2004 invasion of Aedes albopictus (the Asian tiger mosquito) in the Torres Strait Islands (TSI) of Australia. We sampled mosquitoes from 13 TSI villages simultaneously and genotyped 373 mosquitoes at genome-wide single nucleotide polymorphisms (SNPs): 331 from the TSI, 36 from Papua New Guinea (PNG) and four incursive mosquitoes detected in uninvaded regions. Within villages, spatial genetic structure varied substantially but overall displayed isolation by distance and a neighbourhood size of 232-577. Close kin dyads revealed recent movement between islands 31-203 km apart, and deep learning inferences showed incursive Ae. albopictus had travelled to uninvaded regions from both adjacent and nonadjacent islands. Private alleles and a co-ancestry matrix indicated direct gene flow from PNG into nearby islands. Outlier analyses also detected four linked alleles introgressed from PNG, with the alleles surrounding 12 resistance-associated cytochrome P450 genes. By treating dispersal as both an intergenerational process and a set of discrete events, we describe a highly interconnected invasive system.