Genomic survey provides insights into the evolutionary changes that occurred during European expansion of the invasive mosquitofish (Gambusia holbrooki)

Genomic investigations of successful invasions: the picture emerging from recent studies

Invasion biology aims to identify traits and mechanisms that contribute to successful invasions, while also providing general insights into the mechanisms underlying population expansion and adaptation to rapid climate and habitat changes. Certain phenotypic attributes have been linked to successful invasions, and the role of genetics has been critical in understanding adaptation of invasive species. Nevertheless, a comprehensive summary evaluating the most common evolutionary mechanisms associated with successful invasions across species and environments is still lacking. Here we present a systematic review of studies since 2015 that have applied genomic tools to investigate mechanisms of successful invasions across different organisms. We examine demographic patterns such as changes in genomic diversity at the population level, the presence of genetic bottlenecks and gene flow in the invasive range. We review mechanisms of adaptation such as selection from standing genetic variation and de novo mutations, hybridisation and introgression, all of which can have an impact on invasion success. This comprehensive review of recent articles on the genomic diversity of invasive species led to the creation of a searchable database to provide researchers with an accessible resource. Analysis of this database allowed quantitative assessment of demographic and adaptive mechanisms acting in invasive species. A predominant role of admixture in increasing levels of genetic diversity enabling molecular adaptation in novel habitats is the most important finding of our study. The "genetic paradox" of invasive species was not validated in genomic data across species and ecosystems. Even though the presence of genetic drift and bottlenecks is commonly reported upon invasion, a large reduction in genomic diversity is rarely observed. Any decrease in genetic diversity is often relatively mild and almost always restored via gene flow between different invasive populations. The fact that loci under selection are frequently detected suggests that adaptation to novel habitats on a molecular level is not hindered. The above findings are confirmed herein for the first time in a semi-quantitative manner by molecular data. We also point to gaps and potential improvements in the design of studies of mechanisms driving rapid molecular adaptation in invasive populations. These include the scarcity of comprehensive studies that include sampling from multiple native and invasive populations, identification of invasion sources, longitudinal population sampling, and the integration of fitness measures into genomic analyses. We also note that the potential of whole genome studies is often not exploited fully in predicting invasive potential. Comparative genomic studies identifying genome features promoting invasions are underrepresented despite their potential for use as a tool in invasive species control.

Complex patterns shape immune genes diversity during invasion of common raccoon in Europe - Selection in action despite genetic drift

Rapid adaptation is common in invasive populations and is crucial to their long-term success. The primary target of selection in the invasive species' new range is standing genetic variation. Therefore, genetic drift and natural selection acting on existing variation are key evolutionary processes through which invaders will evolve over a short timescale. In this study, we used the case of the raccoon Procyon lotor invasion in Europe to identify the forces shaping the diversity of immune genes during invasion. The genes involved in the defence against infection should be under intense selection pressure in the invasive range where novel pathogens are expected to occur. To disentangle the selective and demographic processes shaping the adaptive immune diversity of its invasive and expanding populations, we have developed species-specific single-nucleotide polymorphism markers located in the coding regions of targeted immune-related genes. We characterised the genetic diversity of 110 functionally important immune genes in two invasive and one native raccoon genetic clusters, each presenting a different demographic history. Despite the strong effect of demographic processes in the invasive clusters, we detected a subset of genes exhibiting the diversity pattern suggestive of selection. The most likely process shaping the variation in those genes was balancing selection. The selected genes belong to toll-like receptors and cytokine-related genes. Our results suggest that the prevalence of selection depends on the level of diversity, that is - less genetically diverse invasive population from the Czech Republic displayed fewer signs of selection. Our results highlight the role of standing genetic variation in adapting to new environment. Understanding the evolutionary mechanisms behind invasion success would enable predicting how populations may respond to environmental change.

Population genomics reveal rapid genetic differentiation in a recently invasive population of Rattus norvegicus

Background

Invasive species bring a serious effect on local biodiversity, ecosystems, and even human health and safety. Although the genetic signatures of historical range expansions have been explored in an array of species, the genetic consequences of contemporary range expansions have received little attention, especially in mammal species. In this study, we used whole-genome sequencing to explore the rapid genetic change and introduction history of a newly invasive brown rat (Rattus norvegicus) population which invaded Xinjiang Province, China in the late 1970s.

Results

Bayesian clustering analysis, principal components analysis, and phylogenetic analysis all showed clear genetic differentiation between newly introduced and native rat populations. Reduced genetic diversity and high linkage disequilibrium suggested a severe population bottleneck in this colonization event. Results of TreeMix analyses revealed that the introduced rats were derived from an adjacent population in geographic region (Northwest China). Demographic analysis indicated that a severe bottleneck occurred in XJ population after the split off from the source population, and the divergence of XJ population might have started before the invasion of XJ. Moreover, we detected 42 protein-coding genes with allele frequency shifts throughout the genome for XJ rats and they were mainly associated with lipid metabolism and immunity, which could be seen as a prelude to future selection analyses in the novel environment of XJ.

Conclusions

This study presents the first genomic evidence on genetic differentiation which developed rapidly, and deepens the understanding of invasion history and evolutionary processes of this newly introduced rat population. This would add to our understanding of how invasive species become established and aid strategies aimed at the management of this notorious pest that have spread around the world with humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00387-z.

Differential immune gene expression associated with contemporary range expansion in two invasive rodents in Senegal

Biological invasions are major anthropogenic changes associated with threats to biodiversity and health. However, what determines the successful establishment and spread of introduced populations remains unclear. Here, we explore several hypotheses linking invasion success and immune phenotype traits, including those based on the evolution of increased competitive ability concept. We compared gene expression profiles between anciently and recently established populations of two major invading species, the house mouse Mus musculus domesticus and the black rat Rattus rattus, in Senegal (West Africa). Transcriptome analyses identified differential expression between anciently and recently established populations for 364 mouse genes and 83 rat genes. All immune-related genes displaying differential expression along the mouse invasion route were overexpressed at three of the four recently invaded sites studied. Complement activation pathway genes were overrepresented among these genes. By contrast, no particular immunological process was found to be overrepresented among the differentially expressed genes of black rat. Changes in transcriptome profiles were thus observed along invasion routes, but with different specific patterns between the two invasive species. These changes may be driven by increases in infection risks at sites recently invaded by the house mouse, and by stochastic events associated with colonization history for the black rat. These results constitute a first step toward the identification of immune eco-evolutionary processes potentially involved in the invasion success of these two rodent species.

Seascape genomics reveals population isolation in the reef-building honeycomb worm, Sabellaria alveolata (L.)

BACKGROUND

Under the threat of climate change populations can disperse, acclimatise or evolve in order to avoid fitness loss. In light of this, it is important to understand neutral gene flow patterns as a measure of dispersal potential, but also adaptive genetic variation as a measure of evolutionary potential. In order to assess genetic variation and how this relates to environment in the honeycomb worm (Sabellaria alveolata (L.)), a reef-building polychaete that supports high biodiversity, we carried out RAD sequencing using individuals from along its complete latitudinal range. Patterns of neutral population genetic structure were compared to larval dispersal as predicted by ocean circulation modelling, and outlier analyses and genotype-environment association tests were used to attempt to identify loci under selection in relation to local temperature data.

RESULTS

We genotyped 482 filtered SNPs, from 68 individuals across nine sites, 27 of which were identified as outliers using BAYESCAN and ARLEQUIN. All outlier loci were potentially under balancing selection, despite previous evidence of local adaptation in the system. Limited gene flow was observed among reef-sites (FST = 0.28 ± 0.10), in line with the low dispersal potential identified by the larval dispersal models. The North Atlantic reef emerged as a distinct population and this was linked to high local larval retention and the effect of the North Atlantic Current on dispersal.

CONCLUSIONS

As an isolated population, with limited potential for natural genetic or demographic augmentation from other reefs, the North Atlantic site warrants conservation attention in order to preserve not only this species, but above all the crucial functional ecological roles that are associated with their bioconstructions. Our study highlights the utility of using seascape genomics to identify populations of conservation concern.

Chromosome-level genome assembly of the female western mosquitofish (Gambusia affinis)

BACKGROUND

The western mosquitofish (Gambusia affinis) is a sexually dimorphic poeciliid fish known for its worldwide biological invasion and therefore an important research model for studying invasion biology. This organism may also be used as a suitable model to explore sex chromosome evolution and reproductive development in terms of differentiation of ZW sex chromosomes, ovoviviparity, and specialization of reproductive organs. However, there is a lack of high-quality genomic data for the female G. affinis; hence, this study aimed to generate a chromosome-level genome assembly for it.

RESULTS

The chromosome-level genome assembly was constructed using Oxford nanopore sequencing, BioNano, and Hi-C technology. G. affinis genomic DNA sequences containing 217 contigs with an N50 length of 12.9 Mb and 125 scaffolds with an N50 length of 26.5 Mb were obtained by Oxford nanopore and BioNano, respectively, and the 113 scaffolds (90.4% of scaffolds containing 97.9% nucleotide bases) were assembled into 24 chromosomes (pseudo-chromosomes) by Hi-C. The Z and W chromosomes of G. affinis were identified by comparative genomic analysis of female and male G. affinis, and the mechanism of differentiation of the Z and W chromosomes was explored. Combined with transcriptome data from 6 tissues, a total of 23,997 protein-coding genes were predicted and 23,737 (98.9%) genes were functionally annotated.

CONCLUSIONS

The high-quality female G. affinis reference genome provides a valuable omics resource for future studies of comparative genomics and functional genomics to explore the evolution of Z and W chromosomes and the reproductive developmental biology of G. affinis.

Optimizing Sample Size for Population Genomic Study in a Global Invasive Lady Beetle, Harmonia Axyridis

Finding optimal sample sizes is critical for the accurate estimation of genetic diversity of large invasive populations. Based on previous studies, we hypothesized that a minimal sample size of 3-8 individuals is sufficient to dissect the population architecture of the harlequin lady beetle, Harmonia axyridis, a biological control agent and an invasive alien species. Here, equipped with a type IIB endonuclease restriction site-associated (2b-RAD) DNA sequencing approach, we identified 13,766 and 13,929 single nucleotide polymorphisms (SNPs), respectively, among native and invasive H. axyridis populations. With this information we simulated populations using a randomly selected 3000 SNPs and a subset of individuals. From this simulation we finally determined that six individuals is the minimum sample size required for the accurate estimation of intra- and inter-population genetic diversity within and across H. axyridis populations. Our findings provide an empirical advantage for population genomic studies of H. axyridis in particular and suggest useful tactics for similar studies on multicellular organisms in general.

Population structure, demographic history and local adaptation of the grass carp

BACKGROUND

Genetic diversity within a species reflects population evolution, ecology, and ability to adapt. Genome-wide population surveys of both natural and introduced populations provide insights into genetic diversity, the evolutionary processes and the genetic basis underlying local adaptation. Grass carp is the most important freshwater foodfish species for food and water weed control. However, there is as yet no overall picture on genetic variations and population structure of this species, which is important for its aquaculture.

RESULTS

We used 43,310 SNPs to infer the population structure, evidence of local adaptation and sources of introduction. The overall genetic differentiation of this species was low. The native populations were differentiated into three genetic clusters, corresponding to the Yangtze, Pearl and Heilongjiang River Systems, respectively. The populations in Malaysia, India and Nepal were introduced from both the Yangtze and Pearl River Systems. Loci and genes involved in putative local selection for native locations were identified. Evidence of both positive and balancing selection was found in the introduced locations. Genes associated with loci under putative selection were involved in many biological functions. Outlier loci were grouped into clusters as genomic islands within some specific genomic regions, which likely agrees with the divergence hitchhiking scenario of divergence-with-gene-flow.

CONCLUSIONS

This study, for the first time, sheds novel insights on the population differentiation of the grass carp, genetics of its strong ability in adaption to diverse environments and sources of some introduced grass carp populations. Our data also suggests that the natural populations of the grass carp have been affected by the aquaculture besides neutral and adaptive forces.

Genomic changes in the biological control agent Cryptolaemus montrouzieri associated with introduction

Biological control is the main purpose of intentionally introducing non-native invertebrate species. The evolutionary changes that occur in the populations of the introduced biological control agents may determine the agent's efficiency and the environmental safety. Here, to explore the pattern and extent of potential genomic changes in the worldwide introduced predatory ladybird beetle Cryptolaemus montrouzieri, we used a reduced-representation sequencing method to analyze the genome-wide differentiation of the samples from two native and five introduced locations. Our analyses based on a total of 53,032 single nucleotide polymorphism loci showed that beetles from the introduced locations in Asia and Europe exhibited significant reductions in genetic diversity and high differentiation compared with the samples from the native Australian range. Each introduced population belonged to a unique genetic cluster, while the beetles from two native locations were much more similar. These genomic patterns were also detected when the dataset was pruned for genomic outlier loci (52,318 SNPs remaining), suggesting that random genetic drift was the main force shaping the genetic diversity and population structure of this biological control agent. Our results provide a genome-wide characterization of polymorphisms in a biological control agent and reveal genomic differences that were influenced by the introduction history. These differences might complicate assessments of the efficiency of biological control and the invasion potential of this species but also indicate the feasibility of selective breeding.

Genomic signatures of local adaptation to the degree of environmental predictability in rotifers

Environmental fluctuations are ubiquitous and thus essential for the study of adaptation. Despite this, genome evolution in response to environmental fluctuations —and more specifically to the degree of environmental predictability– is still unknown. Saline lakes in the Mediterranean region are remarkably diverse in their ecological conditions, which can lead to divergent local adaptation patterns in the inhabiting aquatic organisms. The facultatively sexual rotifer Brachionus plicatilis shows diverging local adaptation in its life-history traits in relation to estimated environmental predictability in its habitats. Here, we used an integrative approach —combining environmental, phenotypic and genomic data for the same populations– to understand the genomic basis of this diverging adaptation. Firstly, a novel draft genome for B. plicatilis was assembled. Then, genome-wide polymorphisms were studied using genotyping by sequencing on 270 clones from nine populations in eastern Spain. As a result, 4,543 high-quality SNPs were identified and genotyped. More than 90 SNPs were found to be putatively under selection with signatures of diversifying and balancing selection. Over 140 SNPs were correlated with environmental or phenotypic variables revealing signatures of local adaptation, including environmental predictability. Putative functions were associated to most of these SNPs, since they were located within annotated genes. Our results reveal associations between genomic variation and the degree of environmental predictability, providing genomic evidence of adaptation to local conditions in natural rotifer populations.

Effects of water pollution and river fragmentation on population genetic structure of invasive mosquitofish

We analyzed variation at the GPI-2 locus and eleven microsatellite loci of eastern mosquitofish Gambusia holbrooki populations introduced to the Ebro River (Spain), sampling above and below a dam (Flix Reservoir) where severe chronic pollution has been well documented. Allele frequency changes at the GPI-2 locus in the sites nearest to the polluted sediments agree with previous results from studies in mercury-exposed populations of this highly invasive fish. Genetic distinction of the mosquitofish collected close to the polluted sediments was detected at the GPI locus but also at the presumptive neutral microsatellite loci. Recent migration rates estimated from microsatellites indicated that around 30% of fish collected in a specific location were immigrants from upstream and downstream sources. Such high migration rates probably contribute to the mosquitofish's invasive success and suggest that the consequences on the mosquitofish regional genetic structured of high levels of water toxicants could be mediated by immigration from other sites, but the effect of pollutants on local diversity might be higher than observed here.

Parallel evolution and adaptation to environmental factors in a marine flatfish: Implications for fisheries and aquaculture management of the turbot (Scophthalmus maximus)

Unraveling adaptive genetic variation represents, in addition to the estimate of population demographic parameters, a cornerstone for the management of aquatic natural living resources, which, in turn, represent the raw material for breeding programs. The turbot (Scophthalmus maximus) is a marine flatfish of high commercial value living on the European continental shelf. While wild populations are declining, aquaculture is flourishing in southern Europe. We evaluated the genetic structure of turbot throughout its natural distribution range (672 individuals; 20 populations) by analyzing allele frequency data from 755 single nucleotide polymorphism discovered and genotyped by double-digest RAD sequencing. The species was structured into four main regions: Baltic Sea, Atlantic Ocean, Adriatic Sea, and Black Sea, with subtle differentiation apparent at the distribution margins of the Atlantic region. Genetic diversity and effective population size estimates were highest in the Atlantic populations, the area of greatest occurrence, while turbot from other regions showed lower levels, reflecting geographical isolation and reduced abundance. Divergent selection was detected within and between the Atlantic Ocean and Baltic Sea regions, and also when comparing these two regions with the Black Sea. Evidence of parallel evolution was detected between the two low salinity regions, the Baltic and Black seas. Correlation between genetic and environmental variation indicated that temperature and salinity were probably the main environmental drivers of selection. Mining around the four genomic regions consistently inferred to be under selection identified candidate genes related to osmoregulation, growth, and resistance to diseases. The new insights are useful for the management of turbot fisheries and aquaculture by providing the baseline for evaluating the consequences of turbot releases from restocking and farming.

Using MiddRAD-seq data to develop polymorphic microsatellite markers for an endangered yew species

Microsatellites are highly polymorphic markers which have been used in a wide range of genetic studies. In recent years, various sources of next-generation sequencing data have been used to develop new microsatellite loci, but compared with the more common shotgun genomic sequencing or transcriptome data, the potential utility of RAD-seq data for microsatellite ascertainment is comparatively under-used. In this study, we employed MiddRAD-seq data to develop polymorphic microsatellite loci for the endangered yew species Taxus florinii. Of 8,823,053 clean reads generated for ten individuals of a population, 94,851 (∼1%) contained microsatellite motifs. These corresponded to 2993 unique loci, of which 526 (∼18%) exhibited polymorphism. Of which, 237 were suitable for designing microsatellite primer pairs, and 128 loci were randomly selected for PCR validation and microsatellite screening. Out of the 128 primer pairs, 16 loci gave clear, reproducible patterns, and were then screened and characterized in 24 individuals from two populations. The total number of alleles per locus ranged from two to ten (mean = 4.875), and within-population expected heterozygosity from zero to 0.789 (mean = 0.530), indicating that these microsatellite loci will be useful for population genetics and speciation studies of T. florinii. This study represents one of few examples to mine polymorphic microsatellite loci from ddRAD data.

Sexual selection on male body size, genital length and heterozygosity: Consistency across habitats and social settings

Spatial and temporal variation in environmental factors and the social setting can help to maintain genetic variation in sexually selected traits if it affects the strength of directional selection. A key social parameter which affects the intensity of, and sometimes predicts the response to, mating competition is the operational sex ratio (OSR; ratio of receptive males to females). How the OSR affects selection for specific male traits is poorly understood. It is also unclear how sexual selection is affected by interactions between the OSR and environmental factors, such as habitat complexity, that alter key male-female interactions such as mate encounter rates. Here, we experimentally manipulated the OSR and habitat complexity and quantified sexual selection on male mosquitofish (Gambusia holbrooki) by directly measuring male reproductive success (i.e. paternity). We show that despite a more equitable sharing of paternity (i.e. higher levels of multiple paternity) under a male-biased OSR, selection on focal male traits was unaffected by the OSR or habitat complexity. Instead, sexual selection consistently, and significantly, favoured smaller bodied males, males with higher genome wide heterozygosity (based on >3,000 SNP markers) and males with a relatively long gonopodium (intromittent organ). Our results show that sexual selection on male body size, relative genital size and heterozygosity in this system is consistent across environments that vary in ecological parameters that are expected to influence mate encounter rates.

Why does inbreeding reduce male paternity? Effects on sexually selected traits

Mating with relatives has often been shown to negatively affect offspring fitness (inbreeding depression). There is considerable evidence for inbreeding depression due to effects on naturally selected traits, particularly those expressed early in life, but there is less evidence of it for sexually selected traits. This is surprising because sexually selected traits are expected to exhibit strong inbreeding depression. Here, we experimentally created inbred and outbred male mosquitofish (Gambusia holbrooki). Inbred males were the offspring of matings between full siblings. We then investigated how inbreeding influenced a number of sexually selected male traits, specifically: attractiveness, sperm number and velocity, as well as sperm competitiveness based on a male's share of paternity. We found no inbreeding depression for male attractiveness or sperm traits. There was, however, evidence that lower heterozygosity decreased paternity due to reduced sperm competitiveness. Our results add to the growing evidence that competitive interactions exacerbate the negative effects of the increased homozygosity that arises when there is inbreeding.

On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes

The first goal of this paper was to overview modern approaches to local adaptation, with a focus on the use of population genomics data to detect signals of natural selection in fishes. Several mechanisms are discussed that may enhance the maintenance of genetic variation and evolutionary potential, which have been overlooked and should be considered in future theoretical development and predictive models: the prevalence of soft sweeps, polygenic basis of adaptation, balancing selection and transient polymorphisms, parallel evolution, as well as epigenetic variation. Research on fish population genomics has provided ample evidence for local adaptation at the genome level. Pervasive adaptive evolution, however, seems to almost never involve the fixation of beneficial alleles. Instead, adaptation apparently proceeds most commonly by soft sweeps entailing shifts in frequencies of alleles being shared between differentially adapted populations. One obvious factor contributing to the maintenance of standing genetic variation in the face of selective pressures is that adaptive phenotypic traits are most often highly polygenic, and consequently the response to selection should derive mostly from allelic co-variances among causative loci rather than pronounced allele frequency changes. Balancing selection in its various forms may also play an important role in maintaining adaptive genetic variation and the evolutionary potential of species to cope with environmental change. A large body of literature on fishes also shows that repeated evolution of adaptive phenotypes is a ubiquitous evolutionary phenomenon that seems to occur most often via different genetic solutions, further adding to the potential options of species to cope with a changing environment. Moreover, a paradox is emerging from recent fish studies whereby populations of highly reduced effective population sizes and impoverished genetic diversity can apparently retain their adaptive potential in some circumstances. Although more empirical support is needed, several recent studies suggest that epigenetic variation could account for this apparent paradox. Therefore, epigenetic variation should be fully integrated with considerations pertaining to role of soft sweeps, polygenic and balancing selection, as well as repeated adaptation involving different genetic basis towards improving models predicting the evolutionary potential of species to cope with a changing world.