Email updates

Keep up to date with the latest news and content from BMC Developmental Biology and BioMed Central.

Open Access Methodology article

Induction of diploid gynogenesis in an evolutionary model organism, the three-spined stickleback (Gasterosteus aculeatus)

Irene E Samonte-Padilla1*, Christophe Eizaguirre12, Jörn P Scharsack3, Tobias L Lenz1 and Manfred Milinski1

Author Affiliations

1 Department of Evolutionary Ecology, Max-Planck Institute of Evolutionary Biology, August-Thienemann Str 2, Ploen 24306, Germany

2 Leibniz Institute for Marine Sciences (IFM GEOMAR), Evolutionary Ecology of Marine Fishes, Dusternbrooker Weg 20, 24105, Kiel, Germany

3 Department of Animal Evolutionary Ecology, Institute for Evolution and Biodiversity, WWU Muenster, Huefferstrasse 1, 48149 Muenster, Germany

For all author emails, please log on.

BMC Developmental Biology 2011, 11:55  doi:10.1186/1471-213X-11-55

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-213X/11/55


Received:20 April 2011
Accepted:12 September 2011
Published:12 September 2011

© 2011 Samonte-Padilla et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Rapid advances in genomics have provided nearly complete genome sequences for many different species. However, no matter how the sequencing technology has improved, natural genetic polymorphism complicates the production of high quality reference genomes. To address this problem, researchers have tried using artificial modes of genome manipulation such as gynogenesis for fast production of inbred lines.

Results

Here, we present the first successful induction of diploid gynogenesis in an evolutionary model system, the three-spined sticklebacks (Gasterosteus aculeatus), using a combination of UV-irradiation of the sperm and heat shock (HS) of the resulting embryo to inhibit the second meiotic division. Optimal UV irradiation of the sperm was established by exposing stickleback sperm to a UV- light source at various times. Heat shock parameters like temperature, duration, and time of initiation were tested by subjecting eggs fertilized with UV inactivated sperm 5, 10, 15, 20, 25, or 30 minutes post fertilization (mpf) to 30°C, 34°C, or 38°C for 2, 4, 6 or 8 minutes. Gynogen yield was highest when stickleback eggs were activated with 2 minutes UV-irradiated sperm and received HS 5 mpf at 34°C for 4 minutes.

Conclusions

Diploid gynogenesis has been successfully performed in three-spined stickleback. This has been confirmed by microsatellite DNA analysis which revealed exclusively maternal inheritance in all gynogenetic fry tested. Ploidy verification by flow cytometry showed that gynogenetic embryos/larvae exhibiting abnormalities were haploids and those that developed normally were diploids, i.e., double haploids that can be raised until adult size.

Background

The genomics community has been enormously revolutionized by first- and next generation sequencing (NGS) technologies [1-3] that generated vast amounts of genetic data ranging from unicellular (bacteria) to multicellular (eukaryotic) genomes. However, no matter how good these technologies are, polymorphisms in a genome complicate the assembly process, results in lower quality, and the contiguity and completeness of assembly is significantly lower than would be expected from a homozygous template [4]. Hence, there is a growing interest for the development of inbred lines, such as haploid and double haploid (DH) lines, which are particularly advantageous in genomics [5,6] because of their homozygosity and their growth potential. Inbred lines are also useful in physical mapping [7] and in genetic mapping [8] which would enable precise positioning of the genes in the genome, and thus facilitate accurate identification of candidate genes. Furthermore, DHs are also useful for mutagenesis and genetic transformation studies [9-14].

An inbred line is normally produced by classical sibling mating which is not only time and space consuming, but also laborious and expensive, especially in organisms with long generation times. Alternatively, it can be produced via androgenesis or gynogenesis which both requires inactivation of the genetic material through chemical means or by ionizing radiations. For species with relatively large eggs (1.5 - 1.7 mm), with thick membrane that hardens upon contact with water and contains numerous oil globules like that of the three-spined stickleback, optimum inactivation of the genetic material of the eggs is difficult and hence gynogenesis should be favored. Gynogenesis is defined by Thorgaard [15] as an all maternal type of inheritance wherein the genetic material of the sperm does not contribute to the resulting embryo. This process occurs naturally in some fish species like in crucian carp and in several species of poecilids [16-18]. Induction of diploid gynogenesis involves egg activation by irradiated homologous or heterologous sperm, and diploidization by retention of the second polar body (meiotic gynogenesis) or suppression of the first mitotic cleavage (mitotic gynogenesis) [19] using strong physical treatments (i.e., shocks). The most commonly used treatments are either low or high temperatures (cold or heat shocks) or hydrostatic pressure. Under the influence of such treatments, the spindle fibers are destroyed and cell division stops resulting in the fusion of the daughter cells, thus forming diploid gynogenetic embryos. Mitotic gynogenesis was regarded as the fastest way to produce DH [completely inbred (homozygous) fish]. It requires only one round of gynogenetic manipulation. Production of DHs via meiotic gynogenesis, on the other hand, requires an additional 2-3 rounds of gynogenetic manipulations. Prior to the boom of the genomics era, artificial gynogenesis has already been applied in many fish species mainly because of its potential value in experimental genetics and aquaculture [20,21]. Furthermore, it has played an important role in fish genetic improvement and control of reproduction [20]. For instance, diploid gynogenetic fish were observed to grow faster and have stronger resistance to disease than haploids [22]. Successful gynogenetic manipulations were done, e.g., in rainbow trout [20], medaka [23], common carp [24], muskellunge [25,26], goldfish [27], Russian sturgeon and starlet [28], mud loach [29], sea lamprey [30], and in Wels- and channel catfish [6,31].

The three-spined stickleback is a teleost that exhibits multiple examples of parallel evolution. It inhabits different marine, estuarine, and freshwater systems all over the Northern hemisphere and is characterized by rapid and repeated ecological and phenotypic divergence from the marine ancestor to the freshwater forms [32,33]. Like other model systems, the linkage and chromosome maps of this fish have already been published [34,35] and its genome sequenced [36]. Genomic information has now been used to understand many fundamental evolutionary problems. For instance, the successful elucidation of the evolution of body armor reduction in three-spined sticklebacks from the phenotypic down to the gene level [37] prompted Gibson [38] to elevate the status of this species to an evolutionary supermodel. The advent of NGS has also advanced studies on the population genomics of this species [39]. However despite all these developments, the problem with polymorphism remains and left a significant part of the genome unassembled. Thus, increasing the need for the production of an inbred three-spined stickleback.

In our knowledge, no inbred line is available in this species partly because the emphasis has always been on the study of natural populations. Information about effective modification of genotype in this species was limited only to ploidy manipulation using temperature shock that resulted in the production of triploid and haploid individuals [40,41]. These haploids would be ideal for genomics but as in many attempts of fish haploid production, these haploids survived only until embryonic development [20,42-46] or right after hatching [10,47,48] and hence do not provide enough template material for downstream processing.

Here, we present a working protocol for the induction of gynogenesis using a combination of sperm UV irradiation and heat shock (HS) treatment. This study aims at producing inbred lines of the three-spined stickleback suited not only for whole genome sequencing, physical and genetic mapping but also for future mutagenesis and experimental studies. Success of gynogenesis is confirmed by microsatellite DNA analysis and the ploidy of the gynogenetic embryos/larvae is verified by flow cytometry.

Results

Screening for the optimum UV irradiation exposure

Optimum UV irradiation of the sperms was done by exposing groups of sperms to the UV light source at various durations. UV irradiation of the sperms worked well for both 2 min and 4 min exposure with fertilization rates of 97% and 93%, respectively (Table 1). All the embryos, however, developed abnormally and did not hatch. They were collected at day 3 or 4 for DNA extraction. The sperms irradiated for 6 and 8 minutes failed to activate the stickleback eggs, thus no embryo was produced. This suggests that long exposure to UV irradiation could already be lethal to three-spined stickleback sperms as evidenced by the presence of non-motile sperms after irradiation in these treatment groups.

Table 1. Embryonic survival of three-spined sticklebacks subjected to UV irradiation at 4 various times resulting to zero fertilization success

The normal controls, on the other hand, were not only successfully fertilized but also have high percentage of survival to the eyed stage (98%) and hatched as expected at day 6 or 7.

Effects of heat shock on chromosome diploidization of the activated stickleback eggs

The effectiveness of heat shock (HS) in promoting diploid gynogenesis, expressed by the number of gynogens produced (total gynogen yield), was clearly shown in two separate experiments. Gynogen yield was computed as the percentage of hatched, surviving fry over the total number of fertilized embryos.

In the first experiment, the optimal initiation time of heat shock was established to be at 5 minutes post fertilization (mpf) in combination with heat shock at 34°C for 4 minutes. This treatment group produced the highest diploid gynogen yield of 22% (Figure 1). A low gynogen yield of 3% was also observed in the treatment group subjected to 30°C applied 20 mpf (Figure 1b). The rest of the treatments produced low yields of diploid gynogens (0-2%).

thumbnailFigure 1. Effects of the initiation time of heat shock (HS) on diploid gynogen production. Eggs activated with UV irradiated sperms at different initiation times were heat-shocked at two different temperatures (30°C and 34°C) for 2 and 4 minutes. Control (C) served as HS control activated with UV irradiated sperms but without HS treatment. Effects of HS initiation time on diploid gynogen yield in terms of HS duration (Figure 1a) and in terms of HS temperature (Figure 1b) were showed. Here, results were presented as means and standard error of the raw data from three independent trials.

Using the previously established HS initiation time of 5 mpf and UV irradiation of 2 min, optimal temperature and duration for applying HS was found to be at 34°C for 4 minutes. This treatment group produced a much improved gynogen yield of 65% (Figure 2) while the remaining groups that were heat shocked at 30°C and 38°C for shorter (2 minutes) and longer (>4 minutes) both yielded low diploid gynogenetic fry (<5%). In addition to HS initiation time (p = 0.023) and HS duration (p = 0.003), HS temperature (p = 0.001) was also found to contribute statistically significant effects on the induction of diploid gynogenesis in three-spined stickleback.

thumbnailFigure 2. Effects of temperature and duration of HS on diploid gynogen yield. Eggs activated with UV irradiated sperm were heat shocked starting at 5 mpf at 30°C and 34°C for 2, 4, 6, and 8 minutes. Control 1 and Control 2 served as HS controls activated with UV irradiated sperms 5 mpf and after all the UV treatments, respectively, but not heat shocked. Here, results were presented as means and standard error of the raw data from three independent trials.

All embryos in the HS controls developed into abnormal embryos/larvae which indicated that UV irradiation for 2 min successfully inactivated the genetic material of the sperms but retained the ability to induce the development of the eggs. Such haploid (n) eggs develop abnormally in the absence of HS treatment but become morphologically normal (2n) when optimal HS condition was applied to them.

Microsatellite DNA analysis of gynogenetic families

Using the here established optimal condition for stickleback diploid gynogenesis induction, additional six gynogenetic families designated Family 1-6 were produced for microsatellite analyses. In families 1 - 4, all gynogenetic embryos/larvae were collected at day 7 and used to assess level of heterozygosity in the gynogenetic individuals. Families 5 and 6 were used to determine the type of gynogenesis that was induced in this experiment, hence only the embryos/larvae that hatched and survived past 120 days were collected and tested. In general, all the gynogenetic embryos/larvae subjected to microsatellite analysis using 9 loci were found to carry only the maternal alleles and as expected all the normal controls had biparental genotypes. Data from the four gynogenetic families produced to assess genetic homogeneity of the resulting gynogens revealed that families 1 to 4 had statistically significant lower heterozygosity index values compared to normal control fish (p<0.05). As shown in Figure 3, the microsatellite data were analyzed in two ways: 1) by using the whole dataset and 2) by excluding individuals homozygous for allele shared by both parents (scored as "0.5"). In both cases, gynogens were found to be more homozygous than the control fish. The degree of individual homozygosity, though, was further intensified when alleles shared by both parents were excluded in the test (Figure 3b).

thumbnailFigure 3. Heterozygosity of gynogenetic families revealed by microsatellite DNA data. Nine microsatellite loci were used to analyse the heterozygosity of four gynogenetic families (Treated Family 1-4) and their corresponding positive controls (Control Family 1-4, no UV irradiation and no HS) were analysed using the whole data set (Figure 3a) and by excluding individuals homozygous for alleles shared by the two parental fish and hence scored as "0.5" (Figure 3b).

Moreover, in families 5 and 6 (Table 2), 97-100% of the samples were heterozygous for the maternal alleles. This suggests that majority of the gynogens are meiotic diploids produced by retention of the second polar body. The remaining 3% found to be homozygous for the maternal alleles could be mitotic gynogens produced by inhibiting first mitotic division. These embryos can be good candidate for a completely homozygous inbred fish.

Table 2. Summary of microsatellite genotyping data of the two gynogenetic families used to assess the mode of gynogenesis occurring in three-spined stickleback

Ploidy level of gynogenetic fry

A total of 145 gynogenetic and 16 normal control embryos/larvae were subjected to flow cytometry for relative DNA quantification. Of the 145 gynogenetic embryos/larvae, 102 exhibited abnormal features while the remaining 43 were morphologically normal. A histogram of the relative DNA content of a haploid (n) embryo is illustrated in Figure 4. It had a dominant peak at mean channel 230 which indicates that most of the cells are in the n(G1) phase and a minor peak at mean channel 460 representing cells in the n(G2/M) phase of the cell cycle. Also shown is a histogram of a diploid embryo characterized by a prominent peak at channel 460. The ploidy of the sample is reflected by the dominant peak suggesting that abnormal embryos are haploid (1n, mean channel 230) and the morphologically normal gynogens are diploid (2n, mean channel 460) with twice as much genetic material as the latter.

thumbnailFigure 4. Representative fluorescent distribution of cells prepared from gynogenetic three-spined stickleback embryos measured by flow cytometry. The mean fluorescence peaks of representative haploid and diploid gynogenetic embryos/larvae have mean fluorescence at mean channel numbers 230 and 460, respectively.

Development of the gynogenetic embryos/larvae

The daily larval development of representative stickleback embryos activated with UV-irradiated sperm and with further HS treatment at 34°C for 4 minutes 5 mpf was monitored for one week (Figure 5). The diploid embryo in Figure 5a followed normal development as described by Kuntz & Radcliffe [49]. However the growth of the haploid embryo started to slow down 48 hours after fertilization. It managed to cover only one half of the circumference of the yolk sac and growth from here on was restricted horizontally increasing the volume of the fish but not its length (Figure 5b). The diploid embryos hatched at day 6 or 7 post fertilization, whereas haploid embryos did not hatch and remained within their enveloping membrane until they die. Death of the haploid embryo could be attributed to the many deformities observed mainly in the head and tail regions, as well as in the abdomen where the yolk sac is attached (Figure 6). Tissue degeneration at the abdominal region prevents absorption of the yolk sac thereby cutting the source of nutrients of the developing embryo leading to the death of the embryo.

thumbnailFigure 5. Embryologic development of representative gynogenetic three-spined stickleback embryos. Development of the embryos leading to a normal (Figure 5a) and abnormal (Figure 5b) gynogenetic larvae after receiving HS treatment 15 mpf for 4 min at 34°C were shown.

thumbnailFigure 6. Commonly observed abnormalities in gynogenetic three-spined stickleback embryos. Abnormal cell proliferations (marked by black arrows) were commonly found in abnormal gynogenetic three-spined sticklebacks. These embryos often died within the first week.

Discussion

Gynogenesis is regarded as a powerful tool in the production of an all female progeny, in studying mechanisms of sex determination, and in the production of completely inbred lines. Inbred lines are useful for whole genome sequencing and for production of mutational and developmental genetic screens. For an evolutionary model organism like the three-spined stickleback, optimization of this method will be immensely important.

In this study, UV irradiation of the sperm and heat shock (HS) of the resulting embryos were used to induce gynogenesis. UV irradiation was used because aside from being effective in deactivating the DNA of the sperm, it is cheap and easy to set-up as compared to ionizing irradiations which are more dangerous and require special containment area [20,21]. Since DNA irradiation is a key step to successful gynogenesis, it is therefore important to determine optimal condition for sperm irradiation. UV irradiation should inactivate the sperm DNA maximally and avoids chromosome fragmentation, because partially inactivated sperm DNA interferes in the embryonic development whereas over-exposure to UV light may already kill the sperm [19]. In this study, 2 to 4 minutes exposure to UV irradiation can be regarded as optimal with fertilization success of 93-97%. All the resulting embryos were morphological abnormal (Figure 5b and 6), confirmed to be haploids by flow cytometry, and had only the maternal genotype. However, UV exposure of more than 4 min in our case was already lethal to the sperms resulting in zero fertilization success. This was also indicated by the observed absence of sperm motility in these treatment groups after UV exposure. Although both 2 and 4 min UV irradiations were successful in inactivating the genetic material of the sperm, we used the former in the succeeding experiments to minimize possible "Hertwig effect" wherein haploid embryos resulting from the total destruction of the genetic material in the sperm by high doses of radiation would survive longer than those exposed to lower UV dosage [50,51].

Application of HS was observed to be effective for the retention of the second polar body in eggs activated with UV irradiated sperms resulting in morphologically normal gynogenetic diploid larvae [25,27,52-54]. For HS to be successful, right combination of the different parameters namely temperature, duration, and initiation time of application is critical. From the review of Komen & Thorgaard [21], it is suggested that optimal temperature for HS is species-specific and could vary from 30°C to up to 42°C. Since HS acts through depolymerization of protein complexes, then the best temperature should be close to the upper limit of tolerance. Likewise, optimal time for applying heat appears to be slightly different [55,56]. Timing of HS is very crucial as this decides the kind of gynogens produced. Generally, when HS is applied early, it targets second meiotic division leading to the retention of the second polar body, thereby producing diploid meiotic gynogens. Late application of heat treatment, on the other hand, results in suppression of the first mitotic cleavage, producing an all homozygote female fry (inbred clone). Here, the best yield of diploid gynogens (65%) was observed in eggs activated with UV irradiated sperm and heat shocked at 34°C 5 minutes post fertilization (mpf) for 4 minutes (Figure 2). Although aiming at different end products, the optimal temperature and duration of HS established here correspond to what was used by Swarup [40,41] in his first attempt to produce haploid and triploid three-spined sticklebacks without UV irradiation. Independently tested, this clearly shows that HS at 34°C for 4 minutes can effectively destabilize and disorganize the microtubules of activated stickleback egg, thereby facilitating diploid gynogenesis. The comparatively low diploid gynogen yields at lower (30°C) and higher (38°C) temperatures indicates that 30°C can still be well tolerated by the developing embryo. The effect of temperature shock at 30°C is very minimal, whereas shock at 38°C could have already interfered in a wider range of developmental mechanisms resulting to further abnormalities that drives the embryos to death. Like in common carp [53], the present study clearly demonstrated that temperature, duration, and initiation time of HS play an important role in the success of diploid gynogenesis in three-spined sticklebacks.

Success of gynogenesis induction in three-spined stickleback is verified by microsatellite analysis. Microsatellite DNA technology is commonly used because of its PCR-based, co-dominant Mendelian inheritance, and highly polymorphic characteristics which facilitate molecular typing of the gynogenetic embryos/larvae [57,58]. Success is measured by the absence of paternal alleles in the resulting gynogenetic individuals. In the six families tested here, it was clear that the paternal alleles are absent and only the maternal alleles are passed on to the gynogenetic embryos/larvae. The predominance of heterozygotes (97%) for the maternal alleles only suggests that the established optimum HS condition promotes retention of the second polar body, thereby producing mostly heterozygous meiotic diploid gynogens which have the expected high degree of homozygosity (Figure 3). With the relatively early timing of heat shock (5 mpf) and due to the mode of egg collection (mainly hand stripping), it could be that the eggs used in this study are still undergoing meiosis when the HS was applied. Flow cytometry remains the best method for ploidy determination in gynogenetic stickleback. This method clearly shows that stickleback gynogenetic embryos/larvae that developed normally were diploids (2n) whereas those that showed abnormalities were haploids (n) with peaks at mean channel numbers 460 and 230, respectively (Figure 4). We observed no mosaic or aneuploid, chromosomal anomalies that are also reliably detected through flow cytometry. With the accuracy of the result obtained here, we find no reason to use more laborious and time consuming classical methods of ploidy determination like chromosome preparation and silver staining of nuclear organizer regions (NORs). Since in the present study, possible effects of doublets on fluorescence were excluded in all samples, the observed predominance of haploid cells at n(G1) phase with a significant 5% of the cells staying at the n(G2/M) is quite striking (Figure 4a). This is comparable to what was observed in muskellunge by Lin and Dabrowski [26] which suggests that haploid cells have an abnormal cell cycle that probably takes longer in time than that of their diploid counterparts. This could also mean production of lower cell numbers, hence the observed stunted growth of haploids.

Conclusions

In summary, a working protocol for the induction of diploid meiotic gynogenesis in three-spined sticklebacks is presented in this study. The optimal condition for the inactivation of the sperm and the HS treatment to effectively suppress the second meiotic division in stickleback eggs were established to be at 2 minutes UV irradiation and HS at 34°C 15 mpf for 4 minutes, respectively. Flow cytometry and microsatellites analysis proved to be an ideal combination of methods to determine the ploidy and the all maternal inheritance in the gynogenetic individuals produced.

Methods

Sticklebacks were collected from Grosse Plöner See (northern Germany) from November 2009 to February 2010. The fish were acclimatized first for 3 months in winter (4°C), another 3 months in spring (12°C), before finally being transferred to individual tanks and placed in summer condition (18°C) until they became reproductively mature. Stickleback males with conspicuously developed red-orange coloration on their throat region and gravid females with extended belly and dilated cloacal opening ready to spawn were used in this study.

Male fish was sacrificed by decapitation. The testes were carefully dissected out, and placed in a 40 μm cell strainer (Falcon, BD Biosciences, Europe) standing inside a small plastic petri-dish containing appropriate amount of Hank's Balanced Salt Solution (HbSS, Sigma). Using a mini pestle, the testes were carefully pressed against the mesh of the strainer to release all the sperms into the HbSS. The motility and relative quantity of the collected sperm was checked under the microscope. Only sperm solutions with active and motile sperms were used in the experiments. Likewise the gravid stickleback female was either sacrificed by decapitation and the ovaries were carefully dissected or hand stripped to collect the eggs. The latter method is advantageous because a female can be used several times.

Sperm Inactivation by UV Irradiation

Successful irradiation means finding the adequate radiation that will result to total genetic inactivation of the sperm without destroying their ability to initiate fertilization. To attain this, we first determined the appropriate duration of UV irradiation. The experiment was done twice and each consisted of two males and five female sticklebacks. Sperm from two male sticklebacks were pooled into a 2-ml sperm suspension which was further divided into five parts: one part was not exposed to UV irradiation and served as control. The other four parts were exposed to a UV light source for to 2, 4, 6, and 8 minutes irradiation, respectively. To ensure optimum irradiation, the sperm suspension was placed on thin watch glass and irradiated at a distance of 30 cm under a 280 ultra-violet (type C) light source with continuous gentle shaking at 20 rpm. The sperm was then used to fertilize groups of pooled eggs (58-93 eggs per group) 5 minutes after dilution and UV irradiation. Fertilization was terminated by addition of 7.5 ml hatchery water, which induces swelling of the egg and closure of the micropyle. Activated eggs were transferred to well-aerated breeding jars and grown at 18°C for 6 to 7 days. The breeding jars were kept in the dark in the first 24 hours to avoid possible photo-repair of the UV-inactivated sperm DNA. The number of fertilized and unfertilized eggs was counted on day 2. Unfertilized eggs were removed immediately to avoid fungal contamination and the embryos were allowed to develop until day 6 or 7. Survival to the eyed stage and hatching was observed and recorded at day 7.

Heat Shock (HS) Treatment to Promote DNA diploidization

Two experiments were conducted to determine the effects of HS on the diploidization of the stickleback egg activated with UV irradiated sperm. In both experiments, stickleback gametes were prepared as described above. Each of the experiments was repeated three times using independent sets of females allowing the effect of each treatment to be tested statistically. Since three-spined sticklebacks normally produce about 150-200 eggs at a time, the pooling of gametes strategy of Lin and Dabrowski [26] was adopted to enable us to investigate a wider range of combinations of HS treatments. All UV irradiations from here on were performed for 2 min.

The first experiment was designed to estimate the effect of initial time of HS on chromosome diploidization. The properly mixed pool of stickleback eggs was divided into 25 groups, each containing around 20-35 eggs. One group served as HS control composed of eggs fertilized with UV-irradiated sperms but not heat-shocked. The other 24 groups were subjected to HS in a full factorial design at 30°C or 34°C regulated water baths for 2 or 4 minutes 5, 10, 15, 20, 25, or 30 mpf. To ensure homogeneous water temperature, the water bath was kept shaking at 30 rpm. A second experiment was intended to further explore the optimal level of temperature and HS duration. The pooled eggs were divided into 14 groups. Two groups that were not heat shocked but were fertilized with irradiated sperms served as controls: HS Control 1 was transferred to the breeding jar immediately 5 mpf whereas HS Control 2 was allowed to continue fertilization until all the HS treatments were finished. The remaining groups were heat-shocked at three different temperatures (30, 34, and 38°C) for 2, 4, 6, or 8 minutes, respectively. Breeding proceeded as previously described. Like in the UV irradiation experiments, the number of fertilized and unfertilized eggs was recorded on day 2. Unfertilized eggs were removed and the fertilized embryos were allowed to grow up to day 6 or 7 after which the embryos were collected and sorted into two groups based on their morphological features: normal embryos that hatched and exhibited normal development and abnormal embryos showing haploid syndromes.

Using the here established optimal condition for stickleback diploid gynogenesis induction, additional six gynogenetic families designated Family 1-6 were produced for microsatellite analyses. In families 1 - 4, all gynogenetic embryos/larvae were collected at day 7 and used to assess level of heterozygosity in the gynogenetic individuals. Families 5 and 6 were used to determine the type of gynogenesis that was induced in this experiment, hence only the embryos/larvae that hatched and survived past 120 days were collected and tested.

Molecular Analyses

Genomic DNA isolation for microsatellite DNA analysis

Genomic DNA from 4-7 day old whole embryos or caudal fin clips from the surviving gynogenetic fry and stickleback parents was extracted using Qiagen DNAeasy Kit (Hilden, Germany). The quantity of the extracted DNA was measured using Nanodrop (PeqLab Biotechnology GmbH, Germany). A total of 9 microsatellites multiplexed in 2 polymerase chain reaction (PCR) protocols [59] were used to confirm the exclusively maternal mode of inheritance of the resulting gynogenetic embryos and fry.

Ploidy determination by flow cytometry

Relative DNA content of embryonic/larval cells was used to infer ploidy of the gynogenetic stickleback embryos/fry. DNA content of cells was measured after DNA staining with propidium iodide by means of flow cytometry (FACS Calibur, BD Biosciences, Belgium) according to the combined protocols used for ploidy determination [60] and for cell cycle analysis of stickleback leukocytes [61] with some modifications. In brief, embryos and fry were collected between 11:00 to 13:00 and individually placed in a 40 μm cell strainer (Falcon, BD Biosciences, Europe) standing inside a small plastic Petri-dish containing 1 ml RPMI 1640 medium (Sigma) with 10% v/v distilled water (R-90). Single cell suspensions from embryos/larvae were produced by pressing it through the nylon mesh of the strainer with the help of a pestle. The resulting cell suspensions were transferred to individual wells of a 96 well deep-well micro titer plate and washed twice (600 × g, 10 min, 4°C) with 0.8 ml cold R-90 medium. After washing, cells were resuspended with 200 μl R-90 and fixed with 800 μl ice cold 99% ethanol. After fixation, plates were centrifuged (600 × g, 10 min, 4°C), the supernatant ethanol was removed and the cells were resuspended with 250 μl Sheath fluid (BD) containing RNAse (0.5 mg L-1) to remove background staining of RNA. The cell suspension was incubated for 10 minutes at room temperature before it was stained with propidium iodide (0.5 mg L-1, Sigma), a red fluorescence dye which intercalates with double stranded DNA. Red fluorescence characteristics of at least 10,000 cells were recorded per sample by means of flow cytometry. For analysis, doublet cells were subtracted from single cells as described by Wersto and co-workers [62]. DNA content of single cells was analyzed by CellQuest (BD Immunocytometry Systems, San Jose, CA).

Data Analysis

Statistical analysis was conducted using R stats v. 2.12.2. Data pertaining to UV and HS optimization were log transformed to meet normal distribution. Analyses of variance (ANOVA) were then conducted on fertilization rate using UV irradiation as independent variable and on diploid gynogen yield with the different HS parameters as independent variables. Values were taken as significant when P < 0.5.

Individual heterozygosity was calculated across all scored loci according to Coulson and co-workers [63] with modification in terms of genotype scoring. In addition to the common genotype scoring of "1" (heterozygous) or "0" (homozygous), an individual homozygous for an allele shared by both parents was scored as "0.5". The difference between control and HS treatment was tested in an additional ANOVA.

Authors' contributions

ISP conceptualized and designed the experiment, carried out all the experiments, and wrote the manuscript. CE provided the markers for microsatellite analysis, conducted most of the statistical analyses, and helped draft the manuscript. JPS helped in the flow cytometry experiment and in drafting the manuscript. TL and MM participated in conceptualizing and designing the experiment, and in drafting the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors are grateful for the technical help extended by Stephanie Bartel, Roswithe Derner, and Ines Schultz. We are also thankful to Arne Traulsen and Mahesh Panchal who gave helpful comments on the manuscript.

References

  1. Mardis ER: Next-generation DNA sequencing methods.

    Annu Rev Genomics. Hum Genet 2008, 9:387-402. PubMed Abstract | Publisher Full Text OpenURL

  2. Mardis ER: The impact of next-generation sequencing technology on genetics.

    Trends Genet 2008, 24:133-141. PubMed Abstract | Publisher Full Text OpenURL

  3. Imelfort M, Duran C, Batley J, Edwards D: Discovering genetic polymorphisms in next-generation sequencing data.

    J Plant Biotechnol 2009, 7:312-317. Publisher Full Text OpenURL

  4. Vinson JP, Jaffe DB, O'Neil K, Karlsson EK, Stange-Thomann N, Anderson S, Mesirov JP, Satou N, Satou Y, Nusbaum C, Birren B, Galagan JE, Lander ES: Assembly of polymorphic genomes: Algorithms and application to Ciona savignyi.

    Genome Res 2005, 8:1127-1135. OpenURL

  5. Talon M, Gmitter FG: Citrus genomics.

    International Journal of Plant Genomics 2008, Article ID 528361. OpenURL

  6. Waldbieser GC, Bosworth BG, Quiniou SM: Production of viable homozygous, double diploid channel catfish (Ictalarus punctatus).

    Mar Biotechnol (NY) 2010, 12(4):380-385. Publisher Full Text OpenURL

  7. Gonzalez VM, Benjak A, Henaff EM, Mir G, Casacuberta JM, Garcia-Mas J, Puigdomenech PS: Sequencing of 6.7 Mb of the melon genome using a BAC pooling strategy.

    BMC Plant Biology 2010, 10e.

    246

    OpenURL

  8. Young WP, Wheeler PA, Coryell VH, Keim P, Thoorgard GH: A detailed linkage map of rainbow trout produced using doubled haploids.

    Genetics 1998, 148(2):839-850. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Green MC: Genetic variants and strains of the laboratory mouse. New York: Gustav Fischer Verlag; 1981. OpenURL

  10. Scheerer PD, Thorgaard GH, Allendorf FW: Genetic analysis of androgenetic rainbow trout.

    J Exp Zool 1991, 260:382-390. PubMed Abstract | Publisher Full Text OpenURL

  11. Hirsch N, Zimmerman LB, Grainger RM: Xenopus the Next Generation: X. tropicalis Genetics and Genomics.

    Dev Dynamics 2002, 225:422-433. Publisher Full Text OpenURL

  12. Noramly S, Zimmerman L, Cox A, Aloise R, Fisher M, Grainger RM: A gynogenetic screen to isolate naturally occurring recessive mutations in Xenopus tropicalis.

    Mech Develop 2005, 122:273-287. Publisher Full Text OpenURL

  13. Pelegri F, Dekens MPS, Schulte-Merker S, Maischein HM, Weiler C, Nusslein-Volhard C: Identification of recessive maternal-effect mutations in the zebrafish using a gynogenesis-based method.

    Dev Dynam 2004, 231:324-325. Publisher Full Text OpenURL

  14. Goda T, Abu-Daya A, Carruthers S, Clark MD, Stemple DL, Zimmerman LB: Genetic screens for mutations affecting development of Xenopus laevis.

    PloS Genetics 2006, 2(6):e91. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Thoorgard GH: Ploidy manipulation and performance.

    Aquaculture 1986, 57:57-64. Publisher Full Text OpenURL

  16. Hubbs C, Hubbs L: Apparent parthenogenesis in nature in a form of fish of hybrid origin.

    Science 1932, 76:628-30. PubMed Abstract | Publisher Full Text OpenURL

  17. Schartl M, Wilde B, Schlupp I, Parzefall J: Evolutionary origin of parthenoform, the amazon molly Poecilia formosa on the basis of a molecular genealogy.

    Evolution 1995, 49:827-835. Publisher Full Text OpenURL

  18. Mezhzherin SV, Kokodiĭ SV: Genetic structure of silver crussian carp Carassius (superspecies auratus) (Linnaeus, 1758) colonies from the Middle Dnieper basin.

    Genetika 2010, 46(6):817-24. PubMed Abstract OpenURL

  19. Pandian TJ, Koteeswaran R: Ploidy induction and sex control in fish.

    Hydrobiologia 1998, 384:167-243. Publisher Full Text OpenURL

  20. Chourrout D: Gynogenesis caused by ultraviolet irradiation of salmonid sperm.

    J Exp Zool 1982, 223(2):175-81. PubMed Abstract | Publisher Full Text OpenURL

  21. Komen H, Thorgaard GH: Androgenesis, gynogenesis, and the production of clones in fishes: A review.

    Aquaculture 2007, 269(1-4):150-173. Publisher Full Text OpenURL

  22. Yuan-Dong S, Chun Z, Shao-Jun L, Min T, Chen Z, Yun L: Induction of gynogenesis in Japanese crucian carp (Carassius cuvieri).

    Acta Genetica Sinica 2006, 33:405-412. PubMed Abstract | Publisher Full Text OpenURL

  23. Naruse K, Iriji K, Shima A, Egami N: The production of cloned fish in the medaka (Oryzias latipes).

    J Exp Zool 1985, 236:335-341. PubMed Abstract | Publisher Full Text OpenURL

  24. Bongers ABJ, Sukkel M, Gort G, Komen J, Richter CJJ: Development and use of genetically uniform strains of common carp in experimental animal research.

    Lab Anim 1998, 32(4):349-363. PubMed Abstract | Publisher Full Text OpenURL

  25. Dabrowski K, Rinchard J, Lin F, Garcia-Abiado MA, Schmidt D: Induction of gynogenesis in muskellunge with irradiated sperm of yellow perch proves diploid muskellunge male homogamety.

    J Exp Zool 2000, 287(1):96-105. PubMed Abstract | Publisher Full Text OpenURL

  26. Lin F, Dabrowski K: Effects of sperm irradiation and heat shock on induction of gynogenesis in muskellunge (Esox masquinongy).

    Can J Fish Aquat Sci 1996, 53:2067-2075. Publisher Full Text OpenURL

  27. Paschos I, Natsis L, Nathanailides C, Kagalou I, Kolettas E: Induction of gynogenesis and androgenesis in goldfish Carassius auratus (var. oranda).

    Reprod Domes Anim 2001, 36(3-4):195-198. OpenURL

  28. Rekubratsii AV, Grunina AS, Barmintsev VA, Golovanova TS, Chudinov OS, Abramova AB, Panchenko NS, Kupchenko SA: Meiotic gynogenesis in the stellate, Russian sturgeon and starlet.

    Ontogenez 2003, 32(2):121-131. OpenURL

  29. Nam YK, Choi GC, Kim DS: An efficient method for blocking the 1st mitotic cleavage of fish zygote using combined thermal treatment, exemplified by mud loach (Misgurnus mizolepis).

    Theriogenology 2004, 61:933-945. PubMed Abstract | Publisher Full Text OpenURL

  30. Rinchard J, Dabrowski K, Garcia-Abiado MA: High efficiency of meotic gynogenesis in sea lamprey Petromyzon marinus.

    J Exp Zool B Mol Dev Evol 2006, 306(6):521-527. PubMed Abstract | Publisher Full Text OpenURL

  31. Fopp-Bayat D: Induction of diploid gynogenesis in Wels catfish (Silurus glanis) using UV-irradiated grass carp (Ctenopharyngodon idella) sperm.

    J Exp Zool A Ecol Genet Physiol 2010, 313(1):22-27. OpenURL

  32. Bell MA, Foster SA: The evolutionary biology of the threespine stickleback. Oxford: Oxford University Press; 1994:571.

  33. Cresko W, McGuigan K, Phillips P, Postlethwait J: Studies of threespine stickleback developmental evolution: progress and promise.

    Genetica 2007, 129:105-126. PubMed Abstract | Publisher Full Text OpenURL

  34. Kingsley DM, Zhu BL, Osoegawa K, De Jong PJ, Schein J, Marra M, Peichel C, Amamiya C, Schluter D, Balabhadra S, Friedlander B, Cha YM, Dickson M, Grimwood J, Schmutz J, Talbot WS, Myers R: New genomic tools for molecular studies of evolutionary change in threespine sticklebacks.

    Behaviour 2004, 141:1331-1344. Publisher Full Text OpenURL

  35. Peichel CL, Nereng KS, Ohgi KA, Cole BLE, Colosimo PF, Buerkle CA, Schluter D, Kingsley DM: The genetic architecture of divergence between threespine stickleback species.

    Nature 2001, 414:901-905. PubMed Abstract | Publisher Full Text OpenURL

  36. Kingsley DM: Sequencing the genome of threespine sticklebacks (Gasterosteus aculeatus). [http:/ / www.genome.gov/ Pages/ Research/ Sequencing/ SeqProposals/ SticklebackSEQ.pdf] webcite

    National Human Genome Research Institute White Paper 2003. OpenURL

  37. Colosimo PF, Hosemann KE, Balabhadra S, Villareal G, Dickson M, Grimwood J, Schmutz J, Myers RM, Schluter D, Kingsley DM: Widespread parallel evolution in sticklebacks by repeated fixation of ectodysplasin alleles.

    Science 2005, 307:1928-1932. PubMed Abstract | Publisher Full Text OpenURL

  38. Gibson G: The synthesis and evolution of a supermodel.

    Science 2005, 307:1890-1891. PubMed Abstract | Publisher Full Text OpenURL

  39. Hohenlohe PA, Bassham S, Etter PD, Stiffler N, Johnson EA, Cresko WA: Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags.

    PLOS Genetics 2010, 6(2):e1000862. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Swarup H: Production of heteroploidy in the three-spined stickleback, Gasterosteus aculeatus (L.).

    Nature 1956, 178:1124-1125. Publisher Full Text OpenURL

  41. Swarup H: Production of haploidy in experimentally treated eggs of Gasterosteus aculeatus.

    Proc Natl Acad Sci India Sect B 1957, 27:274-27842. OpenURL

  42. Purdon CE: Radiation induced gynogenesis and androgenesis in fish.

    Heredity 1969, 24:431-444. PubMed Abstract | Publisher Full Text OpenURL

  43. Chourrout D: Thermal induction of diploid gynogenesis and triploidy on the rainbow trout (Salmo gairdneri Richardson).

    Reprod Nutr Dev 1980, 20:727-733. PubMed Abstract | Publisher Full Text OpenURL

  44. Streisinger G, Walker C, Dower N, Knauber D, Singer F: Production of clones of homozygous diploid zebrafish (Brachydanio rerio).

    Nature 1981, 291:293-296. PubMed Abstract | Publisher Full Text OpenURL

  45. Refstie T, Stoss J, Donaldson EM: Production of all female Coho salmon (Oncorhynchus kisutch) by diploid gynogenesis using irradiated sperm and cold shock.

    Aquaculture 1982, 29:67-82. Publisher Full Text OpenURL

  46. Taniguchi N, Kijima A, Fukai J, Inada Y: Conditions to induce triploid and gynogenetic diploid in Ayu (Plecoglossus altivelis).

    Bull Jpn Soc Sci Fish 1986, 52:49-53. Publisher Full Text OpenURL

  47. Nagy A, Rajki K, Horvath L, Csanyi V: Investigation on carp Cyprinus carpio gynogenesis.

    J Fish Biol 1978, 13:215-224. Publisher Full Text OpenURL

  48. Chourrout D: Use of grayling sperm (Thymallus thymallus) as a marker for the production of gynogenetic rainbow trout (Salmo gairdneri).

    Theor Appl Genet 1986, 72:633-636. Publisher Full Text OpenURL

  49. Kunzt A, Radcliffe L: Notes on the embryology and larval development of twelve teleostean fishes.

    Washington Bull Bur Fish 1915, 35:87-134. OpenURL

  50. Iriji K: Gamma-ray irradiation of the sperm of the fish Oryzias latipes and induction of gynogenesis.

    J Radiat Res 1980, 21:263-270. PubMed Abstract | Publisher Full Text OpenURL

  51. Iriji K, Egami N: Hertwig effect caused by UV irradiation of sperm of Oryzias latipes (teleost) and its photoreactivation.

    Mut Res 1980, 69:241-248. Publisher Full Text OpenURL

  52. Gomelsky B, Mims SD, Onders RJ, Shelton WL, Dabrowski K, Garcia-Abiado MA: Induced gynogenesis in black crappie.

    North American Journal of Aquaculture 2000, 62(1):33-41. Publisher Full Text OpenURL

  53. Khan TA, Bhise MP, Lakra WE: Early heat shock induced diploid gynogenesis in common carp (Cyprinus carpio L.).

    Isr J Aquacult-Bamid 2000, 52:11-20. OpenURL

  54. Chakraborty BK, Miah MI, Mirja MJA, Habib MAB: Induction of gynogenesis in endangered sarpunti, Puntius sarana (Hamilton) and evidence for female homogamety.

    Aquaculture 2006, 258:312-320. Publisher Full Text OpenURL

  55. Palti Y, Li JJ, Thoorgard GH: Improved efficiency of heat and pressure shocks for producing gynogenetic rainbow trout.

    Prog Fish Cult 1997, 59:1-13. Publisher Full Text OpenURL

  56. Palti Y, Shirack A, Cnaani A, Hulata G, Avtalion RR, Ron M: Detection of genes with deleterious alleles in an inbred line of tilapia Oreochromis aureus.

    Aquaculture 2002, 206:151-164. Publisher Full Text OpenURL

  57. Castro J, Bouza C, Sanchez L, Cal RM, Piferrer F, Martinez P: Gynogenesis assessment using microsatellite genetic markers in turbot (Scopthalmus maximus).

    Mar Biotech (NY) 2003, 5(6):584-592. Publisher Full Text OpenURL

  58. Chen SL, Liu YG, Xu MY, Li J: Isolation and characterization of polymorphic microsatellite loci from an EST-library of red sea bream (Chrysophrys major) and cross-species amplification.

    Mol Ecol Notes 2005, 5:215-217. Publisher Full Text OpenURL

  59. Eizaguirre C, Yeates SE, Lenz TL, Kalbe M, Milinski M: MHC-based mate choice combines good genes and maintenance of MHC polymorphism.

    Mol Ecol 2009, 18(15):3316-29.

    Epub 2009 Jun 11

    PubMed Abstract | Publisher Full Text OpenURL

  60. Allen SK: Flow cytometry-assaying experimental polyploid fish and shellfish.

    Aquaculture 1983, 33:317-328. Publisher Full Text OpenURL

  61. Scharsack JP, Koch K, Hammerschmidt K: Who is in control of the stickleback immune system: interactions between Schistocephalus solidus and its specific vertebrate host.

    P Roy Soc Lond B Biol 2007, 274:3151-58. Publisher Full Text OpenURL

  62. Wersto RP, Chrest FJ, Leary JF, Morris C, Stetler-Stevenson MA, Gabrielson E: Doublet discrimination in DNA cell-cycle analysis.

    Cytometry 2001, 46:296-306. PubMed Abstract | Publisher Full Text OpenURL

  63. Coulson TN, Pemberton JM, Albon SD, Beaumont M, Marshall TC, Slate J, Guiness FE, Clutton-Brock TH: Microsatellites reveal heterosis in red deer.

    Proc R Soc Lond B 1998, 265:489-495. Publisher Full Text OpenURL