Concurrent invasions of European starlings in Australia and North America reveal population-specific differentiation in shared genomic regions.

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.

Viral Discovery in the Invasive Australian Cane Toad (Rhinella marina) Using Metatranscriptomic and Genomic Approaches.

Cane toads are a notorious invasive species, inhabiting over 1.2 million km2 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.

Whole-mitogenome analysis unveils previously undescribed genetic diversity in cane toads across their invasion trajectory.

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.

An epigenetic DNA methylation clock for age estimates in Indo-Pacific bottlenose dolphins (Tursiops aduncus).

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.

Building genetic networks using relatedness information: a novel approach for the estimation of dispersal and characterization of group structure in social animals.

Natal dispersal is an important life history trait driving variation in individual fitness, and therefore, a proper understanding of the factors underlying dispersal behaviour is critical to many fields including population dynamics, behavioural ecology and conservation biology. However, individual dispersal patterns remain difficult to quantify despite many years of research using direct and indirect methods. Here, we quantify dispersal in a single intensively studied population of the cooperatively breeding chestnut-crowned babbler (Pomatostomus ruficeps) using genetic networks created from the combination of pairwise relatedness data and social networking methods and compare this to dispersal estimates from re-sighting data. This novel approach not only identifies movements between social groups within our study sites but also provides an estimation of immigration rates of individuals originating outside the study site. Both genetic and re-sighting data indicated that dispersal was strongly female biased, but the magnitude of dispersal estimates was much greater using genetic data. This suggests that many previous studies relying on mark-recapture data may have significantly underestimated dispersal. An analysis of spatial genetic structure within the sampled population also supports the idea that females are more dispersive, with females having no structure beyond the bounds of their own social group, while male genetic structure expands for 750 m from their social group. Although the genetic network approach we have used is an excellent tool for visualizing the social and genetic microstructure of social animals and identifying dispersers, our results also indicate the importance of applying them in parallel with behavioural and life history data.

Mitochondrial DNA offers unique insights into invasion history of the common starling.

Mitochondrial DNA (mtDNA) can be a powerful genetic marker for tracing origins and history of invasive populations. Here, we use mtDNA to address questions relevant to the understanding of invasion pathways of common starlings (Sturnus vulgaris) into Western Australia (WA) and discuss the utility of this marker to provide information useful to invasive species management. Mitochondrial sequence data indicate two geographically restricted genetic groups within Australia. Evidence of dispersal from genetically distinct sources outside the sampled range of starlings in Australia suggests increased vigilance by management agencies may be required to prevent further incursions from widely separated localities. Overall, genetic diversity in Australia was lower than in samples from the native range. Within Australia, genetic diversity was lowest in the most recently colonized area in the west, indicating that demographic bottlenecks have occurred in this area. Evidence of restricted dispersal between localities on the edge of the range expansion (ERE) in WA and other Australian sampling localities suggests that localized control within the ERE may be effective in preventing further range expansion. Signatures of spatial and demographic expansion are present in mismatch analyses from sampling localities located at the ERE, but neutrality indices did not support this finding, suggesting that the former may be more sensitive to recent expansion. Additionally, mismatch analyses support the presence of admixture, which is likely to have occurred pre-introduction. We compare our findings with those from a microsatellite study of the same samples and discuss how the mtDNA analyses used here offer valuable and unique insights into the invasion history of introduced species.

Contrasting mitochondrial diversity of European starlings (Sturnus vulgaris) across three invasive continental distributions.

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.

A comparison of nonlethal sampling methods for amphibian gut microbiome analyses.

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.

Invasive species can't cover their tracks: using microsatellites to assist management of starling (Sturnus vulgaris) populations in Western Australia.

Invasive species are known to cause environmental and economic damage, requiring management by control agencies worldwide. These species often become well established in new environments long before their detection, resulting in a lack of knowledge regarding their history and dynamics. When new invasions are discovered, information regarding the source and pathway of the invasion, and the degree of connectivity with other populations can greatly benefit management strategies. Here we use invasive common starling (Sturnus vulgaris) populations from Australia to demonstrate that genetic techniques can provide this information to aid management, even when applied to highly vagile species over continental scales. Analysis of data from 11 microsatellites in 662 individuals sampled at 17 localities across their introduced range in Australia revealed four populations. One population consisted of all sampling sites from the expansion front in Western Australia, where control efforts are focused. Despite evidence of genetic exchange over both contemporary and historical timescales, gene flow is low between this population and all three more easterly populations. This suggests that localized control of starlings on the expansion front may be an achievable goal and the long-standing practice of targeting select proximal eastern source populations may be ineffective on its own. However, even with low levels of gene flow, successful control of starlings on the expansion front will require vigilance, and genetic monitoring of this population can provide essential information to managers. The techniques used here are broadly applicable to invasive populations worldwide.

Immune and environment-driven gene expression during invasion: An eco-immunological application of RNA-Seq.

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.

Host heterozygosity and genotype rarity affect viral dynamics in an avian subspecies complex.

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.

Expected Shannon Entropy and Shannon Differentiation between Subpopulations for Neutral Genes under the Finite Island Model.

Shannon entropy H and related measures are increasingly used in molecular ecology and population genetics because (1) unlike measures based on heterozygosity or allele number, these measures weigh alleles in proportion to their population fraction, thus capturing a previously-ignored aspect of allele frequency distributions that may be important in many applications; (2) these measures connect directly to the rich predictive mathematics of information theory; (3) Shannon entropy is completely additive and has an explicitly hierarchical nature; and (4) Shannon entropy-based differentiation measures obey strong monotonicity properties that heterozygosity-based measures lack. We derive simple new expressions for the expected values of the Shannon entropy of the equilibrium allele distribution at a neutral locus in a single isolated population under two models of mutation: the infinite allele model and the stepwise mutation model. Surprisingly, this complex stochastic system for each model has an entropy expressable as a simple combination of well-known mathematical functions. Moreover, entropy- and heterozygosity-based measures for each model are linked by simple relationships that are shown by simulations to be approximately valid even far from equilibrium. We also identify a bridge between the two models of mutation. We apply our approach to subdivided populations which follow the finite island model, obtaining the Shannon entropy of the equilibrium allele distributions of the subpopulations and of the total population. We also derive the expected mutual information and normalized mutual information ("Shannon differentiation") between subpopulations at equilibrium, and identify the model parameters that determine them. We apply our measures to data from the common starling (Sturnus vulgaris) in Australia. Our measures provide a test for neutrality that is robust to violations of equilibrium assumptions, as verified on real world data from starlings.

Personality in captivity: more exploratory males reproduce better in an aviary population.

The existence of animal personality is well-established across a wide range of species, with the majority of evidence for this being obtained from individuals held in captivity. However, there has been little work assessing the influence of commonly-measured personality traits on fitness, which is pertinent when the genetic basis of personality is considered. We measured whether the reproductive behaviour and success of zebra finches in a captive mixed-sex aviary environment was influenced by an aspect of their personality, their exploratory behaviour in a single-sex social aviary. We found that more exploratory males made a greater number of breeding attempts and raised more nestlings than less exploratory males. These results were not confounded by extra-pair paternity, which was not related to personality, or by the individuals that did not initiate any reproductive attempts at all. Our work provides evidence that attributes of personality may influence the degree to which individuals cope with, and thrive in a captive environment and this should be accounted for in both experimental design and the interpretation of results. Furthermore, this suggests that there may be selection on these traits as part of the domestication process.

Genetic monogamy despite variable ecological conditions and social environment in the cooperatively breeding apostlebird.

Mating strategies may be context-dependent and may vary across ecological and social contexts, demonstrating the role of these factors in driving the variation in genetic polyandry within and among species. Here, we took a longitudinal approach across 5 years (2006-2010), to study the apostlebird (Struthidea cinerea), an Australian cooperatively breeding bird, whose reproduction is affected by ecological "boom and bust" cycles. Climatic variation drives variation in the social (i.e., group sizes, proportion of males and females) and ecological (i.e., plant and insect abundance) context in which mating occurs. By quantifying variation in both social and ecological factors and characterizing the genetic mating system across multiple years using a molecular parentage analysis, we found that the genetic mating strategy did not vary among years despite significant variation in rainfall, driving primary production, and insect abundance, and corresponding variation in social parameters such as breeding group size. Group sizes in 2010, an ecologically good year, were significantly smaller (mean = 5.8 ± 0.9, n = 16) than in the drought affected years, between 2006 and 2008, (mean = 9.1 ± 0.5, n = 63). Overall, apostlebirds were consistently monogamous with few cases of multiple maternity or paternity (8 of 78 nests) across all years.

The role of the Ord Arid Intrusion in the historical and contemporary genetic division of long-tailed finch subspecies in northern Australia.

The effect of separation by biogeographic features followed by secondary contact can blur taxonomic boundaries and produce complex genetic signatures. We analyzed population structure and gene flow across the range of the long-tailed finch (Poephila acuticauda) in northern Australia (1) to test the hypothesis that Ord Arid Intrusion acted as the causative barrier that led to divergence of P. acuticauda subspecies, (2) to determine whether genetic data support the presence of a gradual cline across the range or a sudden shift, both of which have been suggested based on morphological data, and (3) to estimate levels of contemporary gene flow within this species complex. We collected samples from 302 individuals from 10 localities. Analyses of 12 microsatellite loci and sequence data from 333 base pairs of the mitochondrial control region were used to estimate population structure and gene flow, using analysis of molecular variance (AMOVA), haplotype network analysis, frequency statistics, and clustering methods. Mitochondrial sequence data indicated the presence of three genetic groups (regions) across the range of P. acuticauda. Genetic diversity was highest in the east and lowest in the west. The Ord Arid Intrusion appears to have functioned as a biogeographic barrier in the past, according to mtDNA evidence presented here and evidence from previous studies. The absence of isolation by distance between adjacent regions and the lack of population genetic structure of mtDNA within regions indicates that genetic changes across the range of P. acuticauda subspecies are characterized by discrete breaks between regions. While microsatellite data indicate a complete absence of genetic structure across this species' range, it appears unlikely that this results from high levels of gene flow. Mitochondrial data do not support the presence of contemporary gene flow across the range of this species.

Sex-dependent selection differentially shapes genetic variation on and off the guppy Y chromosome.

Because selection is often sex-dependent, alleles can have positive effects on fitness in one sex and negative effects in the other, resulting in intralocus sexual conflict. Evolutionary theory predicts that intralocus sexual conflict can drive the evolution of sex limitation, sex-linkage, and sex chromosome differentiation. However, evidence that sex-dependent selection results in sex-linkage is limited. Here, we formally partition the contribution of Y-linked and non-Y-linked quantitative genetic variation in coloration, tail, and body size of male guppies (Poecilia reticulata)-traits previously implicated as sexually antagonistic. We show that these traits are strongly genetically correlated, both on and off the Y chromosome, but that these correlations differ in sign and magnitude between both parts of the genome. As predicted, variation in attractiveness was found to be associated with the Y-linked, rather than with the non-Y-linked component of genetic variation in male ornamentation. These findings show how the evolution of Y-linkage may be able to resolve sexual conflict. More generally, they provide unique insight into how sex-specific selection has the potential to differentially shape the genetic architecture of fitness traits across different parts of the genome.