Discordance between nuclear and mitochondrial genomes in sexual and asexual lineages of the freshwater snail Potamopyrgus antipodarum

Genome Evolution and Introgression in the New Zealand mud Snails Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa

We have sequenced, assembled, and analyzed the nuclear and mitochondrial genomes and transcriptomes of Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa, two prosobranch snail species native to New Zealand that together span the continuum from estuary to freshwater. These two species are the closest known relatives of the freshwater species Potamopyrgus antipodarum-a model for studying the evolution of sex, host-parasite coevolution, and biological invasiveness-and thus provide key evolutionary context for understanding its unusual biology. The P. estuarinus and P. kaitunuparaoa genomes are very similar in size and overall gene content. Comparative analyses of genome content indicate that these two species harbor a near-identical set of genes involved in meiosis and sperm functions, including seven genes with meiosis-specific functions. These results are consistent with obligate sexual reproduction in these two species and provide a framework for future analyses of P. antipodarum-a species comprising both obligately sexual and obligately asexual lineages, each separately derived from a sexual ancestor. Genome-wide multigene phylogenetic analyses indicate that P. kaitunuparaoa is likely the closest relative to P. antipodarum. We nevertheless show that there has been considerable introgression between P. estuarinus and P. kaitunuparaoa. That introgression does not extend to the mitochondrial genome, which appears to serve as a barrier to hybridization between P. estuarinus and P. kaitunuparaoa. Nuclear-encoded genes whose products function in joint mitochondrial-nuclear enzyme complexes exhibit similar patterns of nonintrogression, indicating that incompatibilities between the mitochondrial and the nuclear genome may have prevented more extensive gene flow between these two species.

Invasive freshwater snails are less sensitive to population density than native conspecifics

Understanding how and why some species or lineages become invasive is critically important for effectively predicting and mitigating biological invasions. Here, we address an important unanswered question in invasion biology: do key life-history traits of invasive versus native lineages within a species differ in response to key environmental stressors? We focus on the environmental factor of population density, which is a fundamental characteristic of all populations, and investigate how changes in density affect native versus invasive Potamopyrgus antipodarum (New Zealand mudsnail). P. antipodarum has invaded 39 countries and detrimentally affects invaded environments. Previous studies of native and invasive populations and from laboratory experiments have demonstrated that growth and reproduction of P. antipodarum is sensitive to population density, though whether and how this sensitivity varies across native versus invasive lineages remains uncharacterized. We quantified individual growth rate and reproduction in P. antipodarum from multiple distinct native and invasive lineages across three different population density treatments. The growth of native but not invasive lineages decreased as density increased. There was no differential effect of density treatment on embryo production of invasive versus native snails, but a significantly higher proportion of snails were reproductive in high density compared to intermediate density for invasive lineages. In native lineages, there were no significant differences in the relative frequency of reproductive snails across density treatments. While the extent to which these results from our laboratory study can be extrapolated to the more complex natural world remain unclear, our findings are consistent with a scenario where differential sensitivity to population density could help explain why some lineages become successful invaders. Our findings also align with previous studies that show that invasive P. antipodarum lineages exhibit a relatively wide range of tolerance to environmental stressors.

Single-molecule Sequencing of an Animal Mitochondrial Genome Reveals Chloroplast-like Architecture and Repeat-mediated Recombination

Recent advances in long-read sequencing technology have allowed for single-molecule sequencing of entire mitochondrial genomes, opening the door for direct investigation of the mitochondrial genome architecture and recombination. We used PacBio sequencing to reassemble mitochondrial genomes from two species of New Zealand freshwater snails, Potamopyrgus antipodarum and Potamopyrgus estuarinus. These assemblies revealed a ∼1.7 kb structure within the mitochondrial genomes of both species that was previously undetected by an assembly of short reads and likely corresponding to a large noncoding region commonly present in the mitochondrial genomes. The overall architecture of these Potamopyrgus mitochondrial genomes is reminiscent of the chloroplast genomes of land plants, harboring a large single-copy (LSC) region and a small single-copy (SSC) region separated by a pair of inverted repeats (IRa and IRb). Individual sequencing reads that spanned across the Potamopyrgus IRa-SSC-IRb structure revealed the occurrence of a "flip-flop" recombination. We also detected evidence for two distinct IR haplotypes and recombination between them in wild-caught P. estuarinus, as well as extensive intermolecular recombination between single-nucleotide polymorphisms in the LSC region. The chloroplast-like architecture and repeat-mediated mitochondrial recombination we describe here raise fundamental questions regarding the origins and commonness of inverted repeats in cytoplasmic genomes and their role in mitochondrial genome evolution.

Evidence for stronger discrimination between conspecific and heterospecific mating partners in sexual vs. asexual female freshwater snails

Once-useful traits that no longer contribute to fitness tend to decay over time. Here, we address whether the expression of mating-related traits that increase the fitness of sexually reproducing individuals but are likely less useful or even costly to asexual counterparts seems to exhibit decay in the latter. Potamopyrgus antipodarum is a New Zealand freshwater snail characterized by repeated transitions from sexual to asexual reproduction. The frequent coexistence of sexual and asexual lineages makes P. antipodarum an excellent model for the study of mating-related trait loss. Under the presumption (inherent in the Biological Species Concept) that failure to discriminate between conspecific and heterospecific mating partners represents a poor mate choice, we used a mating choice assay including sexual and asexual P. antipodarum females and conspecific (presumed better choice) vs. heterospecific (presumed worse choice) males to evaluate the loss of behavioral traits related to sexual reproduction. We found that sexual females engaged in mating behaviors with conspecific mating partners more frequently and for a greater duration than with heterospecific mating partners. By contrast, asexual females mated at similar frequency and duration as sexual females, but did not mate more often or for longer duration with conspecific vs. heterospecific males. While further confirmation will require inclusion of a more diverse array of sexual and asexual lineages, these results are consistent with a scenario where selection acting to maintain effective mate discrimination in asexual P. antipodarum is weak or ineffective relative to sexual females and, thus, where asexual reproduction is associated with the evolutionary decay of mating-related traits in this system.

Stay in shape: Assessing the adaptive potential of shell morphology and its sensitivity to temperature in the invasive New Zealand mud snail Potamopyrgus antipodarum through phenotypic plasticity and natural selection in Europe

Climate change may force organisms to adapt genetically or plastically to new environmental conditions. Invasive species show remarkable potential for rapid adaptation. The ovoviviparous New Zealand mud snail (NZMS), Potamopyrgus antipodarum, has successfully established across Europe with two clonally reproducing mitochondrial lineages since its arrival in the first half of the 19th century. Its remarkable variation in shell morphology was shown to be fitness relevant. We investigated the effects of temperature on shell morphology across 11 populations from Germany and the Iberian Peninsula in a common garden across three temperatures. We analyzed size and shape using geometric morphometrics. For both, we compared reaction norms and estimated heritabilities. For size, the interaction of temperature and haplotype explained about 50% of the total variance. We also observed more genotype by environment interactions indicating a higher degree of population differentiation than in shape. Across the three temperatures, size followed the expectations of the temperature-size rule, with individuals growing larger in cold environments. Changes in shape may have compensated for changes in size affecting space for brooding embryos. Heritability estimates were relatively high. As indicated by the very low coefficients of variation for clonal repeatability (CV A ), they can probably not be compared in absolute terms. However, they showed some sensitivity to temperature, in haplotype t more so than in z, which was only found in Portugal. The low CV A values indicate that genetic variation among European populations is still restricted with a low potential to react to selection. A considerable fraction of the genetic variation was due to differences between the clonal lineages. The NZMS has apparently not been long enough in Europe to accumulate significant genetic variation relevant for morphological adaptation. As temperature is obviously not the sole factor influencing shell morphology, their interaction will probably not be a factor limiting population persistence under a warming climate in Europe.

Life-history trait variation in native versus invasive asexual New Zealand mud snails

Potamopyrgus antipodarum is a New Zealand freshwater snail that is invasive worldwide. While native P. antipodarum populations are characterized by frequent coexistence between obligately sexual and obligately asexual individuals, only the asexual snails are known to invade other ecosystems. Despite low genetic diversity and the absence of sex, invasive asexual P. antipodarum are highly successful. Here, we quantified variation in three key life-history traits across invasive P. antipodarum lineages and compared this variation to already documented variation in these same traits in asexual native lineages to provide a deeper understanding of why some lineages become invasive. In particular, we evaluated if invasive lineages of P. antipodarum could be successful because they represent life-history variation from native ancestors that could facilitate invasion. We found that invasive snails displayed a non-representative sample of native diversity, with invasive snails growing more slowly and maturing more rapidly than their native counterparts. These results are consistent with expectations of a scenario where invasive lineages represent a subset of native variation that is beneficial in the setting of invasion. Together, these results help illuminate the mechanisms driving the worldwide expansion of invasive populations of these snails.

Patterns of gene expression in ovaries of sexual vs. asexual lineages of a freshwater snail

Why sexual reproduction is so common when asexual reproduction should be much more efficient and less costly remains an open question in evolutionary biology. Comparisons between otherwise similar sexual and asexual taxa allow us to characterize the genetic architecture underlying asexuality, which can, in turn, illuminate how this reproductive mode transition occurred and the mechanisms by which it is maintained or disrupted. Here, we used transcriptome sequencing to compare patterns of ovarian gene expression between actively reproducing obligately sexual and obligately asexual females from multiple lineages of Potamopyrgus antipodarum, a freshwater New Zealand snail characterized by frequent separate transitions to asexuality and coexistence of otherwise similar sexual and asexual lineages. We also used these sequence data to evaluate whether population history accounts for variation in patterns of gene expression. We found that source population was a major source of gene expression variation, and likely more influential than reproductive mode. This outcome for these common garden-raised snails is strikingly similar to earlier results from field-collected snails. While we did not identify a likely set of candidate genes from expression profiles that could plausibly explain how transitions to asexuality occurred, we identified around 1,000 genes with evidence of differential expression between sexual and asexual reproductive modes, and 21 genes that appear to exhibit consistent expression differences between sexuals and asexuals across genetic backgrounds. This second smaller set of genes provides a good starting point for further exploration regarding a potential role in the transition to asexual reproduction. These results mark the first effort to characterize the causes of asexuality in P. antipodarum, demonstrate the apparently high heritability of gene expression patterns in this species, and hint that for P. antipodarum, transitions to asexuality might not necessarily be strongly associated with broad changes in gene expression.

Phylogeography and Cryptic Species Structure of a Locally Adapted Parasite in New Zealand

Phylogeographic patterns of many taxa on New Zealand's South Island are characterised by disjunct distributions that have been attributed to Pleistocene climatic cycles and the formation of the Southern Alps. Pleistocene glaciation has also been implicated in shaping the contemporary genetic differentiation between populations of the aquatic snail Potamopyrgus antipodarum. We investigated whether similar phylogeographic patterns exist for the snail's locally adapted trematode parasite, Atriophallophorus winterbourni. We found evidence for a barrier to gene-flow in sympatry between cryptic, but ecologically divergent species. When focusing on the most common of these species, disjunct geographic distributions are found for mitochondrial lineages that diverged during the Pleistocene. The boundary between these distributions is found in the central part of the South Island and reinforced by low cross-alpine migration. Further support for a vicariant origin of the phylogeographic pattern was found when assessing nuclear multilocus SNP data. Nuclear and mitochondrial population differentiation was concordant in pattern, except for populations in a potential secondary contact zone. Additionally, we found larger than expected differentiation between nuclear- and mitochondrial-based empirical Bayes F_ST estimates (global F_ST : 0.02 versus 0.39 for nuclear and mitochondrial data, respectively). Population subdivision is theoretically expected to be stronger for mitochondrial genomes due to a smaller effective population size, but the strong difference here, together with mito-nuclear discordance in a putative contact zone, is potentially indicative of divergent gene flow of nuclear and mitochondrial genomes.

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits-mode of genomic inheritance-influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S-5.8S-28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA-histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction-with potential functional consequences related to the loss of sexual reproduction.

Phenotypic Variation in Mitochondria-Related Performance Traits Across New Zealand Snail Populations

Mitochondrial function is critical for energy homeostasis and should shape how genetic variation in metabolism is transmitted through levels of biological organization to generate stability in organismal performance. Mitochondrial function is encoded by genes in two distinct and separately inherited genomes-the mitochondrial genome and the nuclear genome-and selection is expected to maintain functional mito-nuclear interactions. The documented high levels of polymorphism in genes involved in these mito-nuclear interactions and wide variation for mitochondrial function demands an explanation for how and why variability in such a fundamental trait is maintained. Potamopyrgus antipodarum is a New Zealand freshwater snail with coexisting sexual and asexual individuals and, accordingly, contrasting systems of separate vs. co-inheritance of nuclear and mitochondrial genomes. As such, this snail provides a powerful means to dissect the evolutionary and functional consequences of mito-nuclear variation. The lakes inhabited by P. antipodarum span wide environmental gradients, with substantial across-lake genetic structure and mito-nuclear discordance. This situation allows us to use comparisons across reproductive modes and lakes to partition variation in cellular respiration across genetic and environmental axes. Here, we integrated cellular, physiological, and behavioral approaches to quantify variation in mitochondrial function across a diverse set of wild P. antipodarum lineages. We found extensive across-lake variation in organismal oxygen consumption and behavioral response to heat stress and differences across sexes in mitochondrial membrane potential but few global effects of reproductive mode. Taken together, our data set the stage for applying this important model system for sexual reproduction and polyploidy to dissecting the complex relationships between mito-nuclear variation, performance, plasticity, and fitness in natural populations.

Mobilizing molluscan models and genomes in biology

Molluscs are among the most ancient, diverse, and important of all animal taxa. Even so, no individual mollusc species has emerged as a broadly applied model system in biology. We here make the case that both perceptual and methodological barriers have played a role in the relative neglect of molluscs as research organisms. We then summarize the current application and potential of molluscs and their genomes to address important questions in animal biology, and the state of the field when it comes to the availability of resources such as genome assemblies, cell lines, and other key elements necessary to mobilising the development of molluscan model systems. We conclude by contending that a cohesive research community that works together to elevate multiple molluscan systems to 'model' status will create new opportunities in addressing basic and applied biological problems, including general features of animal evolution. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Gene Duplication and Gain in the Trematode Atriophallophorus winterbourni Contributes to Adaptation to Parasitism

Gene duplications and novel genes have been shown to play a major role in helminth adaptation to a parasitic lifestyle because they provide the novelty necessary for adaptation to a changing environment, such as living in multiple hosts. Here we present the de novo sequenced and annotated genome of the parasitic trematode Atriophallophorus winterbourni and its comparative genomic analysis to other major parasitic trematodes. First, we reconstructed the species phylogeny, and dated the split of A. winterbourni from the Opisthorchiata suborder to approximately 237.4 Ma (±120.4 Myr). We then addressed the question of which expanded gene families and gained genes are potentially involved in adaptation to parasitism. To do this, we used hierarchical orthologous groups to reconstruct three ancestral genomes on the phylogeny leading to A. winterbourni and performed a GO (Gene Ontology) enrichment analysis of the gene composition of each ancestral genome, allowing us to characterize the subsequent genomic changes. Out of the 11,499 genes in the A. winterbourni genome, as much as 24% have arisen through duplication events since the speciation of A. winterbourni from the Opisthorchiata, and as much as 31.9% appear to be novel, that is, newly acquired. We found 13 gene families in A. winterbourni to have had more than ten genes arising through these recent duplications; all of which have functions potentially relating to host behavioral manipulation, host tissue penetration, and hiding from host immunity through antigen presentation. We identified several families with genes evolving under positive selection. Our results provide a valuable resource for future studies on the genomic basis of adaptation to parasitism and point to specific candidate genes putatively involved in antagonistic host-parasite adaptation.

Intraspecific genetic variation for anesthesia success in a New Zealand freshwater snail

Intraspecific genetic variation can drive phenotypic variation even across very closely related individuals. Here, we demonstrate that genetic differences between snails are a major contributor to wide variation in menthol anesthesia success in an important freshwater snail model system, Potamopyrgus antipodarum. Anesthesia is used to immobilize organisms for experiments and surgical procedures and to humanely mitigate pain. This is the first example of which we are aware of a role for genetic variation in anesthesia success in a mollusk. These findings highlight the fact that using only one strain or lineage for many experiments will not provide a full picture of phenotypic variation, demonstrate the importance of optimizing biomedically relevant techniques and protocols across a variety of genetic backgrounds, illuminate a potential mechanism underlying previously documented challenges in molluscan anesthesia, and set the stage for powerful and humane manipulative experiments in P. antipodarum.

DNA Content Variation and SNP Diversity Within a Single Population of Asexual Snails

A growing body of research suggests that many clonal populations maintain genetic diversity even without occasional sexual reproduction. The purpose of our study was to document variation in single-nucleotide polymorphism (SNP) diversity, DNA content, and pathogen susceptibility in clonal lineages of the New Zealand freshwater snail, Potamopyrgus antipodarum. We studied snails that were collected from multiple field sites around a single lake (Lake Alexandrina), as well as isofemale clonal lineages that had been isolated and maintained in the laboratory. We used the kompetitive allele specific PCR (KASP) method to genotype our samples at 46 nuclear SNP sites, and we used flow cytometry to estimate DNA content. We found high levels of SNP diversity, both in our field samples and in our clonal laboratory lines. We also found evidence of high variation in DNA content among clones, even among clones with identical genotypes across all SNP sites. Controlled pathogen exposures of the laboratory populations revealed variation in susceptibility among distinct clonal genotypes, which was independent of DNA content. Taken together, these results show high levels of diversity among asexual snails, especially for DNA content, and they suggest rapid genome evolution in asexuals.

Skipping sex: A nonrecombinant genomic assemblage of complementary reproductive modules

The unusual occurrence and developmental diversity of asexual eukaryotes remain a puzzle. De novo formation of a functioning asexual genome requires a unique assembly of sets of genes or gene states to disrupt cellular mechanisms of meiosis and gametogenesis, and to affect discrete components of sexuality and produce clonal or hemiclonal offspring. We highlight two usually overlooked but essential conditions to understand the molecular nature of clonal organisms, that is, a nonrecombinant genomic assemblage retaining modifiers of the sexual program, and a complementation between altered reproductive components. These subtle conditions are the basis for physiologically viable and genetically balanced transitions between generations. Genomic and developmental evidence from asexual animals and plants indicates the lack of complementation of molecular changes in the sexual reproductive program is likely the main cause of asexuals' rarity, and can provide an explanatory frame for the developmental diversity and lability of developmental patterns in some asexuals as well as for the discordant time to extinction estimations.

A layover in Europe: Reconstructing the invasion route of asexual lineages of a New Zealand snail to North America

Non-native invasive species are threatening ecosystems and biodiversity worldwide. High genetic variation is thought to be a critical factor for invasion success. Accordingly, the global invasion of a few clonal lineages of the gastropod Potamopyrgus antipodarum is thus both puzzling and has the potential to help illuminate why some invasions succeed while others fail. Here, we used SNP markers and a geographically broad sampling scheme (N = 1617) including native New Zealand populations and invasive North American and European populations to provide the first widescale population genetic assessment of the relationships between and among native and invasive P. antipodarum. We used a combination of traditional and Bayesian molecular analyses to demonstrate that New Zealand populations harbour very high diversity relative to the invasive populations and are the source of the two main European genetic lineages. One of these two European lineages was in turn the source of at least one of the two main North American genetic clusters of invasive P. antipodarum, located in Lake Ontario. The other widespread North American group had a more complex origin that included the other European lineage and two New Zealand clusters. Altogether, our analyses suggest that just a small handful of clonal lineages of P. antipodarum were responsible for invasion across continents. Our findings provide critical information for prevention of additional invasions and control of existing invasive populations and are of broader relevance towards understanding the establishment and evolution of asexual populations and the forces driving biological invasion.

Male phenotypes in a female framework: Evidence for degeneration in sperm produced by male snails from asexual lineages

How changes in selective regimes affect trait evolution is an important open biological question. We take advantage of naturally occurring and repeated transitions from sexual to asexual reproduction in a New Zealand freshwater snail species, Potamopyrgus antipodarum, to address how evolution in an asexual context-including the potential for relaxed selection on male-specific traits-influences sperm morphology. The occasional production of male offspring by the otherwise all-female asexual P. antipodarum lineages affords a unique and powerful opportunity to assess the fate of sperm traits in a context where males are exceedingly rare. These comparisons revealed that the sperm produced by 'asexual' males are markedly distinct from sexual counterparts. We also found that the asexual male sperm harboured markedly higher phenotypic variation and was much more likely to be morphologically abnormal. Together, these data suggest that transitions to asexual reproduction might be irreversible, at least in part because male function is likely to be compromised. These results are also consistent with a scenario where relaxed selection and/or mutation accumulation in the absence of sex translates into rapid trait degeneration.

Regional epigenetic variation in asexual snail populations among urban and rural lakes

Epigenetic variation has the potential to influence environmentally dependent development and contribute to phenotypic responses to local environments. Environmental epigenetic studies of sexual organisms confirm the capacity to respond through epigenetic variation. An epigenetic response could be even more important in a population when genetic variation is lacking. A previous study of an asexual snail, Potamopyrgus antipodarum, demonstrated that different populations derived from a single clonal lineage differed in both shell phenotype and methylation signature when comparing lake versus river populations. Here, we examine methylation variation among lakes that differ in environmental disturbance and pollution histories. Snails were collected from a more pristine rural Lake 1 (Lake Lytle), and two urban lakes, Lake 2 (Capitol Lake) and Lake 3 (Lake Washington) on the Northwest Pacific coast. DNA methylation was assessed for each sample population using methylated DNA immunoprecipitation, MeDIP, followed by next-generation sequencing. The differential DNA methylation regions (DMRs) identified among the different lake comparisons suggested a higher number of DMRs and variation between rural Lake 1 and one urban Lake 2, and between the two urban Lakes 2 and 3, but limited variation between the rural Lake 1 and urban Lake 3. DMR genomic characteristics and gene associations were investigated. The presence of site-specific differences between each of these lake populations suggest an epigenetic response to varied environmental factors. The results do not support an effect of geographic distance in these populations. The role of dispersal distance among lakes, population history, environmental pollution and stably inherited methylation versus environmentally triggered methylation in producing the observed epigenetic variation are discussed. Observations support the proposal that epigenetic alterations may associate with phenotypic variation and environmental factors and history of the different lakes.

Genetic Variation for Mitochondrial Function in the New Zealand Freshwater Snail Potamopyrgus antipodarum

The proteins responsible for mitochondrial function are encoded by 2 different genomes with distinct inheritance regimes, rendering rigorous inference of genotype-phenotype connections intractable for all but a few model systems. Asexual organisms provide a powerful means to address these challenges because offspring produced without recombination inherit both nuclear and mitochondrial genomes from a single parent. As such, these offspring inherit mitonuclear genotypes that are identical to the mitonuclear genotypes of their parents and siblings but different from those of other asexual lineages. Here, we compared mitochondrial function across distinct asexual lineages of Potamopyrgus antipodarum, a New Zealand freshwater snail model for understanding the evolutionary consequences of asexuality. Our analyses revealed substantial phenotypic variation across asexual lineages at 3 levels of biological organization: mitogenomic, organellar, and organismal. These data demonstrate that different asexual lineages have different mitochondrial function phenotypes, likely reflecting heritable variation (i.e., the raw material for evolution) for mitochondrial function in P. antipodarum. The discovery of this variation combined with the methods developed here sets the stage to use P. antipodarum to study central evolutionary questions involving mitochondrial function, including whether mitochondrial mutation accumulation influences the maintenance of sexual reproduction in natural populations.

Adaptive phenotypic plasticity in a clonal invader

Organisms featuring wide trait variability and occurring in a wide range of habitats, such as the ovoviviparous New Zealand freshwater snail Potamopyrgus antipodarum, are ideal models to study adaptation. Since the mid‐19th century, P. antipodarum, characterized by extremely variable shell morphology, has successfully invaded aquatic areas on four continents. Because these obligately and wholly asexual invasive populations harbor low genetic diversity compared to mixed sexual/asexual populations in the native range, we hypothesized that (1) this phenotypic variation in the invasive range might be adaptive with respect to colonization of novel habitats, and (2) that at least some of the variation might be caused by phenotypic plasticity. We surveyed 425 snails from 21 localities across northwest Europe to attempt to disentangle genetic and environmental effects on shell morphology. We analyzed brood size as proxy for fitness and shell geometric morphometrics, while controlling for genetic background. Our survey revealed 10 SNP genotypes nested into two mtDNA haplotypes and indicated that mainly lineage drove variation in shell shape but not size. Physicochemical parameters affected both shell shape and size and the interaction of these traits with brood size. In particular, stronger stream flow rates were associated with larger shells. Our measurements of brood size suggested that relatively larger slender snails with relatively large apertures were better adapted to strong flow than counterparts with broader shells and relatively small apertures. In conclusion, the apparent potential to modify shell morphology plays likely a key role in the invasive success of P. antipodarum; the two main components of shell morphology, namely shape and size, being differentially controlled, the former mainly genetically and the latter predominantly by phenotypic plasticity.

Radical amino acid mutations persist longer in the absence of sex

Harmful mutations are ubiquitous and inevitable, and the rate at which these mutations are removed from populations is a critical determinant of evolutionary fate. Closely related sexual and asexual taxa provide a particularly powerful setting to study deleterious mutation elimination because sexual reproduction should facilitate mutational clearance by reducing selective interference between sites and by allowing the production of offspring with different mutational complements than their parents. Here, we compared the rate of removal of conservative (i.e., similar biochemical properties) and radical (i.e., distinct biochemical properties) nonsynonymous mutations from mitochondrial genomes of sexual versus asexual Potamopyrgus antipodarum, a New Zealand freshwater snail characterized by coexisting and ecologically similar sexual and asexual lineages. Our analyses revealed that radical nonsynonymous mutations are cleared at higher rates than conservative changes and that sexual lineages eliminate radical changes more rapidly than asexual counterparts. These results are consistent with reduced efficacy of purifying selection in asexual lineages allowing harmful mutations to remain polymorphic longer than in sexual lineages. Together, these data illuminate some of the population-level processes contributing to mitochondrial mutation accumulation and suggest that mutation accumulation could influence the outcome of competition between sexual and asexual lineages.

Population genomics and geographical parthenogenesis in Japanese harvestmen (Opiliones, Sclerosomatidae, Leiobunum)

Naturally occurring population variation in reproductive mode presents an opportunity for researchers to test hypotheses regarding the evolution of sex. Asexual reproduction frequently assumes a geographical pattern, in which parthenogenesis-dominated populations are more broadly dispersed than their sexual conspecifics. We evaluate the geographical distribution of genomic signatures associated with parthenogenesis using nuclear and mitochondrial DNA sequence data from two Japanese harvestman sister taxa, Leiobunum manubriatum and Leiobunum globosum. Asexual reproduction is putatively facultative in these species, and female-biased localities are common in habitat margins. Past karyotypic and current cytometric work indicates L. globosum is entirely tetraploid, while L. manubriatum may be either diploid or tetraploid. We estimated species phylogeny, genetic differentiation, diversity, and mitonuclear discordance in females collected across the species range in order to identify range expansion toward marginal habitat, potential for hybrid origin, and persistence of asexual lineages. Our results point to northward expansion of a tetraploid ancestor of L. manubriatum and L. globosum, coupled with support for greater male gene flow in southern L. manubriatum localities. Specimens from localities in the Tohoku and Hokkaido regions were indistinct, particularly those of L. globosum, potentially due to little mitochondrial differentiation or haplotypic variation. Although L. manubriatum overlaps with L. globosum across its entire range, L. globosum was reconstructed as monophyletic with strong support using mtDNA, and marginal support with nuclear loci. Ultimately, we find evidence for continued sexual reproduction in both species and describe opportunities to clarify the rate and mechanism of parthenogenesis.

Phosphorus availability in the source population influences response to dietary phosphorus quantity in a New Zealand freshwater snail

We investigated whether previously documented variation among populations in availability of dietary phosphorus (P) is linked to heterogeneity in growth rate of the New Zealand freshwater snail Potamopyrgus antipodarum on a P-limited diet. We chose this system because P. antipodarum inhabits water bodies that vary in P availability and because P. antipodarum growth rate varies considerably in response to low P. We quantified specific growth rate and alkaline phosphatase (AP) expression in a diverse array of juvenile P. antipodarum fed high vs. low-P diets. We found strong associations between P content of epilithon in the source lake and P. antipodarum growth rate on high vs. low-P diets, with snails collected from lakes with relatively low-P epilithon showing the greatest increase in growth rate on the high-P relative to low-P diet. We also found substantial intraspecific variation in growth response to P limitation. Expression of AP also varied among lineages and was negatively associated with C: P of lake epilithon but did not explain the relationship between C: P in the lake of origin and sensitivity to P limitation. Together, our results demonstrate a strong signature of the P environment in the lake of origin on how this snail responds to P limitation as well as preliminary evidence for intraspecific variation of AP expression in animals.

Signatures of mito-nuclear discordance in Schistosoma turkestanicum indicate a complex evolutionary history of emergence in Europe

High levels of molecular diversity were identified in mitochondrial cytochrome c oxidase (cox1) gene sequences of Schistosoma turkestanicum from Hungary. These cox1 sequences were all specific to Hungary which contrasted with the low levels of diversity seen in the nuclear internal transcribed spacer region (ITS) sequences, the majority of which were shared between China and Iran isolates. Measures of within and between host molecular variation within S. turkestanicum showed there to be substantial differences in molecular diversity, with cox1 being significantly more diverse than the ITS. Measures of haplotype frequencies revealed that each host contained its own subpopulation of genetically unique parasites with significant levels of differentiation. Pairwise mismatch analysis of cox1 sequences indicated S. turkestanicum populations to have a bimodal pairwise difference distribution and to be stable unlike the ITS sequences, which appeared to have undergone a recent population expansion event. Positive selection was also detected in the cox1 sequences, and biochemical modelling of the resulting protein illustrated significant mutational events causing an alteration to the isoelectric point of the cox1 protein, potentially altering metabolism. The evolutionary signature from the cox1 indicates local adaptation and long establishment of S. turkestanicum in Hungary with continual introgression of nuclear genes from Asian isolates. These processes have led to the occurrence of mito-nuclear discordance in a schistosome population.

Genomic evidence for population-specific responses to co-evolving parasites in a New Zealand freshwater snail

Reciprocal co-evolving interactions between hosts and parasites are a primary source of strong selection that can promote rapid and often population- or genotype-specific evolutionary change. These host-parasite interactions are also a major source of disease. Despite their importance, very little is known about the genomic basis of co-evolving host-parasite interactions in natural populations, especially in animals. Here, we use gene expression and sequence evolution approaches to take critical steps towards characterizing the genomic basis of interactions between the freshwater snail Potamopyrgus antipodarum and its co-evolving sterilizing trematode parasite, Microphallus sp., a textbook example of natural coevolution. We found that Microphallus-infected P. antipodarum exhibit systematic downregulation of genes relative to uninfected P. antipodarum. The specific genes involved in parasite response differ markedly across lakes, consistent with a scenario where population-level co-evolution is leading to population-specific host-parasite interactions and evolutionary trajectories. We also used an F_ST -based approach to identify a set of loci that represent promising candidates for targets of parasite-mediated selection across lakes as well as within each lake population. These results constitute the first genomic evidence for population-specific responses to co-evolving infection in the P. antipodarum-Microphallus interaction and provide new insights into the genomic basis of co-evolutionary interactions in nature.

Sensitivity to dietary phosphorus limitation in native vs. invasive lineages of a New Zealand freshwater snail

Why some species and lineages are more likely to be invasive than others is one of the most important unanswered questions in basic and applied biology. In particular, the relative contributions to the invasion process of factors like pre-adaptation to invasiveness in the native range, evolution post-colonization, and random vs. non-random sampling of colonist lineages remain unclear. We use a powerful common garden approach to address the potential for a role for sensitivity to nutrient limitation in determining the invasiveness of particular lineages of Potamopyrgus antipodarum, a New Zealand freshwater snail that has become globally invasive. We quantified specific growth rate (SGR), an important fitness-related trait in this species, under high phosphorus (P) vs. low-P conditions for a diverse set of native and invasive P. antipodarum. This experiment revealed that native-range P. antipodarum experience a more severe decline in SGR in low-P conditions relative to SGR in high-P conditions than their invasive range counterparts. Although these results suggest resilience to P limitation in invasive lineages, the absence of significant absolute differences in SGR between native and invasive lineages indicates that a straightforward connection between response to P limitation and invasiveness in P. antipodarum is unlikely. Regardless, our data demonstrate that invasive vs. native lineages of P. antipodarum exhibit consistently different responses to an important environmental variable that is rarely studied in the context of invasion success. Further studies directed at exploring and disentangling the roles of sampling effects, selection on preexisting variation, and evolution after colonization will be required to provide a comprehensive picture of the role (or lack thereof) of nutrient limitation in the global invasion of P. antipodarum, as well for as other invasive taxa.

Distinct Bacterial Microbiomes in Sexual and Asexual Potamopyrgus antipodarum, a New Zealand Freshwater Snail

Different reproductive strategies and the transition to asexuality can be associated with microbial symbionts. Whether such a link exists within mollusks has never been evaluated. We took the first steps towards addressing this possibility by performing pyrosequencing of bacterial 16S rRNA genes associated with Potamopyrgus antipodarum, a New Zealand freshwater snail. A diverse set of 60 tissue collections from P. antipodarum that were genetically and geographically distinct and either obligately sexual or asexual were included, which allowed us to evaluate whether reproductive mode was associated with a particular bacterial community. 2624 unique operational taxonomic units (OTU, 97% DNA similarity) were detected, which were distributed across ~30 phyla. While alpha diversity metrics varied little among individual samples, significant differences in bacterial community composition and structure were detected between sexual and asexual snails, as well as among snails from different lakes and genetic backgrounds. The mean dissimilarity of the bacterial communities between the sexual and asexual P. antipodarum was 90%, largely driven by the presence of Rickettsiales in sexual snails and Rhodobacter in asexual snails. Our study suggests that there might be a link between reproductive mode and the bacterial microbiome of P. antipodarum, though a causal connection requires additional study.

Effects of polyploidy and reproductive mode on life history trait expression

Ploidy elevation is increasingly recognized as a common and important source of genomic variation. Even so, the consequences and biological significance of polyploidy remain unclear, especially in animals. Here, our goal was to identify potential life history costs and benefits of polyploidy by conducting a large multiyear common garden experiment in Potamopyrgus antipodarum, a New Zealand freshwater snail that is a model system for the study of ploidy variation, sexual reproduction, host-parasite coevolution, and invasion ecology. Sexual diploid and asexual triploid and tetraploid P. antipodarum frequently coexist, allowing for powerful direct comparisons across ploidy levels and reproductive modes. Asexual reproduction and polyploidy are very often associated in animals, allowing us to also use these comparisons to address the maintenance of sex, itself one of the most important unresolved questions in evolutionary biology. Our study revealed that sexual diploid P. antipodarum grow and mature substantially more slowly than their asexual polyploid counterparts. We detected a strong negative correlation between the rate of growth and age at reproductive maturity, suggesting that the relatively early maturation of asexual polyploid P. antipodarum is driven by relatively rapid growth. The absence of evidence for life history differences between triploid and tetraploid asexuals indicates that ploidy elevation is unlikely to underlie the differences in trait values that we detected between sexual and asexual snails. Finally, we found that sexual P. antipodarum did not experience discernable phenotypic variance-related benefits of sex and were more likely to die before achieving reproductive maturity than the asexuals. Taken together, these results suggest that under benign conditions, polyploidy does not impose obvious life history costs in P. antipodarum and that sexual P. antipodarum persist despite substantial life history disadvantages relative to their asexual counterparts.

Response to phosphorus limitation varies among lake populations of the freshwater snail Potamopyrgus antipodarum

Local adaptation – typically recognized as higher values of fitness-related traits for native vs. non-native individuals when measured in the native environment - is common in natural populations because of pervasive spatial variation in the intensity and type of natural selection. Although local adaptation has been primarily studied in the context of biotic interactions, widespread variation in abiotic characteristics of environments suggests that local adaptation in response to abiotic factors should also be common. Potamopyrgus antipodarum, a freshwater New Zealand snail that is an important model system for invasion biology and the maintenance of sexual reproduction, exhibits local adaptation to parasites and rate of water flow. As an initial step to determining whether P. antipodarum are also locally adapted to phosphorus availability, we examined whether populations differ in their responses to phosphorus limitation. We found that field-collected juvenile P. antipodarum grew at a lower rate and reached an important size threshold more slowly when fed a relatively low vs. a relatively high- phosphorus diet. We also detected significant across-population variation in individual growth rate. A marginally significant population-by-dietary phosphorus interaction along with a two-fold difference across populations in the extent of suppression of growth by low phosphorus suggests that populations of P. antipodarum may differ in their response to phosphorus limitation. Local adaptation may explain this variation, with the implication that snails from lakes with relatively low phosphorus availability should be less severely affected by phosphorus limitation than snails from lakes with higher phosphorus availability.

Higher rate of tissue regeneration in polyploid asexual versus diploid sexual freshwater snails

Characterizing phenotypic differences between sexual and asexual organisms is a critical step towards understanding why sexual reproduction is so common. Because asexuals are often polyploid, understanding how ploidy influences phenotype is directly relevant to the study of sex and will provide key insights into the evolution of ploidy-level variation. The well-established association between genome size and cell cycle duration, evidence for a link between genome size and tissue regeneration rate and the growing body of research showing that ploidy influences growth rate and gene expression led us to hypothesize that healing and tissue regeneration might be affected by ploidy-level variation. We evaluated this hypothesis by measuring the rate of regeneration of antenna tissue of Potamopyrgus antipodarum, a New Zealand snail characterized by frequent coexistence between diploid sexuals and polyploid asexuals. Antennae of triploid and presumptive tetraploid asexuals regenerated more rapidly than the antennae of diploid sexuals, but regeneration rate did not differ between triploids and tetraploids. These results suggest either that ploidy elevation has nonlinear positive effects on tissue regeneration and/or that factors associated directly with reproductive mode affect regeneration rate more than ploidy level. The results of this study also indicate that the lower ploidy of sexual P. antipodarum is unlikely to confer advantages associated with more rapid regeneration.