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.

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.

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.

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.

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.

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.

A genetic basis for facultative parthenogenesis in Drosophila.

Facultative parthenogenesis enables sexually reproducing organisms to switch between sexual and asexual parthenogenetic reproduction. To gain insights into this phenomenon, we sequenced the genomes of sexually reproducing and parthenogenetic strains of Drosophila mercatorum and identified differences in the gene expression in their eggs. We then tested whether manipulating the expression of candidate gene homologs identified in Drosophila mercatorum could lead to facultative parthenogenesis in the non-parthenogenetic species Drosophila melanogaster. This identified a polygenic system whereby increased expression of the mitotic protein kinase polo and decreased expression of a desaturase, Desat2, caused facultative parthenogenesis in the non-parthenogenetic species that was enhanced by increased expression of Myc. The genetically induced parthenogenetic Drosophila melanogaster eggs exhibit de novo centrosome formation, fusion of the meiotic products, and the onset of development to generate predominantly triploid offspring. Thus, we demonstrate a genetic basis for sporadic facultative parthenogenesis in an animal.

Transcriptomics and the origin of obligate parthenogenesis.

Despite the presence of obligately parthenogenetic (OP) lineages derived from sexual ancestors in diverse phylogenetic groups, the genetic mechanisms giving rise to the OP lineages remain poorly understood. The freshwater microcrustacean Daphnia pulex typically reproduces via cyclical parthenogenesis. However, some populations of OP D. pulex have emerged due to ancestral hybridization and introgression events between two cyclically parthenogenetic (CP) species D. pulex and D. pulicaria. These OP hybrids produce both subitaneous and resting eggs parthenogenetically, deviating from CP isolates where resting eggs are produced via conventional meiosis and mating. This study examines the genome-wide expression and alternative splicing patterns of early subitaneous versus early resting egg production in OP D. pulex isolates to gain insight into the genes and mechanisms underlying this transition to obligate parthenogenesis. Our differential expression and functional enrichment analyses revealed a downregulation of meiosis and cell cycle genes during early resting egg production, as well as divergent expression patterns of metabolism, biosynthesis, and signaling pathways between the two reproductive modes. These results provide important gene candidates for future experimental verification, including the CDC20 gene that activates the anaphase-promoting complex in meiosis.

No recombination suppression in asexually produced males of Daphnia pulex.

Obligate parthenogenesis (OP) is often thought to evolve by disruption of reductional meiosis and suppression of crossover recombination. In the crustacean Daphnia pulex, OP lineages, which have evolved from cyclical parthenogenetic (CP) ancestors, occasionally produce males that are capable of reductional meiosis. Here, by constructing high-density linkage maps, we find that these males show only slightly and nonsignificantly reduced recombination rates compared to CP males and females. Both meiosis disruption and recombination suppression are therefore sex-limited (or partly so), which speaks against the evolution of OP by disruption of a gene that is essential for meiosis or recombination in both sexes. The findings may be explained by female-limited action of genes that suppress recombination, but previously identified candidate genes are known to be expressed in both sexes. Alternatively, and equally consistent with the data, OP might have evolved through a reuse of the parthenogenesis pathways already present in CP and through their extension to all events of oogenesis. The causal mutations for the CP to OP transition may therefore include mutations in genes involved in oogenesis regulation and may not necessarily be restricted to genes of the "meiosis toolkit." More generally, our study emphasizes that there are many ways to achieve asexuality, and elucidating the possible mechanisms is key to ultimately identify the genes and traits involved.

Evolutionary consequences of loss of sexual reproduction on male-related traits in parthenogenetic lineages of the pea aphid.

Transition from sexual reproduction to parthenogenesis constitutes a major life-history change with deep evolutionary consequences for sex-related traits, which are expected to decay. The pea aphid Acyrthosiphon pisum shows intraspecific reproductive polymorphism, with cold-resistant cyclically parthenogenetic (CP) lineages that alternate sexual and asexual generations and cold-sensitive obligately parthenogenetic (OP) lineages that produce only asexual females but still males. Here, the genotyping of 219 pea aphid lineages collected in cold-winter and mild-winter regions revealed contrasting population structures. Samples from cold-winter regions consisted mostly of distinct multilocus genotypes (MLGs) usually represented by a single sample (101 different MLGs for 111 samples) and were all phenotyped as CP. In contrast, fewer MLGs were found in mild-winter regions (28 MLGs for 108 samples), all but one being OP. Since the males produced by OP lineages are unlikely to pass on their genes (sexual females being rare in mild-winter regions), we tested the hypothesis that their traits could degenerate due to lack of selection by comparing male production and male reproductive success between OP and CP lineages. Male production was indeed reduced in OP lineages, but a less clear pattern was observed for male reproductive success: females mated with OP males laid fewer eggs (fertilized or not) but OP and CP males fertilized the same proportion of eggs. These differences may stem from the type of selective forces: male production may be counter-selected whereas male performances may evolve under the slower process of relaxed selection. The overall effective reproductive capacity of OP males could result from recent sex loss in OP lineages or underestimated reproductive opportunities.

Differentiation of uniparental human embryonic stem cells into granulosa cells reveals a paternal contribution to gonadal development.

Genomic imprinting underlies the mammalian requirement for sexual reproduction. Nonetheless, the relative contribution of the two parental genomes during human development is not fully understood. Specifically, a fascinating question is whether the formation of the gonad, which holds the ability to reproduce, depends on equal contribution from both parental genomes. Here, we differentiated androgenetic and parthenogenetic human pluripotent stem cells (hPSCs) into ovarian granulosa-like cells (GLCs). We show that in contrast to biparental and androgenetic cells, parthenogenetic hPSCs present a reduced capacity to differentiate into GLCs. We further identify the paternally expressed gene IGF2 as the most upregulated imprinted gene upon differentiation. Remarkably, while IGF2 knockout androgenetic cells fail to differentiate into GLCs, the differentiation of parthenogenetic cells supplemented with IGF2 is partly rescued. Thus, our findings unravel a surprising essentiality of genes that are only expressed from the paternal genome to the development of the female reproductive system.

Parthenogenesis in dipterans: a genetic perspective.

Parthenogenesis has been documented in almost every phylum of animals, and yet this phenomenon is largely understudied. It has particular importance in dipterans since some parthenogenetic species are also disease vectors and agricultural pests. Here, we present a catalogue of parthenogenetic dipterans, although it is likely that many more remain to be identified, and we discuss how their developmental biology and interactions with diverse environments may be linked to different types of parthenogenetic reproduction. We discuss how the advances in genetics and genomics have identified chromosomal loci associated with parthenogenesis. In particular, a polygenic cause of facultative parthenogenesis has been uncovered in Drosophila mercatorum, allowing the corresponding genetic variants to be tested for their ability to promote parthenogenesis in another species, Drosophila melanogaster. This study probably identifies just one of many routes that could be followed in the evolution of parthenogenesis. We attempt to account for why the phenomenon has evolved so many times in the dipteran order and why facultative parthenogenesis appears particularly prevalent. We also discuss the significance of coarse genomic changes, including non-disjunction, aneuploidy, and polyploidy and how, together with changes to specific genes, these might relate to both facultative and obligate parthenogenesis in dipterans and other parthenogenetic animals.