Population-genomic analysis identifies a low rate of global adaptive fixation in the proteins of the cyclical parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared to Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction

Uncovering cryptic asexuality in Daphnia magna by RAD sequencing

The breeding systems of many organisms are cryptic and difficult to investigate with observational data, yet they have profound effects on a species’ ecology, evolution, and genome organization. Genomic approaches offer a novel, indirect way to investigate breeding systems, specifically by studying the transmission of genetic information from parents to offspring. Here we exemplify this method through an assessment of self-fertilization vs. automictic parthenogenesis in Daphnia magna. Self-fertilization reduces heterozygosity by 50% compared to the parents, but under automixis, whereby two haploid products from a single meiosis fuse, the expected heterozygosity reduction depends on whether the two meiotic products are separated during meiosis I or II (i.e., central vs. terminal fusion). Reviewing the existing literature and incorporating recombination interference, we derive an interchromosomal and an intrachromosomal prediction of how to distinguish various forms of automixis from self-fertilization using offspring heterozygosity data. We then test these predictions using RAD-sequencing data on presumed automictic diapause offspring of so-called nonmale producing strains and compare them with “self-fertilized” offspring produced by within-clone mating. The results unequivocally show that these offspring were produced by automixis, mostly, but not exclusively, through terminal fusion. However, the results also show that this conclusion was only possible owing to genome-wide heterozygosity data, with phenotypic data as well as data from microsatellite markers yielding inconclusive or even misleading results. Our study thus demonstrates how to use the power of genomic approaches for elucidating breeding systems, and it provides the first demonstration of automictic parthenogenesis in Daphnia

Genetic parallelism between European flat oyster populations at the edge of their natural range

Although all marine ecosystems have experienced global-scale losses, oyster reefs have shown the greatest. Therefore, substantial efforts have been dedicated to restoration of such ecosystems during the last two decades. In Europe, several pilot projects for the restoration of the native European flat oyster, Ostrea edulis, recently begun and recommendations to preserve genetic diversity and to conduct monitoring protocols have been made. In particular, an initial step is to test for genetic differentiation against homogeneity among the oyster populations potentially involved in such programs. Therefore, we conducted a new sampling of wild populations at the European scale and a new genetic analysis with 203 markers to (1) confirm and study in more detail the pattern of genetic differentiation between Atlantic and Mediterranean populations, (2) identify potential translocations that could be due to aquaculture practices and (3) investigate the populations at the fringe of the geographical range, since they seemed related despite their geographic distance. Such information should be useful to enlighten the choice of the animals to be translocated or reproduced in hatcheries for further restocking. After the confirmation of the general geographical pattern of genetic structure and the identification of one potential case of aquaculture transfer at a large scale, we were able to detect genomic islands of differentiation mainly in the form of two groups of linked markers, which could indicate the presence of polymorphic chromosomal rearrangements. Furthermore, we observed a tendency for these two islands and the most differentiated loci to show a parallel pattern of differentiation, grouping the North Sea populations with the Eastern Mediterranean and Black Sea populations, against geography. We discussed the hypothesis that this genetic parallelism could be the sign of a shared evolutionary history of the two groups of populations despite them being at the border of the distribution nowadays

Microplastics Affect Energy Balance and Gametogenesis in the Pearl Oyster Pinctada margaritifera

Plastic pollution in the environment is increasing at global scale. Microplastics (MP) are derived from degradation of larger plastic items or directly produced in microparticles form (<5 mm). Plastics, widely used in structures and equipments of pearl farming, are a source of pollution to the detriment of the lagoon ecosystem. In order to evaluate the impact of MP on the physiology of Pinctada margaritifera, a species of ecological and commercial interests, adult oysters were exposed to polystyrene microbeads (micro-PS of 6 and 10 μm) for 2 months. Three concentrations were tested: 0.25, 2.5, 25 μg L-1 and a control. Ingestion and respiration rate and assimilation efficiency were monitored on a metabolic measurement system to determine the individual energy balance (Scope For Growth, SFG). Effects on reproduction were also assessed. The assimilation efficiency decreased significantly according to micro-PS concentration. The SFG was significantly impacted by a dose-dependent decrease from 0.25 μg L-1 (p < 0.0001). A negative SFG was measured in oysters exposed to 25 μg L-1. Gonads may have provided the missing energy to maintain animals’ metabolism through the production of metabolites derived from germ cells phagocytosis. This study shows that micro-PS significantly impact the assimilation efficiency and more broadly the energy balance of P. margaritifera, with negative repercussions on reproduction

Understanding connectivity of pearl oyster populations within Tuamotu atoll semi-closed lagoons: Cumulative insight from genetics and biophysical modelling approaches

Connectivity affects species demography, (meta)population dynamics, evolution, phylogeny and biogeography. Various methodological approaches are applied to measure connectivity. Biophysical modelling can explore systematically the influence of atmospheric, oceanic and ecological forcing, while genetics measures connectivity patterns within the sampling strategy limit. In the Pacific Ocean pearl farming lagoons, the activity relies on spat collecting of the black lipped pearl oyster Pinctada margaritifera occurring after the larval dispersal phase, which follows spawning from wild or farmed populations. Biophysical 3D modelling and genomic studies have both separately brought insights on within-lagoon connectivity and on the origin of spats. Here, we combined previous genetics results with new realistic biophysical modelling scenarios to elucidate connectivity in Ahe Atoll lagoon. When combined, we identified the weather sequence likely explaining the realized connectivity observations. We discuss the strengths, weaknesses, opportunities and threats of combining these two approaches considering specific pearl farming demographic connectivity questions

Colour plasticity in the shells and pearls of animal graft model Pinctada margaritifera assessed by HSV colour quantification

The bivalve Pinctada margaritifera has the capacity to produce the most varied and colourful pearls in the world. Colour expression in the inner shell is under combined genetic and environmental control and is correlated with the colour of pearls produced when the same individual is used as a graft donor. One major limitation when studying colour phenotypes is grader subjectivity, which leads to inconsistent colour qualification and quantification. Through the use of HSV (Hue Saturation Value) colour space, we created an R package named ‘ImaginR’ to characterise inner shell colour variations in P. margaritifera. Using a machine-learning protocol with a training dataset, ImaginR was able to reassign individual oysters and pearls to predefined human-based phenotype categories. We then tested the package on samples obtained in an experiment testing the effects of donor conditioning depth on the colour of the donor inner shell and colour of the pearls harvested from recipients following grafting and 20 months of culture in situ. These analyses successfully detected donor shell colour modifications due to depth-related plasticity and the maintenance of these modifications through to the harvested pearls. Besides its potential interest for standardization in the pearl industry, this new method is relevant to other research projects using biological models

Evolution of Gene Expression during a Transition from Environmental to Genetic Sex Determination

Genetic sex determination (GSD) can evolve from environmental sex determination (ESD) via an intermediate state in which both coexist in the same population. Such mixed populations are found in the crustacean Daphnia magna, where non-male producers (NMP, genetically determined females) coexist with male producers (MP), in which male production is environmentally inducible and can also artificially be triggered by exposure to juvenile hormone. This makes Daphnia magna a rare model species for the study of evolutionary transitions from ESD to GSD. Although the chromosomal location of the NMP-determining mutation has been mapped, the actual genes and pathways involved in the evolution of GSD from ESD remain unknown. Here, we present a transcriptomic analysis of MP and NMP females under control (female producing) and under hormone exposure conditions. We found similar to 100 differentially expressed genes between MP and NMP under control conditions. Genes in the NMP-determining chromosome region were especially likely to show such constitutive expression differences. Hormone exposure led to expression changes of an additional similar to 100 (MP) to similar to 600 (NMP) genes, with an almost systematic upregulation of those genes in NMP. These observations suggest that the NMP phenotype is not determined by a simple "loss-of-function" mutation. Rather, homeostasis of female offspring production under hormone exposure appears to be an active state, tightly regulated by complex mechanisms involving many genes. In a broader view, this illustrates that the evolution of GSD, while potentially initiated by a single mutation, likely leads to secondary integration involving many genes and pathways

Strong genetic isolation of the black-lipped pearl oyster (Pinctada margaritifera) in the Marquesas archipelago (French Polynesia)

The French Polynesian islands are internationally known for their black pearls, produced by culture of the black lipped pearl oyster Pinctada margaritifera. The ongoing development of hatcheries for P. margaritifera in French Polynesia poses new challenges for the industry, particularly regarding the maintenance of genetic diversity in the hatchery stocks. This emphasizes the necessity to characterize the genetic diversity and differentiation within natural and exploited populations, to carefully select putative parental populations. The present study aimed at validating the phylogenetic status and investigating genetic attributes of populations from the only two non-exploited archipelagos of French Polynesia, the Marquesas archipelago, and the Australes archipelago, never analysed before. We found that individuals from both archipelagos belonged to P. margaritifera species. However, while the Australes population was genetically similar to non-exploited populations of the Tuamotu, the Marquesas populations were highly differentiated from the rest of the populations. This differentiation cannot not be only attributed to geographic distance and aquaculture status, but likely to hydrodynamic barriers allowing vicariant events to take place. Our results add up to other studies describing the Marquesas archipelago as a hotspot for biodiversity and differentiation, with some of the highest levels of endemism and vicariance found among marine species worldwide and provide precious information on available genetic resources for the implementation of P. margaritifera selective breeding and its genetic conservation in French Polynesia

Population Connectivity and Genetic Assessment of Exploited and Natural Populations of Pearl Oysters within a French Polynesian Atoll Lagoon

In French Polynesia, the production and exportation of black pearls through the aquaculture of the black-lip pearl oyster Pinctada margaritifera provide the second largest economic income for the country after tourism. This industry entirely relies on the collection of natural spats from few highly recruiting lagoons. In recent years, pearl oyster producers have experienced variable success rates in spat collection, with significant spatial and temporal variability in spat supply, driving uncertainty in the future of pearl production. This study combines, for the first time in a farmed lagoon, genetic (SNPs), demographic (sex ratio, age), and biophysical data (larval dispersal modelling) to shed new light on population dynamics, connectivity, and spat recruitment in Ahe Atoll, a well-studied pearl farming site. Our results indicate that the geographical structuring of the natural populations and the contribution of both natural and exploited stocks to the production of spats result from the interaction of hydrodynamic features, life history traits and demographic parameters: the northeastern natural populations are older, not well connected to the southwestern natural populations and are not replenished by larvae produced by adjacent exploited populations. Moreover, we observe that the exploited populations did not contribute to larval production during our experiment, despite a sampling period set during the most productive season for spat collection. This is likely the result of a strong male bias in the exploited populations, coupled with a sweepstakes reproductive strategy of the species. Our results warrant further investigations over the future of the northeastern older natural populations and a reflection on the current perliculture techniques

Transition from Environmental to Partial Genetic Sex Determination in Daphnia through the Evolution of a Female-Determining Incipient W Chromosome

Sex chromosomes can evolve during the evolution of genetic sex determination (GSD) from environmental sex determination (ESD). Despite theoretical attention, early mechanisms involved in the transition from ESD to GSD have yet to be studied in nature. No mixed ESD-GSD animal species have been reported, except for some species of Daphnia, small freshwater crustaceans in which sex is usually determined solely by the environment, but in which a dominant female sex-determining locus is present in some populations. This locus follows Mendelian single-locus inheritance, but has otherwise not been characterized genetically. We now show that the sex-determining genomic region maps to the same low-recombining peri-centromeric region of linkage group 3 (LG3) in three highly divergent populations of D. magna, and spans 3.6 Mb. Despite low levels of recombination, the associated region contains signs of historical recombination, suggesting a role for selection acting on several genes thereby maintaining linkage disequilibrium among the 36 associated SNPs. The region carries numerous genes involved in sex differentiation in other taxa, including transformer2 and sox9. Taken together, the region determining the genetic females shows characteristics of a sex-related supergene, suggesting that LG3 is potentially an incipient W chromosome despite the lack of significant additional restriction of recombination between Z and W. The occurrence of the female-determining locus in a pre-existing low recombining region illustrates one possible form of recombination suppression in sex chromosomes. D. magna is a promising model for studying the evolutionary transitions from ESD to GSD and early sex chromosome evolution