Evolutionary origins of genomic adaptations in an invasive copepod

Evolution of ion transporter Na+/K+-ATPase expression in the osmoregulatory maxillary glands of an invasive copepod

While many freshwater invaders originate from saline habitats, the physiological mechanisms involved are poorly understood. We investigated the evolution of ion transporter Na+/K+-ATPase (NKA) protein expression between ancestral saline and freshwater invading populations of the copepod Eurytemora carolleae (Atlantic clade of the E. affinis complex). We compared in situ NKA expression between populations under common-garden conditions at three salinities in the maxillary glands. We found the evolution of reduced NKA expression in the freshwater population under freshwater conditions and reduced plasticity (canalization) across salinities, relative to the saline population. Our results support the hypothesis that maxillary glands are involved in ion reabsorption from excretory fluids at low-salinity conditions in the saline population. However, mechanisms of freshwater adaptation, such as increased ion uptake from the environment, might reduce the need for ion reabsorption in the freshwater population. These patterns of ion transporter expression contribute insights into the evolution of ionic regulation during habitat change.

Local genetic adaptation to habitat in wild chimpanzees

How populations adapt to their environment is a fundamental question in biology. Yet we know surprisingly little about this process, especially for endangered species such as non-human great apes. Chimpanzees, our closest living relatives, are particularly interesting because they inhabit diverse habitats, from rainforest to woodland-savannah. Whether genetic adaptation facilitates such habitat diversity remains unknown, despite having wide implications for evolutionary biology and conservation. Using 828 newly generated exomes from wild chimpanzees, we find evidence of fine-scale genetic adaptation to habitat. Notably, adaptation to malaria in forest chimpanzees is mediated by the same genes underlying adaptation to malaria in humans. This work demonstrates the power of non-invasive samples to reveal genetic adaptations in endangered populations and highlights the importance of adaptive genetic diversity for chimpanzees.

Navigating spatio-temporal microbiome dynamics: Environmental factors and trace elements shape the symbiont community of an invasive marine species

The proliferation of marine invasive species is a mounting concern. While the role of microbial communities in invasive ascidian species is recognized, the role of seasonal shifts in microbiome composition remains largely unexplored. We sampled five individuals of the invasive ascidian Styela plicata quarterly from January 2020 to October 2021 in two harbours, examining gills, tunics, and surrounding water. By analysing Amplicon Sequence Variants (ASVs) and seawater trace elements, we found that compartment (seawater, tunic, or gills) was the primary differentiating factor, followed by harbour. Clear seasonal patterns were evident in seawater bacteria, less so in gills, and absent in tunics. We identified compartment-specific bacteria, as well as seasonal indicator ASVs and ASVs correlated with trace element concentrations. Among these bacteria, we found that Endozoicomonas, Hepatoplasma and Rhodobacteraceae species had reported functions which might be necessary for overcoming seasonality and trace element shifts. This study contributes to understanding microbiome dynamics in invasive holobiont systems, and the patterns found indicate a potential role in adaptation and invasiveness.

When adaptation is slowed down: Genomic analysis of evolutionary stasis in thermal tolerance during biological invasion in a novel climate

Research conducted during the past two decades has demonstrated that biological invasions are excellent models of rapid evolution. Even so, characteristics of invasive populations such as a short time for recombination to assemble optimal combinations of alleles may occasionally limit adaptation to new environments. Here, we investigated such genetic constraints to adaptation in the invasive brown anole (Anolis sagrei)-a tropical ectotherm that was introduced to the southeastern United States, a region with a much colder climate than in its native Caribbean range. We examined thermal physiology for 30 invasive populations and tested for a climatic cline in cold tolerance. Also, we used genomics to identify mechanisms that may limit adaptation. We found no support for a climatic cline, indicating that thermal tolerance did not shift adaptively. Concomitantly, population genomic results were consistent with the occurrence of recombination cold spots that comprise more than half of the genome and maintain long-range associations among alleles in invasive populations. These genomic regions overlap with both candidate thermal tolerance loci that we identified using a standard genome-wide association test. Moreover, we found that recombination cold spots do not have a large contribution to population differentiation in the invasive range, contrary to observations in the native range. We suggest that limited recombination is constraining the contribution of large swaths of the genome to adaptation in invasive brown anoles. Our study provides an example of evolutionary stasis during invasion and highlights the possibility that reduced recombination occasionally slows down adaptation in invasive populations.

Divergence time shapes gene reuse during repeated adaptation

When diverse lineages repeatedly adapt to similar environmental challenges, the extent to which the same genes are involved (gene reuse) varies across systems. We propose that divergence time among lineages is a key factor driving this variability: as lineages diverge, the extent of gene reuse should decrease due to reductions in allele sharing, functional differentiation among genes, and restructuring of genome architecture. Indeed, we show that many genomic studies of repeated adaptation find that more recently diverged lineages exhibit higher gene reuse during repeated adaptation, but the relationship becomes less clear at older divergence time scales. Thus, future research should explore the factors shaping gene reuse and their interplay across broad divergence time scales for a deeper understanding of evolutionary repeatability.

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.

Rapid evolution of phenotypic plasticity in patchy habitats

Phenotypic plasticity may evolve rapidly, enabling a population's persistence in the face of sudden environmental change. Rapid evolution can occur when there is considerable genetic polymorphism at selected loci. We propose that balancing selection could be one of the mechanisms that sustain such polymorphism for plasticity. We use stochastic Monte Carlo simulations and deterministic analysis to investigate the evolution of a plasticity modifier locus in structured populations inhabiting favorable and adverse environments, i.e. patchy habitats. We survey a wide range of parameters including selective pressures on a target (structural) locus, plasticity effects, population sizes, and migration patterns between demes including periodic or continuous bidirectional and source-sink dynamics. We find that polymorphism in phenotypic plasticity can be maintained under a wide range of environmental scenarios in both favorable and adverse environments due to the balancing effect of population structure in patchy habitats. This effect offers a new plausible explanation for the rapid evolution of plasticity in nature: Phenotypic plasticity may rapidly evolve from genetic variation maintained by balancing selection if the population has experienced immigration from populations under different selection regimes.

What Can Evolutionary History Tell Us about the Functioning of Ecological Communities? The ASN Presidential Debate

AbstractIn January 2018, Sharon Strauss, then president of the American Society of Naturalists, organized a debate on the following topic: does evolutionary history inform the current functioning of ecological communities? The debaters-Ives, Lau, Mayfield, and Tobias-presented pro and con arguments, caricatured in standard debating format. Numerous examples show that both recent microevolutionary and longer-term macroevolutionary history are important to the ecological functioning of communities. On the other hand, many other examples illustrate that the evolutionary history of communities or community members does not influence ecological function, or at least not very much. This article aims to provide a provocative discussion of the consistent and conflicting patterns that emerge in the study of contemporary and historical evolutionary influences on community function, as well as to identify questions for further study. It is intended as a thought-provoking exercise to explore this complex field, specifically addressing (1) key assumptions and how they can lead us astray and (2) issues that need additional study. The debaters all agree that evolutionary history can inform us about at least some aspects of community function. The underlying question at the root of the debate, however, is how the fields of ecology and evolution can most profitably collaborate to provide a deeper and broader understanding of ecological communities.

Genome architecture underlying salinity adaptation in the invasive copepod Eurytemora affinis species complex: A review

With climate change, habitat salinity is shifting rapidly throughout the globe. In addition, many destructive freshwater invaders are recent immigrants from saline habitats. Recently, populations of the copepod Eurytemora affinis species complex have invaded freshwater habitats multiple times independently from saline estuaries on three continents. This review discusses features of this species complex that could enhance their evolutionary potential during rapid environmental change. Remarkably, across independent freshwater invasions, natural selection has repeatedly favored the same alleles far more than expected. This high degree of parallelism is surprising, given the expectation of nonparallel evolution for polygenic adaptation. Factors such as population structure and the genome architecture underlying critical traits under selection might help drive rapid adaptation and parallel evolution. Given the preponderance of saline-to-freshwater invasions and climate-induced salinity change, the principles found here could provide invaluable insights into mechanisms operating in other systems and the potential for adaptation in a changing planet.

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.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

A Genomic Quantitative Study on the Contribution of the Ancestral-State Bases Relative to Derived Bases in the Divergence and Local Adaptation of Populus davidiana

Identifying alleles associated with adaptation to new environments will advance our understanding of evolution from the molecular level. Previous studies have found that the Populus davidiana southwest population in East Asia has differentiated from other populations in the range. We aimed to evaluate the contributions of the ancestral-state bases (ASBs) relative to derived bases (DBs) in the local adaptation of P. davidiana in the Yunnan-Guizhou Plateau from a quantitative perspective using whole-genome re-sequencing data from 90 P. davidiana samples from three regions across the species range. Our results showed that the uplift of the Qinghai-Tibet Plateau during the Neogene and associated climate fluctuations during the Middle Pleistocene were likely an important factor in the early divergence of P. davidiana. Highly differentiated genomic regions between populations were inferred to have undergone strong linked natural selection, and ASBs are the chief means by which populations of P. davidiana adapt to novel environmental conditions; however, when adapting to regions with high environmental differences relative to the ancestral range, the proportion of DBs was significantly higher than that of background regions, as ASBs are insufficient to cope with these environments. Finally, a number of genes were identified in the outlier region.

Fluctuating selection and the determinants of genetic variation

Recent studies of cosmopolitan Drosophila populations have found hundreds to thousands of genetic loci with seasonally fluctuating allele frequencies, bringing temporally fluctuating selection to the forefront of the historical debate surrounding the maintenance of genetic variation in natural populations. Numerous mechanisms have been explored in this longstanding area of research, but these exciting empirical findings have prompted several recent theoretical and experimental studies that seek to better understand the drivers, dynamics, and genome-wide influence of fluctuating selection. In this review, we evaluate the latest evidence for multilocus fluctuating selection in Drosophila and other taxa, highlighting the role of potential genetic and ecological mechanisms in maintaining these loci and their impacts on neutral genetic variation.

Evolutionary genomics: Insights from the invasive European starlings

Two fundamental questions for evolutionary studies are the speed at which evolution occurs, and the way that this evolution may present itself within an organism's genome. Evolutionary studies on invasive populations are poised to tackle some of these pressing questions, including understanding the mechanisms behind rapid adaptation, and how it facilitates population persistence within a novel environment. Investigation of these questions are assisted through recent developments in experimental, sequencing, and analytical protocols; in particular, the growing accessibility of next generation sequencing has enabled a broader range of taxa to be characterised. In this perspective, we discuss recent genetic findings within the invasive European starlings in Australia, and outline some critical next steps within this research system. Further, we use discoveries within this study system to guide discussion of pressing future research directions more generally within the fields of population and evolutionary genetics, including the use of historic specimens, phenotypic data, non-SNP genetic variants (e.g., structural variants), and pan-genomes. In particular, we emphasise the need for exploratory genomics studies across a range of invasive taxa so we can begin understanding broad mechanisms that underpin rapid adaptation in these systems. Understanding how genetic diversity arises and is maintained in a population, and how this contributes to adaptability, requires a deep understanding of how evolution functions at the molecular level, and is of fundamental importance for the future studies and preservation of biodiversity across the globe.

Relaxed risk of predation drives parallel evolution of stickleback behavior

The occurrence of similar phenotypes in multiple independent populations derived from common ancestral conditions (viz. parallel evolution) is a testimony of evolution by natural selection. Parallel evolution implies that populations share a common phenotypic response to a common selection pressure associated with habitat similarity. Examples of parallel evolution at genetic and phenotypic levels are fairly common, but the driving selective agents often remain elusive. Similarly, the role of phenotypic plasticity in facilitating early stages of parallel evolution is unclear. We investigated whether the relaxation of predation pressure associated with the colonization of freshwater ponds by nine-spined sticklebacks (Pungitius pungitius) likely explains the divergence in complex behaviors between marine and pond populations, and whether this divergence is parallel. Using laboratory-raised individuals exposed to different levels of perceived predation risk, we calculated vectors of phenotypic divergence for four behavioral traits between habitats and predation risk treatments. We found a significant correlation between the directions of evolutionary divergence and phenotypic plasticity, suggesting that divergence in behavior between habitats is aligned with the response to relaxation of predation pressure. Finally, we show alignment across multiple pairs of populations, and that relaxation of predation pressure has likely driven parallel evolution of behavior in this species.

Mechanisms of Na+ uptake from freshwater habitats in animals

Life in fresh water is osmotically and energetically challenging for living organisms, requiring increases in ion uptake from dilute environments. However, mechanisms of ion uptake from freshwater environments are still poorly understood and controversial, especially in arthropods, for which several hypothetical models have been proposed based on incomplete data. One compelling model involves the proton pump V-type H⁺ ATPase (VHA), which energizes the apical membrane, enabling the uptake of Na⁺ (and other cations) via an unknown Na⁺ transporter (referred to as the "Wieczorek Exchanger" in insects). What evidence exists for this model of ion uptake and what is this mystery exchanger or channel that cooperates with VHA? We present results from studies that explore this question in crustaceans, insects, and teleost fish. We argue that the Na⁺/H⁺ antiporter (NHA) is a likely candidate for the Wieczorek Exchanger in many crustaceans and insects; although, there is no evidence that this is the case for fish. NHA was discovered relatively recently in animals and its functions have not been well characterized. Teleost fish exhibit redundancy of Na⁺ uptake pathways at the gill level, performed by different ion transporter paralogs in diverse cell types, apparently enabling tolerance of low environmental salinity and various pH levels. We argue that much more research is needed on overall mechanisms of ion uptake from freshwater habitats, especially on NHA and other potential Wieczorek Exchangers. Such insights gained would contribute greatly to our general understanding of ionic regulation in diverse species across habitats.

Experimental evolution reveals the synergistic genomic mechanisms of adaptation to ocean warming and acidification in a marine copepod

Metazoan adaptation to global change relies on selection of standing genetic variation. Determining the extent to which this variation exists in natural populations, particularly for responses to simultaneous stressors, is essential to make accurate predictions for persistence in future conditions. Here, we identified the genetic variation enabling the copepod Acartia tonsa to adapt to experimental ocean warming, acidification, and combined ocean warming and acidification (OWA) over 25 generations of continual selection. Replicate populations showed a consistent polygenic response to each condition, targeting an array of adaptive mechanisms including cellular homeostasis, development, and stress response. We used a genome-wide covariance approach to partition the allelic changes into three categories: selection, drift and replicate-specific selection, and laboratory adaptation responses. The majority of allele frequency change in warming (57%) and OWA (63%) was driven by shared selection pressures across replicates, but this effect was weaker under acidification alone (20%). OWA and warming shared 37% of their response to selection but OWA and acidification shared just 1%, indicating that warming is the dominant driver of selection in OWA. Despite the dominance of warming, the interaction with acidification was still critical as the OWA selection response was highly synergistic with 47% of the allelic selection response unique from either individual treatment. These results disentangle how genomic targets of selection differ between single and multiple stressors and demonstrate the complexity that nonadditive multiple stressors will contribute to predictions of adaptation to complex environmental shifts caused by global change.

Genome-wide signatures of synergistic epistasis during parallel adaptation in a Baltic Sea copepod

The role of epistasis in driving adaptation has remained an unresolved problem dating back to the Evolutionary Synthesis. In particular, whether epistatic interactions among genes could promote parallel evolution remains unexplored. To address this problem, we employ an Evolve and Resequence (E&R) experiment, using the copepod Eurytemora affinis, to elucidate the evolutionary genomic response to rapid salinity decline. Rapid declines in coastal salinity at high latitudes are a predicted consequence of global climate change. Based on time-resolved pooled whole-genome sequencing, we uncover a remarkably parallel, polygenic response across ten replicate selection lines, with 79.4% of selected alleles shared between lines by the tenth generation of natural selection. Using extensive computer simulations of our experiment conditions, we find that this polygenic parallelism is consistent with positive synergistic epistasis among alleles, far more so than other mechanisms tested. Our study provides experimental and theoretical support for a novel mechanism promoting repeatable polygenic adaptation, a phenomenon that may be common for selection on complex physiological traits.

Identification and characterization of homeobox gene clusters in harpacticoid and calanoid copepods

In this study, we have identified the entire complement of typical homeobox (Hox) genes (Lab, Pb, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) in harpacticoid and calanoid copepods and compared them with the cyclopoid copepod Paracyclopina nana. The harpacticoid copepods Tigriopus japonicus and Tigriopus kingsejongensis have seven Hox genes (Lab, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) and the Pb and Ftz genes are also present in the cyclopoid copepod P. nana. In the Hox gene cluster of the calanoid copepod Eurytemora affinis, all the Hox genes were present linearly in the genome but the Antp gene was duplicated. Of the three representative copepods, the P. nana Hox gene cluster was the most compact due to its small genome size. The Hox gene expression profile patterns in the three representative copepods were stage-specific. The Lab, Dfd, Scr, Pb, Ftz, and Hox3 genes showed a high expression in early developmental stages but Antp, Ubx, Abd-A, and Abd-B genes were mostly expressed in later developmental stages, implying that these Hox genes may be closely associated with the development of segment identity during early development.

Recognizing salinity threats in the climate crisis

Climate change is causing habitat salinity to transform at unprecedented rates across the globe. While much of the research on climate change has focused on rapid shifts in temperature, far less attention has focused on the effects of changes in environmental salinity. Consequently, predictive studies on the physiological, evolutionary, and migratory responses of organisms and populations to the threats of salinity change are relatively lacking. This omission represents a major oversight, given that salinity is among the most important factors that define biogeographic boundaries in aquatic habitats. In this perspective, we briefly touch on responses of organisms and populations to rapid changes in salinity occurring on contemporary time scales. We then discuss factors that might confer resilience to certain taxa, enabling them to survive rapid salinity shifts. Next, we consider approaches for predicting how geographic distributions will shift in response to salinity change. Finally, we identify additional data that are needed to make better predictions in the future. Future studies on climate change should account for the multiple environmental factors that are rapidly changing, especially habitat salinity.

Novel genome characteristics contribute to the invasiveness of Phragmites australis (common reed)

The rapid invasion of the non‐native Phragmites australis (Poaceae, subfamily Arundinoideae) is a major threat to native wetland ecosystems in North America and elsewhere. We describe the first reference genome for P. australis and compare invasive (ssp. australis) and native (ssp. americanus) genotypes collected from replicated populations across the Laurentian Great Lakes to deduce genomic bases driving its invasive success. Here, we report novel genomic features including a Phragmites lineage‐specific whole genome duplication, followed by gene loss and preferential retention of genes associated with transcription factors and regulatory functions in the remaining duplicates. Comparative transcriptomic analyses revealed that genes associated with biotic stress and defence responses were expressed at a higher basal level in invasive genotypes, but native genotypes showed a stronger induction of defence responses when challenged by a fungal endophyte. The reference genome and transcriptomes, combined with previous ecological and environmental data, add to our understanding of mechanisms leading to invasiveness and support the development of novel, genomics‐assisted management approaches for invasive Phragmites.

Genome of the estuarine oyster provides insights into climate impact and adaptive plasticity

Understanding the roles of genetic divergence and phenotypic plasticity in adaptation is central to evolutionary biology and important for assessing adaptive potential of species under climate change. Analysis of a chromosome-level assembly and resequencing of individuals across wide latitude distribution in the estuarine oyster (Crassostrea ariakensis) revealed unexpectedly low genomic diversity and population structures shaped by historical glaciation, geological events and oceanographic forces. Strong selection signals were detected in genes responding to temperature and salinity stress, especially of the expanded solute carrier families, highlighting the importance of gene expansion in environmental adaptation. Genes exhibiting high plasticity showed strong selection in upstream regulatory regions that modulate transcription, indicating selection favoring plasticity. Our findings suggest that genomic variation and population structure in marine bivalves are heavily influenced by climate history and physical forces, and gene expansion and selection may enhance phenotypic plasticity that is critical for the adaptation to rapidly changing environments.

Breathing space: deoxygenation of aquatic environments can drive differential ecological impacts across biological invasion stages

The influence of climate change on the ecological impacts of invasive alien species (IAS) remains understudied, with deoxygenation of aquatic environments often-overlooked as a consequence of climate change. Here, we therefore assessed how oxygen saturation affects the ecological impact of a predatory invasive fish, the Ponto-Caspian round goby (Neogobius melanostomus), relative to a co-occurring endangered European native analogue, the bullhead (Cottus gobio) experiencing decline in the presence of the IAS. In individual trials and mesocosms, we assessed the effect of high, medium and low (90%, 60% and 30%) oxygen saturation on: (1) functional responses (FRs) of the IAS and native, i.e. per capita feeding rates; (2) the impact on prey populations exerted; and (3) how combined impacts of both fishes change over invasion stages (Pre-invasion, Arrival, Replacement, Proliferation). Both species showed Type II potentially destabilising FRs, but at low oxygen saturation, the invader had a significantly higher feeding rate than the native. Relative Impact Potential, combining fish per capita effects and population abundances, revealed that low oxygen saturation exacerbates the high relative impact of the invader. The Relative Total Impact Potential (RTIP), modelling both consumer species’ impacts on prey populations in a system, was consistently higher at low oxygen saturation and especially high during invader Proliferation. In the mesocosm experiment, low oxygen lowered RTIP where both species were present, but again the IAS retained high relative impact during Replacement and Proliferation stages at low oxygen. We also found evidence of multiple predator effects, principally antagonism. We highlight the threat posed to native communities by IAS alongside climate-related stressors, but note that solutions may be available to remedy hypoxia and potentially mitigate impacts across invasion stages.

Ion Transporter Gene Families as Physiological Targets of Natural Selection During Salinity Transitions in a Copepod

Salinity is a key factor that structures biodiversity on the planet. With anthropogenic change, such as climate change and species invasions, many populations are facing rapid and dramatic changes in salinity throughout the globe. Studies on the copepod Eurytemora affinis species complex have implicated ion transporter gene families as major loci contributing to salinity adaptation during freshwater invasions. Laboratory experiments and population genomic surveys of wild populations have revealed evolutionary shifts in genome-wide gene expression and parallel genomic signatures of natural selection during independent salinity transitions. Our results suggest that balancing selection in the native range and epistatic interactions among specific ion transporter paralogs could contribute to parallel freshwater adaptation. Overall, these studies provide unprecedented insights into evolutionary mechanisms underlying physiological adaptation during rapid salinity change.

Genomic analysis unveils mechanisms of northward invasion and signatures of plateau adaptation in the Asian house rat

The Asian house rat (AHR), Rattus tanezumi, has recently invaded the northern half of China. The AHR is a highly adaptive rat species that has also successfully conquered the Qinghai-Tibet Plateau (QTP) and replaced the brown rat (BR), R. norvegicus, at the edge of the QTP. Here, we assembled a draft genome of the AHR and explored the mechanisms of its northward invasion and the genetic basis underlying plateau adaptation in this species. Population genomic analyses revealed that the northwardly invasive AHRs consisted of two independent and genetically distinct populations which might result from multiple independent primary invasion events. One invasive population exhibited reduced genetic diversity and distinct population structure compared with its source population, while the other displayed preserved genetic polymorphisms and little genetic differentiation from its source population. Genes involved in G-protein coupled receptors and carbohydrate metabolism may contribute to the local adaptation of northern AHRs. In particular, RTN4 was identified as a key gene for AHRs in the QTP that favours adaptation to high-altitude hypoxia. Coincidently, the physiological performance and transcriptome profiles of hypoxia-exposed rats both showed better hypoxia adaptation in AHRs than in BRs that failed to colonize the heart of the QTP, which may have facilitated the replacement of the BR population by the invading AHRs at the edge of the QTP. This study provides profound insights into the multiple origins of the northwardly invasive AHR and the great tolerance to hypoxia in this species.

The salmon louse genome: Copepod features and parasitic adaptations

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification.

Invasion of the body snatchers: the role of parasite introduction in host distribution and response to salinity in invaded estuaries

In dynamic systems, organisms are faced with variable selective forces that may impose trade-offs. In estuaries, salinity is a strong driver of organismal diversity, while parasites shape species distributions and demography. We tested for trade-offs between low-salinity stress and parasitism in an invasive castrating parasite and its mud crab host along salinity gradients of two North Carolina rivers. We performed field surveys every six to eight weeks over 3 years to determine factors influencing parasite prevalence, host abundance, and associated taxa diversity. We also looked for signatures of low-salinity stress in the host by examining its response (time-to-right and gene expression) to salinity. We found salinity and temperature significantly affected parasite prevalence, with low-salinity sites (less than 10 practical salinity units (PSU)) lacking infection, and populations in moderate salinities at warmer temperatures reaching prevalence as high as 60%. Host abundance was negatively associated with parasite prevalence. Host gene expression was plastic to acclimation salinity, but several osmoregulatory and immune-related genes demonstrated source-dependent salinity response. We identified a genetic marker that was strongly associated with salinity against a backdrop of no neutral genetic structure, suggesting possible selection on standing variation. Our study illuminates how selective trade-offs in naturally dynamic systems may shape host evolutionary ecology.

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.

The genome of the European estuarine calanoid copepod Eurytemora affinis: Potential use in molecular ecotoxicology

In this study, we sequenced and assembled the genome of a European estuarine calanoid copepod using Oxford Nanopore PromethION and Illumina HiSeq 2500 platforms. The length of the assembled genome was 776.1 Mb with N50 = 474.9 kb (BUSCO 85.9%), and the genome consisted of 2473 contigs. A total of 18,014 genes were annotated and orthologous gene clusters were analyzed in comparison to other copepods. In addition, genome-wide identification of cytochrome P450s, glutathione S-transferases, and ATP-binding cassette transporters in E. affinis was performed to determine gene repertoire of these detoxification-related gene families. Results revealed the presence of species-specific gene inventories, indicating that these gene families have evolved through species-specific gene loss/expansion processes, possibly due to adaptation to different environmental stressors. Our study provides a new inventory of the European estuarine calanoid copepod E. affinis genome with emphasis on phase I, II, and III detoxification systems.

Rapid genetic adaptation to recently colonized environments is driven by genes underlying life history traits

Background

Uncovering the mechanisms underlying rapid genetic adaptation can provide insight into adaptive evolution and shed light on conservation, invasive species control, and natural resource management. However, it can be difficult to experimentally explore rapid adaptation due to the challenges associated with propagating and maintaining species in captive environments for long periods of time. By contrast, many introduced species have experienced strong selection when colonizing environments that differ substantially from their native range and thus provide a “natural experiment” for studying rapid genetic adaptation. One such example occurred when sea lamprey (Petromyzon marinus), native to the northern Atlantic, naturally migrated into Lake Champlain and expanded their range into the Great Lakes via man-made shipping canals.

Results

Utilizing 368,886 genome-wide single nucleotide polymorphisms (SNPs), we calculated genome-wide levels of genetic diversity (i.e., heterozygosity and π) for sea lamprey collected from native (Connecticut River), native but recently colonized (Lake Champlain), and invasive (Lake Michigan) populations, assessed genetic differentiation between all populations, and identified candidate genes that responded to selection imposed by the novel environments. We observed a 14 and 24% reduction in genetic diversity in Lake Michigan and Lake Champlain populations, respectively, compared to individuals from the Connecticut River, suggesting that sea lamprey populations underwent a genetic bottleneck during colonization. Additionally, we identified 121 and 43 outlier genes in comparisons between Lake Michigan and Connecticut River and between Lake Champlain and Connecticut River, respectively. Six outlier genes that contained synonymous SNPs in their coding regions and two genes that contained nonsynonymous SNPs may underlie the rapid evolution of growth (i.e., GHR), reproduction (i.e., PGR, TTC25, STARD10), and bioenergetics (i.e., OXCT1, PYGL, DIN4, SLC25A15).

Conclusions

By identifying the genomic basis of rapid adaptation to novel environments, we demonstrate that populations of invasive species can be a useful study system for understanding adaptive evolution. Furthermore, the reduction in genome-wide levels of genetic diversity associated with colonization coupled with the identification of outlier genes underlying key life history traits known to have changed in invasive sea lamprey populations (e.g., growth, reproduction) illustrate the utility in applying genomic approaches for the successful management of introduced species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07553-x.

The inflated significance of neutral genetic diversity in conservation genetics

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

Genome Survey Sequencing of In Vivo Mother Plant and In Vitro Plantlets of Mikania cordata

Mikania cordata, the only native congener of the invasive weed Mikania micrantha in China, is an ideal species for comparative study to reveal the invasion mechanism. However, its genome resources are lagging far behind its congener, which limits the comparative genomic analysis. Our goal is to characterize the genome of M. cordata by next-generation sequencing and propose a scheme for long-read genome sequencing. Previous studies have shown that the genomic resources of the host plant would be affected by the endophytic microbial DNA. An aseptic sample of M. cordata will ensure the proper genome in downstream analysis. Because endophytes are ubiquitous in the greenhouse-grown M. cordata, the in vitro culture with cefotaxime or timentin treatment was undertaken to obtain the aseptic plantlets. The in vivo mother plant and in vitro plantlets were used to survey the genome. The microbial contamination in M. cordata was recognized by blast search and eliminated from the raw reads. The decontaminated sequencing reads were used to predict the genome size, heterozygosity, and repetitive rate. The in vivo plant was so contaminated that microbes occupied substantial sequencing resources and misled the scaffold assembly. Compared with cefotaxime, treatment with timentin performed better in cultivating robust in vitro plantlets. The survey result from the in vitro plantlets was more accurate due to low levels of contamination. The genome size was estimated to be 1.80 Gb with 0.50% heterozygosity and 78.35% repetitive rate. Additionally, 289,831 SSRs were identified in the genome. The genome is heavily contaminated and repetitive; therefore, the in vitro culture technique and long-read sequencing technology are recommended to generate a high-quality and highly contiguous genome.

Rapid evolution of genome-wide gene expression and plasticity during saline to freshwater invasions by the copepod Eurytemora affinis species complex

Saline migrants into freshwater habitats constitute among the most destructive invaders in aquatic ecosystems throughout the globe. However, the evolutionary and physiological mechanisms underlying such habitat transitions remain poorly understood. To explore the mechanisms of freshwater adaptation and distinguish between adaptive (evolutionary) and acclimatory (plastic) responses to salinity change, we examined genome-wide patterns of gene expression between ancestral saline and derived freshwater populations of the Eurytemora affinis species complex, reared under two different common-garden conditions (0 versus 15 PSU). We found that evolutionary shifts in gene expression (between saline and freshwater inbred lines) showed far greater changes and were more widespread than acclimatory responses to salinity (0 versus 15 PSU). Most notably, 30-40 genes showing evolutionary shifts in gene expression across the salinity boundary were associated with ion transport function, with inorganic cation transmembrane transport forming the largest Gene Ontology category. Of particular interest was the sodium transporter, the Na⁺ /H⁺ antiporter (NHA) gene family, which was discovered in animals relatively recently. Thirty key ion regulatory genes, such as NHA paralogue #7, demonstrated concordant evolutionary and plastic shifts in gene expression, suggesting the evolution of ion transporter function and plasticity during rapid invasions into novel salinities. Moreover, freshwater invasions were associated with the evolution of reduced plasticity in the freshwater population, again for the same key ion transporters, consistent with the predicted evolution of canalization following adaptation to stressful conditions. Our results have important implications for understanding evolutionary and physiological mechanisms of range expansions by some of the most widespread invaders in aquatic habitats.

Genetic differentiation underlies seasonal variation in thermal tolerance, body size, and plasticity in a short-lived copepod

Organisms experience variation in the thermal environment on several different temporal scales, with seasonality being particularly prominent in temperate regions. For organisms with short generation times, seasonal variation is experienced across, rather than within, generations. How this affects the seasonal evolution of thermal tolerance and phenotypic plasticity is understudied, but has direct implications for the thermal ecology of these organisms. Here we document intra-annual patterns of thermal tolerance in two species of Acartia copepods (Crustacea) from a highly seasonal estuary, showing strong variation across the annual temperature cycle. Common garden, split-brood experiments indicate that this seasonal variation in thermal tolerance, along with seasonal variation in body size and phenotypic plasticity, is likely affected by genetic polymorphism. Our results show that adaptation to seasonal variation is important to consider when predicting how populations may respond to ongoing climate change.

Salinity tolerance and geographical origin predict global alien amphipod invasions

Invasive alien species are driving global biodiversity loss, compromising ecosystem function and service provision, and human, animal and plant health. Habitat characteristics and geographical origin may predict invasion success, and in aquatic environments could be mediated principally by salinity tolerance. Crustacean invaders are causing global problems and we urgently require better predictive power of their invasiveness. Here, we compiled global aquatic gammarid (Crustacea: Amphipoda: Gammaroidea) diversity and examined their salinity tolerances and regions of origin to test whether these factors predict invasion success. Across 918 aquatic species within this superfamily, relatively few gammarids (n = 27, 3%) were reported as aliens, despite extensive invasion opportunities and high numbers of published studies on amphipod invasions. However, reported alien species were disproportionately salt-tolerant (i.e. 32% of brackish-water species), with significantly lower proportions of aliens originating from freshwater and marine environments (both 1%). Alien gammarids also significantly disproportionally originated from the Ponto-Caspian (20% of these taxa) when compared with all ‘other' grouped regions (1%), and principally invaded Eurasian waters, with translocations of salt-tolerant taxa to freshwaters being pervasive. This suggests habitat characteristics, alongside regional contexts, help predict invasibility. In particular, broad environmental tolerances to harsh environments and associated evolutionary history probably promote success of aliens globally.