RESEARCH ARTICLE Post- meiotic mechanism of facultative parthenogenesis in gonochoristic whiptail lizard species David V Ho1,2†, Duncan Tormey3†‡, Aaron Odell1, Aracely A Newton3§, Robert R Schnittker3, Diana P Baumann3, William B Neaves3, Morgan R Schroeder3, Rutendo F Sigauke3, Anthony J Barley4, Peter Baumann1,2,5* 1Department of Biology, Johannes Gutenberg University, Mainz, Germany; 2Institute of Quantitative and Computational Biosciences, Johannes Gutenberg University, Mainz, Germany; 3Stowers Institute for Medical Research, Kansas City, United States; 4School of Mathematical and Natural Sciences, Arizona State University–West Valley Campus, Glendale, United States; 5Institute of Molecular Biology, Mainz, Germany Abstract 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 *For correspondence: high degree of homozygosity at microsatellite loci or low levels of heterozygosity detected in peter@baumannlab.org next- generation sequencing data. This has led to the proposal that second polar body fusion †These authors contributed following the meiotic divisions restores diploidy and thereby mimics fertilization. Here, we show equally to this work that FP occurring in the gonochoristic Aspidoscelis species A. marmoratus and A. arizonae results in Present address: ‡Synthego genome- wide homozygosity, an observation inconsistent with polar body fusion as the underlying Corporation, Redwood City, mechanism of restoration. Instead, a high- quality reference genome for A. marmoratus and anal- United States; §Missouri Western ysis of whole- genome sequencing from multiple FP and control animals reveals that a post- meiotic State University, Saint Joseph, mechanism gives rise to homozygous animals from haploid, unfertilized oocytes. Contrary to the United States widely held belief that females need to be isolated from males to undergo FP, females housed with Competing interest: The authors conspecific and heterospecific males produced unfertilized eggs that underwent spontaneous devel- declare that no competing opment. In addition, offspring arising from both fertilized eggs and parthenogenetic development interests exist. 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 Funding: See page 17 exposes the genetic load and explains the high rate of congenital malformations and embryonic Preprinted: 22 September 2023 mortality associated with FP in many species. Conversely, for animals that develop normally, FP Received: 18 February 2024 could potentially exert strong purifying selection as all lethal recessive alleles are purged in a single Accepted: 17 May 2024 generation. Published: 07 June 2024 Reviewing Editor: Detlef Weigel, Max Planck Institute for Biology Editor's evaluation Tübingen, Germany In this valuable paper, convincing evidence is provided for the production of facultatively parthe- Copyright Ho, Tormey et al. nogenetic whiptail lizards through a gametic duplication. The audience for the work will be broad, This article is distributed under given that parthenogenesis is such a fascinating topic. the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited. Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 1 of 23 Research article Genetics and Genomics Introduction Incidences of facultative parthenogenesis (FP) have been reported to occur in diverse vertebrate clades including bony fish (Lampert et al., 2007), sharks (Chapman et al., 2007; Chapman et al., 2008; Feldheim et al., 2010; Dudgeon et al., 2017), snakes (Dubach et al., 1997; Groot et al., 2003; Booth and Schuett, 2011; Germano and Smith, 2010; Booth et al., 2012; Kinney et al., 2013; Booth et al., 2014; Allen et al., 2018; Card et al., 2021), lizards (Lenk et al., 2005; Watts et al., 2006; Kratochvíl et al., 2020), crocodilians (Booth et al., 2023), and birds (Bartelmez and Riddle, 1924; Olsen and Marsden, 1954; Sarvella, 1973; Schut et al., 2008; Ramachandran and McDaniel, 2018; Parker et al., 2010). The phenomenon was originally mistaken for long- term sperm storage occurring in zoo environments where females were housed without current or recent access to conspecific males. The most parsimonious explanation was, therefore, that the animal had previ- ously been in contact with a male and that stored sperm was responsible for delayed fertilization (Booth and Schuett, 2011; Holt and Lloyd, 2010; Sever and Hamlett, 2002). However, more recent studies involving microsatellite (MS) and/or amplified fragment length polymorphism (AFLP) analyses revealed no paternal contributions, as all alleles detected in the offspring were only found in the maternal ancestors (Groot et al., 2003; Booth and Schuett, 2011; Shibata et al., 2017; Ryder et al., 2021; Levine et al., 2024). Females with no access to males producing solely male (ZW systems) or female (XY systems) offspring that only harbor maternal genetic markers are now considered hallmarks of facultative parthenogenesis. Nevertheless, clear examples of long- term sperm storage have also been documented in the recent literature (Levine et al., 2021), underscoring the need for molecular methods such as MS analysis or sequencing data to elucidate the underlying mechanisms. Originally thought to only occur in captivity, more recent reports indicate that FP occurs in natural populations as well (Booth et al., 2012; Fields et al., 2015). Serious concerns have been raised by conservation biologists, as species with dwindling population densities, including the endangered species Komodo dragon (Watts et al., 2006), small tooth sawfish (Fields et al., 2015), California condor (Ryder et al., 2021), and American crocodile (Booth et al., 2023) are overrepresented among reports of FP. While overrepresentation could be a consequence of an increased likelihood of detection in species that are the subject of intense research and conservation efforts, the observations raise the question if FP is an adaptive trait aiding in the colonization of new areas and mitigating the effects of popula- tion bottlenecks or is simply a neutral trait (Fields et al., 2015). The adaptive trait hypothesis would of course require successful reproduction of FP animals either sexually or parthenogenetically, which to date has only been documented in a few cases (Kratochvíl et al., 2020; Straube et al., 2016). At the same time, the association of FP with increased homozygosity constitutes a concern for conserva- tion biology, as an increase in FP within dwindling populations further accelerates the loss of genetic diversity, exposes deleterious alleles, and could compromise efforts to maintain the existing gene pool in selective breeding programs (Groot et al., 2003; Booth and Schuett, 2016). FP is also being studied as a desirable outcome in the commercial production of poultry. However, examination of tens of thousands of unfertilized eggs from several different avian species and strains has not resulted in economically sustainable hatching rates thus far. One of the highest hatch rates for unfertilized eggs is seen in the Beltsville small white turkey with a rate of 0.88% (Ramachandran and McDaniel, 2018; Olsen, 1975). A better understanding of the triggers and molecular mechanisms underlying FP and the fitness of the resulting offspring are, therefore, needed in a variety of contexts. These include: to understand a fundamental biological mechanism and its significance in vertebrate evolution, to aid in conservation efforts including captive breeding programs, and to possibly harness FP in an agricul- tural context (Ryder et al., 2021). The high level of homozygosity observed in animals produced by FP has been interpreted as evidence for polar body fusion following meiosis II, also known as automixis, leading to the restoration of diploidy in unfertilized eggs (Figure 1A; Chapman et al., 2007; Card et al., 2021; Booth et al., 2023; Booth and Schuett, 2016; Reynolds et al., 2012). If automixis involves the fusion of one of the meiotic products from the first polar body (central automixis), homozygosity will be concentrated near the chromosome ends and heterozygosity will be preferentially retained near the centromeres as premeiotic recombination strongly favors homologs over sister chromatids and homologs segre- gate during the first meiotic division. In contrast, a second polar body fusion (terminal automixis) would reunite sister chromatids, for which heterozygosity is preferentially seen near the chromosome termini (Figure 1B). In many cases, heterozygous and homozygous loci appeared to be inherited in FP Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 2 of 23 Research article Genetics and Genomics A ChromosomePairing Oocyte Recombination Meiosis I Meiosis II 2nd Polar 1st Polar 1st Polar Body Body Bodies B Oocyte 2nd Polar Body 1st Polar Bodies or or 1st Polar Body Fusion Post-meiotic 2nd Polar Body Fusion Duplication C D A. arizonae A. gularis A. marmoratus Figure 1. Overview. (A) Schematic of canonical meiosis. Only one pair of homologous chromosomes is shown using red and blue to distinguish homologs. (B) Schematic of main mechanisms by which a diploid oocyte may be produced in the context of facultative parthenogenesis. First polar body fusion, second polar body fusion, or post- meiotic duplication of chromosomes in the haploid gamete. (C) Photographs of Aspidoscelis arizonae with characteristic blue ventral coloration (top), A. gularis with light spots in dark fields that separate light stripes on dorsum (middle), and A. marmoratus with light and dark reticulated pattern on dorsum (bottom). (D) Microsatellite analysis for the three co-h oused animals and two offspring (ID 8449 and 8450) produced in this enclosure. Alleles are color- coded as follows: A. arizonae male (blue), A. gularis male (green), and A. marmoratus female (red). Differences in shading highlight the two alleles at heterozygous loci. Both offspring are homozygous at all loci with most alleles matching only maternal alleles. For MS16, offspring alleles are not shaded because of this allele being shared between the mother and the A. gularis male. Single nucleotide differences in size are common binning artifacts and, therefore, are not scored as different alleles. The online version of this article includes the following figure supplement(s) for figure 1: Figure supplement 1. Microsatellite analysis of eight loci for the three female Aspidoscelis marmoratus within the enclosure and the two offspring hatched in January 2009. Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 3 of 23 Research article Genetics and Genomics offspring (Groot et al., 2003; Allen et al., 2018; Card et al., 2021; Kratochvíl et al., 2020; Booth et al., 2023), but information as to the genomic location of these loci has been lacking. Central and terminal automixis are also distinguished by the extent of heterozygosity, with lower levels observed in next- generation sequencing data in snakes and crocodiles suggesting terminal automixis as the likely mechanism (Allen et al., 2018; Card et al., 2021; Booth et al., 2023). While facultative parthenogenesis occurs in a wide range of vertebrate species, true obligate parthenogenesis is limited to a few taxa of squamate reptiles including the North American whip- tail lizards of the genus Aspidoscelis (Avise, 2015). Historic hybridization events between distinct gonochoristic species in this clade has given rise to numerous hybrid individuals with the ability to reproduce clonally as all female lineages (Vanzolini, 1993; Reeder et al., 2002). In contrast to the increased homozygosity associated with FP, obligate parthenogenetic species are characterized by the long- term preservation of the high degree of heterozygosity that had its origin in the lineage- founding cross- species hybridization events. Our laboratory has a longstanding interest in the mechanism of obligate parthenogenesis in whiptail lizards (Lutes et al., 2010; Newton et al., 2016). In this context, we are maintaining and propagating individuals of several obligate parthenogenetic as well as gonochoristic species. MS analysis revealed over 20 incidences of FP in the marbled whiptail lizard, A. marmoratus and the Arizona striped whip- tail, A. arizonae. (The taxonomy of the genus Aspidoscelis has undergone frequent revisions and the maternal ancestor of the obligate parthenogenetic species A. neomexicanus was formerly known as the subspecies A. tigris marmoratus and the male ancestor as the subspecies A. inornatus arizonae (Barley et al., 2022a). For the purpose of this manuscript we follow the taxonomic conclusions by Barley et al., 2021). Whiptail lizards have an XX/XY sex determination system (Cole et al., 1969) and all FP offspring are consequently female. The identification of multiple incidences of FP provided us with the opportunity to investigate the mechanism of FP in whiptail lizards through next- generation sequencing. The generation of a genome assembly, in addition to whole-g enome sequencing, allowed us to distinguish between different mechanisms for restoring diploidy in FP animals. To address the question of whether FP is limited to animals in captivity, we examined reduced-r epresentation sequencing (RAD- seq) data of 321 whiptail lizards from 15 gonochoristic species sampled in nature. In aggregate, a combination of MS analysis, next-g eneration sequencing, and cytological analysis allows us to report on both the evidence and mechanism of FP in whiptail lizards and suggest that a baseline incidence of FP may coexist alongside sexual reproduction in some species. Results Identification of FP in A. marmoratus In the context of studying interspecific hybridization among gonochoristic species of whiptail lizards, three female A. marmoratus (ID 122, 4238, 4239) were housed with a male A. arizonae (ID 4272) and a male A. gularis (ID 302) for close to three years. During this period, seven hybrid offspring between A. marmoratus 122 and A. arizonae 4272 were produced and confirmed by MS analysis. These animals will be described in more detail in due course. Surprisingly, two female hatchlings emerged that resembled A. marmoratus rather than the expected products of hybridization with either A. arizonae or A. gularis. Genotyping revealed only a single allele for each of eight MS markers in the two offspring (ID 8449 and 8450, Figure 1D) and identified A. marmoratus (ID 122) as the mother (Figure 1—figure supplement 1). The mother and the male A. arizonae were each heterozygous at five of the eight markers and the A. gularis male at six. Further supporting a uniparental origin of 8449 and 8450, all alleles found in the offspring were also present in the mother (Figure 1D). For seven of the eight markers, neither male shared the allele found in the hatchling lizards, providing strong evidence that neither male fathered the offspring. For the remaining marker MS16, the A. marmoratus mother and the A. gularis male were homozygous for the same allele found in the two offspring, therefore, not allowing a conclusion to be based on this locus. It is important to note that for two of the markers (MS7 and MS15), the two offspring inherited different alleles from the mother, indicating that they are not genetically identical to each other, but have randomly inherited one of the maternal alleles at each locus. Two additional eggs (ID 8394 and 9070) were recovered from the same enclosure and found to contain developing embryos. MS analysis also revealed A. marmoratus 122 as the mother and Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 4 of 23 Research article Genetics and Genomics complete homozygosity at all loci. Given that all of these offspring are female, inherited only maternal alleles, and animal 122 had no history of being housed with a conspecific male during its lifetime, both interspecific hybridization and long- term sperm storage are all but ruled out and FP is strongly supported. FP animals of A. marmoratus presented a unique opportunity to examine the underlying molecular mechanism. The observed homozygosity at all MS loci further promised to aid in the generation of a high- quality genome assembly, as homozygosity circumvents the challenge of collapsing haplotypes into a consensus sequence (Kajitani et al., 2014). To increase homozygosity, inbreeding for 15–20 generations is common practice prior to whole- genome sequencing and genome assembly (Zhang et al., 2019; Fang et al., 2012). However, generation times of more than one year make this a costly and time- consuming strategy for many vertebrate species including A. marmoratus. Genome sequencing and de novo assembly The A. marmoratus genome is distributed over 23 chromosomes as previously demonstrated by meta- phase spread analysis (Lowe and Wright, 1966). We used flow cytometry to compare nuclear DNA content of A. marmoratus erythrocytes from whole blood with cells from three species with well- characterized genome sizes. The nuclear DNA content of A. marmoratus was close to that of Danio rerio (1.4 Gb) and we calculated a haploid genome size for A. marmoratus of 1.55 Gb (Figure 2— figure supplement 1A). Genomic DNA of FP animal 8450 was used to generate short insert paired- end, mate-p air (5 Kb, 8 Kb, 2–15 Kb, 40 Kb), and Chicago (Putnam et al., 2016) libraries for Illumina short-r ead sequencing. The paired- end and mate- pair reads were first assembled with Meraculous (Chapman et al., 2011) yielding an N50 of 1.6 Mb. The subsequent addition of Chicago reads and scaffolding with the HiRise pipeline by Dovetail Genomics produced an assembly of 1,639,530,780 bp distributed over 3826 scaf- folds (Supplementary file 1) and raised the scaffold N50 to 32.22 Mb (Figure 2—figure supplement 1B). With a BUSCO completeness score of 96% the A. marmoratus genome assembly is comparable to other recently released reptilian genome assemblies (Figure 2—figure supplement 1C). Over 98% of the assembled sequences are contained within 90 scaffolds of more than 1 Mb in length, making this assembly highly contiguous. Phylogenetic analysis of shared BUSCO genes with several other reference genomes (Xenopus, zebrafish, medaka, platyfish, tegu, green anole, chicken, mouse, rat, dog, cow, human) confirmed that A. marmoratus is most closely related to the tegu Salvator merianae, another representative of the family Teiidae (Figure 2—figure supplement 2). As transposable elements are a driving force in genome evolution, we examined the repeat content for the A. marmoratus genome. All classes of repeat elements combined amounted to 40.27% of the A. marmoratus assembly, only slightly below that found in other lizards S. merianae and Anolis carolinensis (Figure 2—figure supplement 3A). Strikingly, unclassified repeats make up the largest class of repeat elements in the A. marmoratus genome, an observation that parallels findings in S. merianae. However, a comparison between the unclassified repeats found in A. marmoratus, S. merianae, and A. carolinensis revealed few similar- ities with only around 10% of the unclassified repeats shared between A. marmoratus and S. meri- anae, and no significant overlap between these two species and A. carolinensis (Figure 2—figure supplement 3B). While further characterization of the unclassified repeat elements is needed, it is apparent that an impressive expansion of novel repeat element classes has occurred within this clade. To annotate the A. marmoratus genome, we assembled a total of 119,728 transcripts from RNA- seq data generated from blood and embryo using Trinity (Supplementary file 2). These tran- scripts were subsequently used in the MAKER2 gene annotation pipeline, yielding 25,856 protein- coding genes and 44,461 protein isoforms (Supplementary file 3). To assign a putative function, we used BLASTp to query the UniProtKB/Swiss-P rot database (UniProt Consortium, 2018) and found significant hits for 76% of the putative protein- coding genes. Our assembly and annota- tion pipelines yielded 40 HOX genes and 2 EVX genes in four gene clusters (Figure 2—figure supplement 4). The HOX gene clusters are highly conserved among tetrapods and their complete presence and shared order within each cluster serve as a measure of assembly quality (Kuraku and Meyer, 2009). Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 5 of 23 Research article Genetics and Genomics Assessment of heterozygosity The presence of only one allele for each of the examined MS markers already suggested widespread homozygosity in A. marmoratus produced by FP, consistent with similar observations in other verte- brate species. The highly contiguous genome assembly now afforded us the opportunity to probe the mechanism of FP by searching for regions of heterozygosity and mapping their relative genomic loca- tions. Towards this aim, we performed whole-g enome sequencing for an additional nine animals: four of them produced by FP (ID 12512, 12513, 6993, 9177), two mothers (ID 122, 9721), and three unre- lated control animals (ID 003, 001, S30700; Supplementary file 4). Each mother and the controls were heterozygous at several MS markers confirming their origin through sexual reproduction. Following alignment to the reference genome (Supplementary file 5), we defined heterozygous sites as those covered by an even number of reads with two alleles supported by the same number of reads. Sites covered by an odd number of reads were filtered out for this initial analysis. This stringent requirement was chosen to limit the search to apparent heterozygous sites with strong support, decreasing the chance of false positives. For all ten individuals, the average sequencing coverage ranged between 15.91 and 20.08 (Figure  2—figure supplement 5). For FP animals, the number of heterozygous sites in a 10  Kb sliding window approaches zero for all sites with mean coverage (Figure 2A). For positions with coverage greater than the average, an increase in apparent heterozygosity was observed, due to the collapse of repetitive sequences during the assembly process. Based on this observation, we limited further analysis to positions in the genome where the coverage is equal to the mean sequencing depth (as defined by rounding the mean sequencing coverage value to the next even integer). For example, for animal 003, the average sequencing coverage is 18.31 (Figure 2—figure supplement 5A) and we only considered sites with a coverage of 20 (Figure  2—figure supple- ment 6A). This generated between 30,769 and 53,416 heterozygous sites for which two alleles were equally supported in the sexually produced mothers and control animals (Figure 2—figure supple- ment 6A–E). Far fewer heterozygous sites (between 649 and 928) were observed in the FP animals (Figure 2—figure supplement 6F–J). Plotting the heterozygous sites according to their position in the reference genome illustrates not only their sparsity in the genomes of FP animals, but also reveals their random distribution (Figure 2B). If FP genomes were the product of automixis, regions of homozygosity would be interspersed with regions of heterozygosity. The extent of heterozygosity within the latter would be the same as that observed in the respective mother. The few apparently heterozygous sites identified in FP animals are, therefore, not supporting either form of automixis but are most likely the result of over-a ssembly of repetitive regions (ie. collapsing paralogous loci into a single representative sequence) and a combination of biological (ie. somatic mutations) and technical errors (ie. PCR and sequencing errors). When examining each mother- daughter group, the average number of heterozygous sites per 10 Kb window was greater in the sexually produced mothers across the entire assembly (Figure 2C). For most of the assembly, the average number of heterozygous sites remained close to zero for FP animals. The most notable exception was a region around genomic coordinate 1.0  Gb, but even there the extent of apparent heterozygosity remained below 50% of what is observed in this region for each of the mothers. Examination of the scaffold in question (Scaffold 45) revealed 167 genes annotated as homologous to vomeronasal 2 receptor 26 (Vmn2r26; Figure 2—figure supplement 7, Figure 2—figure supplement 8). Members of this subfamily of receptors are found on the microvilli of the vomeronasal organ, where they are responsible for pheromone detection and play a signifi- cant role in social and environmental responses (Ryba and Tirindelli, 1997; Houck, 2009; Su et al., 2009). Given that this genomic region harbors a large cluster of highly similar genes, the most parsi- monious explanation for the elevated level of apparent heterozygosity is over- assembly. This conclu- sion is further supported by the increase in apparent heterozygosity in this region for the mothers and control animals. In aggregate, our analysis strongly supports genome-w ide homozygosity for FP animals, inconsistent with either central or terminal automixis. Instead, the results favor a post-m eiotic mechanism that restores diploidy by replicating the haploid genome residing in the oocyte following completion of the two meiotic divisions and thereby establishing genome- wide homozygosity in the offspring. Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 6 of 23 Research article Genetics and Genomics A ID 50.0 003 001 Sexually 10.0 122 produced S30700 9721 1.0 6993 9177 8450 FP offspring 12512 0.1 12513 4 8 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 88 92 96 100 Coverage B ID 003 001 122 Sexually produced S30700 9721 6993 9177 8450 FP offspring 12512 12513 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 Genomic Coordinate (Gb) C 3 ID 2 122 1 8450 (FP) 0 3 9721 2 12512 (FP) 1 12513 (FP) 0 Figure 2. Genome- wide homozygosity in animals produced by facultative parthenogenesis. (A) Effect of coverage on the apparent rate of heterozygosity based on evenly split read counts supporting two alleles. Analysis of whole- genome sequencing data for five sexually produced animals (ID 003, 001, 122, S30700, 9721) and five individuals produced by facultative parthenogenetic (FP) animals (ID 6993, 9177, 8450, 12512, 12513) were aligned to the reference genome. In FP animals, the number of heterozygous sites approaches zero for sites with mean coverage (x=̄ 18.37). (B) Scaffolds are ordered from largest to smallest on the x- axis. Red lines indicate borders between ordered scaffolds. Each black dot represents a heterozygous position in the genome defined by having a sequencing coverage equal to the average and equal support for only two alleles. The y-a xis position of each data point is a random value between bounds of area shown to spread the data and better illustrate the density of heterozygous sites. (C) Average heterozygous sites, as defined in (B), per 10 Kb window for mothers (orange) and respective FP daughters (purple and pink). The online version of this article includes the following figure supplement(s) for figure 2: Figure supplement 1. Genome assembly of Aspidoscelis marmoratus. Figure supplement 2. Maximum likelihood tree for 13 vertebrate genomes, based on 1333 single-c opy BUSCOs detected across all species analyzed. Figure supplement 3. Identification of unclassified repeat elements in the A. marmoratus genome. Figure supplement 4. Aspidoscelis marmoratus HOX gene clusters. Figure supplement 5. Distribution of sequence coverage across the genome for each animal. Figure 2 continued on next page Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 7 of 23 Avg. Het. Sites / 10 Kb Het. Sites / 10 Kb (log10 scale) Research article Genetics and Genomics Figure 2 continued Figure supplement 6. Analysis of all positions in the genome at average coverage, for which reads support exactly two alleles. Figure supplement 7. A 1Mb sliding window for Scaffold 45 showing the annotated vomeronasal 2 receptors homologs. Figure supplement 8. Phylogenetic tree showing the evolutionary distance between mouse V2Rs (red branches) and Aspidoscelis marmoratus homologs (black branches). Cryptic FP in A. arizonae Following the identification of several A. marmoratus generated by FP, we genotyped individuals from two other gonochoristic species housed in our laboratory. While no cases of FP were identified among 80 A. gularis produced eggs in captivity, we identified eight incidences of FP among 832 A. arizonae records between October 2007 and July 2018. During the same period, we recorded 15 incidences of FP among 286 A. marmoratus records (Supplementary file 6). Notably, in all cases, eggs undergoing FP development had been laid in enclosures where females were housed with conspecific males or males of a sister species known to mate with the heterospecific females. Isolation from mating part- ners was thus not a significant factor in triggering FP. In one enclosure, three A. arizonae females (ID 12850, 12851, 12852) were housed with a conspecific male (ID 12849; Figure 3A). MS analysis of four hatchlings that originated from a single clutch laid in this enclosure identified animal 12852 as the mother of all four animals. Unexpectedly, two of her offspring were homozygous at all eight loci examined, whereas the two others were heterozygous at all loci, identifying animal 12849 as their father (Figure 3B). Therefore, both fertilized and unfertilized eggs developed alongside each other within the same clutch. Mixoploid erythrocytes and developmental defects Microscopic examination of blood from a newly hatched FP lizard revealed a striking bimodality in the sizes of erythrocyte nuclei when compared to a sexually produced animals (Figure 4A, B). Nuclear size correlates well with DNA content measurements (Walker et al., 1991), suggesting the presence of mixoploidy in the FP animal. Whereas most cells closely resembled those observed in the blood from the sexually produced animal, approximately 10% of red blood cells from the FP animal harbored smaller nuclei, consistent with half the amount of DNA (Figure 4B). In addition, 1.27% contained two small nuclei, indicating that the final cytokinesis during erythrocyte differentiation had failed for some haploid progenitor cells. DNA content analysis by flow cytometry confirmed the presence of haploid cells (Figure 4C). In the blood of sexually produced animals, no haploid or bi- nucleated cells were observed. These observations raise the possibility that the embryonic development of FP animals is initiated with consecutive divisions of a haploid, unfertilized oocyte. At a later stage in development, diploid cells most likely arise via failed cytokinesis. From that point forward, both haploid and diploid cells coexist, and the embryo develops in a mixoploid state. Indicative of a more widespread phenom- enon, mixoploidy was also observed in another FP A. marmoratus and in A. arizonae (Figure 4— figure supplement 1). The observed fraction of haploid cells was closer to 1% in these instances. Genome- wide loss of heterozygosity exposes functionally compromised alleles that were previously covered by intact alleles on the homologous chromosomes. Depending on the extent of this genetic load, one would expect a substantial fraction of oocytes to not develop at all or for defects to mani- fest at various stages of embryonic and post-e mbryonic development. Indeed, of the 23 incidences of FP examined here, only 14 hatched, while the remaining lizards died in ovum (Supplementary file 6). For these nine unhatched eggs, we isolated developed embryos shortly after the expected hatch date and confirmed FP origin by MS analysis. The clutches that produced the 23 confirmed cases of FP contained an additional 24 eggs. For these, development did not initiate or terminated at an earlier stage of development precluding MS analysis. Based on the uncertainty regarding how many of these eggs underwent partial FP development, the incidence of FP may be even higher than reported here. The observation of various malformations in several of the FP embryos and hatchlings further supports that genome- wide homozygosity unmasks deleterious alleles. Notable developmental defects included craniofacial abnormalities such as misaligned jaws, agenesis of eyes, missing limbs, and failed abdominal closure (Figure 4—figure supplement 2). Only six out of 16 FP animals (37.5%) Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 8 of 23 Research article Genetics and Genomics A B 12849 12852 16213 16214 16215 16216 Figure 3. Facultative parthenogenesis is also found in Aspidoscelis arizonae. (A) Microsatellite analysis for the four co-h oused adult animals (ID 12849, 12850, 12851, 12852) and the four hatchlings (ID 16213, 16214, 16215, 16216) produced in this enclosure. Alleles are color-c oded for each potential parent: 12849 male (blue), 12850 female (green), 12851 female (orange), 12852 female (red). Differences in shading highlight the two alleles at heterozygous loci. Offspring 16213 and 16214 are heterozygous at all loci, with most loci having one allele matching 12849 and one allele matching 12852. Offspring 16215 and 16216 are homozygous at all loci, with most alleles matching only the 12852 female. Non-s haded offspring alleles indicate ambiguous inheritance as multiple adult animals share the same allele. Single nucleotide differences in size are common binning artifacts and, therefore, are not scored as different alleles. (B) Pedigree shows the relationship between the four offspring. The single clutch of four contains both sexually (yellow) and facultative parthenogenetically (blue) produced offspring. hatched with no discernable developmental defects (Figure 4—figure supplement 2 A- B). This is in stark contrast to sexually produced animals, where over 98% of hatchlings (n=687) showed no abnormalities. Additionally, most of the defects noted in sexually produced animals were less severe than in FP animals including bulges in tails or truncated digits. While we have not recorded instances of FP animals producing offspring via FP, as described for the whitespotted bamboo shark (Straube et al., 2016), FP A. marmoratus 8450 did produce two eggs while housed in isolation, but these failed to hatch. Analysis of the ovaries of FP animal 8450 as well as germinal vesicles of its FP sister 8449 Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 9 of 23 Research article Genetics and Genomics A B Diploid 89.5% Binucleated 1.27% Haploid 9.22% C 2n 3n 9177 DNA content (Propodium Iodide) Figure 4. Detection of mixoploidy associated with facultative parthenogenesis. (A) Giemsa staining of erythrocytes from a sexually produced Aspidoscelis marmoratus (ID 14744). All cells are diploid (n=601). Scale bar corresponds to 10 µm. (B) Giemsa staining of erythrocytes from a newly hatched facultative parthenogenetic A. marmoratus (ID 9177). Diploid (n=844), smaller haploid (n=87), and binucleated (n=12) cells are evident. Scale bar corresponds to 10 µm. (C) DNA content from erythrocytes determined by propidium iodide staining and detection by flow cytometry. Samples are from a sexually produced A. marmoratus (2n, ID 5358), an obligate triploid parthenogen A. exanguis (3n, ID 4950), and facultative parthenogenesis (FP) A. marmoratus (ID 9177). Number of events scored by flow cytometry were 44,145 (2n), 44,043 (3n), and 44,060 (9177). The FP 9177 sample contained an additional peak to the left of the 2 n peak (90.04%), indicating the presence of haploid cells (9.62%). The online version of this article includes the following figure supplement(s) for figure 4: Figure supplement 1. Mixoploidy detected in both A. marmoratus and A. arizonae. Figure supplement 2. Animals produced by facultative parthenogenesis. Figure supplement 3. Ovaries of Aspidoscelis marmoratus facultative parthenogenesis (FP) animal 8450 and germinal vesicles of FP sister 8449 revealed no differences in structure and anatomy compared to fertile sexually reproducing animals. revealed no differences in structure and anatomy compared to fertile sexually reproducing animals (Figure 4—figure supplement 3). Putative incidences of FP in wild populations of whiptail lizards As FP has been associated with captivity in most species where it has been reported, we examined restriction site- associated DNA sequencing (RAD- seq) data for wild animals across 15 gonochoristic species (Figure 5). Because RAD-s eq is a form of reduced-r epresentation sequencing (Rivera- Colón and Catchen, 2022), we limited further analysis to the 321 individuals that had an average sequencing Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 10 of 23 Research article Genetics and Genomics 0.5 0.4 0.3 0.2 0.1 ep s ae rti us piiu t p ari s tus tus us us lis ki i ris us ris sti c nizo b sta de gu l a a m utt rn sis or at ta ac la at tig u ar co n o i e s c a ine g g in attn r m cid s xl a lin e ma oc se Species Figure 5. Heterozygosity estimates of whiptail lizards collected in nature. Percent heterozygosity estimates from reduced- representation sequencing (RAD- seq) for 321 whiptail lizards from 15 species. All individuals had an average coverage of at least 20. Each point is an individual, and percent heterozygosity was calculated only for sites where the coverage is equal to the average sequencing coverage. Five points with black borders indicate individuals (one angusticeps and four deppii) with low levels of heterozygosity. The heterozygosity of Aspidoscelis deppii ID LDOR30 (marked with arrow) is far less than that observed for individuals of the same species (Rosner’s Test for Outliers within deppii individuals, p<0.001), having only one called heterozygous position. Species (sample size): angusticeps (Booth et al., 2023), arizonae (Lenk et al., 2005), burti (Germano and Smith, 2010), costatus (Shibata et al., 2017), deppii (Ryba and Tirindelli, 1997), gularis (Olsen and Marsden, 1954), guttatus (Reeder et al., 2002), inornatus (Ryder et al., 2021), lineattissimus (Groot et al., 2003), marmoratus (Allen et al., 2018), occidentalis (Dudgeon et al., 2017), sackii (Groot et al., 2003), scalaris (Streisinger et al., 1981), sexlineatus (Germano and Smith, 2010), tigris (Card et al., 2021). The online version of this article includes the following source data for figure 5: Source data 1. Source data of Figure 5. coverage of at least 20. Computational data analysis revealed five animals (one A. angusticeps and four A. deppii) that had very low levels of heterozygosity at positions of average coverage (<0.05% hetero- zygous positions). In contrast, the average level of heterozygosity was 0.261% across the dataset, with the highest heterozygosity value at 0.578%. Of the five low heterozygosity animals, A. deppii (ID LDOR30) showed the most striking level of homozygosity affecting all sites but one (Rosner’s Test for Outliers within A. deppii individuals, log-t ransformation, R=5.127, λ=3.928, p<0.001). This is consis- tent with the pattern of homozygosity observed with the whole-g enome sequencing from FP animals produced in the laboratory. Further fieldwork and analysis will be required to assess the level of FP in natural populations of gonochoristic Aspidoscelis species (and other factors that could influence the observed heterozygosity such as population size, levels of hybridization, and inbreeding). To assess whether all whiptail species can produce viable offspring through FP, larger and broader datasets will be required to compare the incidence of FP between species, especially because animals with developmental defects associated with FP would not have hatched or survived in the wild and would, therefore, not have shown up in the RAD-s eq dataset. Discussion In this study, we report over 20 incidences of facultative parthenogenesis in marbled and Arizona striped whiptail lizards. By sequencing and assembling a highly contiguous A. marmoratus genome, we were able to refute automixis as the underlying mechanism in whiptail lizards. Instead, genome- wide homozygosity raises a possibility of a post- meiotic mechanism involving the activation of embryonic development in unfertilized haploid oocytes. Even though FP whiptail lizards are largely Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 11 of 23 Percent heterozygosity Research article Genetics and Genomics comprised of homozygous diploid cells, a fraction of haploid cells persists through development and is readily detectable in young adults. Such mixoploidy and genome- wide homozygosity come at a price. In eight confirmed and 24 suspected cases of FP, development ceased prior to hatching and most FP animals that hatched showed congenital defects. Nevertheless, FP was observed at a rate of 1% and 5% in A. arizonae and A. marmoratus, respectively, and this occurred in the presence of mating partners. Interestingly, these rates are similar to what has been reported for wild populations of two North American pitviper species (Booth et  al., 2012). Our findings indicate that FP is far more common in some vertebrate species than previously thought. The purifying selection associated with homozygosity may be an important force in generating additional resilience to counteract the effects of population bottlenecks and inbreeding depression. However, support for this hypothesis is predicated on the fitness and reproduction of FP offspring and, therefore, more long-t erm studies on seemingly healthy individuals of FP origin are needed. The start of embryonic development is tightly coupled to fertilization in many vertebrates as the sperm entering the oocyte triggers a signaling cascade that is essential for the completion of female meiosis and initiation of cell division following karyogamy (Sagata, 1996). This process has been mimicked in the laboratory by piercing frog oocytes with a needle to trigger the signal to complete meiosis and initiate replication and division in the haploid oocyte (Le Peuch et al., 1985; Wolf, 1974). In zebrafish, homozygous embryos are routinely generated by fertilization of oocytes with UV- irradi- ated sperm, a treatment that destroys the paternal DNA (Streisinger et al., 1981). An oocyte treated in this manner will replicate the intact maternal genome in the absence of karyogamy. If the haploid oocyte is then subjected to heat shock treatment, cytokinesis is prevented resulting in a pseudodip- loid oocyte that undergoes a second round of DNA replication followed by mitosis. Thus an entirely homozygous diploid embryo starts to develop (Kroeger et al., 2014). In contrast to human interven- tion forcing two consecutive rounds of DNA replication to occur without intervening mitosis at the start of embryonic development, our data indicate that one or multiple rounds of DNA replication and mitosis take place in some haploid oocytes of whiptail lizards prior to a skipped mitosis yielding a homozygous diploid cell followed by mixoploid development. Initiation of development in a haploid state may be conserved in avian species as unfertilized turkey eggs can yield embryos that contain 40% of haploid cells at blastoderm followed by a reduction to 1.3% within the blood of hatched birds (Cassar et al., 1998). Depending on the tissue and ploidy distribution, the presence of haploid cells may contribute to abnormal development of specific tissues reported here and elsewhere (Ito et al., 1991; Tanaka et al., 2004). Successful embryonic development from haploid cells that restore the diploid state by duplication has also been observed in a stick insect species (Pijnacker, 1969). In whiptail lizards, we have not been able to examine post- meiotic oocytes as locating the post- meiotic nucleus within a large yolked egg is inherently difficult. The difficulty is compounded by the unpredict- ability of which eggs will undergo FP development and the need to sacrifice animals to remove eggs. In addition to mixoploidy, genome-w ide homozygosity constitutes another obstacle to normal development as each recessive deleterious allele is exposed in either the hemizygous state (haploid cells) or homozygous state (diploid cells). Indeed, arrested development and abnormal phenotypes are observed in FP whiptails, as well as in FP animals across many other species (Booth et al., 2023; Booth and Schuett, 2016; Adams et al., 2023; Olsen, 1973). It is important to note though that some whiptails of FP origin developed normally, much like their sexually produced counterparts. At the population level, FP leads to a precipitous reduction in genetic diversity as only one set of alleles is inherited in the next generation. While FP could be an adaptive trait to bridge population bottlenecks when mate encounters are infrequent, small populations already rely heavily on inbreeding and FP further reduces the size of the gene pool (Watts et al., 2006; Ryder et al., 2021). While FP can be considered the most extreme example of inbreeding, it is also the most powerful example of genetic purging as it eliminates most deleterious alleles in a single generation. FP in whip- tail lizards and other species could, therefore, be considered a reproductive strategy, akin to mixed- mating systems in plants (Goodwillie et al., 2005). The ability to produce offspring via two different strategies provides a level of reproductive assurance in plants (Busch and Delph, 2012). Indeed, we see parallels to this in our own data in which female whiptails have produced offspring via both sexual reproduction and FP on separate occasions or simultaneously within a single clutch. Within plant species with mixed- mating, there are differences between the rates of selfing and outbreeding between populations, and hypotheses as to why these differences occur include limited pollinator Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 12 of 23 Research article Genetics and Genomics visitation and resource availability (Whitehead et al., 2018). To assess whether the co- occurrence of sexual and FP reproduction in vertebrates can indeed be considered a reproductive strategy rather than biological noise will require further studies to assess the reproductive competence and fecundity of offspring produced by either mode of reproduction. To gain a better understanding of the origin and outcomes of FP in whiptail lizards, it will also be important to identify the triggers. It has been proposed that lack of or limited mate encounters triggers FP, but our data in combination with many other reports (Booth et al., 2012; Booth et al., 2014; Kratochvíl et al., 2020; Ryder et al., 2021; Booth et al., 2011; Feldheim et al., 2023; Larose et al., 2023) rejects the idea that this is the key trigger. Recent work identifying key cell cycle genes inducing FP in two species of Drosophila (Sper- ling et al., 2023) and selection resulting in higher incidences of parthenogenesis in birds (Parker et al., 2010; Olsen, 1975; Olsen et al., 1968) suggest a genetic basis for the initiation of FP. Our study adds two species of whiptail lizards to a growing list of vertebrates capable of FP and establishes that it occurs alongside sexual reproduction in the presence of males. Using whole- genome sequencing, we demonstrate that post- meiotic genome duplication is the underlying mecha- nism. One must now consider the possibility that FP is an adaptive trait and that low rates of successful FP could contribute significantly to genome purification. Sexually mature FP offspring will have a low genetic load and only pass on neutral or mildly deleterious alleles to the next generation. However, a role for FP as an adaptive trait hinges on further studies demonstrating the ability of parthenogens to reproduce themselves either through further FP or sexually. If successive reproduction occurs, FP may reduce the frequency of deleterious alleles within a population, as well as provide reproductive assur- ance when males are scarce. Additional whole-g enome sequencing data for species with documented FP will be needed for a better understanding of the genetic basis, propensity, and evolutionary signif- icance of FP. Materials and methods Key resources table Reagent type (species) or resource Designation Source or reference Identifiers Additional information Biological sample Erythrocytes, tail clippings, (Aspidoscelis spp.) liver isolation This paper, wild populations See Ethics Statement Chemical compound, drug Giemsa stain Sigma GS500 0.40% Chemical compound, drug Schiff’s reagent Fisher Scientific #SS32- 500 Commercial assay or kit Roche gDNA Isolation Kit Roche #11814770001 Commercial assay or kit KAPA HTP kit KAPA KK8234 Nextera Mate- Pair Library Commercial assay or kit Prep Kit Illumina FC- 132–1001 TruSeq RNA Library Prep Kit Commercial assay or kit v2 Illumina RS- 122–2001 Sequence-b ased reagent Oligos for MS analysis This paper See Supplementary file 7 https://jgi.doe.gov/data-and-tools/ Software, algorithm Meraculous 2.0 software-tools/meraculous/ RRID:SCR_010700 Software, algorithm BUSCO 3.0.1 https://busco.ezlab.org/ RRID:SCR_015008 https://cme.h-its.org/exelixis/web/ Software, algorithm RAxML 8.2.11 software/raxml/ RRID:SCR_006086 https://www.repeatmasker.org/ Software, algorithm RepeatModeler 1.0.11 RepeatModeler/ RRID:SCR_015027 Software, algorithm RepeatMasker 4.0.9 https://www.repeatmasker.org/ RRID:SCR_012954 https://blast.ncbi.nlm.nih.gov/Blast. Software, algorithm BLAST 2.6.0 & 2.9.0+ cgi RRID:SCR_004870 Software, algorithm BWA 0.7.15 https://bio-bwa.sourceforge.net/ RRID:SCR_010910 Continued on next page Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 13 of 23 Research article Genetics and Genomics Continued Reagent type (species) or resource Designation Source or reference Identifiers Additional information https://broadinstitute.github.io/ Software, algorithm Picard 1.119 picard/ RRID:SCR_006525 https://gatk.broadinstitute.org/hc/ Software, algorithm GATK 3.5 en-us RRID:SCR_001876 https://github.com/lh3/seqtk Li, 2016 Software, algorithm seqtk 1.2- r94 1.2-r 94 https://github.com/pysam- developers/pysam pysam- developers, 2017 Software, algorithm pysam 0.12.0.1 0.12.0.1 https://github.com/alimanfoo/ pysamstats Miles, 2015 Software, algorithm pysamstats 0.24.3 0.24.3 https://github.com/trinityrnaseq/ trinityrnaseq trinityrnaseq, 2015 Software, algorithm Trinity 2.0.6 https://sourceforge.net/projects/ Software, algorithm seqclean seqclean https://www.yandell-lab.org/software/ Software, algorithm MAKER2 2.31.8 maker.html RRID:SCR_005309 https://interproscan-docs. Software, algorithm Interproscan 5.13–52.0 readthedocs.io/en/latest/ RRID:SCR_005829 https://www.ebi.ac.uk/about/ vertebrate-genomics/software/ Software, algorithm Exonerate 2.4.0 exonerate RRID:SCR_016088 Software, algorithm Geneious 10.1.3 https://www.geneious.com/ RRID:SCR_010519 https://catchenlab.life.illinois.edu/ Software, algorithm Stacks 2.62 stacks/ RRID:SCR_003184 Software, algorithm Micromanager 1.4 https://micro-manager.org/ RRID:SCR_000415 Software, algorithm Flowjo treestar https://www.flowjo.com/ RRID:SCR_008520 Microsatellite analysis DNA was extracted from tail samples for microsatellite genotyping as described in Lutes et al., 2011. PCR products were analyzed by capillary electrophoresis on a 3730 DNA Analyzer and data was analyzed using GeneMapper (v. 4.0). Primer information can be found in Supplementary file 7. Genome size estimation The genome size of A. marmoratus was estimated by fluorescence-a ctivated cell sorting (FACs), in which a standard curve correlating fluorescence intensity of DNA- bound propidium iodide with known genome sizes was generated using cells from fruit flies, zebrafish, and mouse, and then comparing fluorescent intensity with that of erythrocytes from A. marmoratus. Samples were stained using the Sigma PI staining preparation and analyzed on the Influx cytometer. PI fluorescence was collected using the PI Texas red detector with linear amplification and data analysis was performed in FlowJo and Microsoft Excel. DNA isolation, sequencing, and genome assembly for A. marmoratus All genome sequencing libraries generated for the purpose of the A. marmoratus genome assembly were derived from the FP animal 8450. The liver tissue was first dissociated in a 10  mL Dounce Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 14 of 23 Research article Genetics and Genomics homogenizer using the tight- fitting pestle and then processed using the Roche gDNA Isolation Kit (#11814770001, MilliporeSigma, St. Louis, MO, USA). A short insert, high-c overage library was generated using the KAPA HTP kit (KK8234), with 1 µg of gDNA. The resulting library was size selected for fragments between 500–850 bp on a Pippin Prep (Sage Science). Two 40 Kb mate-p air libraries were generated by Lucigen from 1 µg of gDNA using the CviQl and BfaI restriction enzymes, respectively. Each library was sequenced on the Illumina MiSeq using the MiSeq Reagent Kit v2 (500 cycles). An additional three mate- pair libraries were generated, spanning distances of 5 Kb, 8 Kb, and 2–15 Kb, using the Illumina Nextera Mate- Pair Library Prep Kit and 1 µg of gDNA for each. Size selection used the Gel-P lus protocol with Pippin for the 5 and 8  Kb libraries, and the Gel- Free protocol for the 2–15  Kb library. All three libraries were pooled and sequenced on three separate RapidSeq flow cells on an Illumina HiSeq 2500. Chicago libraries were prepared at Dovetail Genomics LLC, Santa Cruz, CA, USA from liver tissue to generate read pairs spanning distances up to 140 Kb and sequenced on an Illumina HiSeq 2500. The combined sequencing data was initially assembled at Dovetail Genomics using Meraculous and their in- house HiRise genome assembly algorithms to generate the A. marmoratus reference genome (AspMar1.0). Assessing assembly completeness In order to assess the completeness of the A. marmoratus reference genome, we used BUSCO (v. 3.0.1) (Simão et  al., 2015) with the vertebrate_od9 dataset containing 2586 genes, with default parameters apart from changing the BLAST cutoff from 1e- 3 to a more stringent value of 1e- 6. We used BUSCO numbers generated in Gao et al., 2017 for Shinisaurus crocodilurus and Alligator missis- sippiensis in Figure 2—figure supplement 1. To perform a phylogenetic analysis, 1333 shared ‘complete’ single-c opy orthologs were identi- fied in the genomes of green anole (Anolis carolinensis, anoCar2), cow (Bos taurus, ARS- UCD1.2), dog (Canis lupus familiaris, CanFam3.1), zebrafish (Danio rerio, danRer10), chicken (Gallus gallus, galGal5), human (Homo sapiens, GRCh38.p13), mouse (Mus musculus, GRCm38.p6), medaka (Oryzias latipes, oryLat2), rat (Rattus norvegicus, Rnor_6.0), Argentine black and white tegu (Salvator meri- anae, HLtupMer3), western clawed frog (Xenopus tropicalis, Xenopus_tropicalis_v9.1), platyfish (Xiphophorus maculatus, X_maculatus-5 .0-m ale). For each amino acid sequence, a multiple sequence alignment was performed with MAFFT (v. 7.305) (Katoh and Standley, 2013). The alignments were concatenated into a supermatrix of 1,112,277 amino acids. Phylogenetic tree topology was estimated using the Maximum Likelihood inference method using the pthreads version of RAxML (v. 8.2.11) and the PROTOGAMMAAUTO model for sequence evolution with 100 bootstrap replicates (Stamatakis, 2014). Repeat identification We quantified and annotated the repetitive DNA content within the A. marmoratus genome assembly by using the RepeatMasker pipeline on A. marmoratus scaffolds greater than or equal to 10 Kb in length. We first generated a de novo list of A. marmoratus repetitive elements using RepeatMod- eler (v. 1.0.11) (Smit and Hubley, 2008). We then used these as input into RepeatMasker (v. 4.0.9) (Smit et al., 2008) using the NCBI/RMBLAST (v. 2.6.0+) search engine. Unclassified repeat element consensus sequences from the RepeatModeler output for each of the three lizards (A. marmoratus, S. merianae, and A. carolinensis) were compared to each other by identifying reciprocal best hits using BLAST (v. 2.9.0+). Whole-genome sequencing, reference genome alignment, and heterozygosity determination Genomic DNA isolated from either liver or tail was prepared for sequencing using the KAPA HTP Library Preparation Kit (KK8234). Stock adapters were used from the Nextflex kit and barcodes were from BioScientific. All libraries were sequenced on the Illumina HiSeq 2500 platform. Whole-g enome sequencing data was aligned to the A. marmoratus reference genome with BWA (v. 0.7.15) (Li and Durbin, 2010) and marked for duplicates with Picard (RRID:SCR_006525; v. 1.119; https://broadinsti- tute.github.io/picard/). Because samples were sequenced over multiple lanes, the alignment files were merged subsequently, and another round of duplication marking was performed. The alignment files were realigned around small insertions and deletions with GATK (v. 3.5) (DePristo et al., 2011). Data Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 15 of 23 Research article Genetics and Genomics corresponding to lizard ID 122’s bam file was down-s ampled to 33% of its original size using seqtk (v. 1.2- r94) to match the expected average genome coverage of the other samples, as this animal was sequenced on one flow cell without multiplexing and, therefore, having much more sequenced reads. The per position nucleotide profiles for each alignment were then generated using a combination of pysam (RRID:SCR_021017, v. 0.12.0.1) (https://github.com/pysam-developers/pysam) and pysam- stats (Miles, 2015. v. 0.24.3) (https://github.com/alimanfoo/pysamstats) to determine the heterozy- gosity at any genomic position. Transcriptome assembly and genome annotation Two poly- A selected stranded RNA-s equencing libraries were generated with the TruSeq RNA Library Prep Kit v2 (RS- 122–2001) and sequenced on an Illumina HiSeq 2500 for the purpose of an A. marmoratus transcriptome assembly. The first library was derived from a blood sample taken from a male animal, and the second library was derived from an embryo incubated at 28 °C and harvested 47–51 days post- egg deposition. Trinity (v. 2.0.6) (Grabherr et al., 2011) was then used to generate an initial transcriptome assembly. The original reads were aligned to this transcriptome assembly using the Trinity companion script align_a nd_ estimate_ abundance. pl. Transcript isoforms with no read support were then filtered out and the remaining assembly was run through seqclean (https://sourceforge.net/projects/seqclean/). Evidence- based annotations for the transcriptome assembly were generated using the MAKER2 pipe- line (v. 2.31.8) (Holt and Yandell, 2011). For MAKER2, the entire UniProtKB/Swiss-P rot database of proteins (UniProt Consortium, 2018) was used and the Repbase data base was used to mask repeats within the MAKER2 framework (Bao et al., 2015). Assigning putative functions to the gene models was performed using BLAST (v. 2.6.0) and Interproscan (v. 5.13–52.0) (Jones et al., 2014). Copy number estimation for the Vomeronasal 2 receptor 26 (Vmn2r26) genes was based on aligning the mouse ortholog (http://www.orthodb.org), to the A. marmoratus reference assembly using Exon- erate (v. 2.4.0) (Slater and Birney, 2005) with a maximum intron size set to 20 Kb. Genes annotated as Vmn2r26 in the MAKER2 annotations were concatenated and aligned using Geneious (v. 10.1.3) (Geneious Prime, 2017) with default settings. The FastTree plugin (v. 1.0) was used to generate the phylogenetic tree from the alignment with default parameters. RAD-sequencing analysis Double digest RAD- sequencing data was derived primarily from previous studies (Barley et al., 2021; Barley et al., 2022b; Barley et al., 2019; Barley et al., 2022a). Fastq files were processed with Stacks (v. 2.62) using process_radtags and then samples from the same species were passed separately using denovo_ map. pl. The script executes all the components of the standard Stacks pipeline. The gstacks program of denovo_ map.p l calls variants for each position in each locus and assigns it either homozy- gous, heterozygous, or unknown. Output files were used to calculate the average coverage for each sample. Samples that have at least a coverage of 20 were considered for subsequent analysis. For each passing sample, the heterozygosity was calculated for sites at average coverage by adding up all heterozygous positions and dividing it by the total. Giemsa staining of erythrocytes Whole blood was collected aseptically from tails using acid citrate dextrose as an anticoagulant. All samples were prepared immediately after collection. A 5 µL aliquot of the diluted blood sample was used to prepare a blood smear on a 25 × 75 × 1 mm microscope slide. Once dry, the slide was placed in 95% ethanol for 5 min. Giemsa stain (0.4%; Sigma, GS500) was applied liberally to cover the slide every 4 min for a total incubation time of 16 min. The prepared slides were imaged on an Axioplan2 imaging microscope equipped with a plan-a pochromat 100 x /1.40 Oil objective and an Axiocam HRc (color) camera (Zeiss). Micromanager (v. 1.4) software was used to acquire the images. The acquired images were then scored visually for the number of haploid, diploid, and binucleated erythrocytes. Feulgen staining of erythrocytes Whole blood was collected aseptically from tails using acid citrate dextrose as an anticoagulant. All samples were prepared immediately after collection. A 5 µL aliquot of the diluted blood sample was used to prepare a blood smear on a 25 × 75 × 1 mm microscope slide. Blood smears were treated Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 16 of 23 Research article Genetics and Genomics with 10% neutral buffered formalin for 5 min at RT, then rinsed twice in distilled water. The slides were immersed in 5 M HCl for 30 min at RT, and then rinsed 2 times in distilled water. Slides were then immersed in Schiff’s reagent (Fisher Scientific #SS32- 500) at RT for 15–30 min until nuclei were stained. The slides were transferred directly to bisulfite water that was prepared by dissolving 2.5 g of potas- sium metabisulfite in 500 mL of water and adjusting the pH to 4.0 by the addition of concentrated HCl. The bisulfite wash was repeated three times with 10–15 sec of agitation. The slide was then washed under running tap water for 2 min and dehydrated by incubating in 70% EtOH for 5 min and then 95% EtOH for 5 min. The preparations were cleared in xylene before mounting. Flow cytometry of erythrocytes Blood collected from animals was treated as previously described for flow cytometry with modifications: ethanol fixation was performed after RNase treatment and propidium iodide staining was performed overnight, followed by sonication to disrupt aggregate cells (Lutes et al., 2011). A minimum number of 50,000 events were collected for each sample. Flowjo (treestar) was used for data analysis. Imaging of ovaries and germinal vesicle Germinal vesicle isolation and acquisition of image stacks by confocal microscopy were performed as described (Lutes et al., 2010). Ovary images were acquired with a Leica M205FA dissection micro- scope with a planar 0.63 X objective using Micromanager (v. 1.4) software. Acknowledgements We are grateful for the husbandry staff at the Stowers Institute for Medical Research (SIMR) espe- cially Rick Kupronis and the team of dedicated reptile technicians, David Jewell, Alex Muensch, Jillian Schieszer, Christina Piraquive, and Kristy Winter, for outstanding husbandry and herpetocultural skills. We thank the members of the New Mexico Department of Game and Fish and the Arizona Game and Fish Department for their support. We thank Veronica Cloud for help with RNA extractions and the molecular biology, flow cytometry, and advanced microscopy core facilities at SIMR for their excellent technical support and the members of the Emergent AI Center at Johannes Gutenberg University Mainz (JGU) for outstanding IT support and helpful discussions. We acknowledge computing time granted on the supercomputer MOGON II at JGU as part of NHR South-W est (n hrsw.d e), as well as the Institute of Molecular Biology Bioinformatics Core for computing resources. This work was funded in part by the Howard Hughes Medical Institute, SIMR, and an Alexander von Humboldt Professor- ship awarded to P.B. at JGU. This work also received support from the National Science Foundation grant DEB-1 754350 and the GenEvo RTG funded by the Deutsche Forschungsgemeinschaft (German Research Foundation) – GRK2526/1 – Projectnr. 40.7023052 and the Institute for Quantitative and Computational Biosciences (IQCB) in Mainz. We thank Charles (Jay) Cole, Uri García-V ázquez, Rick Kupronis, Daniel Lara-T ufiño, Norma Manríquez-M oran, Maximiliano Monroy, Priscilla Neaves, Adrián Nieto Montes de Oca, Harry Taylor, Robert Thomson, and Carol Townsend for assistance with field- work and specimen collection. We also thank the Natural History Museum of Los Angeles County, the UTEP Biodiversity Collections, Jay Cole and the American Museum of Natural History, the Museum of Vertebrate Zoology, Randy Klabacka and the Auburn University Museum of Natural History, and the LSU Museum of Natural Science for providing tissue loans in support of this research. Additional information Funding Funder Grant reference number Author Stowers Institute for Duncan Tormey Medical Research Howard Hughes Medical Peter Baumann Institute Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 17 of 23 Research article Genetics and Genomics Funder Grant reference number Author Alexander von Humboldt- Peter Baumann Stiftung Deutsche GRK2526/1 David V Ho Forschungsgemeinschaft National Science DEB-1754350 Anthony J Barley Foundation The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. Author contributions David V Ho, Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Writing – original draft, Writing – review and editing; Duncan Tormey, Data curation, Software, Investiga- tion, Methodology, Writing – original draft; Aaron Odell, Formal analysis, Validation, Investigation, Visualization, Methodology; Aracely A Newton, Robert R Schnittker, Formal analysis, Investigation, Methodology; Diana P Baumann, Resources, Data curation; William B Neaves, Investigation; Morgan R Schroeder, Rutendo F Sigauke, Software, Formal analysis, Investigation, Methodology; Anthony J Barley, Resources, Investigation, Writing – review and editing; Peter Baumann, Conceptualization, Resources, Supervision, Funding acquisition, Validation, Writing – original draft, Project administra- tion, Writing – review and editing Author ORCIDs David V Ho http://orcid.org/0000-0002-4709-269X Duncan Tormey https://orcid.org/0000-0002-0963-3424 Robert R Schnittker http://orcid.org/0000-0002-2950-8604 Anthony J Barley http://orcid.org/0000-0003-1675-6577 Peter Baumann http://orcid.org/0000-0003-4892-1485 Ethics Animals used in this study were produced in the AAALAC International accredited Stowers Reptile and Aquatics Facility in compliance with protocols approved by the Institutional Animal Care and Use Committee. They descended from breeding stock collected in New Mexico under permit numbers 3199 and 3395 and Arizona under license number SP564133. Tissues used in the RAD- seq analyses were derived from samples collected under permits from the Arizona Game and Fish Department, California Department of Fish and Wildlife, New Mexico Department of Game and Fish, and the Secretariá de Medio Ambiente y Recursos Naturales, Dirección General de Fauna Silvestre of Mexico. Decision letter and Author response Decision letter https://doi.org/10.7554/eLife.97035.sa1 Author response https://doi.org/10.7554/eLife.97035.sa2 Additional files Supplementary files • Supplementary file 1. A. marmoratus genome assembly statistics. • Supplementary file 2. Trinity assembly statistics. • Supplementary file 3. MAKER2 summary. • Supplementary file 4. Information for all A. marmoratus animals sequenced. • Supplementary file 5. Whole genome sequencing alignment statistics. • Supplementary file 6. Animals confirmed by microsatellite analysis to be of facultative parthenogenesis (FP) origin. • Supplementary file 7. Microsatellite primer information. • MDAR checklist Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 18 of 23 Research article Genetics and Genomics Data availability All raw sequencing data pertaining to the AspMar1.0 genome assembly are available at the National Center for Biotechnology Information under project accession number PRJNA360150. All other whole-g enome and RNA sequencing data can be found under PRJNA980964. RAD-s eq data can be located under PRJNA827355, PRJNA707030, PRJNA762930, and PRJNA1016487. Code and raw microsatellite data used for analysis are available at GitHub (copy archived at baumannlab, 2024). The following datasets were generated: Author(s) Year Dataset title Dataset URL Database and Identifier Baumann P 2020 Aspidoscelis marmoratus https://www. ncbi. nlm. NCBI BioProject, isolate:SIMRID8450 nih.g ov/b ioproject/ PRJNA360150 (Marbled whiptail) PRJNA360150 Ho DV, Tormey D, 2024 Sequencing of Aspidoscelis https://www. ncbi.n lm. NCBI BioProject, Odell A, Newton maramoratus nih. gov/ bioproject/ PRJNA980964 AA, Schnittker RR, PRJNA980964 Baumann DP, Neaves WB, Schroeder MR, Sigauke RF, Barley AJ, Baumann P Ho DV, Tormey D, 2024 Post-m eiotic mechanism of https://www.n cbi. nlm. NCBI BioProject, Odell A, Newton facultative parthenogenesis nih. gov/ bioproject/ PRJNA1016487 AA, Schnittker RR, in gonochoristic whiptail PRJNA1016487 Baumann DP, Neaves lizard species WB, Schroeder MR, Sigauke RF, Barley AJ, Baumann P The following previously published datasets were used: Author(s) Year Dataset title Dataset URL Database and Identifier Barley AJ 2021 Evolutionary study https://www.n cbi.n lm. NCBI BioProject, of whiptail lizards nih. gov/ bioproject/ PRJNA707030 (Aspidoscelis) PRJNA707030 Barley AJ 2021 Genetic diversity and the https://www. ncbi.n lm. NCBI BioProject, origins of parthenogenesis nih. gov/b ioproject/ PRJNA762930 in the teiid lizard PRJNA762930 Aspidoscelis laredoensis Barley AJ, de Oca 2022 The evolutionary network https://www. ncbi. nlm. NCBI BioProject, ANM, Normal L, of whiptail lizards reveals nih.g ov/b ioproject/ PRJNA827355 Robert CT predictable outcomes of PRJNA827355 hybridization References Adams L, Lyons K, Monday J, Larkin E, Wyffels J. 2023. Costs of parthenogenesis on growth and longevity in ex situ zebra sharks Stegostoma tigrinum. Endangered Species Research 50:81–91. DOI: https://doi.org/10.3354/ esr01224 Allen L, Sanders KL, Thomson VA. 2018. Molecular evidence for the first records of facultative parthenogenesis in elapid snakes. Royal Society Open Science 5:171901. DOI: https://doi.org/10.1098/rsos.171901, PMID: 29515892 Avise JC. 2015. Evolutionary perspectives on clonal reproduction in vertebrate animals. PNAS 112:8867–8873. DOI: https://doi.org/10.1073/pnas.1501820112, PMID: 26195735 Bao W, Kojima KK, Kohany O. 2015. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mobile DNA 6:11. DOI: https://doi.org/10.1186/s13100-015-0041-9, PMID: 26045719 Barley AJ, Nieto- Montes de Oca A, Reeder TW, Manríquez- Morán NL, Arenas Monroy JC, Hernández- Gallegos O, Thomson RC. 2019. Complex patterns of hybridization and introgression across evolutionary timescales in Mexican whiptail lizards (Aspidoscelis). Molecular Phylogenetics and Evolution 132:284–295. DOI: https://doi.org/10.1016/j.ympev.2018.12.016, PMID: 30562610 Barley AJ, Reeder TW, Nieto- Montes de Oca A, Cole CJ, Thomson RC. 2021. A new diploid parthenogenetic whiptail lizard from sonora, mexico, is the “missing link” in the evolutionary transition to polyploidy. The American Naturalist 198:295–309. DOI: https://doi.org/10.1086/715056, PMID: 34260872 Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 19 of 23 Research article Genetics and Genomics Barley AJ, Cordes JE, Walker JM, Thomson RC. 2022a. Genetic diversity and the origins of parthenogenesis in the teiid lizard Aspidoscelis laredoensis. Molecular Ecology 31:266–278. DOI: https://doi.org/10.1111/mec. 16213, PMID: 34614250 Barley AJ, Nieto- Montes de Oca A, Manríquez- Morán NL, Thomson RC. 2022b. The evolutionary network of whiptail lizards reveals predictable outcomes of hybridization. Science 377:773–777. DOI: https://doi.org/10. 1126/science.abn1593, PMID: 35951680 Bartelmez GW, Riddle O. 1924. On parthenogenetic cleavage and on the rôle of water absorption by the ovum in the formation of the subgerminal cavity in the pigeon’s egg. American Journal of Anatomy 33:57–66. DOI: https://doi.org/10.1002/aja.1000330104 baumannlab. 2024. Aspi_Fp. swh:1:rev:01a8d33519e700d701db67bb159fb5e693933327. Software Heritage. https://archive.softwareheritage.org/swh:1:dir:5cc3a8b9cefd8fb929893207512fb43194ee0a5f;origin=https:// github.com/baumannlab/Aspi_FP;visit=swh:1:snp:25c60a968c6382793a701367d2e6f8e39d0dc415;anchor= swh:1:rev:01a8d33519e700d701db67bb159fb5e693933327 Booth W, Johnson DH, Moore S, Schal C, Vargo EL. 2011. Evidence for viable, non-c lonal but fatherless Boa constrictors. Biology Letters 7:253–256. DOI: https://doi.org/10.1098/rsbl.2010.0793, PMID: 21047849 Booth W, Schuett GW. 2011. Molecular genetic evidence for alternative reproductive strategies in North American pitvipers (Serpentes: Viperidae): long-t erm sperm storage and facultative parthenogenesis. Biological Journal of the Linnean Society 104:934–942. DOI: https://doi.org/10.1111/j.1095-8312.2011.01782.x Booth W, Smith CF, Eskridge PH, Hoss SK, Mendelson JR, Schuett GW. 2012. Facultative parthenogenesis discovered in wild vertebrates. Biology Letters 8:983–985. DOI: https://doi.org/10.1098/rsbl.2012.0666, PMID: 22977071 Booth W, Schuett GW, Ridgway A, Buxton DW, Castoe TA, Bastone G, Bennett C, McMahan W. 2014. New insights on facultative parthenogenesis in pythons. Biological Journal of the Linnean Society 112:461–468. DOI: https://doi.org/10.1111/bij.12286 Booth W, Schuett GW. 2016. The emerging phylogenetic pattern of parthenogenesis in snakes. Biological Journal of the Linnean Society 118:172–186. DOI: https://doi.org/10.1111/bij.12744 Booth W, Levine BA, Corush JB, Davis MA, Dwyer Q, De Plecker R, Schuett GW. 2023. Discovery of facultative parthenogenesis in a new world crocodile. Biology Letters 19:20230129. DOI: https://doi.org/10.1098/rsbl. 2023.0129, PMID: 37282490 Brykczynska U, Tzika AC, Rodriguez I, Milinkovitch MC. 2013. Contrasted evolution of the vomeronasal receptor repertoires in mammals and squamate reptiles. Genome Biology and Evolution 5:389–401. DOI: https://doi. org/10.1093/gbe/evt013, PMID: 23348039 Busch JW, Delph LF. 2012. The relative importance of reproductive assurance and automatic selection as hypotheses for the evolution of self-f ertilization. Annals of Botany 109:553–562. DOI: https://doi.org/10.1093/ aob/mcr219, PMID: 21937484 Card DC, Vonk FJ, Smalbrugge S, Casewell NR, Wüster W, Castoe TA, Schuett GW, Booth W. 2021. Genome- wide data implicate terminal fusion automixis in king cobra facultative parthenogenesis. Scientific Reports 11:7271. DOI: https://doi.org/10.1038/s41598-021-86373-1, PMID: 33790309 Cassar G, John TM, Etches RJ. 1998. Observations on ploidy of cells and on reproductive performance in parthenogenetic turkeys. Poultry Science 77:1457–1462. DOI: https://doi.org/10.1093/ps/77.10.1457, PMID: 9776050 Chapman DD, Shivji MS, Louis E, Sommer J, Fletcher H, Prodöhl PA. 2007. Virgin birth in a hammerhead shark. Biology Letters 3:425–427. DOI: https://doi.org/10.1098/rsbl.2007.0189, PMID: 17519185 Chapman DD, Firchau B, Shivji MS. 2008. Parthenogenesis in a large‐bodied requiem shark, the blacktip Carcharhinus limbatus. Journal of Fish Biology 73:1473–1477. DOI: https://doi.org/10.1111/j.1095-8649.2008. 02018.x Chapman JA, Ho I, Sunkara S, Luo S, Schroth GP, Rokhsar DS. 2011. Meraculous: de novo genome assembly with short paired- end reads. PLOS ONE 6:e23501. DOI: https://doi.org/10.1371/journal.pone.0023501, PMID: 21876754 Cole CJ, Lowe CH, Wright JW. 1969. Sex Chromosomes in Teiid Whiptail Lizards (Genus Cnemidophorus). American Museum of Natural History. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ. 2011. A framework for variation discovery and genotyping using next-g eneration DNA sequencing data. Nature Genetics 43:491–498. DOI: https://doi.org/10.1038/ng.806 Dubach J, Sajewicz A, Pawley R. 1997. Parthenogenesis in the Arafuran filesnake (Acrochordus arafurae). Herpetological Natural History 5:11–18. Dudgeon CL, Coulton L, Bone R, Ovenden JR, Thomas S. 2017. Switch from sexual to parthenogenetic reproduction in a zebra shark. Scientific Reports 7:40537. DOI: https://doi.org/10.1038/srep40537, PMID: 28091617 Fang X, Mou Y, Huang Z, Li Y, Han L, Zhang Y, Feng Y, Chen Y, Jiang X, Zhao W, Sun X, Xiong Z, Yang L, Liu H, Fan D, Mao L, Ren L, Liu C, Wang J, Li K, et al. 2012. The sequence and analysis of a Chinese pig genome. GigaScience 1:16. DOI: https://doi.org/10.1186/2047-217X-1-16, PMID: 23587058 Feldheim KA, Chapman DD, Sweet D, Fitzpatrick S, Prodöhl PA, Shivji MS, Snowden B. 2010. Shark virgin birth produces multiple, viable offspring. The Journal of Heredity 101:374–377. DOI: https://doi.org/10.1093/ jhered/esp129, PMID: 20106913 Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 20 of 23 Research article Genetics and Genomics Feldheim KA, Dubach J, Watson L. 2023. Parthenogenesis in an elasmobranch in the presence of conspecific males. Journal of Fish Biology 102:525–527. DOI: https://doi.org/10.1111/jfb.15268, PMID: 36369968 Fields AT, Feldheim KA, Poulakis GR, Chapman DD. 2015. Facultative parthenogenesis in a critically endangered wild vertebrate. Current Biology 25:R446–R447. DOI: https://doi.org/10.1016/j.cub.2015.04.018, PMID: 26035783 Gao J, Li Q, Wang Z, Zhou Y, Martelli P, Li F, Xiong Z, Wang J, Yang H, Zhang G. 2017. Sequencing, de novo assembling, and annotating the genome of the endangered Chinese crocodile lizard Shinisaurus crocodilurus. GigaScience 6:gix041. DOI: https://doi.org/10.1093/gigascience/gix041 Geneious Prime. 2017. Geneious prime. 0.1. Graphstats. https://www.graphstats.net/geneious-prime Germano DJ, Smith PT. 2010. Molecular evidence for parthenogenesis in the sierra garter snake, thamnophis couchii (Colubridae). The Southwestern Naturalist 55:280–282. DOI: https://doi.org/10.1894/WL-29.1 Goodwillie C, Kalisz S, Eckert CG. 2005. The evolutionary enigma of mixed mating systems in plants: occurrence, theoretical explanations, and empirical evidence. Annual Review of Ecology, Evolution, and Systematics 36:47–79. DOI: https://doi.org/10.1146/annurev.ecolsys.36.091704.175539 Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad- Toh K, Friedman N, et al. 2011. Full-l ength transcriptome assembly from RNA- Seq data without a reference genome. Nature Biotechnology 29:644–652. DOI: https://doi.org/10.1038/nbt.1883, PMID: 21572440 Groot TVM, Bruins E, Breeuwer JAJ. 2003. Molecular genetic evidence for parthenogenesis in the Burmese python, Python molurus bivittatus. Heredity 90:130–135. DOI: https://doi.org/10.1038/sj.hdy.6800210, PMID: 12634818 Holt WV, Lloyd RE. 2010. Sperm storage in the vertebrate female reproductive tract: How does it work so well? Theriogenology 73:713–722. DOI: https://doi.org/10.1016/j.theriogenology.2009.07.002 Holt C, Yandell M. 2011. MAKER2: an annotation pipeline and genome-d atabase management tool for second- generation genome projects. BMC Bioinformatics 12:491. DOI: https://doi.org/10.1186/1471-2105-12-491, PMID: 22192575 Houck LD. 2009. Pheromone communication in amphibians and reptiles. Annual Review of Physiology 71:161– 176. DOI: https://doi.org/10.1146/annurev.physiol.010908.163134, PMID: 18847365 Ito M, Kaneko- Ishino T, Ishino F, Matsuhashi M, Yokoyama M, Katsuki M. 1991. Fate of haploid parthenogenetic cells in mouse chimeras during development. The Journal of Experimental Zoology 257:178–183. DOI: https:// doi.org/10.1002/jez.1402570206, PMID: 1990050 Jones P, Binns D, Chang H- Y, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, Pesseat S, Quinn AF, Sangrador- Vegas A, Scheremetjew M, Yong S- Y, Lopez R, Hunter S. 2014. InterProScan 5: genome- scale protein function classification. Bioinformatics 30:1236–1240. DOI: https://doi.org/10.1093/bioinformatics/ btu031 Kajitani R, Toshimoto K, Noguchi H, Toyoda A, Ogura Y, Okuno M, Yabana M, Harada M, Nagayasu E, Maruyama H, Kohara Y, Fujiyama A, Hayashi T, Itoh T. 2014. Efficient de novo assembly of highly heterozygous genomes from whole- genome shotgun short reads. Genome Research 24:1384–1395. DOI: https://doi.org/10. 1101/gr.170720.113, PMID: 24755901 Katoh K, Standley DM. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology and Evolution 30:772–780. DOI: https://doi.org/10.1093/molbev/ mst010 Kinney ME, Wack RF, Grahn RA, Lyons L. 2013. Parthenogenesis in a Brazilian rainbow boa (Epicrates cenchria cenchria). Zoo Biology 32:172–176. DOI: https://doi.org/10.1002/zoo.21050, PMID: 23086743 Kratochvíl L, Vukić J, Červenka J, Kubička L, Johnson Pokorná M, Kukačková D, Rovatsos M, Piálek L. 2020. Mixed-s ex offspring produced via cryptic parthenogenesis in a lizard. Molecular Ecology 29:4118–4127. DOI: https://doi.org/10.1111/mec.15617, PMID: 32881125 Kroeger PT, Poureetezadi SJ, McKee R, Jou J, Miceli R, Wingert RA. 2014. Production of haploid zebrafish embryos by in vitro fertilization. Journal of Visualized Experiments 01:51708. DOI: https://doi.org/10.3791/ 51708, PMID: 25046024 Kuraku S, Meyer A. 2009. The evolution and maintenance of Hox gene clusters in vertebrates and the teleost- specific genome duplication. The International Journal of Developmental Biology 53:765–773. DOI: https://doi. org/10.1387/ijdb.072533km, PMID: 19557682 Lampert KP, Steinlein C, Schmid M, Fischer P, Schartl M. 2007. A haploid- diploid-t riploid mosaic of the Amazon molly, Poecilia formosa. Cytogenetic and Genome Research 119:131–134. DOI: https://doi.org/10.1159/ 000109629, PMID: 18160792 Larose C, Lavanchy G, Freitas S, Parker DJ, Schwander T. 2023. Facultative parthenogenesis: a transient state in transitions between sex and obligate asexuality in stick insects? Peer Community Journal 3:283. DOI: https:// doi.org/10.24072/pcjournal.283 Le Peuch C j, Picard A, Dorée M. 1985. Parthenogenetic activation decreases the polyphosphoinositide content of frog eggs. FEBS Letters 187:61–64. DOI: https://doi.org/10.1016/0014-5793(85)81214-x, PMID: 2991011 Lenk P, Eidenmueller B, Staudter H, Wicker R, Wink M. 2005. A parthenogenetic Varanus. Amphibia-R eptilia 26:507–514. DOI: https://doi.org/10.1163/156853805774806296 Levine BA, Schuett GW, Booth W. 2021. Exceptional long-t erm sperm storage by a female vertebrate. PLOS ONE 16:e0252049. DOI: https://doi.org/10.1371/journal.pone.0252049, PMID: 34086677 Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 21 of 23 Research article Genetics and Genomics Levine BA, Moresco A, Trout T, Schuett GW, Booth W. 2024. Female long- term sperm storage results in viable offspring in the Himalayan Mountain Pitviper, Ovophis monticola. Zoo Biology 43:183–187. DOI: https://doi. org/10.1002/zoo.21820, PMID: 38234126 Li H, Durbin R. 2010. Fast and accurate long- read alignment with Burrows–Wheeler transform. Bioinformatics 26:589–595. DOI: https://doi.org/10.1093/bioinformatics/btp698 Li H. 2016. seqtk. v1.2- r94. GitHub. https://github.com/lh3/seqtk Lowe CH, Wright JW. 1966. Evolution of Parthenogenetic Species of Cnemidophorus (Whiptail Lizards) in Western North America. Journal of the Arizona Academy of Science 4:81–87. Lutes AA, Neaves WB, Baumann DP, Wiegraebe W, Baumann P. 2010. Sister chromosome pairing maintains heterozygosity in parthenogenetic lizards. Nature 464:283–286. DOI: https://doi.org/10.1038/nature08818, PMID: 20173738 Lutes AA, Baumann DP, Neaves WB, Baumann P. 2011. Laboratory synthesis of an independently reproducing vertebrate species. PNAS 108:9910–9915. DOI: https://doi.org/10.1073/pnas.1102811108, PMID: 21543715 Miles A. 2015. pysamstats. v0.24.3. GitHub. https://github.com/alimanfoo/pysamstats Newton AA, Schnittker RR, Yu Z, Munday SS, Baumann DP, Neaves WB, Baumann P. 2016. Widespread failure to complete meiosis does not impair fecundity in parthenogenetic whiptail lizards. Development 143:4486–4494. DOI: https://doi.org/10.1242/dev.141283, PMID: 27802173 Olsen WW, Marsden SJ. 1954. Natural parthenogenesis in turkey eggs. Science 120:545–546. DOI: https://doi. org/10.1126/science.120.3118.545, PMID: 13195695 Olsen MW, Wilson SP, Marks HL. 1968. Genetic control of parthenogenesis in chickens. The Journal of Heredity 59:41–42. DOI: https://doi.org/10.1093/oxfordjournals.jhered.a107639, PMID: 5656917 Olsen MW. 1973. Longevity and organ weights of Beltsville small white parthenogens and normal turkey males. Poultry Science 52:666–670. DOI: https://doi.org/10.3382/ps.0520666, PMID: 4709417 MW Olsen, Avian parthenogensis. Agricultural Research Service USDA, ARS- NE 65, 1–82, (1975). Parker HM, Kiess AS, Wells JB, Young KM, Rowe D, McDaniel CD. 2010. Genetic selection increases parthenogenesis in Chinese painted quail (Coturnix chinensis). Poultry Science 89:1468–1472. DOI: https://doi. org/10.3382/ps.2009-00388 Pijnacker LP. 1969. Automictic parthenogenesis in the stick insectBacillus rossius Rossi (Cheleutoptera, phasmidae). Genetica 40:393–399. DOI: https://doi.org/10.1007/BF01787364 Putnam NH, O’Connell BL, Stites JC, Rice BJ, Blanchette M, Calef R, Troll CJ, Fields A, Hartley PD, Sugnet CW, Haussler D, Rokhsar DS, Green RE. 2016. Chromosome- scale shotgun assembly using an in vitro method for long- range linkage. Genome Research 26:342–350. DOI: https://doi.org/10.1101/gr.193474.115, PMID: 26848124 pysam-d evelopers. 2017. pysam. v0.12.0.1. GitHub. https://github.com/pysam-developers/pysam Ramachandran R, McDaniel CD. 2018. Parthenogenesis in birds: a review. Reproduction 155:R245–R257. DOI: https://doi.org/10.1530/REP-17-0728, PMID: 29559496 Reeder TW, Cole CJ, Dessauer HC. 2002. Phylogenetic relationships of whiptail lizards of the genus cnemidophorus (Squamata: Teiidae): a test of monophyly, reevaluation of karyotypic evolution, and review of hybrid origins. American Museum Novitates 3365:1–61. DOI: https://doi.org/10.1206/0003-0082(2002) 365<0001:PROWLO>2.0.CO;2 Reynolds RG, Booth W, Schuett GW, Fitzpatrick BM, Burghardt GM. 2012. Successive virgin births of viable male progeny in the checkered gartersnake, Thamnophis marcianus. Biological Journal of the Linnean Society 107:566–572. DOI: https://doi.org/10.1111/j.1095-8312.2012.01954.x Rivera- Colón AG, Catchen J. 2022. Population genomics analysis with RAD, reprised: stacks 2. Methods in Molecular Biology 2498:99–149. DOI: https://doi.org/10.1007/978-1-0716-2313-8_7, PMID: 35727543 Ryba NJ, Tirindelli R. 1997. A new multigene family of putative pheromone receptors. Neuron 19:371–379. DOI: https://doi.org/10.1016/s0896-6273(00)80946-0, PMID: 9292726 Ryder OA, Thomas S, Judson JM, Romanov MN, Dandekar S, Papp JC, Sidak-L oftis LC, Walker K, Stalis IH, Mace M, Steiner CC, Chemnick LG. 2021. Facultative Parthenogenesis in California Condors. The Journal of Heredity 112:569–574. DOI: https://doi.org/10.1093/jhered/esab052, PMID: 34718632 Sagata N. 1996. Meiotic metaphase arrest in animal oocytes: its mechanisms and biological significance. Trends in Cell Biology 6:22–28. DOI: https://doi.org/10.1016/0962-8924(96)81034-8, PMID: 15157528 Sarvella P. 1973. Adult parthenogenetic chickens. Nature 243:171. DOI: https://doi.org/10.1038/243171a0, PMID: 4706288 Schut E, Hemmings N, Birkhead TR. 2008. Parthenogenesis in a passerine bird, the Zebra Finch Taeniopygia guttata. Ibis 150:197–199. DOI: https://doi.org/10.1111/j.1474-919X.2007.00755.x Sever DM, Hamlett WC. 2002. Female sperm storage in reptiles. The Journal of Experimental Zoology 292:187– 199. DOI: https://doi.org/10.1002/jez.1154, PMID: 11754034 Shi P, Zhang J. 2007. Comparative genomic analysis identifies an evolutionary shift of vomeronasal receptor gene repertoires in the vertebrate transition from water to land. Genome Research 17:166–174. DOI: https:// doi.org/10.1101/gr.6040007, PMID: 17210926 Shibata H, Sakata S, Hirano Y, Nitasaka E, Sakabe A. 2017. Facultative parthenogenesis validated by DNA analyses in the green anaconda (Eunectes murinus). PLOS ONE 12:e0189654. DOI: https://doi.org/10.1371/ journal.pone.0189654, PMID: 29236745 Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. 2015. BUSCO: assessing genome assembly and annotation completeness with single- copy orthologs. Bioinformatics 31:3210–3212. DOI: https://doi.org/ 10.1093/bioinformatics/btv351, PMID: 26059717 Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 22 of 23 Research article Genetics and Genomics Slater GSC, Birney E. 2005. Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics 6:31. DOI: https://doi.org/10.1186/1471-2105-6-31, PMID: 15713233 Smit A, Hubley R. 2008. RepeatMasker Open- 4.0. https://www.repeatmasker.org/ [Accessed September 18, 2019]. Smit A, Hubley R, Green R. 2008. Repeatmasker Open- 4.0. 4.0.9. ISB. https://www.repeatmasker.org/ Sperling AL, Fabian DK, Garrison E, Glover DM. 2023. A genetic basis for facultative parthenogenesis in Drosophila. Current Biology 33:3545–3560.. DOI: https://doi.org/10.1016/j.cub.2023.07.006, PMID: 37516115 Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-a nalysis of large phylogenies. Bioinformatics 30:1312–1313. DOI: https://doi.org/10.1093/bioinformatics/btu033 Straube N, Lampert KP, Geiger MF, Weiß JD, Kirchhauser JX. 2016. First record of second- generation facultative parthenogenesis in a vertebrate species, the whitespotted bambooshark Chiloscyllium plagiosum. Journal of Fish Biology 88:668–675. DOI: https://doi.org/10.1111/jfb.12862, PMID: 26727105 Streisinger G, Walker C, Dower N, Knauber D, Singer F. 1981. Production of clones of homozygous diploid zebra fish (Brachydanio rerio). Nature 291:293–296. DOI: https://doi.org/10.1038/291293a0, PMID: 7248006 Su CY, Menuz K, Carlson JR. 2009. Olfactory perception: receptors, cells, and circuits. Cell 139:45–59. DOI: https://doi.org/10.1016/j.cell.2009.09.015, PMID: 19804753 Tanaka M, Yamaha E, Arai K. 2004. Survival capacity of haploid-d iploid goldfish chimeras. Journal of Experimental Zoology. Part A, Comparative Experimental Biology 301:491–501. DOI: https://doi.org/10.1002/ jez.a.48, PMID: 15181643 trinityrnaseq. 2015. trinityrnaseq. v2.0.6. GitHub. https://github.com/trinityrnaseq/trinityrnaseq UniProt Consortium T. 2018. UniProt: the universal protein knowledgebase. Nucleic Acids Research 46:2699. DOI: https://doi.org/10.1093/nar/gky092 Vanzolini PE. 1993. Reptilians Too: Biology of Whiptail Lizards. (Genus Cnemidophorus.) John W. Wright and Laurie J. Vitt, Eds. University of Oklahoma and Oklahoma Museum of Natural History, Norman, 1993. xvi, 417 pp., illus. $29. Herpetologists’ League Special Publication no. 3. Based on a symposium, Norman, OK, Aug. 1984. Science 262:1282–1283. DOI: https://doi.org/10.1126/science.262.5137.1282.b Walker JM, Abuhteba RM, Cordes JE. 1991. Morphological and experimental verification of hybridization between all- female cnemidophorus laredoensis b and gonochoristic cnemidophorus gularis (Squamata: Teiidae). Herpetologica 47:152–164. DOI: https://doi.org/10.2307/3671928 Watts PC, Buley KR, Sanderson S, Boardman W, Ciofi C, Gibson R. 2006. Parthenogenesis in Komodo dragons. Nature 444:1021–1022. DOI: https://doi.org/10.1038/4441021a, PMID: 17183308 Whitehead MR, Lanfear R, Mitchell RJ, Karron JD. 2018. Plant mating systems often vary widely among populations. Frontiers in Ecology and Evolution 6:00038. DOI: https://doi.org/10.3389/fevo.2018.00038 Wolf DP. 1974. The cortical response in Xenopus laevis ova. Developmental Biology 40:102–115. DOI: https:// doi.org/10.1016/0012-1606(74)90112-2, PMID: 4213423 Zhang L, Huang Y, Wang M, Guo Y, Liang J, Yang X, Qi W, Wu Y, Si J, Zhu S, Li Z, Li R, Shi C, Wang S, Zhang Q, Tang Z, Wang L, Li K, Fei J- F, Lan G. 2019. Development and genome sequencing of a laboratory-i nbred miniature pig facilitates study of human diabetic disease. iScience 19:162–176. DOI: https://doi.org/10.1016/j. isci.2019.07.025 Ho, Tormey et al. eLife 2024;0:e97035. DOI: https://doi.org/10.7554/eLife.97035 23 of 23