A novel nucleo-cytoplasmic hybrid clone formed via androgenesis in polyploid gibel carp.

BACKGROUND: Unisexual vertebrates have been demonstrated to reproduce by gynogenesis, hybridogenesis, parthenogenesis, or kleptogenesis, however, it is uncertain how the reproduction mode contributes to the clonal diversity. Recently, polyploid gibel carp has been revealed to possess coexisting dual modes of unisexual gynogenesis and sexual reproduction and to have numerous various clones. Using sexual reproduction mating between clone D female and clone A male and subsequent 7 generation multiplying of unisexual gynogenesis, we have created a novel clone strain with more than several hundred millions of individuals. Here, we attempt to identify genetic background of the novel clone and to explore the significant implication for clonal diversity contribution.
METHODS: Several nuclear genome markers and one cytoplasmic marker, the mitochondrial genome sequence, were used to identify the genetic organization of the randomly sampled individuals from different generations of the novel clone.
RESULTS: Chromosome number, Cot-1 repetitive DNA banded karyotype, microsatellite patterns, AFLP profiles and transferrin alleles uniformly indicated that nuclear genome of the novel clone is identical to that of clone A, and significantly different from that of clone D. However, the cytoplasmic marker, its complete mtDNA genome sequence, is same to that of clone D, and different from that of clone A.
CONCLUSIONS: The present data indicate that the novel clone is a nucleo-cytoplasmic hybrid between the known clones A and D, because it originates from the offspring of gonochoristic sexual reproduction mating between clone D female and clone A male, and contains an entire nuclear genome from the paternal clone A and a mtDNA genome (cytoplasm) from the maternal clone D. It is suggested to arise via androgenesis by a mechanism of ploidy doubling of clone A sperm in clone D ooplasm through inhibiting the first mitotic division. Significantly, the selected nucleo-cytoplasmic hybrid female still maintains its gynogenetic ability. Based on the present and previous findings, we discuss the association of rapid genetic changes and high genetic diversity with various ploidy levels and multiple reproduction modes in several unisexual and sexual complexes of vertebrates and even other invertebrates.

Background

Since Amazon molly Poecilia formosa, the first clonal reproduction vertebrate, was found in 1932[1], about 90 all-female unisexual complexes have been reported in fish, amphibians and reptiles[2]. These unisexual animals have been revealed to reproduce by gynogenesis, hybridogenesis, parthenogenesis, or kleptogenesis[3-6], but it has remained uncertain how the reproduction mode contributes to the clonal diversity. Recently, unique dual modes of unisexual gynogenesis and sexual reproduction have been discovered in polyploid gibel carp[7,8], which gives a significant clue to the formation of clone diversity in unisexual vertebrates.

Gibel carp, also commonly known as silver crucian carp or Prussian carp, which was previously nominated as a subspecies (Carassius auratus gibelio) of goldfish, has been currently suggested as a separate species Carassius gibelio[9]. It was preliminarily uncovered to have unisexual reproduction ability of allogynogenesis early in the last century[10,11]. Allogynogenesis is a form of gynogenesis stimulated by heterologous sperm from other fish species. In allogynogenesis, some supernumerary minichromosomes or chromosomal fragments can be accidently incorporated into the maternal genomes[12,13], but the heterologous sperm does not decondense and fuse with the female pronucleus[8,11]. Diverse karyotypes with 156 or 162 chromosomes[14] and triploid origin relative to goldfish with 100 chromosomes[15] were confirmed in wild gibel carp populations. Moreover, numerous different clones and genetic diversity were discriminated by serum transferrin phenotypes[16], RAPD and SCAR markers[17,18], microsatellite markers[19], transferrin allele polymorphism[20,21], and mtDNA control region sequences[22]. Significantly, an approximate 1%-10% of males identical to the females in genetic background were found in natural habitats, and genetic recombination evidence was obtained in the offspring produced by bisexual mating between different clones[7]. Therefore, genetic evidence indicates that gibel carp have coexisting dual modes of unisexual gynogenesis and gonochoristic sexual reproduction, as it can reproduce all-female offspring through unisexual gynogenesis when the eggs are activated by heterologous sperm from other fish species, and also produces female and male offspring through sexual reproduction when the eggs are inseminated by homologous sperm from the gibel carp males[8].

Theoretically, the dual reproduction modes provide a unique way for creating novel clones in laboratory, because numerous genetic recombination offspring can be obtained by sexual reproduction mating between different clones, and some better and valuable clones can be selected and proliferated by unisexual gynogenesis[8]. To exploit the new breeding potential, we performed numerous sexual mating experiments between different clones since 1996. Clone A and clone D are very diverse among the identified clones. Karyotype of clone A individuals contains 156 chromosomes, which are composed of 36 metacentric(m), 54 submetacentric(sm), 36 subtelocentric(st), 24 acrocentric(t), and six small chromosomes, whereas karyotype of clone D individuals has 162 chromosomes, with 42 m, 54 sm, 36 st, 24 t, and six small chromosomes[14]. Through sexual reproduction mating between clones A and D, we found a few of fast-growing individuals. Significantly, the fast-growing individuals still possess its unisexual reproduction ability of gynogenesis. Thereby, the novel clone A+ was created by originally sexual mating between clone D female and clone A male, and rapidly multiplied up to several hundred millions by subsequent 7 generations of unisexual gynogenesis[8]. Here, we summarize formation process of the novel clone strain and attempt to identify the genetic organization and background by nuclear and cytoplasmic markers including chromosome number count, Cot-1 DNA fluorescent banding karyotype analysis, microsatellite electrophoretic pattern, AFLP profile, transferrin allele identification and mitochondrial genome sequence comparison, because most of them had been proven to be particularly valuable for clone discrimination, diversity evaluation and genealogical relationship analysis in several unisexual vertebrates, such as gynogenetic Amazon molly Poecilia formosa[23,24], gynogenetic Phoxinus eos-neogaeus[25], hybridogenetic Poeciliopsis[26-28], hybridogenetic water frog Rana esculenta [29], kleptogenetic salamanders[30-35], parthenogenetic lizards[36], and the gynogenetic gibel carp[19-22,37]. Based on these studies, we explore and discuss the significant implication for clonal diversity contribution.

Methods

Source of samples

Clone D female and clone A male were used as the maternal and paternal for the propagation experiments of gonochoristic sexual reproduction mating. Control gynogenetic individuals of clone D were inseminated by sperm from red common carp (Cyprinus carpio) to activate the eggs. As described previously, spawning was artificially induced by two intraperitoneal injections with a mixture of acetone-dried carp pituitary, hCG and LRH-A[38]. Ovulated eggs from clone D were divided into two parts and respectively inseminated with sperm from a clone A male and from a red common carp male. The produced offspring were respectively cultured in separate fishponds. After they reached to adults for one year, the phenotype, size and sex were determined. Subsequently, one fast-growing clone A-like individual was selected as the clonal maternal and its unisexual reproduction ability of gynogenesis was demonstrated by Xingguo red common carp (Cyprinus carpio) sperm stimulation, because the offspring are all-females, and identical to the maternal. Moreover, they have been proliferated by 7 successive generations of gynogenesis with Xingguo red common carp (Cyprinus carpio) sperm stimulation, and more than several hundred millions of clonal individuals have been produced. From the fifth generation, the clonal strain as a new aquatic variety has been introduced into more than 10 hatcheries, and their gynogenetic offspring were cultured throughout China. As a stock hatchery, the successively gynogenetic offspring of 7 generations have been maintained in the Guanqiao Experimental Station of the Hydrobiology Institute in Wuhan, Hubei Province, and all the samples for genetic analysis were collected from the Station. Besides the maternal and paternal, generally, ten individuals were randomly sampled from each analyzed clone and generation, and their fin tissues and blood cells were sampled for DNA extraction and DNA content measurement respectively. The fin tissues were stored in 95% ethanol at -20°C and blood cells were fixed by 70% ethanol and stored at 4°C.

Genomic DNA extraction and Cot-1 DNA isolation

Genomic DNA were extracted individually from fin tissues as described previously[17]. Cot-1 DNA was isolated as described previously[39], and labeled with biotin-16-dUTP by a nick translation reaction for FISH.

Chromosome preparation and fluorescence in situ hybridization

Chromosome metaphases were prepared from kidney cells of at least 5 individuals for each clone by the method of kidney cell-phytohemagglutinin (PHA) culture in vivo and then were counterstained with DAPI. Fluorescent in situ hybridization (FISH) of Cot-1 DNA was performed according to the method described previously[15], and images were acquired using a Leica inverted DMIRE2 microscope and a Leica LCS SP2 confocal image system (Leica, Germany).

Genetic marker application

Ten microsatellite primers were applied to amplify microsatellites in genomic DNA and PCR amplification was performed as the described method[19]. As for AFLP analysis, genome DNA was digested with EcoR I and Mse I and then selective amplification was performed using ten primer combinations[37]. Products of microsatellite and AFLP were separated using 6% denaturing polyacrylamide electrophoresis and were visualized by silver staining[40]. Transferrin alleles were amplified using special primer (Tf760 and TF1162) and then were cloned and sequenced as described previously[20,21]. To obtain the complete mitochondrial genome, twenty conserved primers[41] were used to amplify contiguous and overlapping fragments. According to the aligned sequences, another two pairs of primers were designed to amplify the fragments including the varied sites for further validation.

Sperm DNA content measurement

Clone A sperms were sampled from 8 individuals in reproduction season and 20 μl of them were respectively dropped into 1 ml of 70% ethanol. The fixed sperms at 4°C overnight were washed with PBS buffer for 3 times at 1000 rpm for 5 min. The pellets were resuspended in 0.5 ml of 0.5% pepsin in 0.1 M HCl. After the solution were incubated at room temperature for 30 min with gentle shaking, 100 μl of them was stained by adding 2 ml of propidium iodide solution (40 μgPI/ml) containing 4 kU/ml RNase (DNase-free ) at room temperature for 3 h in the dark. After the mixed cells were filtered through a special nylon mesh, the DNA content measurement was performed as described previously[42] by Phoenix Flow Systems (Beckman and Coulter). The blood cells of clone A were also sampled and detected from the same individuals as controls under the same conditions.

Data analysis

Microsatellite and AFLP amplification bands in the gibel carp clones were evaluated based on the electrophoretic patterns and the genetic diversities were calculated by Arlequin 3.11 software[43]. The gene regions of mitochondrial genome were identified with the homologous regions of the complete mitochondrial DNA sequence of triploid hybrids of tetraploid and Japanese crucian carp (GenBank accession number: AY771781). All sequences of transferrin and mitochondrial DNA were aligned and exported with CLUSTAL X [44] and Mega 4.0 software [45].

Results

Formation process of novel clone A+

Figure 1 shows the schematic diagram of a novel clone strain formation through using dual modes of gonochoristic sexual reproduction and unisexual gynogenesis. Firstly, female clone D with 162 chromosomes was sexually mated with male clone A with 156 chromosomes, and the propagated offspring was strictly cultured in separate fishponds. In the gonochoristic sexual reproduction mating, only 8.73% of the fertilized eggs can develop into adults, and the surviving offspring display three different phenotypes. Most of them (over 80%) are identical to the maternal clone D, some (about 15%) differ from clone D and clone A, and a few of individuals (less than 5%) look like the paternal clone A. Significantly, the few individuals grow markedly faster than the maternal clone D and the paternal clone A, and there existed about 61% females and 39% males in them. Then, one individual of them was selected as the clonal maternal and proliferated by 7 successive generations of gynogenesis with Xingguo red common carp (Cyprinus carpio) sperm stimulation. Xingguo red common carp is a regional variety of common carp in Xingguo County, Jiangxi Province and has more than 1300 years of culture history. Owing to the significant growth superiority, the clone has been approved by National Certification Committee for Aquatic Varieties as a new aquatic variety that is suitable for aquaculture in China. From the fifth generation, the novel clonal variety has been introduced into more than 10 hatcheries, and their gynogenetic offspring were cultured throughout China. At present, more than several hundred millions of clonal individuals have been produced in the hatcheries. As a novel clone strain, it has been designated as clone A+, but its genetic organization and the mechanism of clone formation have been still unknown.

Figure 1

A schematic diagram showing formation process of the novel clone A.

Identical chromosome number and banded karyotype to clone A

We firstly performed chromosome number count and Cot-1 DNA fluorescent banding karyotype analysis of the clone A+ metaphases sampled from the fifth generation, and compared with that of clone D and clone A. Figure 2 summarizes the data. In the clone A+ (Figure 2a) and A (Figure 2b), 38% and 36% examined metaphases contain 156 chromosomes, whereas in clone D (Figure 2c), 41% counted metaphases possess 162 chromosomes. Therefore, the modal chromosome number and representative metaphase (Figure 2d, 2e) of clones A+ are the same as those of clone A, which contain 156 chromosomes, whereas the modal chromosome number and representative metaphase of clone D (Figure 2f) possess 162 chromosomes.

Figure 2

Comparisons of chromosome number (a-c), representative metaphases (d-f) and triploid karyotypes of .

Moreover, the Cot-1 repetitive DNA fluorescence banding shows homologous chromosome-specifically patterns, and non-homologous chromosomes have different banded patterns. According to the various fluorescent intensities that locate at centromeric regions and other interstitial regions and the differential chromosome size and shape, we compiled the triploid karyotypes of clones A+, clone A and clone D with three homologous chromosomes. In all the analyzed three karyotypes, there are 9 small chromosomes that are wholly labeled by the Cot-1 repetitive DNA fluorescence, in which three of them are submetacentric chromosomes, and others are six small chromosomes. Also, the short arms of three same size subtelocentric chromosomes are entirely stained by the Cot-1 repetitive DNA fluorescence in all examined three clones. Significantly, clones A+ (Figure 2 g) and clone A (Figure 2 h) possess identical karyotypes that are comprised of 36 m, 54 sm, 36 st, 24 t, and six small chromosomes, whereas the karyotype of clone D (Figure 2i) is different from that of clones A+ and clone A, which contains 42 m, 54 sm, 36 st, 24 t and six small chromosomes.

Identical nuclear genome to clone A

Subsequently, the nuclear genome of clone A+ was further discriminated by three nuclear genome markers, such as microsatellite electrophoretic pattern, AFLP profile, and transferrin allele identification. Firstly, ten pairs of microsatellite primers isolated previously[19] were used to assess genetic differences and similarities of clone A+, clone A and clone D (Table 1). In total, 27 microsatellite alleles were equally detected from clone A+ and clone A, and 25 alleles were recognized from clone D. However, only three alleles were shared by clone A+, clone A and clone D from primer pairs 0039 and 0042. Figure 3a shows three typical microsatellite electrophoretic patterns amplified by primer pair YJ0001, YJ0033 and YJ0039, in which same allele patterns exist in clone A+ and clone A, and various patterns appear in clone D. Additionally, completely identical allele patterns are observed among individuals within each clone because of clonal reproduction of gynogenesis.

Table 1

Primer pair sequences and allele distribution of 10 microsatellites among clone A, clone D and clone A+ of gibel carp*

Locus (YJ)	Primer pair sequence	Annealing tempreture	Allele size (bp)
			Clone A	Clone D	Clone A+
0001	5'-CTGGCATGAAGACTGGCTC-3'5'-CAACAACACATATCAGCTCC-3'	53°C	84, 94, 100	88, 92, 96	84, 94, 100
0003	5'-TGAAGTTATTAGAAAGAGAG-3'5'-CTTGATGATGTCTATGTGTG-3'	53°C	270, 288, 302	260, 290, 292	270, 288, 302
0004	5'-CATAGAGGCGTTTCATAGAG-3'5'-CAGATAAATACAGTAAGCCA-3'	50 - 55°C	210, 220, 224	218, 222, 226	210, 220, 224
0005	5'-TAATAAGGTACATAGTCATAG-3'5'-GTCAGCCTCCACCACGAATC-3'	50 - 55°C	226, 230, 234	222	226, 230, 234
0033	5'-CGGACACAAGAACGCCAAC-3'5'-GGACTGGGCTGAAACTGATG-3'	50 - 55°C	172, 178	176	172, 178
0039	5'-GAAGAATACTTTATGACTGAGG-3'5'-GACCAAGACAGACAGCCCAG-3'	50 - 55°C	138, 150, 156	136, 138, 156	138, 150, 156
0042	5'-GGCCACCTACAGTATATGC-3'5'-GAAAACCAGGACCGACATG-3'	50 - 55°C	110, 114, 116	110, 112, 118	110, 114, 116
0040	5'-CCAGTATTAGGGAGCGTTC-3'5'-GTTTCGTCTTCACAATCAGAA-3'	50 - 55°C	124, 138	132, 148	124, 138
0010	5'-GATGGTTGTGCTGTGAGCT-3'5'-GAGTTCGTTTACATCTGGAC-3'	53°C	150, 156, 166	152, 160, 162	150, 156, 166
0022	5'-CACCAACTTTAGGCACATTTG-3'5'-CCAGACTCCCACGTCATG-3'	53°C	144, 156	140, 162, 170	144, 156
Total	10		27	25	27

*For each clone and each generation, 10 individuals were sampled and analyzed, and clone A+ samples of three generations, such as the third, fifth and seventh generation, were examined. No any individual difference was detected within each clone and among different generations.

Figure 3

Genetic discrimination of clone A(a) Three typical microsatellite electrophoretic patterns amplified by the primer YJ0001, YJ0033 and YJ0039. (b) Two representative AFLP patterns amplified by the primer combinations E2M5 and E3M3. (c) The aligned five transferrin alleles Tf1, Tf2, Tf3, Tf4 and Tf5 identified from three clones.

AFLP profiles produced by 10 primer combinations provided much more abundant genetic information than that of the above microsatellites, and each primer combination amplified an average of 76.5 bands, ranging from 64 to 86 bands (Table 2). As shown in Figure 3b, completely identical AFLP profiles are observed in the sampled individuals of clone A+ and clone A. However, the AFLP profiles of clone D are significantly different from that of clone A+ and clone A. The percentages of the shared bands among three clones range from 66.2% to 82.7%, but 79 clone A+ and clone A-specific bands and 88 clone D- specific bands are respectively detected from the amplified 10 AFLP profiles. The data again indicate that clone A+ nuclear genome is same to clone A, and different from clone D.

Table 2

Comparison of amplification bands by AFLP among clone A, clone D and clone A+ of gibel carp*

		No. of amplified bands					No. of clone-specific bands
No.	AFLP primers	A	A+	D	total	No. of shared bands	A and A+	D
1	E-ACAM-CAG	62	62	67	75	54	8	13
2	E-AGGM-CTT	70	70	69	83	56	14	13
3	E-AACM-CTG	61	61	62	74	49	12	13
4	E-AAGM-CTG	63	63	59	71	51	12	8
5	E-ACTM-CAT	71	71	77	81	67	4	10
6	E-ACGM-CTC	69	69	64	77	56	13	8
7	E-AACM-CTA	66	66	71	75	62	4	9
8	E-AAGM-CAA	81	81	80	85	76	5	4
9	E-AAGM-CTA	59	59	57	64	52	7	5
10	E-ACAM-CTA	75	75	71	80	66	9	5
total		677	677	677	765	589	79	88

*For each clone and each generation, 10 individuals were sampled and analyzed, and clone A+ samples of three generations, such as the third, fifth and seventh generation, were examined. No any individual difference was detected within each clone and among different generations.

Multiple ancient allelic lineages of transferrin had been identified in gibel carp[20,21]. To clone and identify all alleles of trnasferrin from the three clones, a pair of primers Tf760 (CTCCTCAAAGAGCCTCGC) and Tf1162 (CAAGGGCATCTGCTTCCT) was designed according to transferrin cDNA sequences (GenBank accession numbers: AY045574, AF457150, AF457151, AF518744, AF518745, AY323916) of gibel carp because high polymorphism sites had been revealed within this fragment among various transferrin alleles. Using this primer pair, we amplified and identified five different transferrin alleles (named Tf1, Tf2, Tf3, Tf4 and Tf5) from the three gibel carp clones (GenBank accession numbers: JF496199-JF496203). Tf2, Tf4 and Tf5 are 1170 bp in length, while the length of Tf1 and Tf3 are 1184 bp and 1177 bp respectively. From the aligned sequences, various mutations including deletion, transition and transversion are found among the five alleles (Figure 3c). Significantly, Tf1 and Tf2 appear specifically in clone A+ and clone A, whereas Tf4 and Tf5 exist specifically in clone D. Only Tf3 is shared by all of the three clones. The divergence between Tf3 and the other alleles is from 0.080 to 0.087, and it is much higher than that of the other alleles (Table 3).

Table 3

Genetic divergences among the five transferrin alleles

	Tf1	Tf2	Tf3	Tf4	Tf5
Tf1	0.000	0.041	0.087	0.042	0.040
Tf2	0.041	0.000	0.080	0.001	0.011
Tf3	0.087	0.080	0.000	0.081	0.081
Tf4	0.042	0.001	0.081	0.000	0.010
Tf5	0.040	0.011	0.081	0.010	0.000

Same mitochondrial genome to clone D

Because clone A and clone D have been revealed from numerous haplotypes and diverse clones of gibel carp to belong to the same haplotype without any variation in the high polymorphism mtDNA control region sequences[22], therefore, we sequenced the complete mitochondrial genome sequences of clone A+, clone A and clone D using a PCR-based method. The complete mitochondrial genomes of all three clones are 16580 bp long (D-clone, JF496197; A clone JF496198), and contain the identical gene complement (13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and a major non-coding D-loop control region) and gene order as found in most vertebrate mitochondrial genomes[46]. Only 4 single nucleotide polymorphisms are detected among the three complete mitochondrial genomes, and they exist respectively in ND2 gene (5864 nucleotide position), COX|gene (6650 nucleotide position), ATP6 gene (9699 nucleotide position) and a tRNA gene (15265 nucleotide position). Significantly, the 4 variation position nucleotides of clone A+ are all identical to that of clone D, and all different from that of clone A (Figure 4). Moreover, we designed two pairs of variation position-specific primers according to the flanking region sequences at 5864 and 6650 nucleotide variation positions, and used them to amplify and sequence the mtDNA fragments from 10 different individuals of each clone. The sequencing data indicate that both of the variation nucleotides in all individuals of clone A+ are same to that of clone D, and different from that of clone A, which further confirm the genetic consistency of mitochondrial genomes between clone A+ and clone D.

Figure 4

Mitochondrial genome sequence comparison of clone A. The complete genome sequences have been deposited in GenBank and the 4 single nucleotide polymorphisms are shown by the arrows at the corresponding positions.

Genetic uniformity among different generations

To further confirm the above data and to judge the validity, we randomly sampled 10 individuals in different generations including the third generation (F3), the fifth generation (F5) and the seventh generation (F7) from the stock hatchery and comparatively analyzed the genetic profiles with the original maternal clone D and paternal clone A by nuclear markers of microsatellite and AFLP and cytoplasmic marker of mtDNA sequence. As shown in Figure 5, the genetic profiles of microsatellite (Figure 5a) and AFLP (Figure 5b) are highly identical among the samples collected from the three different generations, and are same to that of the paternal clone A, while are significantly different from that of clone D. Simultaneously, no any various nucleotides are detected by the two mtDNA fragment sequences among 30 individuals of the three generations, and their sequences are completely identical to that of clone D (data not shown). The additional data verify the genetic uniformity among different generations of the clone A+ along with the seven generations.

Figure 5

Genetic uniformity among different generations of the clone A. (a) One typical microsatellite electrophoretic pattern of different individuals sampled from the third(F3), fifth (F5) and seventh (F7) generations and the original maternal clone D (D) and paternal clone A (A) that were amplified by the primer YJ0001 in. (b) One representative AFLP pattern of different individuals sampled from the third (F3), fifth (F5) and seventh (F7) generations and the original maternal clone D (D) and paternal clone A (A) that were amplified by the primer combination E3M3.

Clone A sperm DNA content

In comparison with only one peak of its own blood cells, clone A sperm DNA content measurement revealed one main peak and minor peak. As shown in Figure 6, the main peak fluorescence value is about half of the blood cells, whereas the minor peak value is basically same to that of the blood cells. Statistical data of 8 individuals indicate that the minor peak accounts for about 0.91% of the total sperms, implicating the existence of a few proportion of unreduced sperms.

Figure 6

Comparison of the DNA contents of sperms (a) and blood cells (b) of clone A.

Discussion

In this study, we reported a detailed formation process of the novel clone strain, and analyzed and identified the genetic organization and background by nuclear and cytoplasmic genome markers. The current data clearly indicate that the nucleus of clone A+ comes from that of clone A, because its chromosome number, Cot-1 repetitive DNA banded karyotype, and nuclear genome markers including microsatellite patterns, AFLP profiles and transferrin alleles are identical to that of clone A, and significantly different from that of clone D, whereas the cytoplasm comes from that of clone D, as the cytoplasmic marker-mtDNA genome is completely same to that of clone D, and different from that of clone A. Therefore, the major finding in this study is that the novel strain is clonal and represents a novel hybrid genotype with an entire nuclear genome from the paternal clone A and a mtDNA genome (cytoplasm) from the maternal clone D. And, the selected fast-growing individual as clonal maternal might be a nucleo-cytoplasmic hybrid between clone A and clone D, which might be formed via androgenesis of the clone A sperm in the clone D ooplasm.

As far as we know, androgenesis has been observed to occur naturally in interspecific hybrids of Sicilian stick insects[47,48] and triploid Asian clams (genus Corbicula)[49,50]. In fish, natural androgenesis has never been reported, but viable androgenetic fishes have been artificially induced in some ancient polyploid species, such as rainbow trout[51], common carp[52], loach[53,54], and sturgeon[55], and in a artificial allotetraploid hybrid of common carp and red crucian carp[41]. Generally, androgenesis involves loss or inactivation of the female genome and doubling of the paternal genome[56]. In the sexual mating between clone D female and clone A male, most of the fertilized eggs might undergo the fusion of clone D female pronucleus and clone A male pronucleus, and resulted in high mortality during embryo hatching and fingerling cultivation owing to the significant difference of chromosome number and karyotypes between them. The surviving offspring differ from clone D and clone A in morphological phenotype. Indeed, fusion failure of the female and male pronuclei might also occur in a few of the fertilized eggs, in which some developed into clone D-identical individuals, some developed into the clone A-like individuals. The clone D-identical individual might occur through the mechanism of gynogenesis, because the sperm was occasionally recognized as heterologous sperm, just like that from other fish species owing to the significant genetic difference between clones A and D. The clone A-like individual should be suggested to form via androgenesis, because it has an entire nuclear genome from the paternal clone A and a mtDNA genome from the maternal clone D. In contrast to gynogenesis, the entire female nuclear genome of clone D might be extruded, and the offspring might be produced by two possible mechanisms, such as by ploidy doubling of clone A sperm through inhibiting the first mitotic division or by inclusion of clone A unreduced sperm into clone D ooplasm. If the latter is true, it is impossible to differentiate females and males in the androgenetic offspring (Figure 1) according to the XY sex determination mechanism[57], even though the direct evidence for existence of a few proportion of unreduced sperms have been obtained by flow cytometry[42] in the paternal clone A sperms (Figure 6). Therefore, we suggest that the clone A-like androgenetic individuals should be formed by the mechanism of ploidy doubling of clone A sperm through inhibiting the first mitotic division, and thereby the nucleo-cytoplasmic hybrids can differentiate into females and males.

Significantly, the selected nucleo-cytoplasmic hybrid female still maintains its gynogenetic ability, and thereby forms a nucleo-cytoplasmic clone with over several hundred millions of clonal individuals. In the past 20 years, we also obtained one allotetraploid clone with additional genome incorporation[58,59] and one supertriploid clone with subgenomic incorporation[13] in the gibel carp. Actually, unisexual clones were early synthesized in 1973 from Poeciliopsis monacha-lucida[60], and the stably inherited clone with minichromosome incorporation was also found from gynogenetic Amazon molly (Poecilia formosa)[12,61]. However, this androgenetic nucleo-cytoplasmic hybrid clone should be the first case in vertebrates, and might be related to the high polyploidy level and the diverse reproduction modes[8]. In invertebrate freshwater clams of the genus Corbicula that reproduce by androgenesis, ploidy changes including diploid, triploid, and tetraploid exist widely among different species and populations, and the higher ploidy levels have been suggested to be related to the origin of androgenetic reproduction[49,62]. In interspecific hybrids of Sicilian stick insects, various reproduction modes including parthenogenesis, hybridogenesis and androgenesis have been revealed to be associated with hybridization, polyploidy and karyotype re-diploidization, and these modes have been believed to be a noteworthy ability to overcome species-specific reproduction isolation mechanisms[48,63]. Similarly, various ploidy including triploid, supertriploid and tetraploid and different reproduction modes including gynogenesis, androgenesis and sexual reproduction also exist extensively in gibel carp. Therefore, gibel carp and even other unisexual vertebrates may also possess the ability for overcoming species-specific reproduction isolation mechanisms. In unisexual salamanders (genus Ambystoma), Bogart et al. revealed a new reproductive mode of kleptogenesis by which unisexual salamanders steal sperm genome from donors of normally bisexual species[33,34]. Perhaps, various ploidy levels and different reproduction modes are common characteristics for the unisexual animal complexes, regardless of whatever they are invertebrates or vertebrates, and it is these characteristics to result in rapid incorporation of nuclear genome, chromosome, chromosomal segment or mtDNA genome and to cause rapid genetic changes. Moreover, they achieve remarkable ecological success, even though some are rare and localized, and others are abundant [2]. For example, in comparison with the related sexual species Carassius auratus, polyploid gibel carp have wide geographic distribution and occupy more multiple niches from northern Europe to Asia[8].

Indeed, the coexisting dual modes of unisexual gynogenesis and sexual reproduction in gibel carp might be indicative case of a change in reproductive mode from complete unisexual reproduction towards sexual reproduction[8]. In green toads (Bufo viridis)[64,65] and water frogs (Pelophylax esculentus)[66] involving various ploidy levels, sexual reproduction triploids have been observed. Recently, a male-biased mutant family that contains 97.2% males was also found from the gynogenetic progeny of a gibel carp clone[37]. Therefore, various ploidy levels and highly diverse reproductive modes might be the causes for novel clone formation and clonal diversity in unisexual vertebrates, and also make them excellent models for the studies of evolutionary genetics and ecology, both theoretically and empirically. Additionally, as the nucleo-cytoplasmic hybrid clone escapes the genetic and developmental destruction caused by drastic treatments of irradiation and physical shocks in induced androgenesis[55] and nuclear transplantation[67,68], the novel finding will be of great significance for exploiting the genetic breeding approaches in gibel carp[69]. Of course, this finding of nucleo-cytoplasmic hybrid clone can help us to understand some reasons for clonal and genetic diversity production of unisexual animals, but it remains unknown whether similar clones also exist in natural populations of gibel carp. Therefore, further genetic resource survey will be required for elucidating the evolutionary mechanisms of polyploid gibel carp through the nuclear and cytoplasmic markers.

Conclusions

In conclusion, we here identify a novel clone in polyploid gibel carp, and suggest that the novel clone might be a nucleo-cytoplasmic hybrid between the known clones A and D, because it originates from the offspring of gonochoristic sexual reproduction mating between clone D female and clone A male, and contains an entire nuclear genome from the paternal clone A and a mtDNA genome (cytoplasm) from the maternal clone D. It is suggested to arise via androgenesis by a mechanism of ploidy doubling of clone A sperm in clone D ooplasm through inhibiting the first mitotic division. Significantly, the selected nucleo-cytoplasmic hybrid female still maintains its gynogenetic ability, and thereby forms a novel clone with over several hundred millions of clonal individuals proliferated by 7 successive generations of gynogenesis. The androgenetic nucleo-cytoplasmic hybrid clone should be the first case in vertebrates, and might be related to the high polyploidy level and the diverse reproduction modes in gibel carp. Based on the present and previous findings, we discuss the association of rapid genetic changes and high genetic diversity with various ploidy levels and multiple reproduction modes in several unisexual and sexual complexes of vertebrates and even other invertebrates. Therefore, the novel finding will be of great significance not only for exploiting the genetic breeding approaches in gibel carp but also for the studies of evolutionary genetics and ecology in other unisexual animals.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JFG, ZWW and LZ conceived and designed the experiments. HPZ contributed to the chromosome data, DW and ZWW contributed to the AFLP data, FFJ and ZWW contributed to the transferrin allele data, and WG and ZWW contributed to the microsatellite data. ZWW analyzed the data. JFG and ZWW wrote the paper. All authors have read and approved the final manuscript.