The exact distributions of F(IS) under partial asexuality in small finite populations with mutation

GenAPoPop 1.0: A user-friendly software to analyse genetic diversity and structure from partially clonal and selfed autopolyploid organisms

Autopolyploidy is quite common in most clades of eukaryotes. The emergence of sequence-based genotyping methods with individual and marker tags now enables confident allele dosage, overcoming the main obstacle to the democratization of the population genetic approaches when studying ecology and evolution of autopolyploid populations and species. Reproductive modes, including clonality, selfing and allogamy, have deep consequences on the ecology and evolution of population and species. Analysing genetic diversity and its dynamics over generations is one efficient way to infer the relative importance of clonality, selfing and allogamy in populations. GenAPoPop is a user-friendly solution to compute the specific corpus of population genetic indices, including indices about genotypic diversity, needed to analyse partially clonal, selfed and allogamous polysomic populations genotyped with confident allele dosage. It also easily provides the posterior probabilities of quantitative reproductive modes in autopolyploid populations genotyped at two-time steps and a graphical representation of the minimum spanning trees of the genetic distances between polyploid individuals, facilitating the interpretation of the genetic coancestry between individuals in hierarchically structured populations. GenAPoPop complements the previously existing solutions, including SPAGEDI and POLYGENE, to use genotypings to study the ecology and evolution of autopolyploid populations. It was specially developed with a simple graphical interface and workflow, and comes with a simulator to facilitate practical courses and teaching of population genetics for autopolyploid populations.

Parthenogenomics: Insights on mutation rates and nucleotide diversity in parthenogenetic Panagrolaimus nematodes

Asexual reproduction is assumed to lead to the accumulation of deleterious mutations, and reduced heterozygosity due to the absence of recombination. Panagrolaimid nematode species display different modes of reproduction. Sexual reproduction with distinct males and females, asexual reproduction through parthenogenesis in the genus Panagrolaimus, and hermaphroditism in Propanagrolaimus. Here, we compared genomic features of free-living nematodes in populations and species isolated from geographically distant regions to study diversity, and genome-wide differentiation under different modes of reproduction. We firstly estimated genome-wide spontaneous mutation rates in a triploid parthenogenetic Panagrolaimus, and a diploid hermaphroditic Propanagrolaimus via long-term mutation accumulation lines. Secondly, we calculated population genetic parameters including nucleotide diversity, and fixation index (F ST) between populations of asexually and sexually reproducing nematodes. Thirdly, we used phylogenetic network methods on sexually and asexually reproducing Panagrolaimus populations to understand evolutionary relationships between them. The estimated mutation rate was slightly lower for the asexual population, as expected for taxa with this reproductive mode. Natural polyploid asexual populations revealed higher nucleotide diversity. Despite their common ancestor, a gene network revealed a high level of genetic differentiation among asexual populations. The elevated heterozygosity found in the triploid parthenogens could be explained by the third genome copy. Given their tendentially lower mutation rates it can be hypothesized that this is part of the mechanism to evade Muller's ratchet. Our findings in parthenogenetic triploid nematode populations seem to challenge common expectations of evolution under asexuality.

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.

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.

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.

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.

Mixed-Mating Model of Reproduction Revealed in European Phytophthora cactorum by ddRADseq and Effector Gene Sequence Data

A population study of Phytophthora cactorum was performed using ddRADseq sequence variation analysis completed by the analysis of effector genes-RXLR6, RXLR7 and SCR113. The population structure was described by F-statistics, heterozygosity, nucleotide diversity, number of private alleles, number of polymorphic sites, kinship coefficient and structure analysis. The population of P. cactorum in Europe seems to be structured into host-associated groups. The isolates from woody hosts are structured into four groups described previously, while isolates from strawberry form another group. The groups are diverse in effector gene composition and the frequency of outbreeding. When populations from strawberry were analysed, both asexual reproduction and occasional outbreeding confirmed by gene flow among distinct populations were detected. Therefore, distinct P. cactorum populations differ in the level of heterozygosity. The data support the theory of the mixed-mating model for P. cactorum, comprising frequent asexual behaviour and inbreeding alternating with occasional outbreeding. Because P. cactorum is not indigenous to Europe, such variability is probably caused by multiple introductions of different lineages from the area of its original distribution, and the different histories of sexual recombination and host adaptation of particular populations.

Genetic Diversity and Connectivity in Plant Species Differing in Clonality and Dispersal Mechanisms in Wetland Island Habitats

In plants, long-distance dispersal is both attenuated and directed by specific movement vectors, including animals, wind, and/or water. Hence, movement vectors partly shape metapopulation genetic patterns that are, however, also influenced by other life-history traits such as clonal growth. We studied the relationship between area, isolation, plant-species richness, reproduction, and dispersal mechanisms with genetic diversity and divergence in 4 widespread wetland plant-species in a total of 20 island-like kettle-hole habitats surrounded by an intensive agricultural landscape. Our results showed that genetic parameters reflect the reproduction strategies with the highest genetic diversity being observed in the non-clonal, outcrossing Oenanthe aquatica compared to the clonal Lycopus europaeus, Typha latifolia, and Phragmites australis. Lycopus showed a positive relationship between genetic diversity and kettle-hole area, but a negative relationship with the number of neighboring kettle holes (less isolation). Genetic diversity increased with plant-species richness in the clonal species Phragmites and Lycopus; while it decreased in the non-clonal Oenanthe. Finally, genetic divergence and, therefore, connectivity differed between alternative dispersal strategies, where wind-dispersed Typha and Phragmites had a higher gene flow between the analyzed kettle holes compared with the insect-pollinated, hydrochorous Lycopus and Oenanthe. Our study provides information on genetic patterns related to reproduction and dispersal mechanisms of 4 common wetland species contributing to the understanding of the functioning of plant metacommunities occurring in kettle holes embedded in agricultural landscapes.

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.

The Origins of Coca: Museum Genomics Reveals Multiple Independent Domestications from Progenitor Erythroxylum gracilipes

Coca is the natural source of cocaine as well as a sacred and medicinal plant farmed by South American Amerindians and mestizos. The coca crop comprises four closely related varieties classified into two species (Amazonian and Huánuco varieties within Erythroxylum coca Lam., and Colombian and Trujillo varieties within Erythroxylum novogranatense (D. Morris) Hieron.) but our understanding of the domestication and evolutionary history of these taxa is nominal. In this study, we use genomic data from natural history collections to estimate the geographic origins and genetic diversity of this economically and culturally important crop in the context of its wild relatives. Our phylogeographic analyses clearly demonstrate the four varieties of coca comprise two or three exclusive groups nested within the diverse lineages of the widespread, wild species Erythroxylum gracilipes; establishing a new and robust hypothesis of domestication wherein coca originated two or three times from this wild progenitor. The Colombian and Trujillo coca varieties are descended from a single, ancient domestication event in northwestern South America. Huánuco coca was domesticated more recently, possibly in southeastern Peru. Amazonian coca either shares a common domesticated ancestor with Huánuco coca, or it was the product of a third and most recent independent domestication event in the western Amazon basin. This chronology of coca domestication reveals different Holocene peoples in South America were able to independently transform the same natural resource to serve their needs; in this case, a workaday stimulant. [Erythroxylum; Erythroxylaceae; Holocene; Museomics; Neotropics; phylogeography; plant domestication; target-sequence capture.]

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.

Phylogenomic analyses confirm a novel invasive North American Corbicula (Bivalvia: Cyrenidae) lineage

The genus Corbicula consists of estuarine or freshwater clams native to temperate/tropical regions of Asia, Africa, and Australia that collectively encompass both sexual species and clonal (androgenetic) lineages. The latter have become globally invasive in freshwater systems and they represent some of the most successful aquatic invasive lineages. Previous studies have documented four invasive clonal lineages, Forms A, B, C, and Rlc, with varying known distributions. Form A (R in Europe) occurs globally, Form B is found solely in North America, mainly the western United States, Form C (S in Europe) occurs both in European watersheds and in South America, and Rlc is known from Europe. A putative fifth invasive morph, Form D, was recently described in the New World from the Illinois River (Great Lakes watershed), where it occurs in sympatry with Forms A and B. An initial study showed Form D to be conchologically distinct: possessing rust-colored rays and white nacre with purple teeth. However, its genetic distinctiveness using standard molecular markers (mitochondrial cytochrome c oxidase subunit I and nuclear ribosomal 28S RNA) was ambiguous. To resolve this issue, we performed a phylogenomic analysis using 1,699-30,027 nuclear genomic loci collected via the next generation double digested restriction-site associated DNA sequencing method. Our results confirmed Form D to be a distinct invasive New World lineage with a population genomic profile consistent with clonality. A majority (7/9) of the phylogenomic analyses recovered the four New World invasive Corbicula lineages (Forms A, B, C, and D) as members of a clonal clade, sister to the non-clonal Lake Biwa (Japan) endemic, Corbicula sandai. The age of the clonal clade was estimated at 1.49 million years (my; ± 0.401-2.955 my) whereas the estimated ages of the four invasive lineage crown clades ranged from 0.27 to 0.44 my. We recovered very little evidence of nuclear genomic admixture among the four invasive lineages in our study populations. In contrast, 2/6 C. sandai individuals displayed partial nuclear genomic Structure assignments with multiple invasive clonal lineages. These results provide new insights into the origin and maintenance of clonality in this complex system.

A few north Appalachian populations are the source of European black locust

The role of evolution in biological invasion studies is often overlooked. In order to evaluate the evolutionary mechanisms behind invasiveness, it is crucial to identify the source populations of the introduction. Studies in population genetics were carried out on Robinia pseudoacacia L., a North American tree which is now one of the worst invasive tree species in Europe. We realized large-scale sampling in both the invasive and native ranges: 63 populations were sampled and 818 individuals were genotyped using 113 SNPs. We identified clonal genotypes in each population and analyzed between and within range population structure, and then, we compared genetic diversity between ranges, enlarging the number of SNPs to mitigate the ascertainment bias. First, we demonstrated that European black locust was introduced from just a limited number of populations located in the Appalachian Mountains, which is in agreement with the historical documents briefly reviewed in this study. Within America, population structure reflected the effects of long-term processes, whereas in Europe it was largely impacted by human activities. Second, we showed that there is a genetic bottleneck between the ranges with a decrease in allelic richness and total number of alleles in Europe. Lastly, we found more clonality within European populations. Black locust became invasive in Europe despite being introduced from a reduced part of its native distribution. Our results suggest that human activity, such as breeding programs in Europe and the seed trade throughout the introduced range, had a major role in promoting invasion; therefore, the introduction of the missing American genetic cluster to Europe should be avoided.

Intraspecific genetic structure, divergence and high rates of clonality in an amphi-Atlantic starfish

Intraspecific genetic diversity and divergence have a large influence on the adaption and evolutionary potential of species. The widely distributed starfish, Coscinasterias tenuispina, combines sexual reproduction with asexual reproduction via fission. Here we analyse the phylogeography of this starfish to reveal historical and contemporary processes driving its intraspecific genetic divergence. We further consider whether asexual reproduction is the most important method of propagation throughout the distribution range of this species. Our study included 326 individuals from 16 populations, covering most of the species' distribution range. A total of 12 nuclear microsatellite loci and sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene were analysed. COI and microsatellites were clustered in two isolated lineages: one found along the southwestern Atlantic and the other along the northeastern Atlantic and Mediterranean Sea. This suggests the existence of two different evolutionary units. Marine barriers along the European coast would be responsible for population clustering: the Almeria-Oran Front that limits the entrance of migrants from the Atlantic to the Mediterranean, and the Siculo-Tunisian strait that divides the two Mediterranean basins. The presence of identical genotypes was detected in all populations, although two monoclonal populations were found in two sites where annual mean temperatures and minimum values were the lowest. Our results based on microsatellite loci showed that intrapopulation genetic diversity was significantly affected by clonality whereas it had lower effect for the global phylogeography of the species, although still some impact on populations' genetic divergence could be observed between some populations.

Clonal expansion behind a marine diatom bloom

Genetic diversity is what selection acts on, thus shaping the adaptive potential of populations. We studied micro-evolutionary patterns of the key planktonic diatom Pseudo-nitzschia multistriata at a long-term sampling site over 2 consecutive years by genotyping isolates with 22 microsatellite markers. We show that both sex and vegetative growth interplay in shaping intraspecific diversity. We document a brief but massive demographic and clonal expansion driven by strains of the same mating type. The analysis of an extended data set (6 years) indicates that the genetic fingerprint of P. multistriata changed over time with a nonlinear pattern, with intermittent periods of weak and strong diversification related to the temporary predominance of clonal expansions over sexual recombination. These dynamics, rarely documented for phytoplankton, contribute to the understanding of bloom formation and of the mechanisms that drive microevolution in diatoms.

Population differentiation or species formation across the Indian and the Pacific Oceans? An example from the brooding marine hydrozoan Macrorhynchia phoenicea

Assessing population connectivity is necessary to construct effective marine protected areas. This connectivity depends, among other parameters, inherently on species dispersal capacities. Isolation by distance (IBD) is one of the main modes of differentiation in marine species, above all in species presenting low dispersal abilities. This study reports the genetic structuring in the tropical hydrozoan Macrorhynchia phoenicea α (sensu Postaire et al., 2016a), a brooding species, from 30 sampling sites in the Western Indian Ocean and the Tropical Southwestern Pacific, using 15 microsatellite loci. At the local scale, genet dispersal relied on asexual propagation at short distance, which was not found at larger scales. Considering one representative per clone, significant positive F_IS values (from −0.327*** to 0.411***) were found within almost all sites. Gene flow was extremely low at all spatial scales, among sites within islands (<10 km distance) and among islands (100 to >11,000 km distance), with significant pairwise F_ST values (from 0.035*** to 0.645***). A general pattern of IBD was found at the Indo‐Pacific scale, but also within ecoregions in the Western Indian Ocean province. Clustering and network analyses identified each island as a potential independent population, while analysis of molecular variance indicated that population genetic differentiation was significant at small (within island) and intermediate (among islands within province) spatial scales. As shown by this species, a brooding life cycle might be corollary of the high population differentiation found in some coastal marine species, thwarting regular dispersal at distances more than a few kilometers and probably leading to high cryptic diversity, each island housing independent evolutionary lineages.

Generalization of the QST framework in hierarchically structured populations: Impacts of inbreeding and dominance

Q_ST is a differentiation parameter based on the decomposition of the genetic variance of a trait. In the case of additive inheritance and absence of selection, it is analogous to the genic differentiation measured on individual loci, F_ST . Thus, Q_ST -F_ST comparison is used to infer selection: selective divergence when Q_ST  > F_ST , or convergence when Q_ST  < F_ST. The definition of Q-statistics was extended to two-level hierarchical population structures with Hardy-Weinberg equilibrium. Here, we generalize the Q-statistics framework to any hierarchical population structure. First, we developed the analytical definition of hierarchical Q-statistics for populations not at Hardy-Weinberg equilibrium. We show that the Q-statistics values obtained with the Hardy-Weinberg definition are lower than their corresponding F-statistics when F_IS  > 0 (higher when F_IS  < 0). Then, we used an island model simulation approach to investigate the impact of inbreeding and dominance on the Q_ST -F_ST framework in a hierarchical population structure. We show that, while differentiation at the lower hierarchical level (Q_SR ) is a monotonic function of migration, differentiation at the upper level (Q_RT ) is not. In the case of additive inheritance, we show that inbreeding inflates the variance of Q_RT , which can increase the frequency of Q_RT  > F_RT cases. We also show that dominance drastically reduces Q-statistics below F-statistics for any level of the hierarchy. Therefore, high values of Q-statistics are good indicators of selection, but low values are not in the case of dominance.

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).

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

BACKGROUND

Partially clonal organisms are very common in nature, yet the influence of partial asexuality on the temporal dynamics of genetic diversity remains poorly understood. Mathematical models accounting for clonality predict deviations only for extremely rare sex and only towards mean inbreeding coefficient [Formula: see text]. Yet in partially clonal species, both F IS  < 0 and F IS  > 0 are frequently observed also in populations where there is evidence for a significant amount of sexual reproduction. Here, we studied the joint effects of partial clonality, mutation and genetic drift with a state-and-time discrete Markov chain model to describe the dynamics of F IS over time under increasing rates of clonality.

RESULTS

Results of the mathematical model and simulations show that partial clonality slows down the asymptotic convergence to F IS  = 0. Thus, although clonality alone does not lead to departures from Hardy-Weinberg expectations once reached the final equilibrium state, both negative and positive F IS values can arise transiently even at intermediate rates of clonality. More importantly, such "transient" departures from Hardy Weinberg proportions may last long as clonality tunes up the temporal variation of F IS and reduces its rate of change over time, leading to a hyperbolic increase of the maximal time needed to reach the final mean [Formula: see text] value expected at equilibrium.

CONCLUSION

Our results argue for a dynamical interpretation of F IS in clonal populations. Negative values cannot be interpreted as unequivocal evidence for extremely scarce sex but also as intermediate rates of clonality in finite populations. Complementary observations (e.g. frequency distribution of multiloci genotypes, population history) or time series data may help to discriminate between different possible conclusions on the extent of clonality when mean [Formula: see text] values deviating from zero and/or a large variation of F IS over loci are observed.

Interpretation and approximation tools for big, dense Markov chain transition matrices in population genetics

Background

Markov chains are a common framework for individual-based state and time discrete models in evolution. Though they played an important role in the development of basic population genetic theory, the analysis of more complex evolutionary scenarios typically involves approximation with other types of models. As the number of states increases, the big, dense transition matrices involved become increasingly unwieldy. However, advances in computational technology continue to reduce the challenges of “big data”, thus giving new potential to state-rich Markov chains in theoretical population genetics.

Results

Using a population genetic model based on genotype frequencies as an example, we propose a set of methods to assist in the computation and interpretation of big, dense Markov chain transition matrices. With the help of network analysis, we demonstrate how they can be transformed into clear and easily interpretable graphs, providing a new perspective even on the classic case of a randomly mating, finite population with mutation. Moreover, we describe an algorithm to save computer memory by substituting the original matrix with a sparse approximate while preserving its mathematically important properties, including a closely corresponding dominant (normalized) eigenvector. A global sensitivity analysis of the approximation results in our example shows that size reduction of more than 90 % is possible without significantly affecting the basic model results. Sample implementations of our methods are collected in the Python module mamoth.

Conclusion

Our methods help to make stochastic population genetic models involving big, dense transition matrices computationally feasible. Our visualization techniques provide new ways to explore such models and concisely present the results. Thus, our methods will contribute to establish state-rich Markov chains as a valuable supplement to the diversity of population genetic models currently employed, providing interesting new details about evolution e.g. under non-standard reproductive systems such as partial clonality.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-015-0061-5) contains supplementary material, which is available to authorized users.