Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2017; 12(12): e0189654.
Published online 2017 December 13. doi:  10.1371/journal.pone.0189654
PMCID: PMC5728508

Facultative parthenogenesis validated by DNA analyses in the green anaconda (Eunectes murinus)

Hiroki Shibata, Conceptualization, Data curation, Investigation, Methodology, Project administration, Supervision, Validation, Writing – original draft,1,* Shuichi Sakata, Investigation, Resources, Writing – review & editing,2 Yuzo Hirano, Investigation, Resources,2 Eiji Nitasaka, Methodology, Writing – review & editing,3 and Ai Sakabe, Investigation, Resources, Writing – review & editing2
Wolfgang Arthofer, Editor


In reptiles, the mode of reproduction is typically sexual. However, facultative parthenogenesis occurs in some Squamata, such as Komodo dragon (Varanus komodoensis) and Burmese python (Python bivittatus). Here, we report facultative parthenogenesis in the green anaconda (Eunectes murinus). We found two fully developed female neonates and 17 undeveloped eggs in the oviduct of a female anaconda isolated from other individuals for eight years and two months at Ueno Zoo, Japan. To clarify the zygosity of the neonates, we analyzed 18 microsatellite markers of which 16 were informative. We observed only maternal alleles and no paternal alleles for all 16 markers. To examine the possibility of the long-term sperm storage, we estimated allele frequencies in a putative parental stock by genotyping five unrelated founders. If all founders, including the mother, are originated from a single Mendelian population, then the probability that the neonates were produced by sexual reproduction with an unrelated male via long-term sperm storage was infinitesimally small (2.31E-32 per clutch). We also examined samples from two additional offspring that the mother delivered eight years before her death. We consistently observed paternal alleles in these elder offspring, indicating that the mother had switched from sexual reproduction to asexual reproduction during the eight years of isolation. This is the first case of parthenogenesis in Eunectes to be validated by DNA analysis, and suggests that facultative parthenogenesis is widespread in the Boidae.


Among vertebrates, parthenogenesis, the occurrence of unisexual lines whereby reproduction occurs without any involvement of males or their sperm, is most common amongst reptiles. Obligate parthenogenesis (OP) has been reported in more than 20 species of lizards and one species of snake (Indotyphlops braminus) [1]. Facultative parthenogenesis (FP), occasional occurrence of parthenogenesis, in individuals of a species that normally reproduce sexually occurs in at least five families—Boidae, Pythonidae, Viperidae, Acrochordidae and Colubridae [1].

Two modes of FP have been proposed. Systematic facultative parthenogenesis of some species of insect and other invertebrates is thought to be highly adaptive owing to the high hatching success of unfertilized eggs by automixis [2]. The other is accidental FP, referred as Tychoparthenogenesis [2], observed in some insects and vertebrates. The associated very low hatching rate suggests that this form of parthenogenesis is some form of reproduction error [3]. However, accidental FP can be also interpreted as an emergency reproduction after the long period of isolation from mates. The low hatching rate can be explained by effects of multiple deleterious variations in the parthenogens, that are maintained as heterozygotes in sexually reproducing populations.

Two types of accidental FP occur in snakes. Type A, reported from Alethinophidia, Boidae and Pythonidae [48], is characterized by exclusively female offspring. Type B, reported from Caenophidia [915], is characterized by exclusively male offspring. Until recently, all snake species were thought to be heterogametic in females (ZW), so a number of mechanisms for the production of homogametic sons (ZZ) by Type B parthenogenesis appeared plausible. In contrast, mechanisms to produce parthenogenetic daughters by Type A parthenogenesis involved proposing WW individuals, unlikely to be viable except in rare cases. Recently, XY sex chromosomes have been discovered in Boidae and Pythonidae [16] suggesting that the two different parthenogenetic types in snakes simply correspond to the two different modes of sex-determination. The suggestion is that Type A occurs in species with heterogametic males (XY) to yield XX parthenogenetic females, whereas Type B occurs in species with heterogametic females (ZW) to yield ZZ parthenogenetic males.

Three species of Boidae exhibit facultative parthenogenesis that has been validated by DNA evidence, namely, Boa constrictor [5], Epicrates maurus [6] and Epicrates cenchria cenchria [7]. Here we report a case of parthenogenetic reproduction in a long-term isolated female green anaconda Eunectes murinus, validated by DNA evidence using microsatellite markers. Although parthenogenesis has been reported before in E. murinus [17], our findings are the first to definitively report facultative parthenogenesis in Eunectes based on DNA analysis.

Materials and methods

1. Biological samples

A female wild specimen of the green anaconda, the focal mother (Emu-01, Table 1) was caught in Guyana and brought to Ueno Zoo, Tokyo, Japan, on September 7, 2007. One month later, on October 16, she gave birth to 19 neonates, two of which we retained (2007OS-1 (Emu-04) and 2007OS-2 (Emu-05), Table 1). She was kept isolated from males during captivity, and on November 8, 2015, she died of pneumonia. Autopsy revealed two developed fetuses (2015OS-1 (Emu-02) and 2015OS-2 (Emu-03), Table 1) in her oviduct and 17 undeveloped eggs (Fig 1a). Both fetuses were almost fully developed but dead (Fig 1b and 1c). Examination of gonads by dissection confirmed both to be female.

Table 1
Eunectes murinus samples analyzed.
Fig 1
Neonates found in the oviduct of the focal mother (Emu-01) that had been isolated from other snakes for seven years.

Genomic DNA was extracted from muscle tissue of each of the focal mother and the four offspring (Emu-01 to Emu-05,Table 1) using a QIAamp DNA Mini Kit (Qiagen, Hilden, Germany). Genomic DNA was also extracted from shed skin of additional five unrelated specimens (Emu-06 to Emu-10, Table 1) following the protocols of Fetzner (1999) [18] with slight modification.

2. Microsatellite markers

We constructed Illumina NGS library for the focal mother (Emu-01, Table 1) using an Illumina TruSeq DNA PCR-Free LT Sample Prep Kit (Illumina, Tokyo, Japan). We performed 351-bp single read sequencing using MiSeq (Illumina). Filtering the raw reads using PRINSEQ [19] yielded 6,069,144 clean reads with a modal value of 34 for sequence quality. Using a custom perl script, we extracted 15,225 sequences harboring trinucleotide motifs with more than 10 repeat units as candidate sequences for microsatellite markers. We selected 18 of these sequences for further analyses (Genbank ID LC179548-LC179565) (Table 2).

Table 2
Characterization of 18 microsatellite DNA loci designed for Eunectes murinus.

3. Genotyping

We designed PCR primers for the candidate marker region, and added M13 sequence (5’-cacgacgttgtaaaacgac-3’) to the 5’ end of one of the primers for the fluorescent labeling of PCR products. After the initial PCR reaction, we performed second-round PCR for labeling with ROX-labeled M13 primer (5’-cacgacgttgtaaaacgac-3’). The labeled PCR products were examined by fragment analyses on an ABI3130 Sequencer using GeneMapper software version 4.1 (Thermo Fisher Scientific). After the examination of the 18 candidate markers, we could optimize genotyping condition for 16 markers. We genotyped the all ten specimens.

4. Probability calculation

We assumed that the unrelated founders including the focal mother (Emu-01, Emu-06 to Emu-10) are from a Mendelian population in Hardy-Weinberg equilibrium. Dividing the observed number of alleles by number of haplotypes (2n = 12) yielded estimates of allele frequencies for each marker (S1 Table). Under our null hypothesis, the genotype of a putative father of the two neonates is drawn from this allelic distribution.

Then our observed genotypes for the two neonates arise from one of the three cases where p1 and p2 are the frequencies of maternal alleles:

  • Case 1: Heterozygous mother transmitted different alleles to two daughters.
    In this case, the putative father must have had identical heterozygous genotype to that of the focal mother. The frequency of this common genotype in the population is 2p1p2. The probability of the same alleles arising from the parents with the identical heterozygous genotype in each neonate is (1/2)2. Therefore, the probability (Pk) of having two daughters homozygous for different maternal alleles is
    Eight markers, EmuTri001-005, 007, 011 and 015, were consistent with this case.
  • Case 2: Heterozygous mother transmitted the same allele to both daughters.
    In this case, the putative father must be a heterozygote or homozygote for the relevant allele. The frequencies of homozygous and heterozygous males are p12 and 2p1(1 − p1), respectively. The probability that the homozygous father transmits the allele is 1. The probability that the heterozygous father transmits the relevant allele to both daughters is (1/2)2 = 1/4. The probability that the heterozygous mother transmits the relevant allele twice is also 1/4. Therefore, the probability (Pk) of having two daughters homozygous for the same maternal allele in Case 2 is
    Six markers, EmuTri006, 009, 014, 016-018, were consistent with this case.
  • Case 3: Homozygous mother transmitted the same allele to both daughters.
    In this case the putative father must be heterozygote or homozygote for the relevant allele as in Case 2. The frequencies of homozygous and heterozygous males are p12 and 2p1(1 − p1), respectively. The probability that the focal mother transmits the relevant allele is 1. The probability that the homozygous father transmits the allele is also 1. The probability that the heterozygous father transmits the relevant allele twice is 1/4. Therefore, the probability (Pk) of having two daughters homozygous for the same maternal allele in Case 3 is
    Two markers, EmuTri008 and 013 were consistent with this case.

The joint probability is the product of Pk for the 16 makers assuming that they are independent.


We summarized the genotyping results in Table 3. We found the focal mother was heterozygous for 14 out of the 16 microsatellite markers. The two female neonates consistently had only one allele inherited from the focal mother. There were no non-maternal alleles observed in either of the neonates. The lack of paternal alleles with apparent homozygosity of maternal alleles strongly suggests that the neonates were produced by parthenogenesis by the focal mother.

Table 3
Genotypes of the mother, four offspring and five unrelated Eunectes murinus.

The joint probability that two homozygous neonates were produced by sexual reproduction by the focal mother and the putative father, under the assumption that they were drawn from the population characterized by the five unknowns and the mother, and of linkage equilibrium for the markers, was infinitesimally small, 2.31E-32 (Table 4). We therefore excluded the possibility that the two neonates arose through sexual reproduction via long-term sperm storage and confirm our conclusion that the female neonates are parthenogenetic offspring. In contrast, paternal and maternal alleles were present in elder offspring born in 2007 with the same mother– 2007-OS1 (Emu-04) and 2007-OS2 (Emu-05). The focal mother has thus reproduced sexually at least once in the past, and has switched to parthenogenetic reproduction during the 8 years of her isolation.

Table 4
Allele frequencies and probablities estimated from six unrelated Eunectes murinus.


Here, we report the first DNA-validated case of facultative parthenogenesis in the green anaconda, Eunectes murinus. Since facultative parthenogenesis producing only females has been reported in other boid species such as Boa constrictor [5], Epicrates maurus [6] and Epicrates cenchria [7], Type A facultative parthenogenesis is to date the most common form, perhaps the only form, of parthenogenesis in Boidae. Female-biased facultative parthenogenesis has also been reported in the Pythonidae–Python bivittatus [4], Python regius and Malayopython reticulatus [8]–suggesting that it is the ancestral state shared in basal Alethinophidia. Recently, XY sex chromosomes have been discovered in Boa imperator and Python bivittatus [16], and are likely to occur broadly in the Pythonidae and Boidae, including Eunectes murinus. Female offspring produced by Type A parthenogenesis with XX or XO offspring, is now the expectation for this group of snakes subject to further investigation.

Three possible mechanisms for parthenogenesis yielding diploid offspring are recognized: central fusion, terminal fusion and gametic duplication. In species with homogametic females (XX), such as Eunectes murinus, central fusion is expected to produce female offspring only, with the same level of heterozygosity as seen in our focal mother. Since no such case has been ever reported, including the present study, central fusion is highly unlikely to be the mechanism of parthenogenesis in snakes. Terminal fusion produces female-only offspring that are largely homozygous (autozygous), but partially heterozygous owing to recombination during maternal meiosis. Gametic duplication produces female-only offspring that are completely autozygous for maternal chromosomes, even in the region where recombination occurs during maternal meiosis. Limitations of the informative markers used in previous studies of parthenogenesis in Boidae did not allow for a clear distinction between terminal fusion and gametic duplication (Booth et al 2011a using 4 informative markers [5]; Booth et al 2011b using 4 informative markers [6]; Kinney et al 2013 using 7 markers [7]). In the present study, the female offspring were consistently homozygous for all of 16 microsatellite markers, including 14 markers for which the focal mother was heterozygous (Table 3). The exclusive homozygosity (14/14) we observed strongly indicates that the parthenogenesis in Eunectes murinus is produced by gametic duplication, not terminal fusion. This hypothesis could be tested with high-resolution genome-wide analyses such as whole genome sequencing or reduced representational approaches such as ddRAD [20] or DArTSeq [21].

Although highly unlikely for vertebrates, there is a possibility that the parthenogenetic offspring of Eunectes murinus are haploid (XO), developed from unfertilized eggs. Such parthenogenesis by reduction of ploidy has been reported from vertebrates only in whitetip reef shark Triaenodon obesus [22]. To exclude this possibility, it is necessary to examine karyotypes by cytological analyses or by quantitation of nuclear DNA. Unfortunately the parthenogenetic neonates described in the current study were dead upon the autopsy of their mother; we were unable to obtain live cells for cytological analyses or FACS analyses, and so cannot eliminate the slight possibility of haploid offspring.

Since the two parthenogenic offspring were female and dead, there is no direct information on their viability. However, they had no obvious morphological anomalies. A recent suspected case of parthenogenesis in Eunectes murinus has produced viable offspring, but uncertainty remains on the status of this case as no DNA evidence in support of the parthenogenesis was included [17]. DNA-validated parthenogenetic offspring were viable in the closely related Boidae, Boa constrictor, Epicrates maurus and Epicrates cenchria cenchria [57]. Therefore, it is quite likely that the two female pathenogens of Eunectes murinus would have been viable had they not died.

Parthenogenetic events are noticed usually when females are held in long-term isolation, suggesting that the certain period of isolation may be a trigger for parthenogenesis. However, as parthenogenetic events have been reported also in females housed with males [8], the mechanism to trigger parthenogenesis remains unknown. Care must be taken to avoid the reduction of genetic heterogeneity in captive breeding programs in potentially facultative parthenogenetic species by monitoring, with DNA technologies, the occurrence of parthenogenesis.

All of the facultative parthenogenesis reported in Boidae and Pythonidae to date, including the present study, have been exclusively in captivity. Parthenogenesis has been reported in wild populations of other Alethinophidia [13]. However, there has been insufficient effort to detect parthenogenesis in wild boid populations. Hence, parthenogenesis may well also occur in wild boid populations at a low but biologically significant rate. Facultative parthenogenesis via terminal or gametic duplication produces largely autozygous offspring and so results in genetic purging of detrimental variation from the population by enhancing natural selection [23]. Therefore, parthenogenesis may play some role in maintenance of genetic “health” in wild population of parthenogenetic species such as Eunectes murinus. Comparison of heterozygosity in wild populations between parthenogenetic and non-parthenogenetic species may provide opportunity to test the “purging” hypothesis of deleterious alleles.

Supporting information

S1 Table

Estimated allele frequencies of Eunectes murinus markers.



We greatly thank Suma Aqualife Park, Nihondaira Zoo, Noboribetsu Marine Park Nixe, Higashiyama Zoo and Botanical Gardens, Maruyama Zoo and Dr. Shin-ichiro Kawada, National Museum of Nature and Science for providing samples and information for their green anaconda specimen. We also thank Prof. Hidenori Tachida for his technical advice. Prof. Arthur Georges assisted greatly with English expression and grammar on a late draft of the manuscript. All of the authors declare no conflicts of interest associated with this manuscript.

Funding Statement

The authors received no specific funding for this work.

Data Availability

Data Availability

All relevant data are within the paper.


1. Booth W, Schuett GW. 2016. The emerging phylogenetic pattern of parthenogenesis in snakes. Biol J Linn Soc. 118 (2): 172–186.
2. Bell G. 1982. The Masterpiece of Nature: The Evolution and Genetics of Sexuality (Berkeley, CA: University of California Press; ).
3. van der Kooi CJ, Schwander T. 2015. Parthenogenesis: birth of a new lineage or reproductive accident? Curr Biol. 25: R659–61. doi: 10.1016/j.cub.2015.06.055 [PubMed]
4. Groot TVM, Bruins E, Breeuwer JAJ. 2003. Molecular genetic evidence for parthenogenesis in the Burmese python, Python molurus bivittatus. Heredity 90: 130–135. doi: 10.1038/sj.hdy.6800210 [PubMed]
5. Booth W, Johnson DH, Moore S, Schal C, Vargo EL. 2011a. Evidence for viable, non-clonal but fatherless Boa constrictors. Biol Lett. 7: 253–256. [PMC free article] [PubMed]
6. Booth W, Million L, Reynolds RG, Burghardt GM, Vargo EL, Schal C, Tzika AC, Schuett GW. 2011b. Consecutive Virgin Births in the New World Boid Snake, the Colombian Rainbow Boa, Epicrates maurus. J Hered. 102:759–763. [PubMed]
7. Kinney ME, Wack RF, Grahn RA, Lyons L. 2013. Parthenogenesis in a Brazilian rainbow boa (Epicrates cenchria cenchria). Zoo Biol. 32: 172–176. doi: 10.1002/zoo.21050 [PubMed]
8. Booth W, Schuett GW, Ridgway A, Buxton DW, Caston TA, Bastone G, Bennet C, McMahan W. 2014. New insights on facultative parthenogenesis in pythons. Biol J Linn Soc. 112: 461–468.
9. Schuett GW, Fernandez PJ, Gergits WF, Casna NJ, Chiszar D, Smith HM, Mitton JB, Mackessy SP, Odum RA, Demlong MJ. 1997. Production of offspring in the absence of males: evidence for facultative parthenogenesis in bisexual snakes. Herpetol Nat Hist. 5: 1–10.
10. Dubach J, Sajewicz A, Pawley R. 1997. Parthenogenesis in the Arafuran Filesnake (Acrochordus arafurae). Herpetol Nat Hist.5: 11–18.
11. Germano DJ, Smith PT. 2010. Molecular evidence for parthenogenesis in the Sierra Garter Snake, Thamnophis couchii (Colubridae). Southwest Nat. 55: 280–282.
12. Booth W, Schuett GW. 2011. Molecular genetic evidence for alternative reproductive strategies in North American pitvipers (Serpentes: Viperidae): longterm sperm storage and facultative parthenogenesis. Biol J Linn Soc. 104: 934–942.
13. Booth W, Smith C, Eskridge PH, Hoss SK. 2012. Facultative parthenogenesis discovered in wild vertebrates. Biol Lett. 8: 983–985. doi: 10.1098/rsbl.2012.0666 [PMC free article] [PubMed]
14. Reynolds RG, Booth W, Schuett GW, Fitzpatrick BM, Burghardt GM. 2012. Successive virgin births of viable male progeny in the checkered gartersnake, Thamnophis marcianus. Biol J Linn Soc. 107: 566–572.
15. Jordan MA, Perrine-Ripplinger N, Carter ET. 2015. An Independent Observation of Facultative Parthenogenesis in the Copperhead (Agkistrodon contortrix). J Herpetol. 49: 118–121.
16. Gamble T, Castoe TA, Nielsen SV, Banks JL, Card DC, Schield DR, Schuett GW, Booth W. 2017. The discovery of XY sex chromosomes in a Boa and Python. Curr Biol. 27: 2148–2153. doi: 10.1016/j.cub.2017.06.010 [PubMed]
17. O’Shea M., Slater S, Scott R, Smith S, McDonald K, Lawrence B, Kubiak M. 2016. Eunectes murinus (green anaconda) reproduction / facultative parthenogenesis. Herpetol Rev. 47: 73.
18. Fetzner JW Jr. 1999. Extracting high-quality DNA from shed reptile skins: a simplified method. Biotechniques 26: 1052–4. [PubMed]
19. Schmieder R, Edwards R. 2011. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27: 863–864. doi: 10.1093/bioinformatics/btr026 [PMC free article] [PubMed]
20. Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA, Selker EU, Cresko WA, Johnson EA. 2008. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE 3: e3376 doi: 10.1371/journal.pone.0003376 [PMC free article] [PubMed]
21. Sansaloni CP, Petroli CD, Carling J, Hudson CJ, Steane DA, Myburg AA, Grattapaglia D, Vaillancourt RE, Kilian A. 2010. A high-density Diversity Arrays Technology (DArT) microarray for genome-wide genotyping in Eucalyptus. Plant Methods. 6: 16 doi: 10.1186/1746-4811-6-16 [PMC free article] [PubMed]
22. Portnoy DS, Hollenbeck CM, Johnston JS, Casman HM, Gold JR. Parthenogenesis in a whitetip reef shark Triaenodon obesus involves a reduction in ploidy. J Fish Biol. 2014. August;85(2):502–8. doi: 10.1111/jfb.12415 Epub 2014 Jun 6. [PubMed]
23. Hedrick PW. 2007. Virgin birth, genetic variation and inbreeding. Biol Lett. 3: 715–716. doi: 10.1098/rsbl.2007.0293 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science