Contrasting Patterns of Single Nucleotide Polymorphisms and Structural Variation Across Multiple Invasions.

Genetic divergence is the fundamental process that drives evolution and ultimately speciation. Structural variants (SVs) are large-scale genomic differences within a species or population and can cause functionally important phenotypic differences. Characterizing SVs across invasive species will fill knowledge gaps regarding how patterns of genetic diversity and genetic architecture shape rapid adaptation under new selection regimes. Here, we seek to understand patterns in genetic diversity within the globally invasive European starling, Sturnus vulgaris. Using whole genome sequencing of eight native United Kingdom (UK), eight invasive North America (NA), and 33 invasive Australian (AU) starlings, we examine patterns in genome-wide SNPs and SVs between populations and within Australia. Our findings detail the landscape of standing genetic variation across recently diverged continental populations of this invasive avian. We demonstrate that patterns of genetic diversity estimated from SVs do not necessarily reflect relative patterns from SNP data, either when considering patterns of diversity along the length of the organism's chromosomes (owing to enrichment of SVs in subtelomeric repeat regions), or interpopulation diversity patterns (possibly a result of altered selection regimes or introduction history). Finally, we find that levels of balancing selection within the native range differ across SNP and SV of different classes and outlier classifications. Overall, our results demonstrate that the processes that shape allelic diversity within populations is complex and support the need for further investigation of SVs across a range of taxa to better understand correlations between often well-studied SNP diversity and that of SVs.

Genetic and Phenotypic Consequences of Local Transitions between Sexual and Parthenogenetic Reproduction in the Wild.

AbstractTransitions from sexual to asexual reproduction have occurred in numerous lineages, but it remains unclear why asexual populations rarely persist. In facultatively parthenogenetic animals, all-female populations can arise when males are absent or become extinct, and such populations could help to understand the genetic and phenotypic changes that occur in the initial stages of transitions to asexuality. We investigated a naturally occurring spatial mosaic of mixed-sex and all-female populations of the facultatively parthenogenetic Australian phasmid Megacrania batesii. Analysis of single-nucleotide polymorphisms indicated multiple independent transitions between reproductive modes. All-female populations had much lower heterozygosity and allelic diversity than mixed-sex populations, but we found few consistent differences in fitness-related traits between population types. All-female populations exhibited more frequent and severe deformities in their (flight-incapable) wings but did not show higher rates of appendage loss. All-female populations also harbored more ectoparasites in swamp (but not beach) habitats. Reproductive mode explained little variation in female body size, fecundity, or egg hatch rate. Our results suggest that transitions to parthenogenetic reproduction can lead to dramatic genetic changes with little immediate effect on performance. All-female M. batesii populations appear to consist of high-fitness genotypes that might be able to thrive for many generations in relatively constant and benign environments but could be vulnerable to environmental challenges, such as increased parasite abundance.

A biological invasion modifies the dynamics of a host-parasite arms race.

By imposing novel selection pressures on both participants, biological invasions can modify evolutionary 'arms races' between hosts and parasites. A spatially replicated cross-infection experiment reveals strong spatial divergence in the ability of lungworms (Rhabdias pseudosphaerocephala) to infect invasive cane toads (Rhinella marina) in Australia. In areas colonized for longer than 20 years, toads are more resistant to infection by local strains of parasites than by allopatric strains. The situation reverses at the invasion front, where super-infective parasites have evolved. Invasion-induced shifts in genetic diversity and selective pressures may explain why hosts gain advantage over parasites in long-colonized areas, whereas parasites gain advantage at the invasion front.

Maternally derived avian corticosterone affects offspring genome-wide DNA methylation in a passerine species.

Avian embryos develop in an egg composition which reflects both maternal condition and the recent environment of their mother. In birds, yolk corticosterone (CORT) influences development by impacting pre- and postnatal growth, as well as nestling stress responses and development. One possible mechanism through which maternal CORT may affect offspring development is via changes to offspring DNA methylation. We sought to investigate this, for the first time in birds, by quantifying the impact of manipulations to maternal CORT on offspring DNA methylation. We non-invasively manipulated plasma CORT concentrations of egg-laying female zebra finches (Taeniopygia castanotis) with an acute dose of CORT administered around the time of ovulation and collected their eggs. We then assessed DNA methylation in the resulting embryonic tissue and in their associated vitelline membrane blood vessels, during early development (5 days after lay), using two established methods - liquid chromatography-mass spectrometry (LC-MS) and methylation-sensitive amplification fragment length polymorphism (MS-AFLP). LC-MS analysis showed that global DNA methylation was lower in embryos from CORT-treated mothers, compared to control embryos. In contrast, blood vessel DNA from eggs from CORT-treated mothers showed global methylation increases, compared to control samples. There was a higher proportion of global DNA methylation in the embryonic DNA of second clutches, compared to first clutches. Locus-specific analyses using MS-AFLP did not reveal a treatment effect. Our results indicate that an acute elevation of maternal CORT around ovulation impacts DNA methylation patterns in their offspring. This could provide a mechanistic understanding of how a mother's experience can affect her offspring's phenotype.

Brain transcriptome analysis reveals gene expression differences associated with dispersal behaviour between range-front and range-core populations of invasive cane toads in Australia.

Understanding the mechanisms allowing invasive species to adapt to novel environments is a challenge in invasion biology. Many invaders demonstrate rapid evolution of behavioural traits involved in range expansion such as locomotor activity, exploration and risk-taking. However, the molecular mechanisms that underpin these changes are poorly understood. In 86 years, invasive cane toads (Rhinella marina) in Australia have drastically expanded their geographic range westward from coastal Queensland to Western Australia. During their range expansion, toads have undergone extensive phenotypic changes, particularly in behaviours that enhance the toads' dispersal ability. Common-garden experiments have shown that some changes in behavioural traits related to dispersal are heritable. At the molecular level, it is currently unknown whether these changes in dispersal-related behaviour are underlain by small or large differences in gene expression, nor is known the biological function of genes showing differential expression. Here, we used RNA-seq to gain a better understanding of the molecular mechanisms underlying dispersal-related behavioural changes. We compared the brain transcriptomes of toads from the Hawai'ian source population, as well as three distinct populations from across the Australian invasive range. We found markedly different gene expression profiles between the source population and Australian toads. By contrast, toads from across the Australian invasive range had very similar transcriptomic profiles. Yet, key genes with functions putatively related to dispersal behaviour showed differential expression between populations located at each end of the invasive range. These genes could play an important role in the behavioural changes characteristic of range expansion in Australian cane toads.

Historical museum samples enable the examination of divergent and parallel evolution during invasion.

During the Anthropocene, Earth has experienced unprecedented habitat loss, native species decline and global climate change. Concurrently, greater globalization is facilitating species movement, increasing the likelihood of alien species establishment and propagation. There is a great need to understand what influences a species' ability to persist or perish within a new or changing environment. Examining genes that may be associated with a species' invasion success or persistence informs invasive species management, assists with native species preservation and sheds light on important evolutionary mechanisms that occur in novel environments. This approach can be aided by coupling spatial and temporal investigations of evolutionary processes. Here we use the common starling, Sturnus vulgaris, to identify parallel and divergent evolutionary change between contemporary native and invasive range samples and their common ancestral population. To do this, we use reduced-representation sequencing of native samples collected recently in northwestern Europe and invasive samples from Australia, together with museum specimens sampled in the UK during the mid-19th century. We found evidence of parallel selection on both continents, possibly resulting from common global selective forces such as exposure to pollutants. We also identified divergent selection in these populations, which might be related to adaptive changes in response to the novel environment encountered in the introduced Australian range. Interestingly, signatures of selection are equally as common within both invasive and native range contemporary samples. Our results demonstrate the value of including historical samples in genetic studies of invasion and highlight the ongoing and occasionally parallel role of adaptation in both native and invasive ranges.

Captivity induces large and population-dependent brain transcriptomic changes in wild-caught cane toads (Rhinella marina).

Gene expression levels are key molecular phenotypes at the interplay between genotype and environment. Mounting evidence suggests that short-term changes in environmental conditions, such as those encountered in captivity, can substantially affect gene expression levels. Yet, the exact magnitude of this effect, how general it is, and whether it results in parallel changes across populations are not well understood. Here, we take advantage of the well-studied cane toad, Rhinella marina, to examine the effect of short-term captivity on brain gene expression levels, and determine whether effects of captivity differ between long-colonized and vanguard populations of the cane toad's Australian invasion range. We compared the transcriptomes of wild-caught toads immediately assayed with those from toads captured from the same populations but maintained in captivity for seven months. We found large differences in gene expression levels between captive and wild-caught toads from the same population, with an over-representation of processes related to behaviour and the response to stress. Captivity had a much larger effect on both gene expression levels and gene expression variability in toads from vanguard populations compared to toads from long-colonized areas, potentially indicating an increased plasticity in toads at the leading edge of the invasion. Overall, our findings indicate that short-term captivity can induce large and population-specific transcriptomic changes, which has significant implications for studies comparing phenotypic traits of wild-caught organisms from different populations that have been held in captivity.

In an arms race between host and parasite, a lungworm's ability to infect a toad is determined by host susceptibility not parasite preference.

Evolutionary arms races can alter both parasite infectivity and host resistance, and it is difficult to separate the effects of these twin determinants of infection outcomes. We used a co-introduced, invasive host-parasite system (the lungworm Rhabdias pseudosphaerocephala and cane toads Rhinella marina), where rapid adaptation and dispersal have led to population differences in infection resistance. We quantified behavioural responses of parasite larvae to skin-chemical cues of toads from different invasive populations, and rates at which juvenile hosts became infected following standardized exposure to lungworms. Chemical cues from toad skin altered host-seeking behaviour by parasites, similarly among populations. The number of infection attempts (parasite larvae entering the host's body) also did not differ between populations, but rates of successful infection (establishment of adult worm in host lungs) were higher for range-edge toads than for range-core conspecifics. Thus, lower resistance to parasite infection in range-edge juvenile toads appears to be due to less effective immune defences of the host rather than differential behavioural responses of the parasite. In this ongoing host-parasite arms race, changing outcomes appear to be driven by shifts in host immunocompetence.

Signatures of selection in a recent invasion reveal adaptive divergence in a highly vagile invasive species.

A detailed understanding of population genetics in invasive populations helps us to identify drivers of successful alien introductions. Here, we investigate putative signals of selection in Australian populations of invasive common starlings, Sturnus vulgaris, and seek to understand how these have been influenced by introduction history. We used reduced representation sequencing to determine population structure, and identify Single Nucleotide Polymorphisms (SNPs) that are putatively under selection. We found that since their introduction into Australia, starling populations have become genetically differentiated despite the potential for high levels of dispersal, and that starlings have responded to selective pressures imposed by a wide range of environmental conditions across their geographic range. Isolation by distance appears to have played a strong role in determining genetic substructure across the starling's Australian range. Analyses of candidate SNPs that are putatively under selection indicated that aridity, precipitation and temperature may be important factors driving adaptive variation across the starling's invasive range in Australia. However, we also noted that the historic introduction regime may leave footprints on sites flagged as being under adaptive selection, and encourage critical interpretation of selection analyses in non-native populations.

Selection on Mitochondrial Variants Occurs between and within Individuals in an Expanding Invasion.

Mitochondria are critical for life, yet their underlying evolutionary biology is poorly understood. In particular, little is known about interaction between two levels of evolution: between individuals and within individuals (competition between cells, mitochondria or mitochondrial DNA molecules). Rapid evolution is suspected to occur frequently in mitochondrial DNA, whose maternal inheritance predisposes advantageous mutations to sweep rapidly though populations. Rapid evolution is also predicted in response to changed selection regimes after species invasion or removal of pathogens or competitors. Here, using empirical and simulated data from a model invasive bird species, we provide the first demonstration of rapid selection on the mitochondrial genome within individuals in the wild. Further, we show differences in mitochondrial DNA copy number associated with competing genetic variants, which may provide a mechanism for selection. We provide evidence for three rarely documented phenomena: selection associated with mitochondrial DNA abundance, selection on the mitochondrial control region, and contemporary selection during invasion.

Colour polymorphism is likely to be disadvantageous to some populations and species due to genetic architecture and morph interactions.

Polymorphism describes two or more distinct, genetically determined, phenotypes that co-occur in the same population, where the rarest morph is maintained at a frequency above the mutation rate (Ford 1945; Huxley 1955). In a recent opinion piece, we explored a new idea regarding the role of genetic architectures and morph interactions in colour polymorphisms and how this can negatively affect population performance (Bolton et al. 2015). In this issue of Molecular Ecology, Forsman (2016) thoroughly discusses the current evidence for polymorphisms enhancing population performance and critiques the validity of the definitions of polymorphism we use in our original paper. We respond by clarifying that the negative consequences of polymorphisms that we discussed are likely to be most pertinent in species that have a particular set of characteristics, such as strong sexual or social interactions between morphs and discrete genetic architectures. Although it was not our intention to redefine polymorphism, we do believe that there should be further discussion about refining or characterizing balanced polymorphisms with respect to the degree of morph sympatry, discreteness of traits and their underlying genetic architecture, and the types of selection that drive and maintain the variation. The latter describes whether polymorphism is primarily maintained by external factors such as predation pressure or internal factors such as interactions with members of the same species. The contribution of Forsman (2016) is useful to this discussion, and we hope that our exchange of opinions will inspire new empirical and theoretical ideas on the origin and maintenance of colour polymorphisms.

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

The process of biological invasion exposes a species to novel pressures, in terms of both the environments it encounters and the evolutionary consequences of range expansion. Several invaders have been shown to exhibit rapid evolutionary changes in response to those pressures, thus providing robust opportunities to clarify the processes at work during rapid phenotypic transitions. The accelerating pace of invasion of cane toads (Rhinella marina) in tropical Australia during its 80-year history has been well characterized at the phenotypic level, including common-garden experiments that demonstrate heritability of several dispersal-relevant traits. Individuals from the invasion front (and their progeny) show distinctive changes in morphology, physiology and behaviour that, in combination, result in far more rapid dispersal than is true of conspecifics from long-colonized areas. The extensive body of work on cane toad ecology enables us to place into context studies of the genetic basis of these traits. Our analyses of differential gene expression from toads from both ends of this invasion-history transect reveal substantial upregulation of many genes, notably those involved in metabolism and cellular repair. Clearly, then, the dramatically rapid phenotypic evolution of cane toads in Australia has been accompanied by substantial shifts in gene expression, suggesting that this system is well suited to investigating the genetic underpinnings of invasiveness.

The danger within: the role of genetic, behavioural and ecological factors in population persistence of colour polymorphic species.

Polymorphic species have been the focus of important work in evolutionary biology. It has been suggested that colour polymorphic species have specific evolutionary and population dynamics that enable them to persist through environmental changes better than less variable species. We suggest that recent empirical and theoretical work indicates that polymorphic species may be more vulnerable to extinction than previously thought. This vulnerability arises because these species often have a number of correlated sexual, behavioural, life history and ecological traits, which can have a simple genetic underpinning. When exacerbated by environmental change, these alternate strategies can lead to conflict between morphs at the genomic and population levels, which can directly or indirectly affect population and evolutionary dynamics. In this perspective, we identify a number of ways in which the nature of the correlated traits, their underpinning genetic architecture, and the inevitable interactions between colour morphs can result in a reduction in population fitness. The principles illustrated here apply to all kinds of discrete polymorphism (e.g. behavioural syndromes), but we focus primarily on colour polymorphism because they are well studied. We urge further empirical investigation of the genetic architecture and interactions in polymorphic species to elucidate the impact on population fitness.

The genome of a globally invasive passerine, the common myna, Acridotheres tristis.

In an era of global climate change, biodiversity conservation is receiving increased attention. Conservation efforts are greatly aided by genetic tools and approaches, which seek to understand patterns of genetic diversity and how they impact species health and their ability to persist under future climate regimes. Invasive species offer vital model systems in which to investigate questions regarding adaptive potential, with a particular focus on how changes in genetic diversity and effective population size interact with novel selection regimes. The common myna (Acridotheres tristis) is a globally invasive passerine and is an excellent model species for research both into the persistence of low-diversity populations and the mechanisms of biological invasion. To underpin research on the invasion genetics of this species, we present the genome assembly of the common myna. We describe the genomic landscape of this species, including genome wide allelic diversity, methylation, repeats, and recombination rate, as well as an examination of gene family evolution. Finally, we use demographic analysis to identify that some native regions underwent a dramatic population increase between the two most recent periods of glaciation, and reveal artefactual impacts of genetic bottlenecks on demographic analysis.

Is developmental plasticity triggered by DNA methylation changes in the invasive cane toad (Rhinella marina)?

Many organisms can adjust their development according to environmental conditions, including the presence of conspecifics. Although this developmental plasticity is common in amphibians, its underlying molecular mechanisms remain largely unknown. Exposure during development to either 'cannibal cues' from older conspecifics, or 'alarm cues' from injured conspecifics, causes reduced growth and survival in cane toad (Rhinella marina) tadpoles. Epigenetic modifications, such as changes in DNA methylation patterns, are a plausible mechanism underlying these developmental plastic responses. Here we tested this hypothesis, and asked whether cannibal cues and alarm cues trigger the same DNA methylation changes in developing cane toads. We found that exposure to both cannibal cues and alarm cues was associated with local changes in DNA methylation patterns. These DNA methylation changes affected genes putatively involved in developmental processes, but in different genomic regions for different conspecific-derived cues. Genetic background explains most of the epigenetic variation among individuals. Overall, the molecular mechanisms triggered by exposure to cannibal cues seem to differ from those triggered by alarm cues. Studies linking epigenetic modifications to transcriptional activity are needed to clarify the proximate mechanisms that regulate developmental plasticity in cane toads.

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.

Mothers adjust offspring sex to match the quality of the rearing environment.

Theory predicts that mothers should adjust offspring sex ratios when the expected fitness gains or rearing costs differ between sons and daughters. Recent empirical work has linked biased offspring sex ratios to environmental quality via changes in relative maternal condition. It is unclear, however, whether females can manipulate offspring sex ratios in response to environmental quality alone (i.e. independent of maternal condition). We used a balanced within-female experimental design (i.e. females bred on both low- and high-quality diets) to show that female parrot finches (Erythrura trichroa) manipulate primary offspring sex ratios to the quality of the rearing environment, and not to their own body condition and health. Individual females produced an unbiased sex ratio on high-quality diets, but over-produced sons in poor dietary conditions, even though they maintained similar condition between diet treatments. Despite the lack of sexual size dimorphism, such sex ratio adjustment is in line with predictions from sex allocation theory because nutritionally stressed foster sons were healthier, grew faster and were more likely to survive than daughters. These findings suggest that mothers may adaptively adjust offspring sex ratios to optimally match their offspring to the expected quality of the rearing environment.

Genetics and Plasticity Are Responsible for Ecogeographical Patterns in a Recent Invasion.

Patterns of covariation between phenotype and environment are presumed to be reflective of local adaptation, and therefore translate to a meaningful influence on an individual's overall fitness within that specific environment. However, these environmentally driven patterns may be the result of numerous and interacting processes, such as genetic variation, epigenetic variation, or plastic non-heritable variation. Understanding the relative importance of different environmental variables on underlying genetic patterns and resulting phenotypes is fundamental to understanding adaptation. Invasive systems are excellent models for such investigations, given their propensity for rapid evolution. This study uses reduced representation sequencing data paired with phenotypic data to examine whether important phenotypic traits in invasive starlings (Sturnus vulgaris) within Australia appear to be highly heritable (presumably genetic) or appear to vary with environmental gradients despite underlying genetics (presumably non-heritable plasticity). We also sought to determine which environmental variables, if any, play the strongest role shaping genetic and phenotypic patterns. We determined that environmental variables-particularly elevation-play an important role in shaping allelic trends in Australian starlings and may also reinforce neutral genetic patterns resulting from historic introduction regime. We examined a range of phenotypic traits that appear to be heritable (body mass and spleen mass) or negligibly heritable (e.g. beak surface area and wing length) across the starlings' Australian range. Using SNP variants associated with each of these phenotypes, we identify key environmental variables that correlate with genetic patterns, specifically that temperature and precipitation putatively play important roles shaping phenotype in this species. Finally, we determine that overall phenotypic variation is correlated with underlying genetic variation, and that these interact positively with the level of vegetation variation within a region, suggesting that ground cover plays an important role in shaping selection and plasticity of phenotypic traits within the starlings of Australia.

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.

Transcript- and annotation-guided genome assembly of the European starling.

The European starling, Sturnus vulgaris, is an ecologically significant, globally invasive avian species that is also suffering from a major decline in its native range. Here, we present the genome assembly and long-read transcriptome of an Australian-sourced European starling (S. vulgaris vAU), and a second, North American, short-read genome assembly (S. vulgaris vNA), as complementary reference genomes for population genetic and evolutionary characterization. S. vulgaris vAU combined 10× genomics linked-reads, low-coverage Nanopore sequencing, and PacBio Iso-Seq full-length transcript scaffolding to generate a 1050 Mb assembly on 6222 scaffolds (7.6 Mb scaffold N50, 94.6% busco completeness). Further scaffolding against the high-quality zebra finch (Taeniopygia guttata) genome assigned 98.6% of the assembly to 32 putative nuclear chromosome scaffolds. Species-specific transcript mapping and gene annotation revealed good gene-level assembly and high functional completeness. Using S. vulgaris vAU, we demonstrate how the multifunctional use of PacBio Iso-Seq transcript data and complementary homology-based annotation of sequential assembly steps (assessed using a new tool, saaga) can be used to assess, inform, and validate assembly workflow decisions. We also highlight some counterintuitive behaviour in traditional busco metrics, and present buscomp, a complementary tool for assembly comparison designed to be robust to differences in assembly size and base-calling quality. This work expands our knowledge of avian genomes and the available toolkit for assessing and improving genome quality. The new genomic resources presented will facilitate further global genomic and transcriptomic analysis on this ecologically important species.