Search

An understanding of the basic mechanisms responsible for the pathogenesis of liver neoplasms is needed in order to develop better therapeutic strategies. The present study utilized a pharmacogenetic mouse model to assess the role of cytochrome P4501A1 (Cyp1a1) in modulating genetic damage to oncogenic and tumor suppressor loci following in utero exposure to the polycyclic aromatic hydrocarbon, 3-methylcholanthrene (MC). Analysis of the Ha-ras, Ki-ras, INK4a and p53 genes was carried out with lysates from paraffin-embedded liver tissue from transplacentally-treated mice. The lysates were subjected to DNA amplification by the PCR technique followed by allele-specific oligonucleotide hybridization screening and SSCP analysis. All of the 26 neoplasms screened (23 hepatocellular carcinomas, two hepatocellular adenomas and one sarcoma) exhibited a GGC→CGC (GLY13→ARG13) transversion at the Ki-ras gene locus. None of the tumors had Ki-ras mutations at codon 12 of exon 1. Approximately 12% (3/26) of the liver tumors exhibited point mutations in exon 1 of the INK4a gene, with each of the three tumors exhibiting two point mutations. Analysis of exon 2 of the INK4a gene showed the presence of a CCG→CTG (PRO73→LEU73) transition in two of the 26 neoplasms. No mutations were found in exons 1 or 2 of the Ha-ras gene, or in exons 5-8 of the p53 gene. Analysis of tumor RNAs showed overexpression of Ha-ras, cip1 and c-jun in 38% of the liver tumor samples. The results of this study suggest that mutagenic damage to oncogenes and tumor suppressor genes may be critical factors in mediating transplacentally-induced liver tumorigenesis. The fact that Ki-ras mutations were found in all of the tumors suggests that mutation at this gene locus may be an early event in liver tumor pathogenesis, while mutation in tumor suppressor genes may occur later during tumor progression. These combined results are consistent with the pathogenesis of cancer in humans. [ABSTRACT FROM PUBLISHER]

Abstract: The remarkable efficacy of Epstein-Barr virus (EBV)-specific T cells for the treatment of posttransplant lymphomas has not been reproduced for EBV-positive (EBV+) malignancies outside the transplant setting. This is because of, in part, the heterogeneous expression and poor immunogenicity of the viral antigens expressed, namely latent membrane proteins 1 and 2, EBV nuclear antigen 1, and BamHI A rightward reading frame 1 (type-2 [T2] latency). However, EBV lytic cycle proteins are also expressed in certain EBV+ malignancies and, because several EBV lytic cycle proteins are abundantly expressed, have oncogenic activity, and likely contribute to malignancy, we sought and identified viral lytic-cycle transcripts in EBV+ Hodgkin lymphoma biopsies. This provided the rationale for broadening the target antigen-specific repertoire of EBV-specific T cells (EBVSTs) for therapy. We stimulated, peripheral blood mononuclear cells from healthy donors and patients with EBV+ lymphoma with both lytic and latent cycle proteins to produce broad repertoire (BR) EBVSTs. Compared with T2 antigen-specific EBVSTs, BR-EBVSTs more rapidly cleared autologous EBV+ tumors in NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ (NSG) mice and produced higher levels of proinflammatory cytokines that should reactivate the immunosuppressive tumor microenvironment leading to epitope spreading. Our results confirm that lytic cycle antigens are clinically relevant targets for EBV+ lymphoma and underpin the rationale for integrating BR-EBVSTs as a therapeutic approach for relapsed/refractory EBV+ lymphoma (www.clinicaltrials.gov identifiers: #NCT01555892 and #NCT04664179), as well as for other EBV-associated malignancies., (© 2024 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

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., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Invasive species offer insights into rapid adaptation to novel environments. The iconic cane toad ( Rhinella marina ) is an excellent model for studying rapid adaptation during invasion. Previous research using the mitochondrial NADH dehydrogenase 3 ( ND3 ) gene in Hawai'ian and Australian invasive populations found a single haplotype, indicating an extreme genetic bottleneck following introduction. Nuclear genetic diversity also exhibited reductions across the genome in these two populations. Here, we investigated the mitochondrial genomics of cane toads across this invasion trajectory. We created the first reference mitochondrial genome for this species using long-read sequence data. We combined whole-genome resequencing data of 15 toads with published transcriptomic data of 125 individuals to construct nearly complete mitochondrial genomes from the native (French Guiana) and introduced (Hawai'i and Australia) ranges for population genomic analyses. In agreement with previous investigations of these populations, we identified genetic bottlenecks in both Hawai'ian and Australian introduced populations, alongside evidence of population expansion in the invasive ranges. Although mitochondrial genetic diversity in introduced populations was reduced, our results revealed that it had been underestimated: we identified 45 mitochondrial haplotypes in Hawai'ian and Australian samples, none of which were found in the native range. Additionally, we identified two distinct groups of haplotypes from the native range, separated by a minimum of 110 base pairs (0.6%). These findings enhance our understanding of how invasion has shaped the genetic landscape of this species., (© 2024 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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., (© 2024 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

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.

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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.

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

A species' success during the invasion of new areas hinges on an interplay between the demographic processes common to invasions and the specific ecological context of the novel environment. Evolutionary genetic studies of invasive species can investigate how genetic bottlenecks and ecological conditions shape genetic variation in invasions, and our study pairs two invasive populations that are hypothesized to be from the same source population to compare how each population evolved during and after introduction. Invasive European starlings (Sturnus vulgaris) established populations in both Australia and North America in the 19th century. Here, we compare whole-genome sequences among native and independently introduced European starling populations to determine how demographic processes interact with rapid evolution to generate similar genetic patterns in these recent and replicated invasions. Demographic models indicate that both invasive populations experienced genetic bottlenecks as expected based on invasion history, and we find that specific genomic regions have differentiated even on this short evolutionary timescale. Despite genetic bottlenecks, we suggest that genetic drift alone cannot explain differentiation in at least two of these regions. The demographic boom intrinsic to many invasions as well as potential inversions may have led to high population-specific differentiation, although the patterns of genetic variation are also consistent with the hypothesis that this infamous and highly mobile invader adapted to novel selection (e.g., extrinsic factors). We use targeted sampling of replicated invasions to identify and evaluate support for multiple, interacting evolutionary mechanisms that lead to differentiation during the invasion process., (© 2023 John Wiley & Sons Ltd.)

Safely expanding indications for cellular therapies has been challenging given a lack of highly cancer-specific surface markers. Here we explore the hypothesis that tumor cells express cancer-specific surface protein conformations that are invisible to standard target discovery pipelines evaluating gene or protein expression, and these conformations can be identified and immunotherapeutically targeted. We term this strategy integrating cross-linking mass spectrometry with glycoprotein surface capture 'structural surfaceomics'. As a proof of principle, we apply this technology to acute myeloid leukemia (AML), a hematologic malignancy with dismal outcomes and no known optimal immunotherapy target. We identify the activated conformation of integrin β 2 as a structurally defined, widely expressed AML-specific target. We develop and characterize recombinant antibodies to this protein conformation and show that chimeric antigen receptor T cells eliminate AML cells and patient-derived xenografts without notable toxicity toward normal hematopoietic cells. Our findings validate an AML conformation-specific target antigen and demonstrate a tool kit for applying these strategies more broadly., (© 2023. The Author(s).)

The efficacy of therapeutic T-cells is limited by a lack of positive signals and excess inhibitory signaling in tumor microenvironments. We previously showed that a constitutively active IL7 receptor (C7R) enhanced the persistence, expansion, and anti-tumor activity of T-cells expressing chimeric antigen receptors (CARs), and C7R-modified GD2.CAR T-cells are currently undergoing clinical trials. To determine if the C7R could also enhance the activity of T-cells recognizing tumors via their native T-cell receptors (TCRs), we evaluated its effects in Epstein-Barr virus (EBV)-specific T-cells (EBVSTs) that have produced clinical benefits in patients with EBV-associated malignancies. EBVSTs were generated by stimulation of peripheral blood T-cells with overlapping peptide libraries spanning the EBV lymphoma antigens, LMP1, LMP2, and EBNA 1, followed by retroviral vector transduction to express the C7R. The C7R increased STAT5 signaling in EBVSTs and enhanced their expansion over 30 days of culture in the presence or absence of exogenous cytokines. C7R-EBVSTs maintained EBV antigen specificity but were dependent on TCR stimulation for continued expansion. C7R-EBVSTs produced more rapid lymphoma control in a murine xenograft model than unmodified EBVSTs and persisted for longer. The findings have led to a clinical trial, evaluating C7R-EBVSTs for the treatment of refractory or relapsed EBV-positive lymphoma (NCT04664179).

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Stuart, Sherwin, Edwards and Rollins.)

Knowledge of an animal's chronological age is crucial for understanding and predicting population demographics, survival and reproduction, but accurate age determination for many wild animals remains challenging. Previous methods to estimate age require invasive procedures, such as tooth extraction to analyse growth layers, which are difficult to carry out with large, mobile animals such as cetaceans. However, recent advances in epigenetic methods have opened new avenues for precise age determination. These 'epigenetic clocks' present a less invasive alternative and can provide age estimates with unprecedented accuracy. Here, we present a species-specific epigenetic clock based on skin tissue samples for a population of Indo-Pacific bottlenose dolphins ( Tursiops aduncus ) in Shark Bay, Western Australia. We measured methylation levels at 37,492 cytosine-guanine sites (CpG sites) in 165 samples using the mammalian methylation array. Chronological age estimates with an accuracy of ±1 year were available for 68 animals as part of a long-term behavioral study of this population. Using these samples with known age, we built an elastic net model with Leave-One-Out-Cross-Validation, which retained 43 CpG sites, providing an r  = 0.86 and median absolute age error (MAE) = 2.1 years (5% of maximum age). This model was more accurate for our data than the previously published methylation clock based on skin samples of common bottlenose dolphins ( T. truncatus : r  = 0.83, MAE = 2.2) and the multi-species odontocete methylation clock ( r  = 0.68, MAE = 6.8), highlighting that species-specific clocks can have superior performance over those of multi-species assemblages. We further developed an epigenetic sex estimator, predicting sex with 100% accuracy. As age and sex are critical parameters for the study of animal populations, this clock and sex estimator will provide a useful tool for extracting life history information from skin samples rather than long-term observational data for free-ranging Indo-Pacific bottlenose dolphins worldwide., Competing Interests: SH is a founder of the non‐profit Epigenetic Clock Development Foundation, which plans to license several of his patents from his employer UC Regents. The other authors declare no conflicts of interest., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

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., (© 2022 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

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., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Stuart, Sherwin, Cardilini and Rollins.)

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., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

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., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

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.

Cane toads ( Rhinella marina ) are notoriously successful invaders: from 101 individuals brought to Australia in 1935, poisonous toads now cover an area >1.2 million km 2 with adverse effects on native fauna. Despite extensive research on the role of macroparasites in cane toad invasion, viral research is lagging. We compared viral prevalence and diversity between toads in their native range (French Guiana, n =25) and two introduced ranges: Australia ( n =151) and Hawai'i ( n =10) with a metatranscriptomic and metagenomic approach combined with PCR screening. Australian toads almost exclusively harbor one of seven viruses detected globally. Rhimavirus-A ( Picornaviridae ) exhibited low genetic diversity and likely actively infected 9% of sampled Australian toads extending across ~2,000km of Northern Australia and up to the current invasion front. In native range cane toads, we identified multiple phylogenetically distinct viruses ( Iridoviridae , Picornaviridae , Papillomaviridae , and Nackedna-like virus). None of the same viruses was detected in both ranges, suggesting that Australian cane toads have largely escaped the viral infection experienced by their native range counterparts. The novel native range viruses described here are potential biocontrol agents, as Australian toads likely lack prior immunological exposure to these viruses. Overall, our evidence suggests that there may be differences between viruses infecting cane toads in their native vs. introduced ranges, which lays the groundwork for further studies on how these viruses have influenced the toads' invasion history., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Russo, Harding, Yan, Selechnik, Ducatez, DeVore, Zhou, Sarma, Lee, Richardson, Shine, Rollins and White.)

In response to novel environments, invasive populations often evolve rapidly. Standing genetic variation is an important predictor of evolutionary response but epigenetic variation may also play a role. Here, we use an iconic invader, the cane toad ( Rhinella marina ), to investigate how manipulating epigenetic status affects phenotypic traits. We collected wild toads from across Australia, bred them, and experimentally manipulated DNA methylation of the subsequent two generations (G1, G2) through exposure to the DNA methylation inhibitor zebularine and/or conspecific tadpole alarm cues. Direct exposure to alarm cues (an indicator of predation risk) increased the potency of G2 tadpole chemical cues, but this was accompanied by reductions in survival. Exposure to alarm cues during G1 also increased the potency of G2 tadpole cues, indicating intergenerational plasticity in this inducible defence. In addition, the negative effects of alarm cues on tadpole viability (i.e. the costs of producing the inducible defence) were minimized in the second generation. Exposure to zebularine during G1 induced similar intergenerational effects, suggesting a role for alteration in DNA methylation. Accordingly, we identified intergenerational shifts in DNA methylation at some loci in response to alarm cue exposure. Substantial demethylation occurred within the sodium channel epithelial 1 subunit gamma gene ( SCNN1G ) in alarm cue exposed individuals and their offspring. This gene is a key to the regulation of sodium in epithelial cells and may help to maintain the protective epidermal barrier. These data suggest that early life experiences of tadpoles induce intergenerational effects through epigenetic mechanisms, which enhance larval fitness. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

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., (© 2020 John Wiley & Sons Ltd.)

The developmental environment can exert powerful effects on animal phenotype. Recently, epigenetic modifications have emerged as one mechanism that can modulate developmentally plastic responses to environmental variability. For example, the DNA methylation profile at promoters of hormone receptor genes can affect their expression and patterns of hormone release. Across taxonomic groups, epigenetic alterations have been linked to changes in glucocorticoid (GC) physiology. GCs are metabolic hormones that influence growth, development, transitions between life-history stages, and thus fitness. To date, relatively few studies have examined epigenetic effects on phenotypic traits in wild animals, especially in amphibians. Here, we examined the effects of exposure to predation threat (alarm cues) and experimentally manipulated DNA methylation on corticosterone (CORT) levels in tadpoles and metamorphs of the invasive cane toad (Rhinella marina). We included offspring of toads sampled from populations across the species' Australian range. In these animals, exposure to chemical cues from injured conspecifics induces shifts in developmental trajectories, putatively as an adaptive response that lessens vulnerability to predation. We exposed tadpoles to these alarm cues, and measured changes in DNA methylation and CORT levels, both of which are mechanisms that have been implicated in the control of phenotypically plastic responses in tadpoles. To test the idea that DNA methylation drives shifts in GC physiology, we also experimentally manipulated methylation levels with the drug zebularine. We found differentially methylated regions (DMRs) between control tadpoles and their full-siblings exposed to alarm cues, zebularine, or both treatments. However, the effects of these manipulations on methylation patterns were weaker than clutch (e.g., genetic, maternal, etc.) effects. CORT levels were higher in larval cane toads exposed to alarm cues and zebularine. We found little evidence of changes in DNA methylation across the GC receptor gene (NR3C1) promoter region in response to alarm cue or zebularine exposure. In both alarm cue and zebularine-exposed individuals, we found differentially methylated DNA in the suppressor of cytokine signaling 3 gene (SOCS3), which may be involved in predator avoidance behavior. In total, our data reveal that alarm cues have significant impacts on tadpole physiology, but show only weak links between DNA methylation and CORT levels. We also identify genes containing DMRs in tadpoles exposed to alarm cues and zebularine, particularly in range-edge populations, that warrant further investigation., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Integrative and Comparative Biology. All rights reserved. For permissions please email: journals.permissions@oup.com.)

European starlings ( Sturnus vulgaris ) represent one of the most widespread and problematic avian invasive species in the world. Understanding their unique population history and current population dynamics can contribute to conservation efforts and clarify evolutionary processes over short timescales. European starlings were introduced to Central Park, New York in 1890, and from a founding group of about 100 birds, they have expanded across North America with a current population of approximately 200 million. There were also multiple introductions in Australia in the mid-19th century and at least one introduction in South Africa in the late 19th century. Independent introductions on these three continents provide a robust system to investigate invasion genetics. In this study, we compare mitochondrial diversity in European starlings from North America, Australia, and South Africa, and a portion of the native range in the United Kingdom. Of the three invasive ranges, the North American population shows the highest haplotype diversity and evidence of both sudden demographic and spatial expansion. Comparatively, the Australian population shows the lowest haplotype diversity, but also shows evidence for sudden demographic and spatial expansion. South Africa is intermediate to the other invasive populations in genetic diversity but does not show evidence of demographic expansion. In previous studies, population genetic structure was found in Australia, but not in South Africa. Here we find no evidence of population structure in North America. Although all invasive populations share haplotypes with the native range, only one haplotype is shared between invasive populations. This suggests these three invasive populations represent independent subsamples of the native range. The structure of the haplotype network implies that the native-range sampling does not comprehensively characterize the genetic diversity there. This study represents the most geographically widespread analysis of European starling population genetics to date., Competing Interests: None declared., (© 2020 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Noninvasive sampling methods for studying intestinal microbiomes are widely applied in studies of endangered species and in those conducting temporal monitoring during manipulative experiments. Although existing studies show that noninvasive sampling methods among different taxa vary in their accuracy, no studies have yet been published comparing nonlethal sampling methods in adult amphibians. In this study, we compare microbiomes from two noninvasive sample types (faeces and cloacal swabs) to that of the large intestine in adult cane toads, Rhinella marina. We use 16S rRNA gene sequencing to investigate how microbial communities change along the digestive tract and which nonlethal sampling method better represents large intestinal microbiota. We found that cane toads' intestinal microbiota was dominated by Bacteroidetes, Proteobacteria and Firmicutes and, interestingly, we also saw a high proportion of Fusobacteria, which has previously been associated with marine species and changes in frog immunity. The large and small intestine of cane toads had a similar microbial composition, but the large intestine showed higher diversity. Our results indicate that cloacal swabs were more similar to large intestine samples than were faecal samples, and small intestine samples were significantly different from both nonlethal sample types. Our study provides valuable information for future investigations of the cane toad gut microbiome and validates the use of cloacal swabs as a nonlethal method to study changes in the large intestine microbiome. These data provide insights for future studies requiring nonlethal sampling of amphibian gut microbiota., (© 2020 John Wiley & Sons Ltd.)

Purpose: Current protocols for CD19 chimeric antigen receptor-expressing T cells (CD19.CAR-T cells) require recipients to tolerate preinfusion cytoreductive chemotherapy, and the presence of sufficient target antigen on normal or malignant B cells., Patients and Methods: We investigated whether additional stimulation of CD19.CAR-T cells through their native receptors can substitute for cytoreductive chemotherapy, inducing expansion and functional persistence of CD19.CAR-T even in patients in remission of B-cell acute lymphocytic leukemia. We infused a low dose of CD19.CAR-modified virus-specific T cells (CD19.CAR-VST) without prior cytoreductive chemotherapy into 8 patients after allogeneic stem cell transplant., Results: Absent virus reactivation, we saw no CD19.CAR-VST expansion. In contrast, in patients with viral reactivation, up to 30,000-fold expansion of CD19.CAR-VSTs was observed, with depletion of CD19 + B cells. Five patients remain in remission at 42-60+ months., Conclusions: Dual T-cell receptor and CAR stimulation can thus potentiate effector cell expansion and CAR-target cell killing, even when infusing low numbers of effector cells without cytoreduction., (©2019 American Association for Cancer Research.)

Invasive species often evolve rapidly following introduction despite genetic bottlenecks that may result from small numbers of founders; however, some invasions may not fit this "genetic paradox". The invasive cane toad ( Rhinella marina ) displays high phenotypic variation across its introduced Australian range. Here, we used three genome-wide datasets to characterize their population structure and genetic diversity. We found that toads form three genetic clusters: 1) native range toads, 2) toads from the source population in Hawaii and long-established areas near introduction sites in Australia, and 3) toads from more recently established northern Australian sites. Although we find an overall reduction in genetic diversity following introduction, we do not see this reduction in loci putatively under selection, suggesting that genetic diversity may have been maintained at ecologically relevant traits, or that mutation rates were high enough to maintain adaptive potential. Nonetheless, toads encounter novel environmental challenges in Australia, and the transition between genetic clusters occurs at a point along the invasion transect where temperature rises and rainfall decreases. We identify environmentally associated loci known to be involved in resistance to heat and dehydration. This study highlights that natural selection occurs rapidly and plays a vital role in shaping the structure of invasive populations., (Copyright © 2019 Selechnik, Richardson, Shine, DeVore, Ducatez and Rollins.)

Host-pathogen associations change rapidly during a biological invasion and are predicted to impose strong selection on immune function. It has been proposed that the invader may experience an abrupt reduction in pathogen-mediated selection ("enemy release"), thereby favoring decreased investment into "costly" immune responses. Across plants and animals, there is mixed support for this prediction. Pathogens are not the only form of selection imposed on invaders; differences in abiotic environmental conditions between native and introduced ranges are also expected to drive rapid evolution. Here, we use RNA-Seq to assess the expression patterns of immune and environmentally associated genes in the cane toad ( Rhinella marina ) across its invasive Australian range. Transcripts encoding mediators of costly immune responses (inflammation, cytotoxicity) showed a curvilinear relationship with invasion history, with highest expression in toads from oldest and newest colonized areas. This pattern is surprising given theoretical expectations of density dynamics in invasive species and may be because density influences both intraspecific competition and parasite transmission, generating conflicting effects on the strength of immune responses. Alternatively, this expression pattern may be the result of other evolutionary forces, such as spatial sorting and genetic drift, working simultaneously with natural selection. Our findings do not support predictions about immune function based on the enemy release hypothesis and suggest instead that the effects of enemy release are difficult to isolate in wild populations, especially in the absence of information regarding parasite and pathogen infection. Additionally, expression patterns of genes underlying putatively environmentally associated traits are consistent with previous genetic studies, providing further support that Australian cane toads have adapted to novel abiotic challenges., Competing Interests: None declared.

Individual hosts differ extensively in their competence for parasites, but traditional research has discounted this variation, partly because modeling such heterogeneity is difficult. This discounting has diminished as tools have improved and recognition has grown that some hosts, the extremely competent, can have exceptional impacts on disease dynamics. Most prominent among these hosts are the superspreaders, but other forms of extreme competence (EC) exist and others await discovery; each with potentially strong but distinct implications for disease emergence and spread. Here, we propose a framework for the study and discovery of EC, suitable for different host-parasite systems, which we hope enhances our understanding of how parasites circulate and evolve in host communities., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

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.

Cane toads are a notorious invasive species, inhabiting over 1.2 million km 2 of Australia and threatening native biodiversity. The release of pathogenic cane toad viruses is one possible biocontrol strategy yet is currently hindered by the poorly described cane toad virome. Metatranscriptomic analysis of 16 cane toad livers revealed the presence of a novel and full-length picornavirus, Rhimavirus A (RhiV-A), a member of a reptile- and amphibian-specific cluster of the Picornaviridae basal to the Kobuvirus -like group. In the combined liver transcriptome, we also identified a complete genome sequence of a distinct epsilonretrovirus, Rhinella marina endogenous retrovirus (RMERV). The recently sequenced cane toad genome contains 8 complete RMERV proviruses as well as 21 additional truncated insertions. The oldest full-length RMERV provirus was estimated to have inserted 1.9 million years ago (MYA). To screen for these viral sequences in additional toads, we analyzed publicly available transcriptomes from six diverse Australian locations. RhiV-A transcripts were identified in toads sampled from three locations across 1,000 km of Australia, stretching to the current Western Australia (WA) invasion front, while RMERV transcripts were observed at all six sites. Finally, we scanned the cane toad genome for nonretroviral endogenous viral elements, finding three sequences related to small DNA viruses in the family Circoviridae This shows ancestral circoviral infection with subsequent genomic integration. The identification of these current and past viral infections enriches our knowledge of the cane toad virome, an understanding of which will facilitate future work on infection and disease in this important invasive species. IMPORTANCE Cane toads are poisonous amphibians that were introduced to Australia in 1935 for insect control. Since then, their population has increased dramatically, and they now threaten many native Australian species. One potential method to control the population is to release a cane toad virus with high mortality rates, yet few cane toad viruses have been characterized. This study samples cane toads from different Australian locations and uses an RNA sequencing and computational approach to find new viruses. We report novel complete picornavirus and retrovirus sequences that were genetically similar to viruses infecting frogs, reptiles, and fish. Using data generated in other studies, we show that these viral sequences are present in cane toads from distinct Australian locations. Three sequences related to circoviruses were also found in the toad genome. The identification of new viral sequences will aid future studies that investigate their prevalence and potential as agents for biocontrol., (Copyright © 2018 American Society for Microbiology.)

Empirical genetic datasets used for estimating contemporary dispersal in wild populations and to correctly identify dispersers are rarely tested to determine if they are capable of providing accurate results. Here we test whether a genetic dataset provides sufficient information to accurately identify first-generation dispersers. Using microsatellite data from three wild populations of common starlings ( Sturnus vulgaris ), we artificially simulated dispersal of a subset of individuals; we term this 'Simulated Disperser Analysis'. We then ran analyses for diminishing numbers of loci, to assess at which point simulated dispersers could no longer be correctly identified. Not surprisingly, the correct identification of dispersers varied significantly depending on the individual chosen to 'disperse', the number of loci used, whether loci had high or low Polymorphic Information Content and the location to which the dispersers were moved. A review of the literature revealed that studies that have implemented first-generation migrant detection to date have used on average 10 microsatellite loci. Our results suggest at least 27 loci are required to accurately identify dispersers in the study system evaluated here. We suggest that future studies use the approach we describe to determine the appropriate number of markers needed to accurately identify dispersers in their study system; the unique nature of natural systems means that the number of markers required for each study system will vary. Future studies can use Simulated Disperser Analysis on pilot data to test marker panels for robustness to contemporary dispersal identification, providing a powerful tool in the efficient and accurate design of studies using genetic data to estimate dispersal., Competing Interests: Lee A. Rollins is an Academic Editor for PeerJ.

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., (© The Authors 2017. Published by Oxford University Press.)

Genetic diversity at community, population and individual levels is thought to influence the spread of infectious disease. At the individual level, inbreeding and heterozygosity are associated with increased risk of infection and disease severity. Host genotype rarity may also reduce infection risk if pathogens are co-adapted to common or local hosts, but to date, no studies have investigated the relative importance of genotype rarity and heterozygosity for infection in a wild, sexually reproducing vertebrate. With beak and feather disease virus (BFDV) infection in a wild parrot (Platycercus elegans), we show that both heterozygosity and genotype rarity of individual hosts predicted infection, but in contrasting ways. Heterozygosity was negatively associated with probability of infection, but not with infection load. In contrast, increased host genotype rarity was associated with lower viral load in infected individuals, but did not predict infection probability. These effects were largely consistent across subspecies, but were not evident at the population level. Subspecies and age were also strongly associated with infection. Our study provides novel insights into infection dynamics by quantifying rarity and diversity simultaneously. We elucidate roles that host genetic diversity can play in infection dynamics, with implications for understanding population divergence, intraspecific diversity and conservation.

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., Competing Interests: Lee A. Rollins is an Academic Editor for PeerJ.

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., Competing Interests: Competing Interests: The authors have declared that no competing interest exists.

In socially monogamous species, individuals can use extra-pair paternity and offspring sex allocation as adaptive strategies to ameliorate costs of genetic incompatibility with their partner. Previous studies on domesticated Gouldian finches (Erythrura gouldiae) demonstrated a genetic incompatibility between head colour morphs, the effects of which are more severe in female offspring. Domesticated females use differential sex allocation, and extra-pair paternity with males of compatible head colour, to reduce fitness costs associated with incompatibility in mixed-morph pairings. However, laboratory studies are an oversimplification of the complex ecological factors experienced in the wild and may only reflect the biology of a domesticated species. This study aimed to examine the patterns of parentage and sex ratio bias with respect to colour pairing combinations in a wild population of the Gouldian finch. We utilized a novel PCR assay that allowed us to genotype the morph of offspring before the morph phenotype develops and to explore bias in morph paternity and selection at the nest. Contrary to previous findings in the laboratory, we found no effect of pairing combinations on patterns of extra-pair paternity, offspring sex ratio or selection on morphs in nestlings. In the wild, the effect of morph incompatibility is likely much smaller, or absent, than was observed in the domesticated birds. Furthermore, the previously studied domesticated population is genetically differentiated from the wild population, consistent with the effects of domestication. It is possible that the domestication process fostered the emergence (or enhancement) of incompatibility between colour morphs previously demonstrated in the laboratory., (© 2016 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2016 European Society For Evolutionary Biology.)