- Research article
- Open Access
Oligo-painting and GISH reveal meiotic chromosome biases and increased meiotic stability in synthetic allotetraploid Cucumis ×hytivus with dysploid parental karyotypes
BMC Plant Biology volume 19, Article number: 471 (2019)
Meiosis of newly formed allopolyploids frequently encounter perturbations induced by the merging of divergent and hybridizable genomes. However, to date, the meiotic properties of allopolyploids with dysploid parental karyotypes have not been studied in detail. The allotetraploid Cucumis ×hytivus (HHCC, 2n = 38) was obtained from interspecific hybridization between C. sativus (CC, 2n = 14) and C. hystrix (HH, 2n = 24) followed by chromosome doubling. The results of this study thus offer an excellent opportunity to explore the meiotic properties of allopolyploids with dysploid parental karyotypes.
In this report, we describe the meiotic properties of five chromosomes (C5, C7, H1, H9 and H10) and two genomes in interspecific hybrids and C. ×hytivus (the 4th and 14th inbred family) through oligo-painting and genomic in situ hybridization (GISH). We show that 1) only two translocations carrying C5-oligo signals were detected on the chromosomes C2 and C4 of one 14th individual by the karyotyping of eight 4th and 36 14th plants based on C5- and C7-oligo painting, and possible cytological evidence was observed in meiosis of the 4th generation; 2) individual chromosome have biases for homoeologous pairing and univalent formation in F1 hybrids and allotetraploids; 3) extensive H-chromosome autosyndetic pairings (e.g., H-H, 25.5% PMCs) were observed in interspecific F1 hybrid, whereas no C-chromosome autosyndetic pairings were observed (e.g. C-C); 4) the meiotic properties of two subgenomes have significant biases in allotetraploids: H-subgenome exhibits higher univalent and chromosome lagging frequencies than C-subgenome; and 5) increased meiotic stability in the S14 generation compared with the S4 generation, including synchronous meiosis behavior, reduced incidents of univalent and chromosome lagging.
These results suggest that the meiotic behavior of two subgenomes has dramatic biases in response to interspecific hybridization and allopolyploidization, and the meiotic behavior harmony of subgenomes is a key subject of meiosis evolution in C. ×hytivus. This study helps to elucidate the meiotic properties and evolution of nascent allopolyploids with the dysploid parental karyotypes.
Interspecific hybridization and allopolyploidization frequently result in a ‘genomic shock’ that causes rapid genetic and epigenetic changes, due to the merging of two or more divergent and hybridizable genomes [1,2,3]. The meiosis abnormality of newly formed allopolyploids, as an immediate consequence of allopolyploidization, results in extensive abnormal chromosome pairing, imbalanced chromosome segregation and karyotype variations [4,5,6]. Cytogenetic studies have shown that different subgenomes have distinct meiotic behavior stabilities, and individual chromosomes have biases for chromosome loss and/or gain [7,8,9]. The extent of meiotic instability, genome structural and genetic/epigenetic changes may vary considerably in different polyploid species, suggesting that these changes depend on the origin and evolutionary differences between parental species . To date, most studies on genomic variations and meiosis evolution have been conducted in nascent allopolyploids with analogous subgenomic karyotypes, such as same or close subgenomic chromosome numbers . However, less attention has been devoted to synthetic allopolyploids with dysploid parental karyotypes.
Molecular cytogenetics is indispensable for studying the evolution of the polyploid genome, and it can intuitively visualize the dynamics of the genome. Most significantly, chromosome painting based on fluorescence in situ hybridization (FISH) has been verified as a powerful tool for identifying chromosomes and investigating chromosome rearrangements during evolution [12, 13]. Unfortunately, the issue concerning chromosome identification and tracking meiotic behavior based on chromosome painting has been less thoroughly explored due to complicated chromosome synteny and the absence of suitable probes, especially in interspecific hybrids and allopolyploids. New strategies for developing chromosome painting probes in plants have been inspired by the development of plant whole genome sequencing projects and technical advances in DNA synthesis. Single-copy sequence-based chromosome painting has been developed and applied to identify individual chromosomes, track chromosome pairing and detect rearrangements in cucumber, potato, maize and poplar [14,15,16,17,18,19]. This approach greatly facilitates the identification of homoeologous chromosomes, diagnosis of chromosome abnormalities, and tracking of pairing behavior in interspecific hybrids and allopolyploids.
The allotetraploid C. ×hytivus J. F. Chen & J. H. Kirkbr. (HHCC, 2n = 4x = 38) was synthesized through interspecific hybridization between cucumber (C. sativus L. ‘BejingJietou’, CC, 2n = 2x = 14) and its sister species C. hystrix (HH, 2n = 2x = 24) followed by chromosome doubling [20, 21] (Fig. 1). Like other newly formed allopolyploids, the allotetraploid C. ×hytivus exhibits rapid genetic changes and extensive meiotic instability [11, 22, 23]. The cucumber and C. hystrix were derived from a common ancestral species via decreasing dysploidy (n = 12 to n = 7) approximately 4.6 million years ago [24, 25]. Different evolutionary fates have produced differences in karyotypes, genetics, cytological characteristics and specific traits between cucumber and C. hystrix. The interspecific hybrid and allotetraploid C. ×hytivus with distinctive subgenomic karyotypes provided a new means of revealing the genetic relationships between parental species compared to neo-allopolyploids with analogous subgenomic karyotypes, determining whether the synthetic allotetraploid C. ×hytivus exhibits novel meiotic properties and elucidating how to achieve stable meiosis.
We developed two oligo-probe pools from cucumber chromosomes 5 (C5) and 7 (C7) to identify homoeologous chromosomes, detect chromosomal rearrangements and track individual chromosome meiotic pairing in F1 hybrid and its derived allotetraploid. We analyzed the meiotic properties of the F1 hybrids and two different generations of C. ×hytivus. The meiotic behavior of individual chromosomes was significantly different in homoeologous pairing and univalent formation. The meiotic behavior of these chromosomes (C5, H9, H10 and C7, H1) in S14 generation was improved in the absence of extensive chromosome reshuffling located on the given chromosomes. GISH experiments were performed to investigate the possible differences of two subgenomic meiotic behaviors in the S4 and S14 inbred families. The meiotic behavior biases between two subgenomes were observed in S4 generation, especially univalent formation, chromosome lagging and asynchronous meiotic rhythm. After 10 generations of self-pollination, we document the increased meiotic stability in S14 generation, including synchronous meiosis, increased normal bivalents and reduced lagging chromosome. In addition, many autosyndetic pairings of H-chromosomes were observed in the F1 hybrids. Given the clear parental genetic background, the synthetic allotetraploid C. ×hytivus could serve as a uniquely traceable system to explore allopolyploidization and meiotic evolution of allopolyploids with distinctive parental karyotypes.
Identification of individual chromosomes using oligo-painting
To study the meiotic behavior of individual chromosomes in interspecific hybrid and allotetraploid C. ×hytivus, we developed two oligo-painting probe pools of C5 and C7 based on the cucumber genome (see Materials and Methods). The average oligo densities of the two painting probes were 0.86 (C5) and 1.32 (C7) oligos per kilobase, respectively (Additional file 1: Figure S2). As expected, both probes produced bright and nearly uniform FISH signals on mitotic metaphase and meiotic pachytene chromosomes of cucumber (Fig. 2a, b). No extra oligo-FISH signals were detected on any other chromosomes (Fig. 2c), indicating that these two probes were highly specific to C5 and C7. There was an unambiguous and large signal gap on the pachytene chromosomes of C5 and C7 (Fig. 2d), which contain filtered centromeric regions and 45 s rDNA loci.
The homoeologous chromosomes of C5 and C7 were identified by C5- and C7-oligo probes, Type III and 45 s rDNA probes at mitosis metaphase of C. hystrix, F1 hybrid and allotetraploid C. ×hytivus (Fig. 3). Oligo-FISH signals from C5-oligo probes were unambiguously detected on two C. hystrix chromosomes, 9 (H9) and 10 (H10), which were distinguished by oligo-FISH intensity and 45 s rDNA signal located on the pericentromeric region of H10 (Fig. 3a, b, c) [11, 25]. The C. hystrix chromosome 1 (H1) was painted by C7-oligo probes, which were distinguished by cucumber-specific Type III centromere probes (Fig. 3e, f, g). Then, the two oligo-probes were hybridized to meiotic pachytene chromosomes of S14 generation of C. ×hytivus (Fig. 3d, h). No unambiguous signal gaps were observed on meiotic pachytene chromosomes. Consequently, these chromosome painting patterns can rapid identify individual chromosomes of C5, C7, and H1, H9, H10, diagnose karyotypic variation and trace meiotic pairing of these chromosomes in F1 hybrids and allotetraploid C. ×hytivus.
Homoeologous pairing and autosyndetic pairing in F1 hybrid
Interspecific hybrids provide a new route for studying the genetic relationships among parental chromosomes through the assessment of meiotic chromosome pairing. We investigated the meiotic pairing of C5 and C7 and their homoeologous chromosomes after interspecific hybridization based on the oligo-FISH patterns described above (Fig. 4). We analyzed 132 and 106 well-resolved PMCs at meiotic metaphase I (MI) taken from four F1 plants to trace chromosome pairing using C5 and C7-oligo probes, respectively (Table 1). The C5, H9 and H10 were completely unpaired as univalents in 74 (56.1%) PMCs (Fig. 4a). We found that 30 (22.7%) PMCs harbored intersubgenomic bivalents of C5 and H9 (Fig. 4b), and 17 (12.9%) PMCs harbored bivalents of C5 and H10 (Fig. 4c). The C7 and H1 were completely unpaired as univalents in 87 (82.1%) PMCs (Fig. 4e). The bivalents of C7 and H1 were found in 19 (17.9%) PMCs (Fig. 4f, white arrow). Similarly, one to three intergenomic bivalents were detected in 45 (42.5%) PMCs (Fig. 4g, white arrows). The quadrivalent and trivalents were also detected in one and eight PMCs, respectively (Fig. 4f, g, red arrows).
Interestingly, we found that 10 (7.6%) PMCs contained the autosyndetic pairings of H9 and H10 (Fig. 4d, inset). The pairing of H9 and another chromosome was detected in one PMC (Fig. 4b, box), whereas the pairing of H10 and another chromosome was found in two PMCs (Fig. 4b, red arrow). Two H-chromosomes (chromosomes from the H-genome) were autosyndetic pairing as a bivalent configuration in 27 (25.5%) PMCs (Fig. 4h). However, autosyndetic pairing of C-chromosomes was not observed in the 106 PMCs investigated. Considering each subgenome as a whole, the H-subgenome showed a higher autosyndetic pairing frequency. These results suggest that the meiotic behavior of the H-genome is more susceptible to the shock of interspecific hybridization than is the C-genome.
Two translocations and possible cytological evidence detected based on oligo-painting
In our previous research, aneuploids or large-scale chromosomal rearrangements were not detected in 15 individuals from the S13 generation of allotetraploid C. ×hytivus based on the fosmid-FISH results . This result might be attributed to few the individuals investigated, and the fosmid clones did not cover potential chromosomal rearrangement regions. We analyzed chromosomal variations via C5- and C7-oligo-painting in eight S4 and 36 S14 allotetraploid individuals. No visible chromosomal rearrangement events were observed in any individuals, except one S14 plant, in which there were two large translocations detected (Fig. 5, Additional file 1: Figure S3). Extra C5-oligo signals were detected on the long arm and near the centromeric region of two pairs of chromosomes with large 45 s rDNA locus (Fig. 5a, b). According to the karyotype of C. ×hytivus , we selected two fosmid clones (chr2–41 and chr4–37) to identify the two chromosomes with extra C5-oligo signals (Fig. 5c). We did not detect extra C7-oligo signals localized on other chromosomes, indicating that the chromosomes C7 and/or H1 do not experience nonhomologous rearrangements with other chromosomes in the analyzed plants.
Notably, mispairing of C5 and other chromosome was observed in two PMCs from the S4 generation (Fig. 6d, e). We performed successive FISH experiments using Type III probes to identify the two chromosomes. The FISH results indicated that the chromosome with strong 45 s rDNA loci belonged to the C-subgenome (Fig. 6d, inset), and the other chromosome belonged to the H-subgenome (Fig. 6e, inset). These abnormal pairings may cause potential chromosomal rearrangements that could be stably inherited into the subsequent generations. Three strong 45 s rDNA loci were located on three C-subgenome chromosomes C1, C2 and C4 in allotetraploid C. ×hytivus . One of these chromosomes likely formed a mispairing configuration with C5 and caused chromosomal translocations, which provided possible evidence for the two translocations in one S14 individual (Fig. 5).
Meiotic behavior biases of individual chromosomes
To explore possible meiotic pairing differences among chromosomes at meiosis I of allotetraploid C. ×hytivus, we performed oligo-painting on meiotic chromosomes at meiotic pachytene, MI and anaphase I (AI). We analyzed chromosome pairing in a total of 76 PMCs at the meiotic pachytene stage from the S4 generation (Additional file 1: Figure S4). The C5 and H9/H10 partially paired in 24 (31.6%) PMCs (Additional file 1: Figure S4a, b, c). Completely unpaired C7 chromosome strings were detected only in three PMCs (Additional file 1: Figure S4d), and the long arm of C7 and H1 were unpaired in five (6.6%) PMCs (Additional file 1: Figure S4e, g) and six (8.7%) PMCs (Additional file 1: Figure S4f), respectively. However, the C7 and H1 partially paired or fully paired were not observed in all PMCs, although C7 and H1 are highly syntenic. In addition, we observed a considerable amount of unpaired chromosome fragments and single chromosome strings in the PMCs of S4 generation (Additional file 1: Figure S4, red arrows), which may result in potential chromosome rearrangement events and univalents.
We examined the paring behaviors of five chromosomes (C5, H9, H10 and C7, H1) at MI and AI of meiosis in eight S4 individuals (Fig. 6). A total of 123 and 112 well-resolved MI PMCs was analyzed to trace the meiotic paring behavior of C5, H9, H10 and C7, H1, respectively (Table 1). The results showed that 61% of PMCs harbored exclusive homologous bivalents of the C5, H9 and H10 (Fig. 6a), and 90.2% of PMCs contained exclusive homologous bivalents of C7 and H1 (Fig. 6g). Then, we calculated the frequency of univalents and homoeologous pairings of five chromosomes at MI (Table 1). We found that these chromosomes showed clear differences in their propensities to be in univalent state and homoeologous pairing (Fig. 8f and Table 1). Specifically, among the five chromosomes, H10 showed the highest frequency of univalent formation, whereas both C7 and H1 showed the lowest frequency of univalent formation. Homoeologous pairing of C5 and H10 was observed in 11 (8.94%) PMCs (Fig. 6c). However, the homoeologous pairing of C5/H9 and C7/H1 as bivalents or multivalents was not observed in all PMCs. Therefore, these chromosomes are biases for homoeologous pairing and univalent formation. Unequal segregation of H9 and H10 was detected at AI (Fig. 6f). Similarly, lagging H-chromosomes were also detected at AI (Fig. 6i).
Meiotic chromosome biases of the subgenome cause the low fertility of allotetraploid C. ×hytivus
We investigated possible differences in the meiotic behavior between two subgenomes in MI and AI PMCs taken from eight S4 and ten S14 generations through GISH experiments (Table 2). We observed many PMCs of S4 generation with clear abnormalities, including the presence of asynchronous meiosis at MI (Fig. 7a, b), univalents (Fig. 7c), intergenomic pairings (Fig. 7d, Additional file 1: Figure S5) and lagging chromosomes at AI/TI (Fig. 7e, f). As expected, most of the chromosomes still maintained bivalent configurations in each PMC. Notably, most C-bivalents could be precisely and tightly aligned on the equatorial plate, whereas the H-bivalents seemed irregularly dispersed or lagged behind the C-bivalents at MI (Fig. 6g, Fig. 7a). Some of the H-bivalents were just arrived at the equatorial plate when the C-bivalents had already begun to segregate (Fig. 7b). Furthermore, 73.4% of the MI PMCs showed lagged H-bivalents, which was significantly higher than the 15.4% of PMCs with lagged C-bivalents in S4 generation (Fig. 8g). Similarly, the frequency of lagged H-bivalents was also significantly higher than that of C-bivalents in the S14 generation (Fig. 8g). These results indicated that the meiosis of two subgenomes was significantly asynchronous; that is, the H-subgenome required more time to complete the meiosis process than the C-subgenome in one nucleus.
The H-subgenome showed significantly higher univalent formation and lagging frequencies than the C-subgenome in both the S4 and S14 generations (Table 2). The results indicated that the propensities of the two subgenomes for meiotic behavior were significantly different. The meiotic behavior of the H-subgenome is more responsive to allopolyploidization, which eventually contributed to extensive H-chromosome lagging at subsequent stages (Table 2, Fig. 7e, f). Lagged chromosomes caused male gametes with incomplete chromosome complements, which resulted in the low pollen fertility of allotetraploid C. ×hytivus. In conclusion, meiotic behavior biases of two subgenomes are the primary factor for the low fertility of allotetraploid C. ×hytivus.
Increased meiotic stability by harmony between two subgenomes
We examined the meiotic behavior of C5, H9, H10 and C7, H1 in 137 and 104 well-resolved PMCs, respectively, taken from ten S14 plants with no nonhomologous chromosomal rearrangement based on C5/C7 oligo-painting (Table 1). Compared with the S4 generation, the frequency of homologous bivalents of the five chromosomes was significantly improved in the S14 generation (Table 1, Fig. 8a, b). Correspondingly, the univalent frequencies of H1, C5, H9 and H10 were significantly lower than those of S4 generation (Fig. 8f). The chromosome H10 still had the highest univalent frequency among the five chromosomes (Fig. 8f). The univalent frequencies of C7 between the S4 and S14 generations were statistically insignificant (Fig. 8f). These results indicated that the meiotic stability of these chromosomes was improved in the absence of extensive chromosome reshuffling located on the given chromosomes.
In analyzed 131 well-resolved MI PMCs from the S14 generation, 98 (74.8%) of PMCs contained 19 exclusive homologous bivalents, which was significantly higher than 49 (34.3%) detected in the S4 generation (Table 2). The frequencies of meiotic abnormalities in the S14 generation were remarkably reduced compared with the S4 generation (Table 2, Fig. 8c), including univalents, intergenomic pairing and chromosome lagging. Importantly, the frequency of MI PMCs containing asynchronous meiosis was significantly reduced in the S14 generation (Table 2). In MI PMCs with asynchronous meiosis, only a few (one to three) H-bivalents have not yet reached the equatorial plate in S14 generation (Fig. 8d, e), which was significantly reduced compared with S4 generation with one-seven lagged H-bivalents (Additional file 1: Figure S6). The results showed that two subgenomes have achieved a certain degree of synchronous meiosis at MI. Predictably, the pollen fertility of S14 plants was prominently higher compared to S4 plants, but still remained at a relatively low level (less than 50% on mean) (Additional file 1: Figure S1). These results demonstrated that meiosis stability can be increased through continuous selfing, and harmony between subgenomes.
Meiotic affinities of homoeologous chromosomes
Evaluating homoeologous pairing in interspecific hybrids can provide important information for the affinities among homoeologous chromosomes . The pairing preferences between chromosomes and their potential partners pay less attention in interspecific hybrids and allopolyploids due to the absence of FISH technology to identify individual chromosomes. The cucumber chromosome C5 originated from fusion of inferred chromosomes AK9 and AK10 (ancestor karyotype) similar to H9 and H10 of C. hystrix. Cucumber chromosome C7 and C. hystrix chromosome H1 are evolved from the ancestor chromosome AK1. These chromosomes had experienced strikingly different evolutionary fates [25, 27]. In the analyzed PMCs of F1 hybrids, the homoeologous pairing of C5 and H9 (22.7%) had a higher frequency than pairing of C5 and H10 (12.9%), which may be explained by the collinearity structure of C5 subtelomeric regions with H9 . When an individual chromosome has multiple homoeologous chromosomes, these chromosomes may be sorted by preferential chromosome pairing, because the recognition among chromosomes occurs preferentially between subtelomeric regions [28, 29]. The homoeologous pairing of C7 and H1 was observed only in 19 (17.9%) MI PMCs of the F1 hybrid, even though C7 is highly conserved with H1 of C. hystrix. Therefore, the pairing affinities among homoeologous chromosomes may be attributable to accumulated changes in the chromosomal structure during evolution.
Interestingly, 7.6% of PMCs containing autosyndetic pairing of H9 and H10 were observed in the F1 hybrid. H9 and H10 paired with other H-chromosomes only in one and two PMCs, respectively. Pairing of H1 and other H-chromosomes were not observed in all PMCs. The autosyndetic pairings within one genome implied that segmental homology and/or common origin between the involved chromosomes . The results indicated that H9 and H10 have a certain meiotic affinity and that H1 has a lower meiotic affinity with other H-chromosomes. Meiotic behavior of these chromosomes is dramatically different in response to interspecific hybridization, which appears to be related to the evolutionary fates of these chromosomes. The AK9 and AK10 have two distinct evolutionary fates: evolving into one chromosome (C5) after the fusion event or forming present the H9 and H10 . The meiotic affinity between H9 and H10 may increase our understanding of the fusion event of AK9 and AK10.
Among the five chromosomes analyzed in two different generations of allotetraploid C. ×hytivus, H10 showed a higher univalent frequency than other chromosomes. H10 with an asymmetric centromere position (acrocentric) may form an unstable bivalent, because a certain arm length is required for stable bivalent formation . Unpaired H10 may disturb the pairing of C5, causing C5 to have a relatively high univalent frequency and pairing with H10. Our results showed that the univalent frequencies of C7, H1 and H9 are relatively lower with no C7-H1 and C5-H9 pairings. These results suggest that C7, H1 and H9 can normally recognize homologous chromosomes and form stable bivalents. Homoeologous pairing and recombination are attributed to not only the chromosome structure but also the interaction of genetic factors, such as ph 1 in allohexaploid wheat, PrBn in allotetraploid Brassica napus, and BYS in allotetraploid Arabidopsis suecica . However, the pairing control locus may be insufficient to ensure the exclusive diploid-like meiotic chromosome pairing of newly formed allopolyploids, and additional modifications are still needed in meiotic stabilization . In addition, the timing of chromosome condensation also affects chromosome pairing .
Subgenome biases in meiotic chromosome behavior
Interspecific hybrids and allopolyploids face many challenges caused by large-scale conflict among divergent parental genomes, commonly referred to as ‘genomic shock’ associated with a de novo merger of two or more divergent parental genomes. Because unique evolutionary fates shaped each parent, producing species-specific karyotypes, genetic and epigenetic differences . In interspecific hybrids and allopolyploids, subgenome dominance appears to be a common result that responds differentially to genomic shock, such as subgenome biases in epigenetic regulation, gene expression and homoeologous exchanges . Similarly, the meiotic chromosome behavior of subgenomes also has biases; for example, a subgenome exhibits higher univalent and chromosome lost frequencies or maintains relatively stable meiotic behavior.
The autosyndetic pairing, which is involved in chromosomes from the same genome, was observed in several Paphiopedilum and Brassica interspecific hybrids [33, 34]. Similarly, many PMCs (25.5% of 106 PMCs) containing H-chromosome autosyndetic pairings were observed in the F1 hybrid. However, autosyndetic pairing was not observed between any two C-chromosomes. These results suggest that the meiotic pairing of two genomes has biases in the F1 hybrid. The cucumber and C. hystrix have derived from an ancestral species via the mechanisms of decreasing dysploidy (n = 12 to n = 7), leading to a high level of evolutionary differences between two genomes [24, 25]. The magnitude of genomic changes in response to genomic shock appears to be correlated with the degree of subgenome evolutionary differences . The autosyndetic pairing of H-chromosomes indicates that the meiotic behavior of the H-genome may be more susceptible by interspecific hybridization than the C-genome in the F1 hybrid. Therefore, the different responses of the H- and C-genomes to genomic shock or interspecific hybridization may be attributed to the evolutionary differences between the two parental species.
Extensive cytological studies have shown that homologous bivalents are the dominant pairing configuration in the meiosis of allopolyploids, and individual chromosome/subgenomes have biases in terms of meiotic behavior in some allopolyploids [5, 7, 9]. The meiotic behavior biases among subgenomes were observed in synthesized allohexaploid wheat and Brassica. The distinct subgenome stability was B > A > C in Brassica allohexaploid and D > A > B in synthesized wheat allohexaploid; that is, the C-subgenome of Brassica allohexaploid and the B-subgenome of wheat allohexaploid showed the highest frequencies of univalent formation and chromosome loss [7,8,9]. Indeed, the meiotic behavior biases between two subgenomes were also observed in allotetraploid C. ×hytivus, including asynchronous meiosis, univalent formation, and lagging chromosomes. The meiotic stability is C-subgenome > H-subgenome in C. ×hytivus, the H-subgenome showed higher univalents and chromosome lagging frequencies. However, no obvious asynchronous rhythm was observed in some newly formed allopolyploids, such as synthetic wheat allotetraploid and allohexaploid [9, 35], Brassica allotetraploid and allohexaploid [7, 34]. The asynchronous meiotic rhythm between subgenomes, originating from differences in progression through the meiotic cell cycle between parental subgenomes [36, 37], may be attributed to distinctive karyotypes and evolutionary differences of parental species. The H-subgenome requires a long duration of the cell cycle to complete meiosis, and is a major cause of the high frequency of H-chromosome lagging in Cucumis allotetraploid . Only 9.9% of MI PMCs and 9.7% of AI PMCs contained intergenomic pairing and C-chromosome lagging in S14 generation, respectively. Conversely, 55.8% of AI PMCs contained H-chromosome lagging in S14 generation, indicating that the relatively unstable meiotic behavior of the H-subgenome is primarily responsible for the low fertility of allotetraploid C. ×hytivus. The meiotic behavior of the two subgenomes is dramatically different both in the F1 hybrid and C. ×hytivus, likely associated with the differences of karyotypes, inherent cytological characteristics and evolutionary history between two parental species.
Meiotic harmony of two subgenomes is a key subject of meiotic stability in synthetic Cucumis allotetraploid
Extensive chromosome reshuffling and aneuploidy were observed in wheat allopolyploids [8, 35], Brassica allopolyploids  and Tragopogon neo-allopolyploids [40, 41]. However, no aneuploidy was detected in the analyzed allotetraploid C. ×hytivus, possibly because that the zygotes containing aneuploidy cannot develop into full shape seeds, and then were neglected when collecting seeds or selecting for germination. Notably, two large translocations and possible evidence were detected using oligo-based chromosome painting. However, the rearrangement events of C7/H1 and other chromosomes were not detected in this study. Our results imply that the potentials of individual chromosomes for chromosomal rearrangements and/or structural variations may be different during the process of polyploid evolution. Recently, a chimeric gene, responsible for delayed leaf maturation in the allotetraploid C. ×hytivus, was cloned, which originated from homoeologous recombination . This result implies that the Cucumis allotetraploid experienced homoeologous exchanges or recombination. We also observed that many PMCs contained intergenomic pairing, which was the cytological basis of homoeologous exchanges in Cucumis allotetraploid. It is possible that the rearranged fragments and homoeologous exchanges containing beneficial genes from wild relatives, which can be traced by oligo-FISH, were transferred into crop cultivars through introgression breeding.
Given the subgenome biases for meiotic chromosome behavior, how to coordinate the meiotic behavior between two subgenomes has become a key subject of meiotic stabilization in allotetraploid C. ×hytivus. We analyzed the behavior of five chromosomes (C5, H9, H10 and C7, H1) in two generations, and found that the abnormal meiotic behavior of these chromosomes was improved in the absence of these chromosome reshuffling. Meiotic stability is significantly improved through continuous selfing process, including reduced univalent formation, intergenomic pairing and chromosome lagging. The most obvious evidence is that the two subgenomes achieved a certain degree of synchronization at MI of meiosis in the S14 generation (Additional file 1: Figure S6). The number of H-bivalents that did not reach the equatorial plate was significantly reduced in S14 generation compared with S4 generation. The meiotic stabilization mechanisms possibly accelerated in several early generations of C. ×hytivus, but still require many generations to stabilize. Similar to the situation in other neo-allopolyploid plant taxa, the synthesized allotetraploid C. ×hytivus also experienced rapid genetic and epigenetic changes [22, 23, 43, 44]. These changes laid the foundation for meiotic harmony between two subgenomes. The results of this study will provide an interesting case to explore the meiotic evolution in allopolyploids with distinctive subgenomic karyotypes.
Taken together, the results of this study indicate that oligo-FISH is a powerful and efficient technique for chromosome identification, chromosome variation diagnosis, and tracing of meiotic behavior, especially in polyploids with complex genomes. In interspecific F1 hybrids and allotetraploids, the meiotic biases of five chromosomes (C5, H9, H10 and C7, H1) for homoeologous pairing and univalent formation may be attributed to these chromosomes pairing affinities and structure. The difference in the meiotic properties of the two genomes for autosyndetic pairing, asynchronous meiosis, univalent formation and chromosome lagging indicates that the two genomes differentially respond to interspecific hybridization and allopolyploidization. These differential responses are associated with differences in karyotypes, inherent cytological characteristics and evolutionary history between two parental species. Thus, the harmony of difference meiotic behavior between two subgenomes has become a key subject of meiotic stabilization in allotetraploid C. ×hytivus. Indeed, the meiotic stability increased in S14 generation, including synchronous meiosis, reduced univalent formation and chromosome lagging. Our analysis provided new insights into the meiotic properties and meiosis stabilization of nascent allotetraploids with dysploid parental karyotypes.
The plant materials used for this study included two diploid parents (C. sativus ‘CC3’ and C. hystrix), their interspecific F1 hybrid and synthetic allotetraploid C. ×hytivus (Fig. 1) [20, 21], provided by the state key lab of Cucurbit Genetics and Germplasm Enhancement of Nanjing Agricultural College. All different generations of C. ×hytivus were obtained from the same 3th (S3) inbred family, which was obtained from a single S0 plant of C. ×hytivus by selfing. The seeds (> 20) were randomly selected from each generation seed set for breeding the next generation by selfing. We currently obtain 14th generation through continuous selfing. Four different generations (eight individuals of 4th (S4), five individuals of 8th (S8), six individuals of 11th (S11) and 36 individuals of 14th (S14)) of C. ×hytivus were randomly chosen from last generation seed sets for pollen viability, karyotyping and meiosis analysis. The pollen viability analysis showed that the pollen fertility of S14 generation was significantly higher than that of S4 generation, and had a statistically highly significant difference (P = 0.0009). Therefore, we selected all S4 and S14 plants to conduct mitotic analysis based on C5 and C7-oligo probes. Four F1, eight S4 and ten S14 (randomly selected from 36 S14 plants) individuals were used for meiotic analysis based on oligo-FISH and GISH (Additional file 2 Table S1). All the materials were grown in a greenhouse at Baima Teaching and Research Base Modern Agricultural Science and Technology Zone of Nanjing Agricultural University, Nanjing, China.
Root tips and young male flower buds of all materials were collected and fixed in Carnoy’s solution at 4 °C for at least 1 day. Specifically, the root tips of 8 S4 and 36 S14 individuals were separately collected to detect chromosomal variations via chromosome painting. The young male flower buds of four F1 individuals were collected and mixed together multiple times during flowering. The eight S4 individuals and ten S14 individuals were divided into three groups. The young male flower buds of each group were respectively collected multiple times. One anther from each flower bud was dissected to examine the developmental stage of PMCs. The remaining anthers from the same flower bud were collected into the corresponding groups for chromosome preparation. The procedure of chromosome preparations was performed as described previously [11, 18] with some modifications. The fixed root tips were digested with an enzyme mixture containing 4% cellulose R-10 (Yakult), 2% pectinase (Sigma-Aldrich) and 0.1% pectolase (Yakult) in 0.01 M citrate buffer (pH = 4.8), at 37 °C for 40–60 min. The anthers were collected and digested using enzyme mixtures, including 4% cellulose R-10 (Yakult), 4% pectinase (Sigma-Aldrich) and 2% pectolase (Yakult) at 37 °C for 50–70 min (meiotic pachytene) and 2–3 h (meiotic metaphase and anaphase). Finally, these digested root tips and anthers were smeared onto slides as described previously . The slides with well-spread chromosomes will be prepared for FISH and GISH experiments.
The cytoplasm was removed by pepsin treatment to facilitate penetration of the probes. The slides with well-spread chromosomes were treated as described protocol  with some modifications: washing slides two times in 2× SSC in a Coplin jar for 3 min, followed by treating with 0.1 mg/mL pepsin (Sigma-Aldrich) in 10 mM HCl at 37 °C for 40 s - 1 min, and washing two times in 2× SSC at room temperature for 3 min. Finally, postfix slides with 4% formaldehyde in 2× SSC for 10 min, followed by washing two times in 2× SSC for 5 min, dehydrating in 70, 90 and 100% ethanol for 5 min each, and being left to air-dry.
Probe preparation and oligo-FISH
Previous studies indicated that C7 is highly conserved and preserves a complete synteny with C. hystrix chromosome H1, whereas C5 corresponds to two C. hystrix chromosomes H9 and H10 [11, 25]. The oligo-probes of C5 and C7 were developed using the oligo selection software Chorus (https://github.com/forrestzhang/Chorus) . Briefly, the repetitive sequences in the cucumber ‘Chinese Long’ genome (ftp://cucurbitgenomics.org/pub/cucurbit/genome/cucumber/Chinese_long/v2/, v2 Genome) were filtered using RepeatMasker (http://www.repeatmasker.org/). The filtered C5 and C7 sequences were divided into oligos (45 nt) with a step size of 5 nt. Each oligo was aligned to the cucumber genome to filter out those with duplicates in the genome (> 75% similarity over all 45 nt). Oligo with dTm > 10 (dTm = Tm – hairpin Tm) were kept to build an oligo probe database. We adjusted the number of oligos across the chromosomes to ensure that the oligo probes can cover the entire chromosomes. Specifically, oligos targeting 100 kb per 300 kb were chosen for chromosomes C5 and C7 (Additional file 3: Table S2). A total of 27,392 oligos per oligo pool were synthesized de novo in parallel by Mycroarray (Ann Arbor, MI), and were labeled following published protocols [14, 19].
The two oligo probes and two satellite DNA sequences, Type III and 45 s rDNA, were used for identifying chromosomes at mitosis and meiosis. Genomic DNA was extracted from cucumber and C. hystrix using the CTAB method , which were then labeled as GISH probes for distinguishing two subgenomes at meiosis of allotetraploid C. ×hytivus . All the experimental procedures for FISH were performed as previously described [14, 18]. The final images contrast was processed using ADOBE PHOTOSHOP CC (Adobe, http://www.adobe.com). The pachytene chromosomes in Fig. 2 were straightened using ImageJ software (https://imagej.nih.gov/ij/).
Five biological replicates were prepared for pollen viability. Fifteen male flowers were randomly collected from each generation of allotetraploid C. ×hytivus for each biological replicate. Pollen grains were collected and stained with modified Carbol-fuchsin solution, and more than 2000 pollen grains per biological replicate were counted under a stereomicroscope. The percentage of plump pollen grains was calculated to represent the pollen viability of allotetraploid C. ×hytivus.
Meiotic behaviors of five chromosomes and two subgenomes were counted in Additional file 2: Table S1. Statistical tests for each comparison and graphical analysis were executed in GraphPad Prism 7 (https://www.graphpad.com). An F-test was used to test for differences in the ranges of SD, and the pairwise Student’s t-test was used for comparisons pollen viability of four different generations and meiotic chromosome behavior of two different generations.
Availability of data and materials
All relevant supporting data can be found within the Additional files accompanying this article. Additional file 3: Table S2 contains all information about the number and locations of selected oligos of cucumber chromosomes 5 and 7. The Chorus software used for oligo-FISH probe design is freely available (https://github.com/forrestzhang/Chorus).
Anaphase I of meiosis
Fluorescence in situ hybridization
Genomic in situ hybridization
Metaphase I of meiosis
Pollen mother cell
Bird KA, VanBuren R, Puzey JR, Edger PP. The causes and consequences of subgenome dominance in hybrids and recent polyploids. New Phytol. 2018;220(1):87–93.
Adams KL, Wendel JF. Polyploidy and genome evolution in plants. Curr Opin Plant Biol. 2005;8(2):135–41.
Doyle JJ, Flagel LE, Paterson AH, Rapp RA, Soltis DE, Soltis PS, Wendel JF. Evolutionary genetics of genome merger and doubling in plants. Annu Rev Genet. 2008;42(1):443–61.
Cifuentes M, Grandont L, Moore G, Chèvre AM, Jenczewski E. Genetic regulation of meiosis in polyploid species: new insights into an old question. New Phytol. 2010;186(1):29–36.
Hollister JD. Polyploidy: adaptation to the genomic environment. New Phytol. 2015;205(3):1034–9.
Bomblies K, Jones G, Franklin C, Zickler D, Kleckner N. The challenge of evolving stable polyploidy: could an increase in “crossover interference distance” play a central role? Chromosoma. 2016;125(2):287–300.
Zhou J, Tan C, Cui C, Ge X, Li Z. Distinct subgenome stabilities in synthesized Brassica allohexaploids. Theor Appl Genet. 2016;129(7):1257–71.
Zhang H, Bian Y, Gou X, Zhu B, Xu C, Qi B, Li N, Rustgi S, Zhou H, Han F. Persistent whole-chromosome aneuploidy is generally associated with nascent allohexaploid wheat. Proc Natl Acad Sci U S A. 2013;110(9):3447–52.
Bian Y, Yang C, Ou X, Zhang Z, Wang B, Ma W, Gong L, Zhang H, Liu B. Meiotic chromosome stability of a newly formed allohexaploid wheat is facilitated by selection under abiotic stress as a spandrel. New Phytol. 2018;220(1):262–77.
Soltis DE, Visger CJ, Marchant DB, Soltis PS. Polyploidy: pitfalls and paths to a paradigm. Am J Bot. 2016;103(7):1146–66.
Wang Y, Zhao Q, Qin X, Yang S, Li Z, Li J, Lou Q, Chen J. Identification of all homoeologous chromosomes of newly synthetic allotetraploid Cucumis× hytivus and its wild parent reveals stable subgenome structure. Chromosoma. 2017;126(6):713–28.
Rens W, Fu B, O'Brien PC, Ferguson-Smith M. Cross-species chromosome painting. Nat Protoc. 2006;1(2):783.
Nie W, Wang J, Su W, Wang D, Tanomtong A, Perelman PL, Graphodatsky AS, Yang F. Chromosomal rearrangements and karyotype evolution in carnivores revealed by chromosome painting. Heredity. 2012;108(1):17–27.
Han Y, Zhang T, Thammapichai P, Weng Y, Jiang J. Chromosome-specific painting in Cucumis species using bulked oligonucleotides. Genetics. 2015;200(3):771–9.
Braz GT, He L, Zhao H, Zhang T, Semrau K, Rouillard JM, Torres GA, Jiang J. Comparative oligo-FISH mapping: an efficient and powerful methodology to reveal karyotypic and chromosomal evolution. Genetics. 2018;208(2):513–23.
Xin H, Zhang T, Han Y, Wu Y, Shi J, Xi M, Jiang J. Chromosome painting and comparative physical mapping of the sex chromosomes in Populus tomentosa and Populus deltoides. Chromosoma. 2018;127(3):313–21.
He L, Braz GT, Torres GA, Jiang J. Chromosome painting in meiosis reveals pairing of specific chromosomes in polyploid Solanum species. Chromosoma. 2018;127(4):505–13.
Lou Q, Zhang Y, He Y, Li J, Jia L, Cheng C, Guan W, Yang S, Chen J. Single-copy gene-based chromosome painting in cucumber and its application for chromosome rearrangement analysis in Cucumis. Plant J. 2014;78(1):169–79.
Albert PS, Zhang T, Semrau K, Rouillard J-M, Kao Y-H, Wang C-JR, Danilova TV, Jiang J, Birchler JA. Whole-chromosome paints in maize reveal rearrangements, nuclear domains, and chromosomal relationships. Proc Natl Acad Sci U S A. 2019;116(5):1679–85.
Chen JF, Kirkbride JH Jr. A new synthetic species of Cucumis (Cucurbitaceae) from interspecific hybridization and chromosome doubling. Brittonia. 2000;52(4):315–9.
J-f C, Staub JE, Tashiro Y, Isshiki S, Miyazaki S. Successful interspecific hybridization between Cucumis sativus L. and C. hystrix Chakr. Euphytica. 1997;96(3):413–9.
Chen L, Lou Q, Zhuang Y, Chen J, Zhang X, Wolukau JN. Cytological diploidization and rapid genome changes of the newly synthesized allotetraploids Cucumis × hytivus. Planta. 2007;225(3):603–14.
Jiang B, Lou Q, Wu Z, Zhang W, Wang D, Mbira KG, Weng Y, Chen J. Retrotransposon- and microsatellite sequence-associated genomic changes in early generations of a newly synthesized allotetraploid Cucumis × hytivus Chen & Kirkbride. Plant Mol Biol. 2011;77(3):225.
Sebastian P, Schaefer H, Telford IR, Renner SS. Cucumber (Cucumis sativus) and melon (C. melo) have numerous wild relatives in Asia and Australia, and the sister species of melon is from Australia. Proc Natl Acad Sci U S A. 2010;107(32):14269–73.
Yang L, Koo DH, Li D, Zhang T, Jiang J, Luan F, Renner SS, Henaff E, Sanseverino W, Garcia-Mas J. Next-generation sequencing, FISH mapping and synteny-based modeling reveal mechanisms of decreasing dysploidy in Cucumis. Plant J. 2014;77(1):16–30.
Naranjo T. The use of homoeologous pairing in the identification of homoeologous relationships in Triticeae. Hereditas. 1992;116(3):219–23.
Yang L, Koo DH, Li Y, Zhang X, Luan F, Havey MJ, Jiang J, Weng Y. Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly. Plant J. 2012;71(6):895–906.
Corredor E, Lukaszewski AJ, Pachón P, Allen DC, Naranjo T. Terminal regions of wheat chromosomes select their pairing partners in meiosis. Genetics. 2007;177(2):699–706.
Grusz AL, Sigel EM, Witherup C. Homoeologous chromosome pairing across the eukaryote phylogeny. Mol Phylogenet Evol. 2017;117:83–94.
Truco MJ, Hu J, Sadowski J, Quiros CF. Inter- and intra-genomic homology of the Brassica genomes: implications for their origin and evolution. Theor Appl Genet. 1996;93(8):1225–33.
Schubert I, Vu GT. Genome stability and evolution: attempting a holistic view. Trends Plant Sci. 2016;21(9):749–57.
Knight E, Greer E, Draeger T, Thole V, Reader S, Shaw P, Moore G. Inducing chromosome pairing through premature condensation: analysis of wheat interspecific hybrids. Funct Integr Genomics. 2010;10(4):603–8.
Lee Y-I, Chang F-C, Chung M-C. Chromosome pairing affinities in interspecific hybrids reflect phylogenetic distances among lady's slipper orchids (Paphiopedilum). Ann Bot. 2011;108(1):113–21.
Cui C, Ge X, Gautam M, Kang L, Li Z. Cytoplasmic and genomic effects on meiotic pairing in Brassica hybrids and allotetraploids from pair crosses of three cultivated diploids. Genetics. 2012;191(3):725–38.
Gou X, Bian Y, Zhang A, Zhang H, Wang B, Lv R, Li J, Zhu B, Gong L, Liu B. Transgenerationally precipitated meiotic chromosome instability fuels rapid karyotypic evolution and phenotypic diversity in an artificially constructed allotetraploid wheat (AADD). Mol Biol Evol. 2018;35(5):1078–91.
Simon WLC. In a battle between parental chromosomes, a failure to reload. Proc Natl Acad Sci U S A. 2011;108(33):13361–2.
Storme ND, Mason A. Plant speciation through chromosome instability and ploidy change: cellular mechanisms, molecular factors and evolutionary relevance. Curr Plant Biol. 2014;1(1):10–33.
Han Y, Pan J, Thammapichai P, Li Z, Weng Y. Asynchronous meiosis in Cucumis hystrix–cucumber synthetic tetraploids resulting in low male fertility. Crop J. 2016;4(4):275–9.
Xiong Z, Gaeta RT, Pires JC. Homoeologous shuffling and chromosome compensation maintain genome balance in resynthesized allopolyploid Brassica napus. Proc Natl Acad Sci U S A. 2011;108(19):7908–13.
Chester M, Gallagher JP, Symonds VV, Cruz da Silva AV, Mavrodiev EV, Leitch AR, Soltis PS, Soltis DE. Extensive chromosomal variation in a recently formed natural allopolyploid species, Tragopogon miscellus (Asteraceae). Proc Natl Acad Sci U S A. 2012;109(4):1176–81.
Chester M, Riley RK, Soltis PS, Soltis DE. Patterns of chromosomal variation in natural populations of the neoallotetraploid Tragopogon mirus (Asteraceae). Heredity. 2015;114(3):309–17.
Yu X, Wang X, Hyldgaard B, Zhu Z, Zhou R, Kjaer KH, Ouzounis T, Lou Q, Li J, Cai Q, et al. Allopolyploidization in Cucumis contributes to delayed leaf maturation with repression of redundant homoeologous genes. Plant J. 2018;94(2):393–404.
Zhuang Y, Chen J-F, Jahn M. Expression and sequence variation of the cucumber por gene in the synthesized allotetraploid Cucumis × hytivus. Mol Biol Rep. 2008;36(7):1725.
Zhuang Y, Chen J-F. Changes of gene expression in early generations of the synthetic allotetraploid Cucumis × hytivus Chen et Kirkbride. Genet Resour Crop Evol. 2009;56(8):1071.
Lysak MA, Mandáková T. Analysis of plant meiotic chromosomes by chromosome painting. Methods Mol Biol. 2013;990(2):13–24.
Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8(19):4321–6.
We thank Dr. Mengli Xi (Nanjing Forestry University, China) for technical assistance. We are very grateful to the editor and reviewers for critically evaluating the manuscript and providing constructive comments for its improvement.
This research was supported by National Natural Science Foundation of China (Key Program, No.31430075) and National Natural Science Foundation of China (No. 31471872). The funding agencies provided funding to the research projects, but played no role in the design of study, collection and analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Comparison of pollen viability of four generations of C. ×hytivus allotetraploid (S4, S8, S11 and S14). Figure S2. Locations and density of 27,392 oligos along the sequence map of cucumber chromosomes 5 (a) and 7 (b). Figure S3. Two translocations carrying C5-oligo signals were detected in other PMCs from the same individual as Fig. 5. Figure S4. Tracing chromosome pairing at meiotic pachytene of S4 generation. Figure S5. Three representative PMCs with multivalents and one lagged C-bivalent at metaphase I in S4 generation. Figure S6. A heatmap depicting the frequency of different lagged bivalent numbers in the asynchronous meiotic PMCs of the S4 and S14 generation.
Meiotic analysis of F1, S4 and S14 plants based on oligo-painting and GISH.
Location and sequence information of C5- and C7-oligos.
About this article
Cite this article
Zhao, Q., Wang, Y., Bi, Y. et al. Oligo-painting and GISH reveal meiotic chromosome biases and increased meiotic stability in synthetic allotetraploid Cucumis ×hytivus with dysploid parental karyotypes. BMC Plant Biol 19, 471 (2019). https://doi.org/10.1186/s12870-019-2060-z
- Interspecies hybridization
- Meiotic instability
- Karyotype variation