How clonal are clones? A quest for loss of heterozygosity during asexual reproduction in Daphnia magna

Due to the lack of recombination, asexual organisms are predicted to accumulate mutations and show high levels of within‐individual allelic divergence (heterozygosity); however, empirical evidence for this prediction is largely missing. Instead, evidence of genome homogenization during asexual reproduction is accumulating. Ameiotic crossover recombination is a mechanism that could lead to long genomic stretches of loss of heterozygosity (LOH) and unmasking of mutations that have little or no effect in heterozygous state. Therefore, LOH might be an important force for inducing variation among asexual offspring and may contribute to the limited longevity of asexual lineages. To investigate the genetic consequences of asexuality, here we used high‐throughput sequencing of Daphnia magna for assessing the rate of LOH over a single generation of asexual reproduction. Comparing parthenogenetic daughters with their mothers at several thousand genetic markers generated by restriction site‐associated DNA (RAD) sequencing resulted in high LOH rate estimation that largely overlapped with our estimates for the error rate. To distinguish these two, we Sanger re‐sequenced the top 17 candidate RAD‐loci for LOH, and all of them proved to be false positives. Hence, even though we cannot exclude the possibility that short stretches of LOH occur in genomic regions not covered by our markers, we conclude that LOH does not occur frequently during asexual reproduction in D. magna and ameiotic crossovers are very rare or absent. This finding suggests that clonal lineages of D. magna will remain genetically homogeneous at least over time periods typically relevant for experimental work.     

Minimum sample sizes for invasion genomics: Empirical investigation in an invasive whitefly

Analysis of population genetics provides insights into the evolutionary processes, among which the sample size choice is per se a crucial issue in the analysis. Genome‐wide high‐throughput techniques based on RADseq have been increasingly used in studies on the population genomics of invasive species. However, there is little information available regarding optimal sample sizes for analyzing population genomics of invasive species. In this study, we first use type IIB endonucleases restriction site‐associated DNA (2b‐RAD) to mine thousands of single nucleotide polymorphisms (SNPs) for native and introduced populations in Q1 clade (SPB and 17JN) and Q2 clade (ISQ and UAS0601) of the whitefly, Bemisia tabaci (Gennadius) MED (also known as B. tabaci biotype Q). Then, we used resampling techniques to create simulated populations with a random subset of individuals and 3,000 SNPs to determine how many individuals should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We calculated the intrapopulation genetic diversity parameters (unbiased expected heterozygosity, observed heterozygosity, and the number of effect alleles) and pairwise genetic differentiation F_ST; finally, an ad hoc statistic, ΔK, was used to determine the optimal value. Our results showed that a sample size greater than four individuals (n ≥ 4) has little impact on estimates of genetic diversity within whitefly populations; moreover, precise estimate of F_ST can be easily achieved at a very small simple size (n = 3 or 4). Our results will provide in‐depth understanding of the optimization of sampling schemes in population genomics of invasive species.     

The Exact Distributions of F IS under Partial Asexuality in Small Finite Populations with Mutation

Reproductive systems like partial asexuality participate to shape the evolution of genetic diversity within populations, which is often quantified by the inbreeding coefficient F _IS. Understanding how those mating systems impact the possible distributions of F _IS values in theoretical populations helps to unravel forces shaping the evolution of real populations. We proposed a population genetics model based on genotypic states in a finite population with mutation. For populations with less than 400 individuals, we assessed the impact of the rates of asexuality on the full exact distributions of F _IS, the probabilities of positive and negative F _IS, the probabilities of fixation and the probabilities to observe changes in the sign of F _IS over one generation. After an infinite number of generations, we distinguished three main patterns of effects of the rates of asexuality on genetic diversity that also varied according to the interactions of mutation and genetic drift. Even rare asexual events in mainly sexual populations impacted the balance between negative and positive F _IS and the occurrence of extreme values. It also drastically modified the probability to change the sign of F _IS value at one locus over one generation. When mutation prevailed over genetic drift, increasing rates of asexuality continuously increased the variance of F _IS that reached its highest value in fully asexual populations. In consequence, even ancient asexual populations showed the entire F _IS spectrum, including strong positive F _IS. The prevalence of heterozygous loci only occurred in full asexual populations when genetic drift dominated.     

RIBOSOMAL RNA GENE DIVERSITY,EFFECTIVE POPULATION SIZE,AND EVOLUTIONARY LONGEVITY IN ASEXUAL GLOMEROMYCOTA

Arbuscular mycorrhizal fungi (phylum Glomeromycota) are among the oldest and most successful symbionts of land plants. With no evidence of sexual reproduction, their evolutionary success is inconsistent with the prediction that asexual taxa are vulnerable to extinction due to accumulation of deleterious mutations. To explore why Glomeromycota defy this prediction, we studied ribosomal RNA (rRNA) gene evolution in the Claroideoglomus lineage and estimated effective population size, N_e, in C. etunicatum. We found that rRNA genes of these fungi exhibit unusual and complex patterns of molecular evolution. In C. etunicatum, these patterns can be collectively explained by an unexpectedly large N_e combined with imperfect genome‐wide and population‐level rRNA gene repeat homogenization. The mutations accumulated in rRNA gene sequences indicate that natural selection is effective at purging deleterious mutations in the Claroideoglomus lineage, which is also consistent with the large N_e of C. etunicatum. We propose that in the near absence of recombination, asexual reproduction involving massively multinucleate spores typical for Glomeromycota is responsible for the improved efficacy of selection relative to drift. We postulate that large effective population sizes contribute to the evolutionary longevity of Glomeromycota.     

Mating system and environmental variation drive patterns of adaptation in Boechera spatifolia (Brassicaceae)

Determining the relative contribution of population genetic processes to the distribution of natural variation is a major goal of evolutionary biology. Here, we take advantage of variation in mating system to test the hypothesis that local adaptation is constrained by asexual reproduction. We explored patterns of variation in ecological traits and genome‐wide molecular markers in Boechera spatifolia (Brassicaceae), a species that contains both apomictic (asexual) and sexual individuals. Using a combination of quantitative genetics, neutral genetic (SSR) and genome‐wide single nucleotide polymorphism, we assessed the hypothesis that asexual lineages should have reduced signatures of adaptation relative to sexual conspecifics. All three measures (traits, SSRs, SNPs) demonstrated that apomicts are genetically distinct from sexuals, regardless of population location. Additionally, phylogenetic clustering revealed that the apomictic group shared a single common ancestor. Across the landscape, sexual genome‐wide SNP variation was strongly associated with latitude (r² > 0.9), indicating that sexual populations have differentiated across an environmental gradient. Furthermore, flowering time and growth rate, as assessed in a common garden, strongly covary with the elevation and latitude of the source population. Despite a wide geographic distribution that largely overlaps with sexual populations, there was little evidence for differentiation in molecular markers or quantitative characters among apomictic populations. Combined, these data indicated that, in contrast to asexual populations, sexual populations show evidence of local adaptation.     

Single‐nucleotide polymorphism discovery in Leptographium longiclavatum,a mountain pine beetle‐associated symbiotic fungus,using whole‐genome resequencing

Single‐nucleotide polymorphisms (SNPs) are rapidly becoming the standard markers in population genomics studies; however, their use in nonmodel organisms is limited due to the lack of cost‐effective approaches to uncover genome‐wide variation, and the large number of individuals needed in the screening process to reduce ascertainment bias. To discover SNPs for population genomics studies in the fungal symbionts of the mountain pine beetle (MPB), we developed a road map to discover SNPs and to produce a genotyping platform. We undertook a whole‐genome sequencing approach of Leptographium longiclavatum in combination with available genomics resources of another MPB symbiont, Grosmannia clavigera. We sequenced 71 individuals pooled into four groups using the Illumina sequencing technology. We generated between 27 and 30 million reads of 75 bp that resulted in a total of 1, 181 contigs longer than 2 kb and an assembled genome size of 28.9 Mb (N₅₀ = 48 kb, average depth = 125x). A total of 9052 proteins were annotated, and between 9531 and 17 266 SNPs were identified in the four pools. A subset of 206 genes (containing 574 SNPs, 11% false positives) was used to develop a genotyping platform for this species. Using this roadmap, we developed a genotyping assay with a total of 147 SNPs located in 121 genes using the Illumina^® Sequenom iPLEX Gold. Our preliminary genotyping (success rate = 85%) of 304 individuals from 36 populations supports the utility of this approach for population genomics studies in other MPB fungal symbionts and other fungal nonmodel species.     

EVOLUTION OF ASEXUALITY VIA DIFFERENT MECHANISMS IN GRASS THRIPS (THYSANOPTERA: Aptinothrips)

Asexual lineages can derive from sexual ancestors via different mechanisms and at variable rates, which affects the diversity of the asexual population and thereby its ecological success. We investigated the variation and evolution of reproductive systems in Aptinothrips, a genus of grass thrips comprising four species. Extensive population surveys and breeding experiments indicated sexual reproduction in A. elegans, asexuality in A. stylifer and A. karnyi, and both sexual and asexual lineages in A. rufus. Asexuality in A. stylifer and A. rufus coincides with a worldwide distribution, with sexual A. rufus lineages confined to a limited area. Inference of molecular phylogenies and antibiotic treatment revealed different causes of asexuality in different species. Asexuality in A. stylifer and A. karnyi has most likely genetic causes, while it is induced by endosymbionts in A. rufus. Endosymbiont‐community characterization revealed presence of Wolbachia, and lack of other bacteria known to manipulate host reproduction. However, only 69% asexual A. rufus females are Wolbachia‐infected, indicating that either an undescribed endosymbiont causes asexuality in this species or that Wolbachia was lost in several lineages that remained asexual. These results open new perspectives for studies on the maintenance of mixed sexual and asexual reproduction in natural populations.     

Origins of asexuality in <Emphasis Type="Italic">Bryobia</Emphasis>mites (Acari: Tetranychidae)

Background  

Obligate asexual reproduction is rare in the animal kingdom. Generally, asexuals are considered evolutionary dead ends that are unable to radiate. The phytophagous mite genus Bryobia contains a large number of asexual species. In this study, we investigate the origin and evolution of asexuality using samples from 111 populations in Europe, South Africa and the United States, belonging to eleven Bryobia species. We also examine intraspecific clonal diversity for one species, B. kissophila, by genotyping individuals from 61 different populations. Knowledge on the origin of asexuality and on clonal diversity can contribute to our understanding of the paradox of sex.     

Genetic variation and asexual reproduction in the facultatively parthenogenetic cockroach Nauphoeta cinerea: implications for the evolution of sex

Asexual reproduction could offer up to a two‐fold fitness advantage over sexual reproduction, yet higher organisms usually reproduce sexually. Even in facultatively parthenogenetic species, where both sexual and asexual reproduction is sometimes possible, asexual reproduction is rare. Thus, the debate over the evolution of sex has focused on ecological and mutation‐elimination advantages of sex. An alternative explanation for the predominance of sex is that it is difficult for an organism to accomplish asexual reproduction once sexual reproduction has evolved. Difficulty in returning to asexuality could reflect developmental or genetic constraints. Here, we investigate the role of genetic factors in limiting asexual reproduction in Nauphoeta cinerea, an African cockroach with facultative parthenogenesis that nearly always reproduces sexually. We show that when N. cinerea females do reproduce asexually, offspring are genetically identical to their mothers. However, asexual reproduction is limited to a nonrandom subset of the genotypes in the population. Only females that have a high level of heterozygosity are capable of parthenogenetic reproduction and there is a strong familial influence on the ability to reproduce parthenogenetically. Although the mechanism by which genetic variation facilitates asexual reproduction is unknown, we suggest that heterosis may facilitate the switch from producing haploid meiotic eggs to diploid, essentially mitotic, eggs.     

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide‐toxic habitats of a neotropical fish (Poecilia mexicana)

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H₂S) and compared two population pairs of sulphide‐adapted and ancestral fish by sequencing population pools of >200 individuals (Pool‐Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection‐mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide‐adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of J_corr = 0.0032 at the level of SNPs, divergent genome regions (J_corr = 0.0061) and genes therein (J_corr = 0.0091). At the level of metabolic pathways, the overlap was J_corr = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.     

A chromosome‐scale genome assembly of Antheraea pernyi (Saturniidae,Lepidoptera)

Antheraea pernyi is a semi‐domesticated lepidopteran insect species valuable to the silk industry, human health, and ecological tourism. Owing to its economic influence and developmental properties, it serves as an ideal model for investigating divergence of the Bombycoidea super family. However, studies on the karyotype evolution and functional genomics of A. pernyi are limited by scarce genomic resource. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first high‐quality A. pernyi genome from a single male individual. The genome is 720.67 Mb long with 49 chromosomes and a 13.77‐Mb scaffold N50. Approximately 441.75 Mb, accounting for 60.74% of the genome, was identified as repeats. The genome comprises 21,431 protein‐coding genes, 85.22% of which were functionally annotated. Comparative genomics analysis suggested that A. pernyi diverged from its common ancestor with A. yamamai ~30.3 million years ago, and that chromosome fission contributed to the increased chromosome number. The genome assembled in this work will not only facilitate future research on A. pernyi and related species but also help to progress comparative genomics analyses in Lepidoptera.     

LARGE POPULATION SIZE PREDICTS THE DISTRIBUTION OF ASEXUALITY IN SCALE INSECTS

Understanding why some organisms reproduce by sexual reproduction while others can reproduce asexually remains an important unsolved problem in evolutionary biology. Simple demography suggests that asexuals should outcompete sexually reproducing organisms, because of their higher intrinsic rate of increase. However, the majority of multicellular organisms have sexual reproduction. The widely accepted explanation for this apparent contradiction is that asexual lineages have a higher extinction rate. A number of models have indicated that population size might play a crucial role in the evolution of asexuality. The strength of processes that lead to extinction of asexual species is reduced when population sizes get very large, so that the long‐term advantage of sexual over asexual reproduction may become negligible. Here, we use a comparative approach using scale insects (Coccoidea, Hemiptera) to show that asexuality is indeed more common in species with larger population density and geographic distribution and we also show that asexual species tend to be more polyphagous. We discuss the implication of our findings for previously observed patterns of asexuality in agricultural pests.     

Survival in the Tropics despite isolation,inbreeding and asexual reproduction: insights from the genome of the world’s southernmost poplar (Populus ilicifolia)

Species are becoming extinct at unprecedented rates as a consequence of human activity. Hence it is important to understand the evolutionary dynamics of species with already small population sizes. Populus ilicifolia is a vulnerable poplar species that is isolated from other poplar species and is uniquely adapted to the Tropics. It has a very limited size, reproduces partly clonally and is therefore an excellent case study for conservation genomics. We present here the first annotated draft genome of P. ilicifolia, characterize genome‐wide patterns of polymorphisms and compare those to other poplar species with larger natural ranges. P. ilicifolia experienced a more prolonged and severe decline of effective population size (N_e) and signs of genetic erosion than any other poplar species with which it was compared. At present, the species has the lowest genome‐wide genetic diversity, the highest abundance of long runs of homozygosity, high inbreeding levels as well as a high overall accumulation of deleterious variants. However, more effective purging of severely deleterious variants and adaptation to the Tropics may have contributed to its survival. Hence, in spite of its limited genetic variation, it is certainly worth pursuing the conservation efforts of this unique species.     

Automixis in Artemia: solving a century‐old controversy

Parthenogenesis (reproduction through unfertilized eggs) encompasses a variety of reproduction modes with (automixis) or without (apomixis) meiosis. Different modes of automixis have very different genetic and evolutionary consequences but can be particularly difficult to tease apart. In this study, we propose a new method to discriminate different types of automixis from population‐level genetic data. We apply this method to diploid Artemia parthenogenetica, a crustacean whose reproductive mode remains controversial despite a century of intensive cytogenetic observations. We focus on A. parthenogenetica from two western Mediterranean populations. We show that they are diploid and that markers remain heterozygous in cultures maintained up to ~36 generations in the laboratory. Moreover, parallel patterns of population‐wide heterozygosity levels between the two natural populations strongly support the conclusion that diploid A. parthenogenetica reproduce by automictic parthenogenesis with central fusion and low, but nonzero recombination. This settles a century‐old controversy on Artemia, and, more generally, suggests that many automictic organisms harbour steep within‐chromosome gradients of heterozygosity due to a transition from clonal transmission in centromere‐proximal regions to a form of inbreeding similar to self‐fertilization in centromere‐distal regions. Such systems therefore offer a new avenue for contrasting the genomic consequences of asexuality and inbreeding.     

Hybrid asexuality as a primary postzygotic barrier between nascent species: On the interconnection between asexuality,hybridization and speciation

Although sexual reproduction is ubiquitous throughout nature, the molecular machinery behind it has been repeatedly disrupted during evolution, leading to the emergence of asexual lineages in all eukaryotic phyla. Despite intensive research, little is known about what causes the switch from sexual reproduction to asexuality. Interspecific hybridization is one of the candidate explanations, but the reasons for the apparent association between hybridization and asexuality remain unclear. In this study, we combined cross‐breeding experiments with population genetic and phylogenomic approaches to reveal the history of speciation and asexuality evolution in European spined loaches (Cobitis). Contemporary species readily hybridize in hybrid zones, but produce infertile males and fertile but clonally reproducing females that cannot mediate introgressions. However, our analysis of exome data indicates that intensive gene flow between species has occurred in the past. Crossings among species with various genetic distances showed that, while distantly related species produced asexual females and sterile males, closely related species produce sexually reproducing hybrids of both sexes. Our results suggest that hybridization leads to sexual hybrids at the initial stages of speciation, but as the species diverge further, the gradual accumulation of reproductive incompatibilities between species could distort their gametogenesis towards asexuality. Interestingly, comparative analysis of published data revealed that hybrid asexuality generally evolves at lower genetic divergences than hybrid sterility or inviability. Given that hybrid asexuality effectively restricts gene flow, it may establish a primary reproductive barrier earlier during diversification than other “classical” forms of postzygotic incompatibilities. Hybrid asexuality may thus indirectly contribute to the speciation process.     

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

The presence and extent of mitonuclear discordance in coexisting sexual and asexual lineages provides insight into 1) how and when asexual lineages emerged, and 2) the spatial and temporal scales at which the ecological and evolutionary processes influencing the evolution of sexual and asexual reproduction occur. Here, we used nuclear single‐nucleotide polymorphism (SNP) markers and a mitochondrial gene to characterize phylogeographic structure and the extent of mitonuclear discordance in Potamopyrgus antipodarum. This New Zealand freshwater snail is often used to study the evolution and maintenance of sex because obligately sexual and obligately asexual individuals often coexist. While our data indicate that sexual and asexual P. antipodarum sampled from the same lake population are often genetically similar, suggesting recent origin of these asexuals from sympatric sexual P. antipodarum, we also found significantly more population structure in sexuals vs. asexuals. This latter result suggests that some asexual lineages originated in other lakes and/or in the relatively distant past. When comparing mitochondrial and nuclear population genetic structure, we discovered that one mitochondrial haplotype (‘1A’) was rare in sexuals, but common and widespread in asexuals. Haplotype 1A frequency and nuclear genetic diversity were not associated, suggesting that the commonness of this haplotype cannot be attributed entirely to genetic drift and pointing instead to a role for selection.     

Fine‐scale spatiotemporal patterns of genetic variation reflect budding dispersal coupled with strong natal philopatry in a cooperatively breeding mammal

The relatedness structure of animal populations is thought to be a critically important factor underlying the evolution of mating systems and social behaviours. While previous work has shown that population structure is shaped by many biological processes, few studies have investigated how these factors vary over time. Consequently, we explored the fine‐scale spatiotemporal genetic structure of an intensively studied population of cooperatively breeding banded mongooses (Mungos mungo) over a 10‐year period. Overall population structure was strong (average F_ST = 0.129) but groups with spatially overlapping territories were not more genetically similar to one another than noncontiguous groups. Instead, genetic differentiation was associated with historical group‐fission (budding) events, with new groups diverging from their parent groups over time. Within groups, relatedness was high within but not between the sexes, although the latter increased over time since group formation due to group founders being replaced by philopatric young. This trend was not mirrored by a decrease in average offspring heterozygosity over time, suggesting that close inbreeding may often be avoided, even when immigration into established groups is virtually absent and opportunities for extra‐group matings are rare. Fine‐scale spatiotemporal population structure could have important implications in social species, where relatedness between interacting individuals is a vital component in the evolution of patterns of inbreeding avoidance, reproductive skew and kin‐selected helping and harming.     

Strikingly high levels of heterozygosity despite 20 years of inbreeding in a clonal honey bee

Inbreeding (the mating between closely related individuals) often has detrimental effects that are associated with loss of heterozygosity at overdominant loci, and the expression of deleterious recessive alleles. However, determining which loci are detrimental when homozygous, and the extent of their phenotypic effects, remains poorly understood. Here, we utilize a unique inbred population of clonal (thelytokous) honey bees, Apis mellifera capensis, to determine which loci reduce individual fitness when homozygous. This asexual population arose from a single worker ancestor approximately 20 years ago and has persisted for at least 100 generations. Thelytokous parthenogenesis results in a 1/3 of loss of heterozygosity with each generation. Yet, this population retains heterozygosity throughout its genome due to selection against homozygotes. Deep sequencing of one bee from each of the three known sub‐lineages of the population revealed that 3,766 of 10,884 genes (34%) have retained heterozygosity across all sub‐lineages, suggesting that these genes have heterozygote advantage. The maintenance of heterozygosity in the same genes and genomic regions in all three sub‐lineages suggests that nearly every chromosome carries genes that show sufficient heterozygote advantage to be selectively detrimental when homozygous.