INTRODUCTION
The genus Strongyloides comprises over 50 species of nematodes that parasitize mammals, amphibians, reptiles and birds (Viney and Lok, Reference Viney and Lok2007). Strongyloides nematodes have complex but interesting life cycles. Infection by Strongyloides begins when the infective third stage larvae (iL3) attach to and penetrate the host skin. Once inside the host, they moult twice into parasitic adults and settle in the small intestine of the host. Then, the parasitic adults, all female, produce eggs by parthenogenesis (Fig. 1A). Once the eggs or hatched larvae are excreted from the host, they develop via the homogonic route into iL3 forms or the heterogonic route into free-living stages that reproduce sexually outside the host (Fig. 1A); sexual reproduction in the free-living generation was shown for two species, Strongyloides ratti (Viney, Reference Viney1996) and Strongyloides papillosus (Eberhardt et al. Reference Eberhardt, Mayer and Streit2007). The progeny of free-living adults are entirely female and develop into iL3 forms in most Strongyloides species, including S. ratti and Strongyloides stercoralis, with a few exceptions (Yamada et al. Reference Yamada, Matsuda, Nakazawa and Arizono1991; Streit, Reference Streit2008).
Parthenogenetic reproduction by parasitic females was found to be mitotic in cytological studies for S. papillosus, Stronglyoides ransomi (Triantaphyllou and Moncol, Reference Triantaphyllou and Moncol1977) and for S. ratti (Chitwood and Graham, Reference Chitwood and Graham1940), and using molecular markers for S. ratti (Viney, Reference Viney1994) and S. papillosus (Nemetschke et al. Reference Nemetschke, Eberhardt, Hertzberg and Streit2010). In these species, the progeny from parasitic females can develop into one of three distinct morphologies: free-living females, free-living males (both via the heterogonic route) or iL3s (homogonic route). Although the ratio of these morphologies is influenced by environmental factors such as host immune response, temperature and pH (Arizono, Reference Arizono1976; Moncol and Triantaphyllou, Reference Moncol and Triantaphyllou1978; Nwaorgu, Reference Nwaorgu1983; Viney, Reference Viney1996; Harvey et al. Reference Harvey, Gemmill, Read and Viney2000; Minato et al. Reference Minato, Kimura, Shintoku and Uga2008; Sakamoto and Uga, Reference Sakamoto and Uga2013), the detailed mechanisms which determine their route of development remain unclear. Moreover, S. stercoralis and S. ratti harbour three pairs of chromosomes, and one of them is a sex (X) chromosome. Females of these species have an XX and free-living males have an XO karyotype (Streit, Reference Streit2008). Genetic material homologous to two chromosomes in S. ratti, namely I and X, appears combined into one chromosome in S. papillosus. Additionally, an XX/XO karyotype appears to be functionally restored in males by a chromatin diminution event (Nemetschke et al. Reference Nemetschke, Eberhardt, Hertzberg and Streit2010).
Strongyloides venezuelensis parasitizes rodent species including rats, mice and Mongolian gerbil Meriones unguiculatus and is distributed worldwide (Brumpt, Reference Brumpt1934; Wertheim and Lengy, Reference Wertheim and Lengy1964; Little, Reference Little1966; Hasegawa et al. Reference Hasegawa, Orido, Sato and Otsuru1988). Strongyloides venezuelensis, alongside S. ratti, is one of the most widely used laboratory models to study Strongyloides infection and mucosal immunity (Sato and Toma, Reference Sato and Toma1990; El-Malky et al. Reference El-Malky, Maruyama, Al-Harthi, El-Beshbishi and Ohta2013). Although both S. ratti and S. venezuelensis have rodents as hosts, the individual mechanisms by which they establish parasitism are thought to differ from each other. Phylogenetic analysis with ribosomal RNA gene suggested that S. venezuelensis is not as closely related to S. ratti as to S. papillosus and the primate parasite Strongyloides fuelleborni (Dorris et al. Reference Dorris, Viney and Blaxter2002). The migration patterns of the two rodent parasites in the host (Takamure, Reference Takamure1995), as well as mechanisms of host immunomodulation for successful parasitism, differ from each other as well (Wilkes et al. Reference Wilkes, Bleay, Paterson and Viney2007; Matsumoto et al. Reference Matsumoto, Sasaki, Yasuda, Takai, Muramatsu, Yoshimoto and Nakanishi2013). More interestingly, observation of free-living adults is less common in S. venezuelensis as compared with S. ratti (Hasegawa et al. Reference Hasegawa, Orido, Sato and Otsuru1988; Harvey et al. Reference Harvey, Gemmill, Read and Viney2000), suggesting that reproduction strategies used by the two species may also be different. Thus, the parasitic abilities of the two species may have evolved independently from each other, and a comparison between the two may provide invaluable insights into understanding the parasitic mechanism and its evolution in Strongyloides nematodes.
Though S. venezuelensis is a widely used laboratory model, some aspects of its basic biology remain unclear. In this report, we re-examined the phylogenetic position of S. venezuelensis within the genus Strongyloides, using 18S rRNA and the D3 expansion segment of 28S rRNA genes. We also examined the developmental routes of progeny from the parasitic females, chromosome behaviour in the germ cells and early embryos in S. venezuelensis. Our analyses suggest that acquisition of rodent parasitism occurred independently in S. venezuelensis and S. ratti.
MATERIALS AND METHODS
Culturing and handling nematodes
Strongyloides venezuelensis HH1 and S. ratti TDI isolates were used in this study. They were isolated from Okinawa, Japan (Hasegawa et al. Reference Hasegawa, Orido, Sato and Otsuru1988) and Tokyo, Japan, respectively, and maintained in the Parasitology laboratory of the University of Miyazaki, using male Wistar rats. Infectious aliquots were prepared by faecal culture using filter paper at 27 °C for 2 days and 5 days for S. venezuelensis and S. ratti, respectively (Sato and Toma, Reference Sato and Toma1990). The nematodes were washed three times in distilled water, and administered by subcutaneous injection. Strongyloides papillosus iL3 was supplied from Dr Ayako Yoshida of University of Miyazaki. Parastrongyloides trichosuri DNA (strain KNP from Warwick Grant's lab) was supplied by Berriman lab from Wellcome Trust Sanger Institute, UK. Strongyloides stercoralis was collected in Yangon, Myanmar.
PCR conditions and DNA sequencing
Nearly full-length 18S ribosomal DNA (rDNA) was amplified with primers 988F–1912R and 1813F–2646R from a lysate of a single nematode as described previously (Holterman et al. Reference Holterman, Van Der Wurff, Van Den Elsen, Van Megen, Bongers, Holovachov, Bakker and Helder2006). These primers amplify 1652 bp out of 1754 bp of full-length 18S rDNA of C. elegans (NR_000053·1). The D3 region of 28S rDNA was amplified using primers D3A-D3B (Nunn et al. Reference Nunn, Theisen, Christensen and Arctander1996). PCR amplifications were carried out in 30 μL reaction mixtures containing 15 μL GoTaq Green Master Mix (Promega), 0·5 μ m of each primer, and 1 μL of appropriately diluted nematode lysate under thermal-cycling conditions of 94 °C for 1 min, followed by 30 cycles of 94 °C for 30 s, 53 °C for 30 s and 72 °C for 1 min. PCR products were purified before sequencing using a MinElute 96 UF PCR purification plate (QIAGEN). DNA sequencing was performed using the BigDye Terminator 3.1 kit and ABI PRISM 3700 or 3130 Genetic Analyzer (Applied Biosystems).
Phylogenetic analyses
Nearly full-length 18S rRNA gene and D3 expansion segments of 28S rRNA gene of Strongyloides species were used for phylogenetic reconstruction. Parastrongyloides trichosuri was used as an outgroup taxon. The newly obtained and published sequences for each gene were aligned using MAFFT (Katoh et al. Reference Katoh, Misawa, Kuma and Miyata2002) with iterative refinement method (FFT-NS-i). Both datasets were concatenated using Concatenator v.1.1.0. (Pina-Martins and Paulo, Reference Pina-Martins and Paulo2008). Phylogenetic analyses of the sequence dataset were performed with maximum likelihood (ML) using PAUP*4b10 (Wilgenbusch and Swofford, Reference Wilgenbusch and Swofford2003) and Bayesian inference (BI) using MrBayes 3.1.2 (Huelsenbeck and Ronquist, Reference Huelsenbeck and Ronquist2001). The best-fitted model of DNA evolution was obtained using ModelTest v.2 (Darriba et al. Reference Darriba, Taboada, Doallo and Posada2012) with the Akaike Information Criterion (AIC). The Akaike-supported model, base frequency, proportion of invariable sites, gamma distribution shape parameters and substitution rates in the AIC were then used in phylogenetic analyses. BI analysis under the TVM+G model for D3 expansion segment of 28S and the TVM+I+G model for 18S rDNA was initiated with a random starting tree and run with four Metropolis-coupled Markov chain Monte Carlo (MCMC) for 1×106 generations in the combined dataset. The MCMC were sampled every 100 generations. Two runs were performed for each analysis. After discarding burn-in samples and evaluating convergence, the remaining samples were retained for further analyses. The topologies were used to generate a 50% majority rule consensus tree. Trees were visualized using TreeView (Page, Reference Page1996). In ML analysis the estimation of the support for each node was obtained by bootstrap analysis with 1000 heuristic search replicates using the previously obtained models. Posterior probabilities and bootstrap support are given on respective clades.
Free-living occurrence observation
Three rats were infected with 500 iL3s of S. venezuelensis or S. ratti. Faeces samples from each rat were collected at 8, 12, 14, 16 and 24 days post-infection (d.p.i.) for S. venezuelensis and at 8 d.p.i. for S. ratti. Faeces (approx. 0·5 g) were cultured on a 2% (w/v) agar plate at various temperatures (19, 25, or 30 °C) for 3 days and free-living male and female nematodes were counted. Three agar plates for each day and each condition were used. A portion of the faecal samples (c. 1·5 g) was diluted in distilled water to determine the eggs/larvae per gram (epg/lpg) of faeces. All statistical analyses were performed using the software package R version 2.15.2 (http://www.r-project.org).
Early embryogenesis
Parasitic females were obtained from rat intestine at 5–10 d.p.i. Fresh parasites were transferred to a pre-warmed (37 °C) agarose pad (4% (w/v) agarose) prepared on a microscope slide (Shaham, Reference Shaham2006) and covered with a silicon grease-rimmed cover slip and viewed under a Nomarski microscope (IX71, Olympus, Japan). All procedures involving the microscopic observations were performed at 37 °C.
DAPI staining/microscopy
Parasitic females were fixed with ice-cold methanol and stained with DAPI as previously described (Nemetschke et al. Reference Nemetschke, Eberhardt, Hertzberg and Streit2010). To collect the eggs, parasitic females were allowed to lay eggs in PBS, the eggs were then squashed and stained as previously described with the exception of using methanol instead of acetic acid for fixation (Albertson et al. Reference Albertson, Nwaorgu and Sulston1979). Microscopic observations were carried out using a confocal laser scanning microscope (LSM700, Zeiss).
RESULTS
Phylogenetic position of S. venezuelensis
In order to determine the phylogenetic position of S. venezuelensis in the genus, we used nearly full-length 18S rRNA and D3 expansion segments of 28S rRNA genes. We sequenced these regions in S. venezuelensis, S. papillosus, S. stercoralis, S. ratti and P. trichosuri for phylogenetic analysis. 18S sequences from other Strongyloides species obtained from the public database were also included in the analysis. Our phylogenetic tree using P. trichosuri as an outgroup indicated that the Strongyloides species can be divided into two clades (Fig. 2): one clade comprising of S. papillosus, S. fuelleborni, S. venezuelensis, Strongyloides callosciureus, Strongyloides robustus and Strongyloides cebus, and the other including S. stercoralis, Strongyloides procyonis and S. ratti. Within the former clade, S. venezuelensis was further clustered into a subclade with S. papillosus and S. fuelleborni.
Occurrences of free-living adults in S. venezuelensis
As the ratio of free-living females/males vs females developing homogonically into iL3s is influenced by several factors, including host immunity and other environmental conditions (temperature, pH etc.), in other Strongyloides species (Arizono, Reference Arizono1976; Moncol and Triantaphyllou, Reference Moncol and Triantaphyllou1978; Nwaorgu, Reference Nwaorgu1983; Viney, Reference Viney1996; Harvey et al. Reference Harvey, Gemmill, Read and Viney2000; Minato et al. Reference Minato, Kimura, Shintoku and Uga2008; Sakamoto and Uga, Reference Sakamoto and Uga2013), we used two parameters (days post infection and faeces incubation temperature) to examine the occurrence of free-living nematodes (FLNs) in S. venezuelensis.
In a total of 162 attempts, comprising of six d.p.i.s (8, 12, 14, 16, 20 and 24) and three temperature conditions (19, 25 and 30 °C), each with three rats and triplicate, we observed FLNs of S. venezuelensis in only 19 cases. As this frequency is much lower than that reported in S. ratti (Harvey et al. Reference Harvey, Gemmill, Read and Viney2000; Minato et al. Reference Minato, Kimura, Shintoku and Uga2008; Sakamoto and Uga, Reference Sakamoto and Uga2013), we validated our method using S. ratti in place of S. venezuelensis. In the S. ratti experiment we observed FLNs in all culture plates tested (Table S1) and the numbers of FLNs (8·32 to 17·4 FLNs per 1000 eggs) were similar to those reported previously (Minato et al. Reference Minato, Kimura, Shintoku and Uga2008). Therefore we confirmed that this low frequency (19 out of 162 cases) was not due to errors in handling of the samples. Of note is that in some culture plates (52 cases) we did not observe the homogonically developed nematodes (iL3s) although the faeces contained a sufficient number of eggs, which may suggest that the development of S. venezuelensis is more sensitive to environmental conditions than that of S. ratti. Furthermore, free-living males were never detected in S. venezuelensis (Table 1) while approximately half of FLNs in the S. ratti experiment were males (Table S1). Consequently, no hatching was observed in eggs derived from the free-living females (data not shown).
Fifty-four culture plates were used at each temperature.
Numbers in parentheses represent numbers of incidents of free-living observation (i.e. 17(7/54) indicates a total number of 17 free-living nematodes found in 7 out of 54 culture plates).
In the 19 plates with S. venezuelensis FLNs, 7 plates (17 FLNs in total) were from 25 °C culture and 12 plates (46 FLNs) were from 30 °C culture (Table 1). We did not observe any FLNs in 19 °C cultures (Table 1).
We observed S. venezuelensis FLNs in 1, 3, 4, 3, 2, and 2 cases at 8, 12, 14, 16, 20 and 24 d.p.i., respectively (Fig. 3) and the number of FLNs on each d.p.i. was 3, 6, 6, 43, 4, and 4, respectively. Although the number of eggs per gram of faeces (epg) decreased as the d.p.i. increased (Fig. 3), the number of instances that FLNs were observed did not change significantly (Generalized Linear Model with binomial error distribution: df = 52, P>0·9). Then, we sought to investigate the relationship between epg and occurrence of FLNs. Figure S1 shows a plot of epg and number of FLNs per gram of faeces. The trend observed was that fewer epg showed more numbers of FLNs per egg (Fig. S1).
The highest number of FLNs was obtained with faeces collected at 16 d.p.i. and cultured at 30 °C (34 FLNs, 6·58 per 1000 eggs). However this was an outlier as the other two plates under the same conditions (same rat, d.p.i. and temperature) showed much lower numbers of FLNs (0·24 or 0 per 1000 eggs).
Germ cells in the parasitic female
Strongyloides venezuelensis parasitic females have two elongated (didelphic) gonads extending from the vulva to the head or the tail directions, reflexed around the beginning of the intestine or anus, continuing around the vulva (Little, Reference Little1966) (Fig. 4A, B). Egg development and embryogenesis progress as stages move along the tract from the germinal tissue to the vulva. Both gonads spiral about themselves around the intestine (Fig. 4A, B, E). As with many nematodes the distal part of the germ line was a syncytium; nuclei were located at the inner surface of the germ line (Fig. 4C–E). Chromosomes in all nuclei in the distal region dispersed peripherally, and nuclei were condensed in the ‘germinal zone’ (Fig. 4F). As nuclei moved away from the germinal zone they were packed one by one into a cell in the ‘growth zone’ and became oocytes (Fig. 4G). Germinal vesicles appeared and chromosomes became condensed during the oocyte movement towards the proximal region, and they seemed to be arrested at this stage until they passed through the oviduct (Fig. 4G). We detected four chromosomes and this state could be prometaphase of the ‘maturation division’ (see below). We did not identify any sperm or sperm nuclei in the gonad in either bright-field or DAPI-stained observations. From these results, we conclude that parasitic females of S. venezuelensis reproduce parthenogenetically.
Early embryogenesis of eggs from the parasitic female
To see if the embryo starts development without fertilization, we observed oocyte maturation and early embryogenesis in the parasitic female sequentially with light microscopy. After the mature oocyte passed through the oviduct into the uterus, the eggshell was formed (Fig. 5A). The embryo shape was oval and the long axis of the embryo in the uterus was in parallel to the anterior-posterior axis of the mother. The germinal vesicle disappeared, and then only one pronucleus appeared at the lagging side-pole of the embryo (Fig. 5A, B). At this stage a protruded polar body-like structure was observed adjacent to the pronucleus (Fig. 5B). Pseudocleavage furrow and cytoplasmic streaming were observed (Fig. 5C), the pronucleus moved towards the middle of the cell (Fig. 5D) and its membrane broke down (Fig. 5E, F). Subsequently the embryo divided to form the two-cell stage (Fig. 5G, H) and the four-cell stage (Fig. 5I). All embryos (n = 7) observed in these experiments developed to first-stage larvae and hatched successfully (not shown), suggesting that the process observed was normal development.
Chromosome structure and behaviour during early embryogenesis
To clarify the mechanism by which ploidy is maintained, we observed chromosome structure and behaviour during early embryogenesis by DAPI staining. Developmental stages of the fixed and stained embryos were decided by comparison to the living embryos observed above. After passing through the oviduct into the uterus, the oocyte resumed mitotic cell division (Fig. 6A). This division always occurred asymmetrically at the lagging side-pole of the embryo (Fig. 6A) to produce one polar body. Polar bodies stained clearly as blue dots, but often disappeared as embryogenesis progressed. One pronucleus was reconstructed with vesicle (Fig. 6B), moved to the centre, and chromosomes were visible in the prophase stage (Fig. 6C). Chromatids fused and aligned along the metaphase plate (Fig. 6D) and migrated to each pole (Fig. 6E) subsequently forming two-cell stage embryos (Fig. 6F). The cleavage occurred asymmetrically with a bigger blastomere in the leading pole than the lagging one, which is contrary to the observations in Caenorhabditis elegans (Wallenfang and Seydoux, Reference Wallenfang and Seydoux2000) and other Caenorhabditis nematodes (Brauchle et al. Reference Brauchle, Kiontke, Macmenamin, Fitch and Piano2009). Most eggs were laid as two-cell stage embryos. Then, we collected the laid eggs, squashed and stained with DAPI. Four chromosomes, two longer than the others, were observed in prophase cells in the four-cell stage embryos (Fig. 6G, H), suggesting that the karyotype of S. venezuelensis is 2n = 4. We confirmed in nine eggs, all of which had four chromosomes, and no specimens showed different compositions of chromosomes.
DISCUSSION
Strongyloides spp. are unique among parasitic nematodes in having both parasitic and free-living stages in their life cycle. Since the parasitic females produce eggs by parthenogenesis, the free-living stage is the only occasion for them to exchange genetic materials with other individuals. Eggs that are parthenogenetically produced by the parasitic females can develop into three morphs: free-living males, free-living females or iL3. In S. ratti the frequency of free-living development is as high as 60% (Minato et al. Reference Minato, Kimura, Shintoku and Uga2008), which suggests frequent exchanges of genetic materials take place in this species. In this study we observed a much lower frequency of free-living development in S. venezuelensis than S. ratti. Furthermore, we did not find any free-living males. This suggests that S. venezuelensis rarely, or possibly never, exchanges genetic materials with other individuals in its life cycle. It is still possible, however, that there are unknown triggers to stimulate the nematode to have more free-living males and enable them to perform genetic exchanges. This could be an artificial loss of ability due to the stress of prolonged maintenance in laboratory, though absence of free-living males was reported previously in a recent field isolate of S. venezuelensis (Hasegawa et al. Reference Hasegawa, Orido, Sato and Otsuru1988). It would be of value to investigate free-living occurrence rates in other strains or wild isolates to confirm if this is a general characteristic of the species or the observation in this study is an exceptional case.
We noted that the number of free-living females increased with days after infection, although it was not clear which parameter was more influential, days post infection (i.e. ageing of worms and effect of host immunity) or dispersed density of eggs in the faeces. We also observed a higher number of free-living females in 30 and 25 °C cultures than at 19 °C. A similar trend for temperature conditions was also observed with other Strongyloides species (Arizono, Reference Arizono1976; Nwaorgu, Reference Nwaorgu1983; Minato et al. Reference Minato, Kimura, Shintoku and Uga2008; Sakamoto and Uga, Reference Sakamoto and Uga2013).
We observed four chromosomes in S. venezuelensis somatic cells, suggesting S. venezuelensis has a karyotype of 2n = 4. It is known that S. ratti has three pairs of chromosomes (2n = 6) and one of them is a sex chromosome. On the other hand, S. papillosus females have two pairs of chromosomes (2n = 4). The lengths of the two are significantly different and the longer one is suggested to be the result of a fusion of a sex chromosome with one of the autosomes of 2n = 6 (Triantaphyllou and Moncol, Reference Triantaphyllou and Moncol1977). The number of chromosomes in S. venezuelensis could also result from such a fusion. However, the lengths of the two pairs of the four S. venezuelensis chromosomes did not differ from each other to the extent seen in S. papillosus (Nemetschke et al. Reference Nemetschke, Eberhardt, Hertzberg and Streit2010). It would be very interesting to correlate the evolutionary history of reproduction in Strongyloides species with their chromosome structures.
Nemetschke et al. showed that S. papillosus employs chromatin diminution to exclude chromosome regions corresponding to the S. ratti sex chromosome in males (Albertson et al. Reference Albertson, Nwaorgu and Sulston1979; Nemetschke et al. Reference Nemetschke, Eberhardt, Hertzberg and Streit2010). Although S. venezuelensis has the same number of chromosomes as S. papillosus and chromosome diminution could generate males in this species, we did not observe such events or any different chromosome composition in this study. This correlates with the absence of males in S. venezuelensis, although the basis for the lack of chromatin diminution is unclear: it could be that we were unable to generate S. venezuelensis males via this procedure or that males were not observed because S. venezuelensis is not capable of chromatin diminution.
Because embryogenesis occurred without fertilization and the chromosome number in mature oocytes, which are located proximal to the oviduct, was the same as that of the somatic cells it is likely the S. venezuelensis reproduces via mitotic parthenogenesis. Even without sperm stimulation, an embryo of S. venezuelensis had a polarity and produced the polar body on the lagging pole of the embryo. This orientation is opposite from C. elegans (Wallenfang and Seydoux, Reference Wallenfang and Seydoux2000). Parasitic stages of Strongyloides species including S. ratti or S. papillosus have been reported to have only female sex and reproduce by mitotic parthenogenesis (Zaffagnini, Reference Zaffagnini1973; Triantaphyllou and Moncol, Reference Triantaphyllou and Moncol1977; Viney, Reference Viney1994). Our observation is consistent with these reports. In order to confirm this, molecular studies would be useful by checking the progenies have the same genotype (clones) or have genetic variations due to cross-overs of sister chromosomes.
Our phylogenetic analysis using genes for 18S and D3 of 28S ribosomal RNA showed two well-supported clusters in Strongyloides species, namely a group including S. ratti and S. stercoralis and another including S. papillosus and S. venezuelensis (Fig. 2). This is consistent with the gonad morphologies of these species, wherein members of the S. papillosus group have spiral morphology, while those of the S. ratti group have straight morphology (Bartlett, Reference Bartlett1995; Little, Reference Little1966; Sato et al. Reference Sato, Torii, Une and Ooi2007). The number of chromosomes observed for Strongyloides species was also consistent with the phylogeny; karyotypes known thus far are 2n = 6 for S. ratti and S. stercoralis and 2n = 4 for S. papillosus and S. venezuelensis. Strongyloides ratti and S. venezuelensis are both parasites of rodents. However they are clearly different in terms of rRNA phylogeny, chromosome number and gonad morphology. Therefore, it is likely that they acquired their ability to parasitize rodents independently from each other.
SUPPLEMENTARY MATERIAL
To view supplementary material for this article, please visit http://dx.doi.org/S0031182014001036.
ACKNOWLEDGEMENTS
The authors thank Mark Viney and Yasunobu Maeda for reading and commenting on the manuscript, Nancy Holroyd, Matthew Berriman, Warwick Grant, Ayako Yoshida and Yoichiro Horii for supplying DNA or other biological materials, and Valerie Williamson for valuable suggestions about embryo observation. This study was partly supported by JSPS KAKENHI [grant numbers 24780044, 24659190, 26292178 and 26460510]. AH was supported by the Association for Preventive Medicine of Japan.
DATA DEPOSITION
Sequence data has been deposited in GenBank/EMBL/DDBJ under accession numbers AB923880–923889.