Single-step generation of homozygous knockout/knock-in individuals in an extremotolerant parthenogenetic tardigrade using DIPA-CRISPR.

Tardigrades are small aquatic invertebrates known for their remarkable tolerance to diverse extreme stresses. To elucidate the in vivo mechanisms underlying this extraordinary resilience, methods for genetically manipulating tardigrades have long been desired. Despite our prior success in somatic cell gene editing by microinjecting Cas9 ribonucleoproteins (RNPs) into the body cavity of tardigrades, the generation of gene-edited individuals remained elusive. In this study, employing an extremotolerant parthenogenetic tardigrade species, Ramazzottius varieornatus, we established conditions that led to the generation of gene-edited tardigrade individuals. Drawing inspiration from the direct parental CRISPR (DIPA-CRISPR) technique employed in several insects, we simply injected a concentrated Cas9 RNP solution into the body cavity of parental females shortly before their initial oviposition. This approach yielded gene-edited G0 progeny. Notably, only a single allele was predominantly detected at the target locus for each G0 individual, indicative of homozygous mutations. By co-injecting single-stranded oligodeoxynucleotides (ssODNs) with Cas9 RNPs, we achieved the generation of homozygously knocked-in G0 progeny, and these edited alleles were inherited by G1/G2 progeny. This is the first example of heritable gene editing in the entire phylum of Tardigrada. This establishment of a straightforward method for generating homozygous knockout/knock-in individuals not only facilitates in vivo analyses of the molecular mechanisms underpinning extreme tolerance, but also opens up avenues for exploring various topics, including Evo-Devo, in tardigrades.

Post-meiotic mechanism of facultative parthenogenesis in gonochoristic whiptail lizard species.

Facultative parthenogenesis (FP) has historically been regarded as rare in vertebrates, but in recent years incidences have been reported in a growing list of fish, reptile, and bird species. Despite the increasing interest in the phenomenon, the underlying mechanism and evolutionary implications have remained unclear. A common finding across many incidences of FP is either a high degree of homozygosity at microsatellite loci or low levels of heterozygosity detected in next-generation sequencing data. This has led to the proposal that second polar body fusion following the meiotic divisions restores diploidy and thereby mimics fertilization. Here, we show that FP occurring in the gonochoristic Aspidoscelis species A. marmoratus and A. arizonae results in genome-wide homozygosity, an observation inconsistent with polar body fusion as the underlying mechanism of restoration. Instead, a high-quality reference genome for A. marmoratus and analysis of whole-genome sequencing from multiple FP and control animals reveals that a post-meiotic mechanism gives rise to homozygous animals from haploid, unfertilized oocytes. Contrary to the widely held belief that females need to be isolated from males to undergo FP, females housed with conspecific and heterospecific males produced unfertilized eggs that underwent spontaneous development. In addition, offspring arising from both fertilized eggs and parthenogenetic development were observed to arise from a single clutch. Strikingly, our data support a mechanism for facultative parthenogenesis that removes all heterozygosity in a single generation. Complete homozygosity exposes the genetic load and explains the high rate of congenital malformations and embryonic mortality associated with FP in many species. Conversely, for animals that develop normally, FP could potentially exert strong purifying selection as all lethal recessive alleles are purged in a single generation.

Wolbachia symbionts control sex in a parasitoid wasp using a horizontally acquired gene.

Host reproduction can be manipulated by bacterial symbionts in various ways. Parthenogenesis induction is the most effective type of reproduction manipulation by symbionts for their transmission. Insect sex is determined by regulation of doublesex (dsx) splicing through transformer2 (tra2) and transformer (tra) interaction. Although parthenogenesis induction by symbionts has been studied since the 1970s, its underlying molecular mechanism is unknown. Here we identify a Wolbachia parthenogenesis-induction feminization factor gene (piff) that targets sex-determining genes and causes female-producing parthenogenesis in the haplodiploid parasitoid Encarsia formosa. We found that Wolbachia elimination repressed expression of female-specific dsx and enhanced expression of male-specific dsx, which led to the production of wasp haploid male offspring. Furthermore, we found that E. formosa tra is truncated and non-functional, and Wolbachia has a functional tra homolog, termed piff, with an insect origin. Wolbachia PIFF can colocalize and interact with wasp TRA2. Moreover, Wolbachia piff has coordinated expression with tra2 and dsx of E. formosa. Our results demonstrate the bacterial symbiont Wolbachia has acquired an insect gene to manipulate the host sex determination cascade and induce parthenogenesis in wasps. This study reveals insect-to-bacteria horizontal gene transfer drives the evolution of animal sex determination systems, elucidating a striking mechanism of insect-microbe symbiosis.

Identification of gene isoforms and their switching events between male and female embryos of the parthenogenetic crustacean Daphnia magna.

The cladoceran crustacean Daphnia exhibits phenotypic plasticity, a phenomenon that leads to diverse phenotypes from one genome. Alternative usage of gene isoforms has been considered a key gene regulation mechanism for controlling different phenotypes. However, to understand the phenotypic plasticity of Daphnia, gene isoforms have not been comprehensively analyzed. Here we identified 25,654 transcripts derived from the 9710 genes expressed during environmental sex determination of Daphnia magna using the long-read RNA-Seq with PacBio Iso-Seq. We found that 14,924 transcripts were previously unidentified and 5713 genes produced two or more isoforms. By a combination of Illumina short-read RNA-Seq, we detected 824 genes that implemented switching of the highest expressed isoform between females and males. Among the 824 genes, we found isoform switching of an ortholog of CREB-regulated transcription coactivator, a major regulator of carbohydrate metabolism in animals, and a correlation of this switching event with the sexually dimorphic expression of carbohydrate metabolic genes. These results suggest that a comprehensive catalog of isoforms may lead to understanding the molecular basis for environmental sex determination of Daphnia. We also infer the applicability of the full-length isoform analyses to the elucidation of phenotypic plasticity in Daphnia.

Plant kleptomaniacs: geographical genetic patterns in the amphi-apomictic Rubus ser. Glandulosi (Rosaceae) reveal complex reticulate evolution of Eurasian brambles.

BACKGROUND AND AIMS: Rubus ser. Glandulosi provides a unique model of geographical parthenogenesis on a homoploid (2n = 4x) level. We aim to characterize evolutionary and phylogeographical patterns in this taxon and shed light on the geographical differentiation of apomicts and sexuals. Ultimately, we aim to evaluate the importance of phylogeography in the formation of geographical parthenogenesis. METHODS: Rubus ser. Glandulosi was sampled across its Eurasian range together with other co-occurring Rubus taxa (587 individuals in total). Double-digest restriction site-associated DNA sequencing (ddRADseq) and modelling of suitable climate were used for evolutionary inferences. KEY RESULTS: Six ancestral species were identified that contributed to the contemporary gene pool of R. ser. Glandulosi. Sexuals were introgressed from Rubus dolichocarpus and Rubus moschus in West Asia and from Rubus ulmifolius agg., Rubus canescens and Rubus incanescens in Europe, whereas apomicts were characterized by alleles of Rubus subsect. Rubus. Gene flow between sexuals and apomicts was also detected, as was occasional hybridization with other taxa. CONCLUSIONS: We hypothesize that sexuals survived the last glacial period in several large southern refugia, whereas apomicts were mostly restricted to southern France, whence they quickly recolonized Central and Western Europe. The secondary contact of sexuals and apomicts was probably the principal factor that established geographical parthenogenesis in R. ser. Glandulosi. Sexual populations are not impoverished in genetic diversity along their borderline with apomicts, and maladaptive population genetic processes probably did not shape the geographical patterns.

Identification of Parthenogenesis-Inducing Effector Proteins in Wolbachia.

Bacteria in the genus Wolbachia have evolved numerous strategies to manipulate arthropod sex, including the conversion of would-be male offspring to asexually reproducing females. This so-called "parthenogenesis induction" phenotype can be found in a number of Wolbachia strains that infect arthropods with haplodiploid sex determination systems, including parasitoid wasps. Despite the discovery of microbe-mediated parthenogenesis more than 30 yr ago, the underlying genetic mechanisms have remained elusive. We used a suite of genomic, computational, and molecular tools to identify and characterize two proteins that are uniquely found in parthenogenesis-inducing Wolbachia and have strong signatures of host-associated bacterial effector proteins. These putative parthenogenesis-inducing proteins have structural homology to eukaryotic protein domains including nucleoporins, the key insect sex determining factor Transformer, and a eukaryotic-like serine-threonine kinase with leucine-rich repeats. Furthermore, these proteins significantly impact eukaryotic cell biology in the model Saccharomyces cerevisiae. We suggest that these proteins are parthenogenesis-inducing factors and our results indicate that this would be made possible by a novel mechanism of bacterial-host interaction.

Variants of a major DNA satellite discriminate parental subgenomes in a hybrid parthenogenetic lizard Darevskia unisexualis (Darevsky, 1966).

Hybrid parthenogenetic animals are an exceptionally interesting model for studying the mechanisms and evolution of sexual and asexual reproduction. A diploid parthenogenetic lizard Darevskia unisexualis is a result of an ancestral cross between a maternal species Darevskia raddei nairensis and a paternal species Darevskia valentini and presents a unique opportunity for a cytogenetic and computational analysis of a hybrid karyotype. Our previous results demonstrated a significant divergence between the pericentromeric DNA sequences of the parental Darevskia species; however, an in-depth comparative study of their pericentromeres is still lacking. Here, using target sequencing of microdissected pericentromeric regions, we reveal and compare the repertoires of the pericentromeric tandem repeats of the parental Darevskia lizards. We found species-specific sequences of the major pericentromeric tandem repeat CLsat, which allowed computational prediction and experimental validation of fluorescent DNA probes discriminating parental chromosomes within the hybrid karyotype of D. unisexualis. Moreover, we have implemented a generalizable computational method, based on the optimization of the Levenshtein distance between tandem repeat monomers, for finding species-specific fluorescent probes for pericentromere staining. In total, we anticipate that our comparative analysis of Darevskia pericentromeric repeats, the species-specific fluorescent probes that we found and the pipeline that we developed will form a basis for the future detailed cytogenomic studies of a wide range of natural and laboratory hybrids.

Transcriptomic basis of sex loss in the pea aphid.

BACKGROUND: Transitions from sexual to asexual reproduction are common in eukaryotes, but the underlying mechanisms remain poorly known. The pea aphid-Acyrthosiphon pisum-exhibits reproductive polymorphism, with cyclical parthenogenetic and obligate parthenogenetic lineages, offering an opportunity to decipher the genetic basis of sex loss. Previous work on this species identified a single 840 kb region controlling reproductive polymorphism and carrying 32 genes. With the aim of identifying the gene(s) responsible for sex loss and the resulting consequences on the genetic programs controlling sexual or asexual embryogenesis, we compared the transcriptomic response to photoperiod shortening-the main sex-inducing cue-of a sexual and an obligate asexual lineage of the pea aphid, focusing on heads (where the photoperiodic cue is detected) and embryos (the final target of the cue). RESULTS: Our analyses revealed that four genes (one expressed in the head, and three in the embryos) of the region responded differently to photoperiod in the two lineages. We also found that the downstream genetic programs expressed during embryonic development of a future sexual female encompass ∼1600 genes, among which miRNAs, piRNAs and histone modification pathways are overrepresented. These genes mainly co-localize in two genomic regions enriched in transposable elements (TEs). CONCLUSIONS: Our results suggest that the causal polymorphism(s) in the 840 kb region somehow impair downstream epigenetic and post-transcriptional regulations in obligate asexual lineages, thereby sustaining asexual reproduction.

Distinct geographic parthenogenesis in spite of niche conservatism and a single ploidy level: A case of Rubus ser. Glandulosi (Rosaceae).

Asexual organisms often differ in their geographic distributions from their sexual relatives. This phenomenon, termed geographic parthenogenesis, has long been known, but the underlying factors behind its diverse patterns have been under dispute. Particularly problematic is an association between asexuality and polyploidy in most taxa. Here, we present a new system of geographic parthenogenesis on the tetraploid level, promising new insights into this complex topic. We used flow cytometric seed screen and microsatellite genotyping to characterise the patterns of distribution of sexuals and apomicts and genotypic distributions in Rubus ser. Glandulosi across its range. Ecological modelling and local-scale vegetation and soil analyses were used to test for niche differentiation between the reproductive groups. Apomicts were detected only in North-western Europe, sexuals in the rest of the range in Europe and West Asia, with a sharp borderline stretched across Central Europe. Despite that, we found no significant differences in ecological niches. Genotypic richness distributions suggested independence of the reproductive groups and a secondary contact. We argue that unless a niche differentiation (resulting from polyploidy and/or hybridity) evolves, the main factors behind the patterns of geographic parthenogenesis in plants are phylogeographic history and neutral microevolutionary processes, such as clonal turnover.

Genetic and Phenotypic Consequences of Local Transitions between Sexual and Parthenogenetic Reproduction in the Wild.

AbstractTransitions from sexual to asexual reproduction have occurred in numerous lineages, but it remains unclear why asexual populations rarely persist. In facultatively parthenogenetic animals, all-female populations can arise when males are absent or become extinct, and such populations could help to understand the genetic and phenotypic changes that occur in the initial stages of transitions to asexuality. We investigated a naturally occurring spatial mosaic of mixed-sex and all-female populations of the facultatively parthenogenetic Australian phasmid Megacrania batesii. Analysis of single-nucleotide polymorphisms indicated multiple independent transitions between reproductive modes. All-female populations had much lower heterozygosity and allelic diversity than mixed-sex populations, but we found few consistent differences in fitness-related traits between population types. All-female populations exhibited more frequent and severe deformities in their (flight-incapable) wings but did not show higher rates of appendage loss. All-female populations also harbored more ectoparasites in swamp (but not beach) habitats. Reproductive mode explained little variation in female body size, fecundity, or egg hatch rate. Our results suggest that transitions to parthenogenetic reproduction can lead to dramatic genetic changes with little immediate effect on performance. All-female M. batesii populations appear to consist of high-fitness genotypes that might be able to thrive for many generations in relatively constant and benign environments but could be vulnerable to environmental challenges, such as increased parasite abundance.

Aneuploidy during development in facultative parthenogenetic Drosophila.

From concatenated chromosomes to polyploidization, large-scale genome changes are known to occur in parthenogenetic animals. Here, we report mosaic aneuploidy in larval brains of facultatively parthenogenetic Drosophila. We identified a background of aneuploidy in D. mercatorum strains and found increased levels of aneuploidy in the larval brain tissue of animals arising parthenogenetically versus those arising from sexual reproduction. There is also intra-individual variation in germline-derived aneuploidy within the same strain. To determine if this is a general feature of facultative parthenogenesis in drosophilids, we compared sexually reproduced and parthenogenetic offspring from an engineered facultative parthenogenetic strain of D. melanogaster. In addition to germline-derived aneuploidy, this revealed somatic aneuploidy that increased by up to fourfold in parthenogens compared to sexually reproduced offspring. Therefore, the genetic combination identified in D. mercatorum that causes facultative parthenogenesis in D. melanogaster results in aneuploidy, which indicates that the loss of mitotic control resulting in parthenogenesis causes subsequent genome variation within the parthenogenetic offspring. Our findings challenge the assumption that parthenogenetic offspring are near genetic replicas of their mothers.

Temperature and hypoxia trigger developmental phenotypic plasticity of cardiorespiratory physiology and growth in the parthenogenetic marbled crayfish, Procambarus virginalis Lyko, 2017.

Attempting to differentiate phenotypic variation caused by environmentally-induced alterations in gene expression from that caused by actual allelic differences can be experimentally difficult. Environmental variables must be carefully controlled and then interindividual genetic differences ruled out as sources of phenotypic variation. We investigated phenotypic variability of cardiorespiratory physiology as well as biometric traits in the parthenogenetically-reproducing marbled crayfish Procambarus virginalis Lyko, 2017, all offspring being genetically identical clones. Populations of P. virginalis were reared from eggs tank-bred at four different temperatures (16, 19, 22 and 25 °C) or two different oxygen levels (9.5 and 20 kPa). Then, at Stage 3 and 4 juvenile stages, physiological (heart rate, oxygen consumption) and morphological (carapace length, body mass) variables were measured. Heart rate and oxygen consumption measured at 23 °C showed only small effects of rearing temperature in Stage 3 juveniles, with larger effects evident in older, Stage 4 juveniles. Additionally, coefficients of variation were calculated to compare our data to previously published data on P. virginalis as well as sexually-reproducing crayfish. Comparison revealed that carapace length, body mass and heart rate (but not oxygen consumption) indeed showed lower, yet notable coefficients of variation in clonal crayfish. Yet, despite being genetically identical, significant variation in their morphology and physiology in response to different rearing conditions nonetheless occurred in marbled crayfish. This suggests that epigenetically induced phenotypic variation might play a significant role in asexual but also sexually reproducing species.

Examining Wolbachia-Induced Parthenogenesis in Hymenoptera.

The maternally transmitted reproductive manipulator Wolbachia can impact sex ratios of its arthropod host by different mechanisms, ultimately promoting the spread of infection across a population. One of these reproductive phenotypes, parthenogenesis induction (PI), is characterized by the asexual production of female offspring, which in many cases results in an entirely female population. Cases of Wolbachia-mediated PI are most common in the order Hymenoptera, specifically in parasitoid wasps. The complex sex determination pathways of hymenopterans, their diverse life histories, the multiple cytogenetic mechanisms of PI, and the lack of males make functional studies of parthenogenesis induction challenging. Here, we describe the mechanisms of PI, outline methods to recognize and cure PI-Wolbachia infection, and note possible complications when working with PI-Wolbachia strains and their parthenogenetic hosts.

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes.

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Single hybrid population but multiple parental individuals at the origin of parthenogenetic rock lizards Darevskia sapphirina and D. bendimahiensis Schmidtler, & Eiselt Darevsky (1994) endemic to the area of Lake Van in East Turkey.

Among vertebrates, obligate parthenogenesis is only found in Squamata, where it always has a hybrid origin and a few lizard genera contain most of the known hybridogenous parthenogenetic taxa. Parthenogenesis thus seems to be pre-conditioned at the genus level, but it is not clear how often the encounter between two parental sexually reproducing species can result in the parthenogenetic offspring, nor whether the success of such hybridization event requires certain conditions or the specific time frame. To address this question, we studied the rock lizards of genus Darevskia, where a pair of parental species, D. valentini and D. raddei, as well as their parthenogenetic daughter species D. bendimahiensis and D. sapphirina, are found in close proximity NE of the Lake Van in East Anatolia. Using ddRAD-seq genotyping on 19 parental and 18 hybrid individuals, we found that (i) all parthenogenetic individuals from both D. bendimahiensis and D. sapphirina have a monophyletic origin tracing back to a single initial hybrid population, but their current genetic variation is geographically structured; (ii) unlike the most probable paternal ancestor, the genetically closest extant population of the maternal ancestor is not the geographically nearest one; and (iii) in the parthenogens, about 1% of loci carry multiple haplotypes, frequently differentiated by multiple substitutions. This pattern, in addition to biases in the relative frequency of haplotypes of maternal and paternal origin, does not appear compatible with a scenario of the entire parthenogenic clonal population having descended from a single pair of parental individuals. Instead, the data suggest that multiple parental individual ancestries still persist in the parthenogenetic gene pool. This supports the notion that although hybridization leading to parthenogenesis is generally rare at the level of species, it may be more common at the individual/population level once the right conditions are met.

Reproductive biology: A genetic recipe for parthenogenesis.

New work reveals differences in oogenic gene expression between parthenogenetic and sexually reproducing Drosophila mercatorum strains. Recapitulating those changes in D. melanogaster oocytes induced parthenogenesis in this normally sexually reproducing species, providing molecular insight into how these reproductive modes arise.

Contrasting Evolutionary Patterns Between Sexual and Asexual Lineages in a Genomic Region Linked to Reproductive Mode Variation in the pea aphid.

Although asexual lineages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles.

Demonstration of Parthenogenetic Reproduction in a Pet Ball Python (Python regius) through Analysis of Early-Stage Embryos.

Parthenogenesis is an asexual form of reproduction, normally present in various animal and plant species, in which an embryo is generated from a single gamete. Currently, there are some species for which parthenogenesis is supposed but not confirmed, and the mechanisms that activate it are not well understood. A 10-year-old, wild-caught female ball python (Python regius) laid four eggs without any prior contact with a male. The eggs were not incubated and, after 3 days, were submitted to the University of Parma for analysis due to the suspicion of potential embryo presence. Examination of the egg content revealed residual blood vessels and a small red spot, indicative of an early-stage embryo. DNA was extracted from the three deceased embryos and from the mother's blood, five microsatellites were analyzed to ascertain the origin of the embryos. The captive history data, together with the genetic microsatellite analysis approach, demonstrated the parthenogenetic origin of all three embryos. The embryos were homozygous for each of the maternal microsatellites, suggesting a terminal fusion automixis mode of development.

Soft selection reduces loss of heterozygosity in asexual reproduction.

The adaptive value of sexual reproduction is still debated in evolutionary theory. It has been proposed that the advantage of sexual reproduction over asexual reproduction is to promote genetic diversity, to prevent the accumulation of harmful mutations or to preserve heterozygosity. Since these hypothetical advantages depend on the type of asexual reproduction, understanding how selection affects the taxonomic distribution of each type could help us discriminate between existing hypotheses. Here, I argue that soft selection, competition among embryos or offspring in selection arenas prior to the hard selection of the adult phase, reduces loss of heterozygosity in certain types of asexual reproduction. Since loss of heterozygosity leads to the unmasking of recessive deleterious mutations in the progeny of asexual individuals, soft selection facilitates the evolution of these types of asexual reproduction. Using a population genetics model, I calculate how loss of heterozygosity affects fitness for different types of apomixis and automixis, and I show that soft selection significantly reduces loss of heterozygosity, hence increases fitness, in apomixis with suppression of the first meiotic division and in automixis with central fusion, the most common types of asexual reproduction. Therefore, if sexual reproduction evolved to preserve heterozygosity, soft selection should be associated with these types of asexual reproduction. I discuss the evidence for this prediction and how this and other observations on the distribution of different types of asexual reproduction in nature is consistent with the heterozygosity hypothesis.