Rare sex or out of reach equilibrium? The dynamics of F IS in partially clonal organisms

Bumblebees mediate landscape effects on a forest herb's population genetic structure in European agricultural landscapes

Spatially isolated plant populations in agricultural landscapes exhibit genetic responses not only to habitat fragmentation per se but also to the composition of the landscape matrix between habitat patches. These responses can only be understood by examining how the landscape matrix influences among-habitat movements of pollinators and seed vectors, which act as genetic linkers among populations. We studied the forest herb Polygonatum multiflorum and its associated pollinator and genetic linker, the bumblebee Bombus pascuorum, in three European agricultural landscapes. We aimed to identify which landscape features affect the movement activity of B. pascuorum between forest patches and to assess the relative importance of these features in explaining the forest herb's population genetic structure. We applied microsatellite markers to estimate the movement activity of the bumblebee as well as the population genetic structure of the forest herb. We modelled the movement activity as a function of various landscape metrics. Those metrics found to explain the movement activity best were then used to explain the population genetic structure of the forest herb. The bumblebee movement activity was affected by the cover of maize fields and semi-natural grasslands on a larger spatial scale and by landscape heterogeneity on a smaller spatial scale. For some measures of the forest herb's population genetic structure, that is, allelic richness, observed heterozygosity and the F-value, the combinations of landscape metrics, which explained the linker movement activity best, yielded lower AICc values than 95% of the models including all possible combinations of landscape metrics. Synthesis: The genetic linker, B. pascuorum, mediates landscape effects on the population genetic structure of the forest herb P. multiflorum. Our study indicates, that the movement of the genetic linker among forest patches, and thus the pollen driven gene flow of the herb, depends on the relative value of floral resources in the specific landscape setting. Noteworthy, the population genetic structure of the long-lived, clonal forest herb species correlated with recent land-use types such as maize, which have been existing for not more than a few decades within these landscapes. This underscores the short time in which land-use changes can influence the evolutionary potential of long-lived wild plants.

Approximate Bayesian Computation applied to time series of population genetic data disentangles rapid genetic changes and demographic variations in a pathogen population

Adaptation can occur at remarkably short timescales in natural populations, leading to drastic changes in phenotypes and genotype frequencies over a few generations only. The inference of demographic parameters can allow understanding how evolutionary forces interact and shape the genetic trajectories of populations during rapid adaptation. Here we propose a new Approximate Bayesian Computation (ABC) framework that couples a forward and individual-based model with temporal genetic data to disentangle genetic changes and demographic variations in a case of rapid adaptation. We test the accuracy of our inferential framework and evaluate the benefit of considering a dense versus sparse sampling. Theoretical investigations demonstrate high accuracy in both model and parameter estimations, even if a strong thinning is applied to time series data. Then, we apply our ABC inferential framework to empirical data describing the population genetic changes of the poplar rust pathogen following a major event of resistance overcoming. We successfully estimate key demographic and genetic parameters, including the proportion of resistant hosts deployed in the landscape and the level of standing genetic variation from which selection occurred. Inferred values are in accordance with our empirical knowledge of this biological system. This new inferential framework, which contrasts with coalescent-based ABC analyses, is promising for a better understanding of evolutionary trajectories of populations subjected to rapid adaptation.

Ice holes microrefugia harbor genetically and functionally distinct populations of Vaccinium vitis-idaea (Ericaceae)

In the mountain terrain, ice holes are little depressions between rock boulders that are characterized by the exit of cold air able to cool down the rock surface even in summer. This cold air creates cold microrefugia in warmer surroundings that preserve plant species probably over thousands of years under extra-zonal climatic conditions. We hypothesized that ice hole populations of the model species Vaccinium vitis-idaea (Ericaceae) show genetic differentiation from nearby zonal subalpine populations, and high functional trait distinctiveness, in agreement with genetic patterns. We genotyped almost 30,000 single nucleotide polymorphisms using restriction site-associated DNA sequencing and measured eight functional traits indicative of individual performance and ecological strategies. Genetic results showed high differentiation among the six populations suggesting isolation. On siliceous bedrock, ice hole individuals exhibited higher levels of admixture than those from subalpine populations which could have experienced more bottlenecks during demographic fluctuations related to glacial cycles. Ice hole and subalpine calcareous populations clearly separated from siliceous populations, indicating a possible effect of bedrock in shaping genetic patterns. Trait analysis reflected the bedrock effect on populations' differentiation. The significant correlation between trait and genetic distances suggests the genetic contribution in shaping intraspecific functional differentiation. In conclusion, extra-zonal populations reveal a prominent genetic and phenotypic differentiation determined by history and ecological contingency. Therefore, microrefugia populations can contribute to the overall variability of the species and lead to intraspecific-driven responses to upcoming environmental changes.

Clones on the run: The genomics of a recently expanded partially clonal species

Why species that in their core areas mainly reproduce sexually become enriched with clones in marginal populations ("geographic parthenogenesis") remains unclear. Earlier hypotheses have emphasized that selection might promote clonality because it protects locally adapted genotypes. On the other hand, it also hampers recombination and adaptation to changing conditions. The aim of the present study was to investigate the early stages of range expansion in a partially clonal species and what drives an increase in cloning during such expansion. We used genome-wide sequencing to investigate the origin and evolution of large clones formed in a macroalgal species (Fucus vesiculosus) during a recent expansion into the postglacial Baltic Sea. We found low but persistent clonality in core populations, while at range margins, large dominant clonal lineages had evolved repeatedly from different sexual populations. A range expansion model showed that even when asexual recruitment is less favourable than sexual recruitment in core populations, repeated bottlenecks at the expansion front can establish a genetically eroded clonal wave that spreads ahead of a sexual wave into the new area. Genetic variation decreases by drift following repeated bottlenecks at the expansion front. This results in the emerging clones having low expected heterozygosity, which corroborated our empirical observations. We conclude that Baker's Law (clones being favoured by uniparental reproductive assurance in new areas) can play an important role during range expansion in partially clonal species, resulting in a complex spatiotemporal mosaic of clonal and sexual lineages that might persist during thousands of generations.

Long range gene flow beyond predictions from oceanographic transport in a tropical marine foundation species

The transport of passively dispersed organisms across tropical margins remains poorly understood. Hypotheses of oceanographic transportation potential lack testing with large scale empirical data. To address this gap, we used the seagrass species, Halodule wrightii, which is unique in spanning the entire tropical Atlantic. We tested the hypothesis that genetic differentiation estimated across its large-scale biogeographic range can be predicted by simulated oceanographic transport. The alternative hypothesis posits that dispersal is independent of ocean currents, such as transport by grazers. We compared empirical genetic estimates and modelled predictions of dispersal along the distribution of H. wrightii. We genotyped eight microsatellite loci on 19 populations distributed across Atlantic Africa, Gulf of Mexico, Caribbean, Brazil and developed a biophysical model with high-resolution ocean currents. Genetic data revealed low gene flow and highest differentiation between (1) the Gulf of Mexico and two other regions: (2) Caribbean-Brazil and (3) Atlantic Africa. These two were more genetically similar despite separation by an ocean. The biophysical model indicated low or no probability of passive dispersal among populations and did not match the empirical genetic data. The results support the alternative hypothesis of a role for active dispersal vectors like grazers.

Genetic variation and clonal diversity in floating aquatic plants: Comparative genomic analysis of water hyacinth species in their native range

Many eukaryotic organisms reproduce by sexual and asexual reproduction. Genetic diversity in populations can be strongly dependent on the relative importance of these two reproductive modes. Here, we compare the amounts and patterns of genetic diversity in related water hyacinths that differ in their propensity for clonal propagation - highly clonal Eichhornia crassipes and moderately clonal E. azurea (Pontederiaceae). Our comparisons involved genotype-by-sequencing (GBS) of 137 E. crassipes ramets from 60 locations (193,495 nucleotide sites) and 118 E. azurea ramets from 53 locations (198,343 nucleotide sites) among six hydrological basins in central South America, the native range of both species. We predicted that because of more prolific clonal propagation, E. crassipes would exhibit lower clonal diversity than E. azurea. This prediction was supported by all measures of clonal diversity that we examined. Eichhornia crassipes also had a larger excess of heterozygotes at variant sites, another signature of clonality. However, genome-wide heterozygosity was not significantly different between the species. Eichhornia crassipes had weaker spatial genetic structure and lower levels of differentiation among hydrological basins than E. azurea, probably because of higher clonality and more extensive dispersal of its free-floating life form. Our findings for E. crassipes contrast with earlier studies from the invasive range which have reported very low levels of clonal diversity and extensive geographic areas of genetic uniformity.

Population Genetics of the Blueberry Gall Midge, Dasineura oxycoccana (Diptera: Cecidomyiidae), on Blueberry and Cranberry and Testing Invasion Scenarios

We compared the population genetic structure between populations of the blueberry gall midge-Dasineura oxycoccana (Johnson) (Diptera: Cecidomyiidae)-from blueberry and cranberry and determined the genetic relationships among geographical subgroups by genotyping 632 individuals from 31 different populations from their native USA regions (New Jersey, Michigan, and Georgia) and from invaded Korean regions using 12 microsatellite loci. Our population genetic analyses showed a clear separation between the two host-associated D. oxycoccana populations from blueberry and cranberry. Using data from only the blueberry-associated D. oxycoccana populations, we identified five genetically isolated subgroups. An analysis of the approximate Bayesian computation suggests that the invasive D. oxycoccana population from Korea appears to have been introduced from an unsampled source population rather than directly from its native range. Our findings will allow for an easier identification of the source of D. oxycoccana into newly invaded regions, as well as to determine their association with blueberry and cranberry, which based on our results can be considered as two distinct species.

Pollinator movement activity influences genetic diversity and differentiation of spatially isolated populations of clonal forest herbs

In agricultural landscapes, forest herbs live in small, spatially isolated forest patches. For their long-term survival, their populations depend on animals as genetic linkers that provide pollen- or seed-mediated gene flow among different forest patches. However, whether insect pollinators serve as genetic linkers among spatially isolated forest herb populations in agricultural landscapes remains to be shown. Here, we used population genetic methods to analyze: (A) the genetic diversity and genetic differentiation of populations of two common, slow-colonizing temperate forest herb species [Polygonatum multiflorum (L.) All. and Anemone nemorosa L.] in spatially isolated populations within three agricultural landscapes in Germany and Sweden and (B) the movement activity of their most relevant associated pollinator species, i.e., the bumblebee Bombus pascuorum (Scopoli, 1,763) and the hoverfly Melanostoma scalare (Fabricus, 1,794), respectively, which differ in their mobility. We tested whether the indicated pollinator movement activity affected the genetic diversity and genetic differentiation of the forest herb populations. Bumblebee movement indicators that solely indicated movement activity between the forest patches affected both genetic diversity and genetic differentiation of the associated forest herb P. multiflorum in a way that can be explained by pollen-mediated gene flow among the forest herb populations. In contrast, movement indicators reflecting the total movement activity at a forest patch (including within-forest patch movement activity) showed unexpected effects for both plant-pollinator pairs that might be explained by accelerated genetic drift due to enhanced sexual reproduction. Our integrated approach revealed that bumblebees serve as genetic linkers of associated forest herb populations, even if they are more than 2 km apart from each other. No such evidence was found for the forest associated hoverfly species which showed significant genetic differentiation among forest patches itself. Our approach also indicated that a higher within-forest patch movement activity of both pollinator species might enhance sexual recruitment and thus diminishes the temporal buffer that clonal growth provides against habitat fragmentation effects.

Species and population genomic differentiation in Pocillopora corals (Cnidaria, Hexacorallia)

Correctly delimiting species and populations is a prerequisite for studies of connectivity, adaptation and conservation. Genomic data are particularly useful to test species differentiation for organisms with few informative morphological characters or low discrimination of cytoplasmic markers, as in Scleractinians. Here we applied Restriction site Associated DNA sequencing (RAD-sequencing) to the study of species differentiation and genetic structure in populations of Pocillopora spp. from Oman and French Polynesia, with the objectives to test species hypotheses, and to study the genetic structure among sampling sites within species. We focused here on coral colonies morphologically similar to P. acuta (damicornis type β). We tested the impact of different filtering strategies on the stability of the results. The main genetic differentiation was observed between samples from Oman and French Polynesia. These samples corresponded to different previously defined primary species hypotheses (PSH), i.e., PSHs 12 and 13 in Oman, and PSH 5 in French Polynesia. In Oman, we did not observe any clear differentiation between the two putative species PSH 12 and 13, nor between sampling sites. In French Polynesia, where a single species hypothesis was studied, there was no differentiation between sites. Our analyses allowed the identification of clonal lineages in Oman and French Polynesia. The impact of clonality on genetic diversity is discussed in light of individual-based simulations.

Improving sustainable crop protection using population genetics concepts

Growing genetically resistant plants allows pathogen populations to be controlled and reduces the use of pesticides. However, pathogens can quickly overcome such resistance. In this context, how can we achieve sustainable crop protection? This crucial question has remained largely unanswered despite decades of intense debate and research effort. In this study, we used a bibliographic analysis to show that the research field of resistance durability has evolved into three subfields: (i) 'plant breeding' (generating new genetic material), (ii) 'molecular interactions' (exploring the molecular dialogue governing plant-pathogen interactions) and (iii) 'epidemiology and evolution' (explaining and forecasting of pathogen population dynamics resulting from selection pressure(s) exerted by resistant plants). We argue that this triple split of the field impedes integrated research progress and ultimately compromises the sustainable management of genetic resistance. After identifying a gap among the three subfields, we argue that the theoretical framework of population genetics could bridge this gap. Indeed, population genetics formally explains the evolution of all heritable traits, and allows genetic changes to be tracked along with variation in population dynamics. This provides an integrated view of pathogen adaptation, in particular via evolutionary-epidemiological feedbacks. In this Opinion Note, we detail examples illustrating how such a framework can better inform best practices for developing and managing genetically resistant cultivars.

The Role of Genetic Factors in the Differential Invasion Success of Two Spartina Species in China

Biological invasions have become one of the greatest threats to global biodiversity and ecosystem conservation. Most previous studies have revealed how successful invasive species adapt to new environments and climate change through phenotypic and genetic evolution. Some researchers suggested that understanding unsuccessful or less successful biological invasions might be important for understanding the relationships between invasion adaptability and climate factors. We compared the sexual reproduction ability, genetic diversity, and gene × environment interaction in two intentionally introduced alien species in China (Spartina anglica and Spartina alterniflora) based on restriction site-associated DNA (RAD) sequencing. After more than 50 years, the distribution of S. alterniflora has rapidly expanded, while S. anglica has experienced extreme dieback. A total of 212,939 single nucleotide polymorphisms (SNPs) for the two Spartina species were used for analysis. The multilocus genotype (MLG) analysis revealed that clonal reproduction was the prevalent mode of reproduction in both species, indicating that a change in the mode of reproduction was not the key factor enabling successful invasion by Spartina. All genetic diversity indicators (He, Ho, π) in S. alterniflora populations were at least two times higher than those in S. anglica populations, respectively (p < 0.001). Furthermore, the population genetic structure and stronger patterns of climate-associated loci provided support for rapid adaptive evolution in the populations of S. alterniflora in China. Altogether, our results highlight the importance of genetic diversity and local adaptation, which were driven by multiple source populations, in increasing the invasiveness of S. alterniflora.

Cloning and Selfing Affect Population Genetic Variation in Simulations of Outcrossing, Sexual Sea Stars

AbstractMany sea stars are well known for facultative or obligate asexual reproduction in both the adult and larval life-cycle stages. Some species and lineages are also capable of facultative or obligate hermaphroditic reproduction with self-fertilization. However, models of population genetic variation and empirical analyses of genetic data typically assume only sexual reproduction and outcrossing. A recent reanalysis of previously published empirical data (microsatellite genotypes) from two studies of one of the most well-known sea star species (the crown-of-thorns sea star; Acanthaster sp.) concluded that cloning and self-fertilization in that species are rare and contribute little to patterns of population genetic variation. Here we reconsider that conclusion by simulating the contribution of cloning and selfing to genetic variation in a series of models of sea star demography. Simulated variation in two simple models (analogous to previous analyses of empirical data) was consistent with high rates of cloning or selfing or both. More realistic scenarios that characterize population flux in sea stars of ecological significance, including outbreaks of crown-of-thorns sea stars that devastate coral reefs, invasions by Asterias amurensis, and epizootics of sea star wasting disease that kill Pisaster ochraceus, also showed significant but smaller effects of cloning and selfing on variation within subpopulations and differentiation between subpopulations. Future models or analyses of genetic variation in similar study systems might benefit from simulation modeling to characterize possible contributions of cloning or selfing to genetic variation in population samples or to understand the limits on inferring the effects of cloning or selfing in nature.

Detection of Cryptic Sex in Automictic Populations: Theoretical Expectations and a Case Study in Cataglyphis Desert Ants

Reproductive strategies are diverse and a whole continuum of mixed systems lies between strict sexuality and strict clonality (apomixis), including automixis, a parthenogenetic mode of reproduction involving a meiosis and increasing homozygosity over generations. These various systems impact the genetic structure of populations, which can therefore be used to infer reproductive strategies in natural populations. Here, we first develop a mathematical model, validated by simulations, to predict heterozygosity and inbreeding in mixed sexual-automictic populations. It highlights the predominant role of the rate of heterozygosity loss experienced during automixis (γ), which is locus dependent. When γ is low, mixed populations behave like purely sexual ones until sex becomes rare. In contrast, when γ is high, the erosion of genetic diversity is tightly correlated to the rate of sex, so that the individual inbreeding coefficient can inform on the ratio of sexual/asexual reproduction. In the second part of this study, we used our model to test the presence of cryptic sex in a hybridogenetic Cataglyphis ant where new queens are produced parthenogenetically, leaving males with an apparent null fitness while they are essential to colony development as sperm is required to produce workers. Occasional sexual production of queens could resolve this paradox by providing males some fertile progeny. To determine whether this occurs in natural populations, we simulated genotypic datasets in a population under various regimes of sexual vs. asexual reproduction for queen production and compared the distribution of inbreeding, expected heterozygosity and inter-individual relatedness coefficients with those observed in a natural population of Cataglyphis mauritanica using microsatellites. Our simulations show that the distribution of inter-individual relatedness coefficients was particularly informative to assess the relative rate of sexual/asexual reproduction, and our dataset was compatible with pure parthenogenesis but also with up to 2% sexual reproduction. Our approach, implemented in an R script, should be useful to assess reproductive strategies in other biological models.

Genomic signatures of clonality in the deep water kelp Laminaria rodriguezii

The development of population genomic approaches in non-model species allows for renewed studies of the impact of reproductive systems and genetic drift on population diversity. Here, we investigate the genomic signatures of partial clonality in the deep water kelp Laminaria rodriguezii, known to reproduce by both sexual and asexual means. We compared these results with the species Laminaria digitata, a closely related species that differs by different traits, in particular its reproductive mode (no clonal reproduction). We analysed genome-wide variation with dd-RAD sequencing using 4,077 SNPs in L. rodriguezii and 7,364 SNPs in L. digitata. As predicted for partially clonal populations, we show that the distribution of FIS within populations of L. rodriguezii is shifted toward negative values, with a high number of loci showing heterozygote excess. This finding is the opposite of what we observed within sexual populations of L. digitata, characterized by a generalized deficit in heterozygotes. Furthermore, we observed distinct distributions of FIS among populations of L. rodriguezii, which is congruent with the predictions of theoretical models for different levels of clonality and genetic drift. These findings highlight that the empirical distribution of FIS is a promising feature for the genomic study of asexuality in natural populations. Our results also show that the populations of L. rodriguezii analysed here are genetically differentiated and probably isolated. Our study provides a conceptual framework to investigate partial clonality on the basis of RAD-sequencing SNPs. These results could be obtained without any reference genome, and are therefore of interest for various non-model species.

The Combined Effect of Haplodiplonty and Partial Clonality on Genotypic and Genetic Diversity in a Finite Mutating Population

Partial clonality is known to affect the genetic composition and evolutionary trajectory of diplontic (single, free-living diploid stage) populations. However, many partially clonal eukaryotes exhibit life cycles in which somatic development occurs in both haploid and diploid individuals (haplodiplontic life cycles). Here, we studied how haplodiplontic life cycles and partial clonality structurally constrain, as immutable parameters, the reshuffling of genetic diversity and its dynamics in populations over generations. We assessed the distribution of common population genetic indices at different proportions of haploids, rates of clonality, mutation rates, and sampling efforts. Our results showed that haplodiplontic life cycles alone in finite populations affect effective population sizes and the ranges of distributions of population genetic indices. With nonoverlapping generations, haplodiplonty allowed the evolution of 2 temporal genetic pools that may diverge in sympatry due to genetic drift under full sexuality and clonality. Partial clonality in these life cycles acted as a homogenizing force between those 2 pools. Moreover, the combined effects of proportion of haploids, rate of clonality, and the relative strength of mutation versus genetic drift impacts the distributions of population genetics indices, rendering it difficult to transpose and use knowledge accumulated from diplontic or haplontic species. Finally, we conclude by providing recommendations for sampling and analyzing the population genetics of partially clonal haplodiplontic taxa.

Exploring the Genetic Consequences of Clonality in Haplodiplontic Taxa

Partially clonality is an incredibly common reproductive mode found across all the major eukaryotic lineages. Yet, population genetic theory is based on exclusive sexuality or exclusive asexuality, and partial clonality is often ignored. This is particularly true in haplodiplontic eukaryotes, including algae, ferns, mosses, and fungi, where somatic development occurs in both the haploid and diploid stages. Haplodiplontic life cycles are predicted to be correlated with asexuality, but tests of this prediction are rare. Moreover, there are unique consequences of having long-lived haploid and diploid stages in the same life cycle. For example, clonal processes uncouple the life cycle such that the repetition of the diploid stage via clonality leads to the loss of the haploid stage. Here, we surveyed the literature to find studies that had genotyped both haploid and diploid stages and recalculated population genetic summary metrics for seven red algae, one green alga, three brown algae, and three mosses. We compared these data to recent simulations that explicitly addressed the population genetic consequences of partial clonality in haplodiplontic life cycles. Not only was partial clonality found to act as a homogenizing force, but the combined effects of proportion of haploids, rate of clonality, and the relative strength of mutation versus genetic drift impacts the distributions of population genetic indices. We found remarkably similar patterns across commonly used population genetic metrics between our empirical and recent theoretical expectations. To facilitate future studies, we provide some recommendations for sampling and analyzing population genetic parameters for haplodiplontic taxa.

Complex evolution in Aphis gossypii group (Hemiptera: Aphididae), evidence of primary host shift and hybridization between sympatric species

Aphids provide a good model system to understand the ecological speciation concept, since the majority of the species are host-specific, and they spend their entire lifecycle on certain groups of host plants. Aphid species that apparently have wide host plant ranges have often turned out to be complexes of host-specialized biotypes. Here we investigated the various host-associated populations of the two recently diverged species, Aphis gossypii and A. rhamnicola, having multiple primary hosts, to understand the complex evolution with host-associated speciation. Using mitochondrial DNA marker and nine microsatellite loci, we reconstructed the haplotype network, and analyzed the genetic structure and relationships. Approximate Bayesian computation was also used to infer the ancestral primary host and host-associated divergence, which resulted in Rhamnus being the most ancestral host for A. gossypii and A. rhamnicola. As a result, Aphis gossypii and A. rhamnicola do not randomly use their primary and secondary host plants; rather, certain biotypes use only some secondary and specific primary hosts. Some biotypes are possibly in a diverging state through specialization to specific primary hosts. Our results also indicate that a new heteroecious race can commonly be derived from the heteroecious ancestor, showing strong evidence of ecological specialization through a primary host shift in both A. gossypii and A. rhamnicola. Interestingly, A. gossypii and A. rhamnicola shared COI haplotypes with each other, thus there is a possibility of introgression by hybridization between them by cross-sharing same primary hosts. Our results contribute to a new perspective in the study of aphid evolution by identifying complex evolutionary trends in the gossypii sensu lato complex.

The Contribution of Clonality to Population Genetic Structure in the Sea Anemone, Diadumene lineata

Ecological and evolutionary processes differ depending on how genetic diversity is organized in space. For clonal organisms, the organization of both genetic and genotypic diversity can influence the fitness effects of competition, the mating system, and reproductive mode, which are key drivers of life cycle evolution. Understanding how individual reproductive behavior contributes to population genetic structure is essential for disentangling these forces, particularly in species with complex and plastic life cycles. The widespread sea anemone, Diadumene lineata, exhibits temperature-dependent fission, which contributes to predictable variation in clonal rate along the Atlantic coast of the United States, part of its non-native range. Because warmer conditions lead to higher rates of clonality, we expected to find lower genotypic and genetic diversity in lower versus higher latitude populations. We developed primers for 11 microsatellite loci and genotyped 207 anemones collected from 8 sites ranging from Florida to Massachusetts. We found clonal influence at all sites, and as predicted, the largest clones were found at lower latitude sites. We also found genetic signatures of sex in the parts of the range where gametogenesis is most common. Evidence of sex outside the native range is novel for this species and provides insights into the dynamics of this successful invader. Our findings also illustrate challenges that partially clonal taxa pose for eco-evolutionary studies, such as difficulty sampling statistically robust numbers of genets and interpretating common population genetic metrics. For example, we found high among-locus variation in FIS, which makes the meaning of mean multilocus FIS unclear.

The discernible and hidden effects of clonality on the genotypic and genetic states of populations: Improving our estimation of clonal rates

Partial clonality is widespread across the tree of life, but most population genetic models are designed for exclusively clonal or sexual organisms. This gap hampers our understanding of the influence of clonality on evolutionary trajectories and the interpretation of population genetic data. We performed forward simulations of diploid populations at increasing rates of clonality (c), analysed their relationships with genotypic (clonal richness, R, and distribution of clonal sizes, Pareto β) and genetic (F_IS and linkage disequilibrium) indices, and tested predictions of c from population genetic data through supervised machine learning. Two complementary behaviours emerged from the probability distributions of genotypic and genetic indices with increasing c. While the impact of c on R and Pareto β was easily described by simple mathematical equations, its effects on genetic indices were noticeable only at the highest levels (c > 0.95). Consequently, genotypic indices allowed reliable estimates of c, while genetic descriptors led to poorer performances when c < 0.95. These results provide clear baseline expectations for genotypic and genetic diversity and dynamics under partial clonality. Worryingly, however, the use of realistic sample sizes to acquire empirical data systematically led to gross underestimates (often of one to two orders of magnitude) of c, suggesting that many interpretations hitherto proposed in the literature, mostly based on genotypic richness, should be reappraised. We propose future avenues to derive realistic confidence intervals for c and show that, although still approximate, a supervised learning method would greatly improve the estimation of c from population genetic data.

After a catastrophe, a little bit of sex is better than nothing: Genetic consequences of a major earthquake on asexual and sexual populations

Catastrophic events can have profound effects on the demography of a population and consequently on genetic diversity. The dynamics of postcatastrophic recovery and the role of sexual versus asexual reproduction in buffering the effects of massive perturbations remain poorly understood, in part because the opportunity to document genetic diversity before and after such events is rare. Six natural (purely sexual) and seven cultivated (mainly clonal due to farming practices) populations of the red alga Agarophyton chilense were surveyed along the Chilean coast before, in the days after and 2 years after the 8.8 magnitude earthquake in 2010. The genetic diversity of sexual populations appeared sensitive to this massive perturbation, notably through the loss of rare alleles immediately after the earthquake. By 2012, the levels of diversity returned to those observed before the catastrophe, probably due to migration. In contrast, enhanced rates of clonality in cultivated populations conferred a surprising ability to buffer the instantaneous loss of diversity. After the earthquake, farmers increased the already high rate of clonality to maintain the few surviving beds, but most of them collapsed rapidly. Contrasting fates between sexual and clonal populations suggest that betting on strict clonality to sustain production is risky, probably because this extreme strategy hampered adaptation to the brutal environmental perturbation induced by the catastrophe.

New insights into the population genetics of partially clonal organisms: When seagrass data meet theoretical expectations

Seagrass meadows are among the most important coastal ecosystems in terms of both spatial extent and ecosystem services, but they are also declining worldwide. Understanding the drivers of seagrass meadow dynamics is essential for designing sound management, conservation and restoration strategies. However, poor knowledge of the effect of clonality on the population genetics of natural populations severely limits our understanding of the dynamics and connectivity of meadows. Recent modelling approaches have described the expected distributions of genotypic and genetic descriptors under increasing clonal rates, which may help us better understand and interpret population genetics data obtained for partial asexuals. Here, in the light of these recent theoretical developments, we revisited population genetics data for 165 meadows of four seagrass species. Contrasting shoot lifespan and rhizome turnover led to the prediction that the influence of asexual reproduction would increase along a gradient from Zostera noltii to Zostera marina, Cymodocea nodosa and Posidonia oceanica, with increasing departure from Hardy-Weinberg equilibrium (F_is ), mostly towards heterozygote excess, and decreasing genotypic richness (R). This meta-analysis provides a nested validation of this hypothesis at both the species and meadow scales through a significant relationship between F_is and R within each species. By empirically demonstrating the theoretical expectations derived from recent modelling approaches, this work calls for the use of Hardy-Weinberg equilibrium (F_is ) rather than only the strongly sampling-sensitive R to assess the importance of clonal reproduction (c), at least when the impact of selfing on F_is can be neglected. The results also emphasize the need to revise our appraisal of the extent of clonality and its influence on the dynamics, connectivity and evolutionary trajectory of partial asexuals in general, including in seagrass meadows, to develop the most accurate management strategies.

Microsatellite Genotyping Corroborated Loss of Genetic Diversity in Clarias batrachus as a Result of Lack of Regulatory Reforms in Aquaculture

In India, over the past 50 years, aquaculture practices of species such as those used for Clarias batrachus were developed without adequate regulatory oversight. In these situations, it is important to consider the influence that genetic factors can have on such vulnerable aquaculture species. Population genetic structure can be evaluated through the use of neutral molecular markers, and this can aid in predicting the risk of the demise of populations and for framing management strategies to conserve remaining populations. The study presented here reports on the genetic status of C. batrachus populations through the analysis of data collected using 22 microsatellite markers from seven natural and one hatchery population. The mean values for observed heterozygosity across loci within populations ranged from 0.242 to 0.485. Measures of genetic differentiation were low overall, with mean values for F_ST of 0.270, F_IS of 0.113 and F_IT of 0.353. An AMOVA analysis revealed that percentages of variation among and within populations were 27.16 and 6.86, respectively, and Bayesian clustering analyses showed a population subdivision consisting of five clusters with admixture of haplotypes from other populations leading to genetic bottleneck. We also examined how hatchery management factors leading to excessive exchanges of fish between river systems through could impact the structure of the C. batrachus populations. Overall, this study shows how the systematic use of molecular markers can facilitate the development of management policies for these populations and for the development of a comprehensive set of rules for hatcheries and aquaculture practices, including avoidance of excessive homozygosity by avoiding repeated use of feral broodstock and their interrogation.

Neutral processes forming large clones during colonization of new areas

In species reproducing both sexually and asexually clones are often more common in recently established populations. Earlier studies have suggested that this pattern arises due to natural selection favouring generally or locally successful genotypes in new environments. Alternatively, as we show here, this pattern may result from neutral processes during species' range expansions. We model a dioecious species expanding into a new area in which all individuals are capable of both sexual and asexual reproduction, and all individuals have equal survival rates and dispersal distances. Even under conditions that favour sexual recruitment in the long run, colonization starts with an asexual wave. After colonization is completed, a sexual wave erodes clonal dominance. If individuals reproduce more than one season, and with only local dispersal, a few large clones typically dominate for thousands of reproductive seasons. Adding occasional long-distance dispersal, more dominant clones emerge, but they persist for a shorter period of time. The general mechanism involved is simple: edge effects at the expansion front favour asexual (uniparental) recruitment where potential mates are rare. Specifically, our model shows that neutral processes (with respect to genotype fitness) during the population expansion, such as random dispersal and demographic stochasticity, produce genotype patterns that differ from the patterns arising in a selection model. The comparison with empirical data from a post-glacially established seaweed species (Fucus radicans) shows that in this case, a neutral mechanism is strongly supported.

Causes and consequences of large clonal assemblies in a poplar hybrid zone

Asexual reproduction is a common and fundamental mode of reproduction in plants. Although persistence in adverse conditions underlies most known cases of clonal dominance, proximal genetic drivers remain unclear, in particular for populations dominated by a few large clones. In this study, we studied a clonal population of the riparian tree Populus alba in the Douro river basin (northwestern Iberian Peninsula) where it hybridizes with Populus tremula, a species that grows in highly contrasted ecological conditions. We used 73 nuclear microsatellites to test whether genomic background (species ancestry) is a relevant cause of clonal success, and to assess the evolutionary consequences of clonal dominance by a few genets. Additional genotyping-by-sequencing data were produced to estimate the age of the largest clones. We found that a few ancient (over a few thousand years old) and widespread genets dominate the population, both in terms of clone size and number of sexual offspring produced. Interestingly, large clones possessed two genomic regions introgressed from P. tremula, which may have favoured their spread under stressful environmental conditions. At the population level, the spread of large genets was accompanied by an overall ancient (>0.1 Myr) but soft decline of effective population size. Despite this decrease, and the high clonality and dominance of sexual reproduction by large clones, the Douro hybrid zone still displays considerable genetic diversity and low inbreeding. This suggests that even in extreme cases as in the Douro, asexual and sexual dominance of a few large, geographically extended individuals does not threaten population survival.

Divergence within and among Seaweed Siblings (Fucus vesiculosus and F. radicans) in the Baltic Sea

Closely related taxa provide significant case studies for understanding evolution of new species but may simultaneously challenge species identification and definition. In the Baltic Sea, two dominant and perennial brown algae share a very recent ancestry. Fucus vesiculosus invaded this recently formed postglacial sea 8000 years ago and shortly thereafter Fucus radicans diverged from this lineage as an endemic species. In the Baltic Sea both species reproduce sexually but also recruit fully fertile new individuals by asexual fragmentation. Earlier studies have shown local differences in morphology and genetics between the two taxa in the northern and western Bothnian Sea, and around the island of Saaremaa in Estonia, but geographic patterns seem in conflict with a single origin of F. radicans. To investigate the relationship between northern and Estonian distributions, we analysed the genetic variation using 9 microsatellite loci in populations from eastern Bothnian Sea, Archipelago Sea and the Gulf of Finland. These populations are located in between earlier studied populations. However, instead of bridging the disparate genetic gap between N-W Bothnian Sea and Estonia, as expected from a simple isolation-by-distance model, the new populations substantially increased overall genetic diversity and showed to be strongly divergent from the two earlier analysed regions, showing signs of additional distinct populations. Contrasting earlier findings of increased asexual recruitment in low salinity in the Bothnian Sea, we found high levels of sexual reproduction in some of the Gulf of Finland populations that inhabit extremely low salinity. The new data generated in this study supports the earlier conclusion of two reproductively isolated but very closely related species. However, the new results also add considerable genetic and morphological complexity within species. This makes species separation at geographic scales more demanding and suggests a need for more comprehensive approaches to further disentangle the intriguing relationship and history of the Baltic Sea fucoids.

Effects of complex life cycles on genetic diversity: cyclical parthenogenesis

Neutral patterns of population genetic diversity in species with complex life cycles are difficult to anticipate. Cyclical parthenogenesis (CP), in which organisms undergo several rounds of clonal reproduction followed by a sexual event, is one such life cycle. Many species, including crop pests (aphids), human parasites (trematodes) or models used in evolutionary science (Daphnia), are cyclical parthenogens. It is therefore crucial to understand the impact of such a life cycle on neutral genetic diversity. In this paper, we describe distributions of genetic diversity under conditions of CP with various clonal phase lengths. Using a Markov chain model of CP for a single locus and individual-based simulations for two loci, our analysis first demonstrates that strong departures from full sexuality are observed after only a few generations of clonality. The convergence towards predictions made under conditions of full clonality during the clonal phase depends on the balance between mutations and genetic drift. Second, the sexual event of CP usually resets the genetic diversity at a single locus towards predictions made under full sexuality. However, this single recombination event is insufficient to reshuffle gametic phases towards full-sexuality predictions. Finally, for similar levels of clonality, CP and acyclic partial clonality (wherein a fixed proportion of individuals are clonally produced within each generation) differentially affect the distribution of genetic diversity. Overall, this work provides solid predictions of neutral genetic diversity that may serve as a null model in detecting the action of common evolutionary or demographic processes in cyclical parthenogens (for example, selection or bottlenecks).