A genetic perspective on rapid evolution in cane toads (Rhinella marina)

Rhinella marina oocytes: a suitable alternative expression system for functional characterization of aquaglyceroporins

Amphibian oocytes have been extensively used for heterologous expression of membrane proteins for studying their biochemical and biophysical properties. So far, Xenopus laevis is the main amphibian used as oocytes source to express aquaglyceroporins in order to assess water and solutes permeability. However, this well-established amphibian model represents a threat to the biodiversity in many countries, especially in those from tropical regions. For that reason, the import of Xenopus laevis is subjected to strict control, which essentially has restricted its use in these regions. Therefore, a wider variety of expression systems for aquaglyceroporins is needed. Rhinella marina is extensively distributed in the Americas and its native range spreads from South America to Texas, US. Here we report the use of Rhinella marina oocytes as an alternative expression system for aquaglyceroporins and demonstrated its suitability to determine the permeability to water and non-ionic solutes. Rhinella marina oocytes were able to functionally express channels from human and the protozoan pathogen Trypanosoma brucei, two very distant organisms on the evolutionary scale. Permeability values obtained from Rhinella marina oocytes expressing members of aquaporin family were similar and comparable to those values reported in the literature for the same channels expressed in Xenopus laevis oocytes.

Improving amphibian genomic resources: a multitissue reference transcriptome of an iconic invader

Background

Cane toads (Rhinella marina) are an iconic invasive species introduced to 4 continents and well utilized for studies of rapid evolution in introduced environments. Despite the long introduction history of this species, its profound ecological impacts, and its utility for demonstrating evolutionary principles, genetic information is sparse. Here we produce a de novo transcriptome spanning multiple tissues and life stages to enable investigation of the genetic basis of previously identified rapid phenotypic change over the introduced range.

Findings

Using approximately 1.9 billion reads from developing tadpoles and 6 adult tissue-specific cDNA libraries, as well as a transcriptome assembly pipeline encompassing 100 separate de novo assemblies, we constructed 62 202 transcripts, of which we functionally annotated ∼50%. Our transcriptome assembly exhibits 90% full-length completeness of the Benchmarking Universal Single-Copy Orthologs data set. Robust assembly metrics and comparisons with several available anuran transcriptomes and genomes indicate that our cane toad assembly is one of the most complete anuran genomic resources available.

Conclusions

This comprehensive anuran transcriptome will provide a valuable resource for investigation of genes under selection during invasion in cane toads, but will also greatly expand our general knowledge of anuran genomes, which are underrepresented in the literature. The data set is publically available in NCBI and GigaDB to serve as a resource for other researchers.

Draft genome assembly of the invasive cane toad, Rhinella marina

Background

The cane toad (Rhinella marina formerly Bufo marinus) is a species native to Central and South America that has spread across many regions of the globe. Cane toads are known for their rapid adaptation and deleterious impacts on native fauna in invaded regions. However, despite an iconic status, there are major gaps in our understanding of cane toad genetics. The availability of a genome would help to close these gaps and accelerate cane toad research.

Findings

We report a draft genome assembly for R. marina, the first of its kind for the Bufonidae family. We used a combination of long-read Pacific Biosciences RS II and short-read Illumina HiSeq X sequencing to generate 359.5 Gb of raw sequence data. The final hybrid assembly of 31,392 scaffolds was 2.55 Gb in length with a scaffold N50 of 168 kb. BUSCO analysis revealed that the assembly included full length or partial fragments of 90.6% of tetrapod universal single-copy orthologs (n = 3950), illustrating that the gene-containing regions have been well assembled. Annotation predicted 25,846 protein coding genes with similarity to known proteins in Swiss-Prot. Repeat sequences were estimated to account for 63.9% of the assembly.

Conclusions

The R. marina draft genome assembly will be an invaluable resource that can be used to further probe the biology of this invasive species. Future analysis of the genome will provide insights into cane toad evolution and enrich our understanding of their interplay with the ecosystem at large.

Spinal arthritis in cane toads across the Australian landscape

Loss of fitness can be a consequence of selection for rapid dispersal ability in invasive species. Increased prevalence of spinal arthritis may occur in cane toad populations at the invasion front as a cost of increased invasiveness, but our knowledge of the ecological drivers of this condition is lacking. We aimed to determine the factors explaining the prevalence of spinal arthritis in populations across the Australian landscape. We studied populations across a gradient of invasion histories. We collected 2415 toads over five years and determined the presence and size of spondylosis for each individual. We examined the effect of host size, leg length and invasion history on the prevalence of spondylosis. Host size was a significant predictor of spondylosis across populations. Contrary to our expectation, the overall prevalence of spondylosis was not positively related to invasion history and did not correlate with toad relative leg length. Rather than invasion age, the latitude at which populations were sampled provided an alternate explanation for the prevalence of spondylosis in cane toad populations and suggested that the incidence of this condition did not increase as a physiological cost of invasion, but is instead related to physical variables, such as climate.

Benchmarking differential expression analysis tools for RNA-Seq: normalization-based vs. log-ratio transformation-based methods

BACKGROUND

Count data generated by next-generation sequencing assays do not measure absolute transcript abundances. Instead, the data are constrained to an arbitrary "library size" by the sequencing depth of the assay, and typically must be normalized prior to statistical analysis. The constrained nature of these data means one could alternatively use a log-ratio transformation in lieu of normalization, as often done when testing for differential abundance (DA) of operational taxonomic units (OTUs) in 16S rRNA data. Therefore, we benchmark how well the ALDEx2 package, a transformation-based DA tool, detects differential expression in high-throughput RNA-sequencing data (RNA-Seq), compared to conventional RNA-Seq methods such as edgeR and DESeq2.

RESULTS

To evaluate the performance of log-ratio transformation-based tools, we apply the ALDEx2 package to two simulated, and two real, RNA-Seq data sets. One of the latter was previously used to benchmark dozens of conventional RNA-Seq differential expression methods, enabling us to directly compare transformation-based approaches. We show that ALDEx2, widely used in meta-genomics research, identifies differentially expressed genes (and transcripts) from RNA-Seq data with high precision and, given sufficient sample sizes, high recall too (regardless of the alignment and quantification procedure used). Although we show that the choice in log-ratio transformation can affect performance, ALDEx2 has high precision (i.e., few false positives) across all transformations. Finally, we present a novel, iterative log-ratio transformation (now implemented in ALDEx2) that further improves performance in simulations.

CONCLUSIONS

Our results suggest that log-ratio transformation-based methods can work to measure differential expression from RNA-Seq data, provided that certain assumptions are met. Moreover, these methods have very high precision (i.e., few false positives) in simulations and perform well on real data too. With previously demonstrated applicability to 16S rRNA data, ALDEx2 can thus serve as a single tool for data from multiple sequencing modalities.

Testing for Short- and Long-Term Thermal Plasticity in Corticosterone Responses of an Ectothermic Vertebrate

Phenotypic plasticity, broadly defined as the capacity of one genotype to produce more than one phenotype, is a key mechanism for how animals adapt to environmental (including thermal) variation. Vertebrate glucocorticoid hormones exert broad-scale regulation of physiological, behavioral, and morphological traits that influence fitness under many life-history or environmental contexts. Yet the capacity for vertebrates to demonstrate different types of thermal plasticity, including rapid compensation or longer acclimation in glucocorticoid hormone function, when subject to different environmental temperature regimes remains poorly addressed. Here, we explore whether patterns of urinary corticosterone metabolites respond (i.e., evidence of acclimation) to repeated short-term and sustained long-term temperature exposures in an amphibian, the cane toad (Rhinella marina). In response to three repeated short (30-min) high-temperature (37°C) exposures (at 10-d intervals), toads produced urinary corticosterone metabolite responses of sequentially greater magnitude, relative to controls. However, toads subjected to 4 wk of acclimation to either cool (18°C)- or warm (30°C)-temperature environments did not differ significantly in their urinary corticosterone metabolite responses during exposure to a thermal ramp (18°-36°C). Together, these results indicate that adult toads had different, including limited, capacities for their glucocorticoid responses to demonstrate plasticity to different regimes of environmental temperature variation. We advocate further research as necessary to identify plasticity, or lack thereof, in glucocorticoid physiology, to better understand how vertebrates can regulate organismal responses to environmental variation.

Rapid adaptive response to a Mediterranean environment reduces phenotypic mismatch in a recent amphibian invader

Invasive species frequently cope with ecological conditions that are different from those to which they adapted, presenting an opportunity to investigate how phenotypes change across short time scales. In 2000, the guttural toad Sclerophrys gutturalis was first detected in a peri-urban area of Cape Town, where it is now invasive. The ability of the species to invade Cape Town is surprising as the area is characterized by a Mediterranean climate significantly drier and colder than that of the native source area. We measured field hydration state of guttural toads from the invasive Cape Town population and a native source population from Durban. We also obtained from laboratory trials: rates of evaporative water loss and water uptake, sensitivity of locomotor endurance to hydration state, critical thermal minimum (CT_min) and sensitivity of CT_min to hydration state. Field hydration state of invasive toads was significantly lower than that of native toads. Although the two populations had similar rates of water loss and uptake, invasive toads were more efficient in minimizing water loss through postural adjustments. In locomotor trials, invasive individuals noticeably outperformed native individuals when dehydrated but not when fully hydrated. CT_min was lower in invasive individuals than in native individuals, independent of hydration state. Our results indicate that an invasive population that is only 20 years old shows adaptive responses that reduce phenotypic mismatch with the novel environment. The invasion potential of the species in Cape Town is higher than we could infer from its characteristics in the native source population.

New Weapons in the Toad Toolkit: A Review of Methods to Control and Mitigate the Biodiversity Impacts of Invasive Cane Toads (Rhinella Marina)

Our best hope of developing innovative methods to combat invasive species is likely to come from the study of high-profile invaders that have attracted intensive research not only into control, but also basic biology. Here we illustrate that point by reviewing current thinking about novel ways to control one of the world’s most well-studied invasions: that of the cane toad in Australia. Recently developed methods for population suppression include more effective traps based on the toad’s acoustic and pheromonal biology. New tools for containing spread include surveillance technologies (e.g., eDNA sampling and automated call detectors), as well as landscape-level barriers that exploit the toad’s vulnerability to desiccation—a strategy that could be significantly enhanced through the introduction of sedentary, range-core genotypes ahead of the invasion front. New methods to reduce the ecological impacts of toads include conditioned taste aversion in free-ranging predators, gene banking, and targeted gene flow. Lastly, recent advances in gene editing and gene drive technology hold the promise of modifying toad phenotypes in ways that may facilitate control or buffer impact. Synergies between these approaches hold great promise for novel and more effective means to combat the toad invasion and its consequent impacts on biodiversity.

The thermal dependency of locomotor performance evolves rapidly within an invasive species

Biological invasions can stimulate rapid shifts in organismal performance, via both plasticity and adaptation. We can distinguish between these two proximate mechanisms by rearing offspring from populations under identical conditions and measuring their locomotor abilities in standardized trials. We collected adult cane toads (Rhinella marina) from invasive populations that inhabit regions of Australia with different climatic conditions. We bred those toads and raised their offspring under common‐garden conditions before testing their locomotor performance. At high (but not low) temperatures, offspring of individuals from a hotter location (northwestern Australia) outperformed offspring of conspecifics from a cooler location (northeastern Australia). This disparity indicates that, within less than 100 years, thermal performance in cane toads has adapted to the novel abiotic challenges that cane toads have encountered during their invasion of tropical Australia.

Physiological plasticity in a successful invader: rapid acclimation to cold occurs only in cool-climate populations of cane toads (Rhinella marina)

Physiological plasticity may facilitate invasion of novel habitats; but is such plasticity present in all populations of the invader or is it elicited only by specific climatic challenges? In cold-climate areas of Australia, invasive cane toads (Rhinella marina) can rapidly acclimate to cool conditions. To investigate whether this physiological plasticity is found in all invasive cane toads or is only seen in cool climates, we measured the acclimation ability of toads from across Australia and the island of Hawai'i. We collected toads from the field and placed them at either 12 or 24°C for 12 h before measuring their righting response as a proxy for critical thermal minimum (CTmin). Toads from the coolest Australian region (New South Wales) demonstrated plasticity (as previously reported), with exposure to 12°C (vs. 24°C) decreasing CTmin by 2°C. In toads from other Australian populations, CTmins were unaffected by our thermal treatments. Hawai'ian toads from a cool, wet site also rapidly acclimated to cool conditions, whereas those from warmer and drier Hawai'ian sites did not. Thermal plasticity has diverged among populations of invasive cane toads, with rapid acclimation manifested only in two cool-climate populations from widely separated sites. Predictions about the potential range of invasive species thus must consider the possibility of geographic (intraspecific) heterogeneity in thermal plasticity; data from other parts of the species' range may fail to predict levels of plasticity elicited by thermal challenges.

Competition and pesticide exposure affect development of invasive (Rhinella marina) and native (Fejervarya vittigera) rice paddy amphibian larvae

Increased pesticide use in rice agricultural ecosystems may alter competitive interactions between invasive and native amphibian species. We conducted an experiment with two rice paddy amphibians found in Luzon, Philippines, the invasive cane toad (Rhinella marina) and the endemic Luzon wart frog (Fejervarya vittigera), to determine whether exposure to a common herbicide, butachlor, drives competitive interactions in favor of the invasive amphibian. Our results revealed that competition had a strong effect on the development of both species, but in opposing directions; Luzon wart frog tadpoles were smaller and developed slower than when raised alone, whereas cane toad tadpoles were larger and developed faster. Contrary to our predictions, development and survival of endemic wart frog tadpoles was not affected by butachlor, whereas invasive cane toad tadpoles were affected across several endpoints including gene expression, body size, and survival. We also observed an interaction between pesticide exposure and competition for the cane toad, where survival declined but body size and expression of thyroid sensitive genes increased. Taken together, our findings indicate that the success of the cane toad larvae in rice fields may be best explained by increased rates of development and larger body sizes of tadpoles in response to competition with native Luzon wart frog tadpoles rather than lower sensitivity to a common pesticide. Our results for the cane toad also provide evidence that butachlor can disrupt thyroid hormone mediated development in amphibians, and further demonstrate that important species interactions such as competition can be affected by pesticide exposure in aquatic ecosystems.

propr: An R-package for Identifying Proportionally Abundant Features Using Compositional Data Analysis

In the life sciences, many assays measure only the relative abundances of components in each sample. Such data, called compositional data, require special treatment to avoid misleading conclusions. Awareness of the need for caution in analyzing compositional data is growing, including the understanding that correlation is not appropriate for relative data. Recently, researchers have proposed proportionality as a valid alternative to correlation for calculating pairwise association in relative data. Although the question of how to best measure proportionality remains open, we present here a computationally efficient R package that implements three measures of proportionality. In an effort to advance the understanding and application of proportionality analysis, we review the mathematics behind proportionality, demonstrate its application to genomic data, and discuss some ongoing challenges in the analysis of relative abundance data.

propr: An R-package for Identifying Proportionally Abundant Features Using Compositional Data Analysis

In the life sciences, many assays measure only the relative abundances of components in each sample. Such data, called compositional data, require special treatment to avoid misleading conclusions. Awareness of the need for caution in analyzing compositional data is growing, including the understanding that correlation is not appropriate for relative data. Recently, researchers have proposed proportionality as a valid alternative to correlation for calculating pairwise association in relative data. Although the question of how to best measure proportionality remains open, we present here a computationally efficient R package that implements three measures of proportionality. In an effort to advance the understanding and application of proportionality analysis, we review the mathematics behind proportionality, demonstrate its application to genomic data, and discuss some ongoing challenges in the analysis of relative abundance data.

Effects of invasion history on physiological responses to immune system activation in invasive Australian cane toads

The cane toad (Rhinella marina) has undergone rapid evolution during its invasion of tropical Australia. Toads from invasion front populations (in Western Australia) have been reported to exhibit a stronger baseline phagocytic immune response than do conspecifics from range core populations (in Queensland). To explore this difference, we injected wild-caught toads from both areas with the experimental antigen lipopolysaccharide (LPS, to mimic bacterial infection) and measured whole-blood phagocytosis. Because the hypothalamic-pituitary-adrenal axis is stimulated by infection (and may influence immune responses), we measured glucocorticoid response through urinary corticosterone levels. Relative to injection of a control (phosphate-buffered saline), LPS injection increased both phagocytosis and the proportion of neutrophils in the blood. However, responses were similar in toads from both populations. This null result may reflect the ubiquity of bacterial risks across the toad’s invaded range; utilization of this immune pathway may not have altered during the process of invasion. LPS injection also induced a reduction in urinary corticosterone levels, perhaps as a result of chronic stress.

Genetics of dispersal

Dispersal is a process of central importance for the ecological and evolutionary dynamics of populations and communities, because of its diverse consequences for gene flow and demography. It is subject to evolutionary change, which begs the question, what is the genetic basis of this potentially complex trait? To address this question, we (i) review the empirical literature on the genetic basis of dispersal, (ii) explore how theoretical investigations of the evolution of dispersal have represented the genetics of dispersal, and (iii) discuss how the genetic basis of dispersal influences theoretical predictions of the evolution of dispersal and potential consequences.

Dispersal has a detectable genetic basis in many organisms, from bacteria to plants and animals. Generally, there is evidence for significant genetic variation for dispersal or dispersal‐related phenotypes or evidence for the micro‐evolution of dispersal in natural populations. Dispersal is typically the outcome of several interacting traits, and this complexity is reflected in its genetic architecture: while some genes of moderate to large effect can influence certain aspects of dispersal, dispersal traits are typically polygenic. Correlations among dispersal traits as well as between dispersal traits and other traits under selection are common, and the genetic basis of dispersal can be highly environment‐dependent.

By contrast, models have historically considered a highly simplified genetic architecture of dispersal. It is only recently that models have started to consider multiple loci influencing dispersal, as well as non‐additive effects such as dominance and epistasis, showing that the genetic basis of dispersal can influence evolutionary rates and outcomes, especially under non‐equilibrium conditions. For example, the number of loci controlling dispersal can influence projected rates of dispersal evolution during range shifts and corresponding demographic impacts. Incorporating more realism in the genetic architecture of dispersal is thus necessary to enable models to move beyond the purely theoretical towards making more useful predictions of evolutionary and ecological dynamics under current and future environmental conditions. To inform these advances, empirical studies need to answer outstanding questions concerning whether specific genes underlie dispersal variation, the genetic architecture of context‐dependent dispersal phenotypes and behaviours, and correlations among dispersal and other traits.

ORIGINS AND GENETIC DIVERSITY OF INTRODUCED POPULATIONS OF THE PUERTO RICAN RED-EYED COQUÍ, ELEUTHERODACTYLUS ANTILLENSIS, IN SAINT CROIX (U.S. VIRGIN ISLANDS) AND PANAMÁ

The Red-eyed Coquí, Eleutherodactylus antillensis, is a terrestrial frog endemic to the Puerto Rican Bank (Puerto Rico and numerous islands and cays off its eastern coast), in the eastern Caribbean Sea. The species was likely introduced in Saint Croix, an island c. 100 km southeast of Puerto Rico, in the late 1930s, and in Panamá City, Panamá, in the late 1950s or early 1960s, but the source(s) of these introductions are unknown. We analyzed sequence data from one mtDNA locus and four nuDNA introns to infer the origin(s) of the Saint Croix and Panamá City populations and quantify their genetic diversity. Saint Croix and Panamanian populations do not share any haplotypes, and they cluster with different native populations, suggesting that they are derived from separate sources in the Puerto Rican Bank. Patterns of population structure trace the probable sources of E. antillensis in Saint Croix to islands off Puerto Rico's eastern coast, which include Vieques, Culebra, Saint Thomas, Saint John, Tortola, and Virgin Gorda, and possibly to eastern Puerto Rico as well. In contrast, Panamá City E. antillensis probably originated from either western or eastern Puerto Rico. Genetic diversity in the introduced populations is similar to or lower than in populations in the species' native range, indicating that genetic diversity has not increased in the alien frogs. Our findings may facilitate the development of preventive measures to minimize introductions of non-native amphibians in the Caribbean and Central America.

De Novo Assembly of the Liver Transcriptome of the European Starling, Sturnus vulgaris

The European starling, Sturnus vulgaris, is a prolific and worldwide invasive species that also has served as an important model for avian ecological and invasion research. Although the genome sequence recently has become available, no transcriptome data have been published for this species. Here, we have sequenced and assembled the S. vulgaris liver transcriptome, which will provide a foundational resource for further annotation and validation of the draft genome. Moreover, it will be important for ecological and evolutionary studies investigating the genetic factors underlying rapid evolution and invasion success in this global invader.

Costs of immunity and their role in the range expansion of the house sparrow in Kenya

There are at least two reasons to study traits that mediate successful range expansions. First, dispersers will found new populations and thus impact the distribution and evolution of species. Second, organisms moving into new areas will influence the fate of resident communities, directly competing with or indirectly affecting residents by spreading non-native or spilling-back native parasites. The success of invaders in new areas is likely mediated by a counterbalancing of costly traits. In new areas where threats are comparatively rare, individuals that grow rapidly and breed prolifically should be at an advantage. High investment in defenses should thus be disfavored. In the present study, we compared the energetic, nutritional and collateral damage costs of an inflammatory response among Kenyan house sparrow (Passer domesticus) populations of different ages, asking whether costs were related to traits of individuals from three different capture sites. Kenya is among the world's most recent range expansions for this species, and we recently found that the expression of Toll-like receptors (TLRs), leukocyte receptors that instigate inflammatory responses when bound to microbial elements, was related to the range expansion across the country. Here, we found (contrary to our expectations) that energetic and nutritional costs of inflammation were higher, but damage costs were lower, in range-edge compared to core birds. Moreover, at the individual level, TLR-4 expression was negatively related to commodity costs (energy and a critical amino acid) of inflammation. Our data thus suggest that costs of inflammation, perhaps mediated by TLR expression, might mitigate successful range expansions.

The things they carried: The pathogenic effects of old and new parasites following the intercontinental invasion of the Australian cane toad (Rhinella marina)

Brought to Australia in 1935 to control agricultural pests (from French Guiana, via Martinique, Barbados, Jamaica, Puerto Rico and Hawai'i), repeated stepwise translocations of small numbers of founders enabled the cane toad (Rhinella marina) to escape many parasites and pathogens from its native range. However, the infective organisms that survived the journey continue to affect the dynamics of the toad in its new environment. In Australia, the native-range lungworm Rhabdias pseudosphaerocephala decreases its host's cardiac capacity, as well as growth and survival, but not rate of dispersal. The lungworm is most prevalent in long-colonised areas within the toads' Australian range, and absent from the invasion front. Several parasites and pathogens of Australian taxa have host-shifted to cane toads in Australia; for example, invasion-front toads are susceptible to spinal arthritis caused by the soil bacterium, Ochrobactrum anthropi. The pentastome Raillietiella frenata has host-shifted to toads and may thereby expand its Australian range due to the continued range expansion of the invasive toads. Spill-over and spill-back of parasites may be detrimental to other host species; however, toads may also reduce parasite loads in native taxa by acting as terminal hosts. We review the impact of the toad's parasites and pathogens on the invasive anuran's biology in Australia, as well as collateral effects of toad-borne parasites and pathogens on other host species in Australia. Both novel and co-evolved pathogens and parasites may have played significant roles in shaping the rapid evolution of immune system responses in cane toads within their invaded range.

•

Invasive cane toads have lost many parasites due to serial translocations.

•

One native lungworm (Rhabdias pseudosphaerocephala) has been retained.

•

Toads have also acquired novel parasites and pathogens from Australian hosts.

•

Toads either amplify parasite numbers or act as a parasite sink.

•

Differences in immune function exist between toad populations within Australia.

It is lonely at the front: contrasting evolutionary trajectories in male and female invaders

Invasive species often exhibit rapid evolutionary changes, and can provide powerful insights into the selective forces shaping phenotypic traits that influence dispersal rates and/or sexual interactions. Invasions also may modify sexual dimorphism. We measured relative lengths of forelimbs and hindlimbs of more than 3000 field-caught adult cane toads (Rhinella marina) from 67 sites in Hawai'i and Australia (1-80 years post-colonization), along with 489 captive-bred individuals from multiple Australian sites raised in a 'common garden' (to examine heritability and reduce environmental influences on morphology). As cane toads spread from east to west across Australia, the ancestral condition (long limbs, especially in males) was modified. Limb length relative to body size was first reduced (perhaps owing to natural selection on locomotor ability), but then increased again (perhaps owing to spatial sorting) in the invasion vanguard. In contrast, the sex disparity in relative limb length has progressively decreased during the toads' Australian invasion. Offspring reared in a common environment exhibited similar geographical divergences in morphology as did wild-caught animals, suggesting a genetic basis to the changes. Limb dimensions showed significant heritability (2-17%), consistent with the possibility of an evolved response. Cane toad populations thus have undergone a major shift in sexual dimorphism in relative limb lengths during their brief (81 years) spread through tropical Australia.

Population structure and dispersal routes of an invasive parasite, Fascioloides magna, in North America and Europe

Background

Fascioloides magna (Trematoda: Fasciolidae) is an important liver parasite of a wide range of free-living and domestic ruminants; it represents a remarkable species due to its large spatial distribution, invasive character, and potential to colonize new territories. The present study provides patterns of population genetic structure and admixture in F. magna across all enzootic regions in North America and natural foci in Europe, and infers migratory routes of the parasite on both continents.

Methods

In total, 432 individuals from five North American enzootic regions and three European foci were analysed by 11 microsatellite loci. Genetic data were evaluated by several statistical approaches: (i) the population genetic structure of F. magna was inferred using program STRUCTURE; (ii) the genetic interrelationships between populations were analysed by PRINCIPAL COORDINATES ANALYSIS; and (iii) historical dispersal routes in North America and recent invasion routes in Europe were explored using MIGRATE.

Results

The analysis of dispersal routes of the parasite in North America revealed west-east and south-north lineages that partially overlapped in the central part of the continent, where different host populations historically met. The exact origin of European populations of F. magna and their potential translocation routes were determined. Flukes from the first European focus, Italy, were related to F. magna from northern Pacific coast, while parasites from the Czech focus originated from south-eastern USA, particularly South Carolina. The Danube floodplain forests (third and still expanding focus) did not display relationship with any North American population; instead the Czech origin of the Danube population was indicated. A serial dilution of genetic diversity along the dispersion route across central and eastern Europe was observed. The results of microsatellite analyses were compared to previously acquired outputs from mitochondrial haplotype data and correlated with past human-directed translocations and natural migration of the final cervid hosts of F. magna.

Conclusions

The present study revealed a complex picture of the population genetic structure and interrelationships of North American and European populations, global distribution and migratory routes of F. magna and an origin of European foci.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1811-z) contains supplementary material, which is available to authorized users.

Multiple dispersal vectors drive range expansion in an invasive marine species

The establishment and subsequent spread of invasive species is widely recognized as one of the most threatening processes contributing to global biodiversity loss. This is especially true for marine and estuarine ecosystems, which have experienced significant increases in the number of invasive species with the increase in global maritime trade. Understanding the rate and mechanisms of range expansion is therefore of significant interest to ecologists and conservation managers alike. Using a combination of population genetic surveys, environmental DNA (eDNA) plankton sampling and hydrodynamic modelling, we examined the patterns of introduction of the predatory Northern Pacific seastar (Asterias amurensis) and pathways of secondary spread within southeast Australia. Genetic surveys across the invasive range reveal some genetic divergence between the two main invasive regions and no evidence of ongoing gene flow, a pattern that is consistent with the establishment of the second invasive region via a human-mediated translocation event. In contrast, hydrodynamic modelling combined with eDNA plankton sampling demonstrated that the establishment of range expansion populations within a region is consistent with natural larval dispersal and recruitment. Our results suggest that both anthropogenic and natural dispersal vectors have played an important role in the range expansion of this species in Australia. The multiple modes of spread combined with high levels of fecundity and a long larval duration in A. amurensis suggests it is likely to continue its range expansion and significantly impact Australian marine ecosystems.

Mixed population genomics support for the central marginal hypothesis across the invasive range of the cane toad (Rhinella marina) in Australia

Understanding factors that cause species' geographic range limits is a major focus in ecology and evolution. The central marginal hypothesis (CMH) predicts that species cannot adapt to conditions beyond current geographic range edges because genetic diversity decreases from core to edge due to smaller, more isolated edge populations. We employed a population genomics framework using 24 235-33 112 SNP loci to test major predictions of the CMH in the ongoing invasion of the cane toad (Rhinella marina) in Australia. Cane toad tissue samples were collected along broad-scale, core-to-edge transects across their invasive range. Geographic and ecological core areas were identified using GIS and habitat suitability indices from ecological niche modelling. Bayesian clustering analyses revealed three genetic clusters, in the northwest invasion-front region, northeast precipitation-limited region and southeast cold temperature-limited region. Core-to-edge patterns of genetic diversity and differentiation were consistent with the CMH in the southeast, but were not supported in the northeast and showed mixed support in the northwest. Results suggest cold temperatures are a likely contributor to southeastern range limits, consistent with CMH predictions. In the northeast and northwest, ecological processes consisting of a steep physiological barrier and ongoing invasion dynamics, respectively, are more likely explanations for population genomic patterns than the CMH.

Patterns of Genetic Variability in Island Populations of the Cane Toad (Rhinella marina) from the Mouth of the Amazon

The Amazonian coast has several unique geological characteristics resulting from the interaction between drainage pattern of the Amazon River and the Atlantic Ocean. It is one of the most extensive and sedimentologically dynamic regions of the world, with a large number of continental islands mostly formed less than 10,000 years ago. The natural distribution of the cane toad (Rhinella marina), one of the world's most successful invasive species, in this complex Amazonian system provides an intriguing model for the investigation of the effects of isolation or the combined effects of isolation and habitat dynamic changes on patterns of genetic variability and population differentiation. We used nine fast-evolving microsatellite loci to contrast patterns of genetic variability in six coastal (three mainlands and three islands) populations of the cane toad near the mouth of the Amazon River. Results from Bayesian multilocus clustering approach and Discriminant Analyses of Principal Component were congruent in showing that each island population was genetically differentiated from the mainland populations. All FST values obtained from all pairwise comparisons were significant, ranging from 0.048 to 0.186. Estimates of both recent and historical gene flow were not significantly different from zero across all population pairs, except the two mainland populations inhabiting continuous habitats. Patterns of population differentiation, with a high level of population substructure and absence/restricted gene flow, suggested that island populations of R. marina are likely isolated since the Holocene sea-level rise. However, considering the similar levels of genetic variability found in both island and mainland populations, it is reliable to assume that they were also isolated for longer periods. Given the genetic uniqueness of each cane toad population, together with the high natural vulnerability of the coastal regions and intense human pressures, we suggest that these populations should be treated as discrete units for conservation management purposes.

De Novo Assembly and Characterization of the Invasive Northern Pacific Seastar Transcriptome

Invasive species are a major threat to global biodiversity but can also serve as valuable model systems to examine important evolutionary processes. While the ecological aspects of invasions have been well documented, the genetic basis of adaptive change during the invasion process has been hampered by a lack of genomic resources for the majority of invasive species. Here we report the first larval transcriptomic resource for the Northern Pacific Seastar, Asterias amurensis, an invasive marine predator in Australia. Approximately 117.5 million 100 base-pair (bp) paired-end reads were sequenced from a single RNA-Seq library from a pooled set of full-sibling A. amurensis bipinnaria larvae. We evaluated the efficacy of a pre-assembly error correction pipeline on subsequent de novo assembly. Error correction resulted in small but important improvements to the final assembly in terms of mapping statistics and core eukaryotic genes representation. The error-corrected de novo assembly resulted in 115,654 contigs after redundancy clustering. 41,667 assembled contigs were homologous to sequences from NCBI’s non-redundant protein and UniProt databases. We assigned Gene Ontology, KEGG Orthology, Pfam protein domain terms and predicted protein-coding sequences to > 36,000 contigs. The final transcriptome dataset generated here provides functional information for 18,319 unique proteins, comprising at least 11,355 expressed genes. Furthermore, we identified 9,739 orthologs to P. miniata proteins, evaluated our annotation pipeline and generated a list of 150 candidate genes for responses to several environmental stressors that may be important for adaptation of A. amurensis in the invasive range. Our study has produced a large set of A. amurensis RNA contigs with functional annotations that can serve as a resource for future comparisons to other echinoderm transcriptomes and gene expression studies. Our data can be used to study the genetic basis of adaptive change and other important evolutionary processes during a successful invasion.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

An ideal weed: plasticity and invasiveness in Polygonum cespitosum

The introduced Asian plant Polygonum cespitosum has only recently become invasive in northeastern North America, spreading into sunny as well as shaded habitats. We present findings from a multiyear case study of this ongoing species invasion, drawing on field environmental measurements, glasshouse plasticity and resurrection experiments, and molecular genetic (microsatellite) data. We focus in particular on patterns of individual phenotypic plasticity (norms of reaction), their diversity within and among populations in the species' introduced range, and their contribution to its potential to evolve even greater invasiveness. Genotypes from introduced-range P. cespitosum populations have recently evolved to express greater adaptive plasticity to full sun and/or dry conditions without any loss of fitness in shade. Evidently, this species may evolve the sort of "general-purpose genotypes" hypothesized by Herbert Baker to characterize an "ideal weed." Indeed, we identified certain genotypes capable of extremely high reproductive output across contrasting conditions, including sunny, shaded, moist, and dry. Populations containing these high-performance genotypes had consistently higher fitness in all glasshouse habitats; there was no evidence for local adaptive differentiation among populations from sunny, shaded, moist, or dry sites. Norm of reaction data may provide valuable insights to invasion biology: the presence of broadly adaptive, high-performance genotypes can promote a species' ecological spread while providing the fuel for increased invasiveness to evolve.

What we still don't know about invasion genetics

Publication of The Genetics of Colonizing Species in 1965 launched the field of invasion genetics and highlighted the value of biological invasions as natural ecological and evolutionary experiments. Here, we review the past 50 years of invasion genetics to assess what we have learned and what we still don't know, focusing on the genetic changes associated with invasive lineages and the evolutionary processes driving these changes. We also suggest potential studies to address still-unanswered questions. We now know, for example, that rapid adaptation of invaders is common and generally not limited by genetic variation. On the other hand, and contrary to prevailing opinion 50 years ago, the balance of evidence indicates that population bottlenecks and genetic drift typically have negative effects on invasion success, despite their potential to increase additive genetic variation and the frequency of peak shifts. Numerous unknowns remain, such as the sources of genetic variation, the role of so-called expansion load and the relative importance of propagule pressure vs. genetic diversity for successful establishment. While many such unknowns can be resolved by genomic studies, other questions may require manipulative experiments in model organisms. Such studies complement classical reciprocal transplant and field-based selection experiments, which are needed to link trait variation with components of fitness and population growth rates. We conclude by discussing the potential for studies of invasion genetics to reveal the limits to evolution and to stimulate the development of practical strategies to either minimize or maximize evolutionary responses to environmental change.