The DcPS1 RNAi line produces diploid male gametes in carnation
To explore DcPS1 function in carnation, an RNAi construct was generated utilizing a 500 bp region, and this fragment was inserted into an RNAi vector. 35S:RNAi-DcPS1 was introduced into carnation calli using Agrobacterium-mediated T-DNA transfer. One independent DcPS1 RNAi line was recovered, and transgene integration was verified by PCR (Supplementary Figs. 2A, 3). The results showed that the expression of DcPS1 was significantly decreased (t test, P < 0.05) in transgenic plants compared with the wild type (Supplementary Fig. 2B).
At the tetrad stage of pollen mother cells (PMCs) in the WT and DcPS1 RNAi lines, a total of 3302 and 2489 PMCs were observed, respectively. Among them, the vast majority (89.95%) were tetrads in WT (Fig. 1A). In addition, we found 0.18% triads (Fig. 1B), 9.18% dyads (Fig. 1C), 0.67% monads (Fig. 1D) and 0.03% polyads (Fig. 1E), as well as a few abnormal dyads (Fig. 1F) and tetrads (Fig. 1G, H) in WT. The dyads and triads eventually produced two 2n pollen grains and one 2n pollen grain (Fig. 1I), respectively. In the DcPS1 RNAi line, however, we found 71.03% tetrads, 0.24% triads, 28.6% dyads, and 0.12% monads. The production of dyads in the DcPS1 RNAi line was observed at a higher frequency than that in wild type (Fig. 1J, K). These data strongly supported the conclusion that DcPS1 is involved in diplogamete production in carnation.
Both normal and abnormal meiosis were observed during microsporogenesis in the WT and DcPS1 RNAi lines. In normal meiosis, the metaphase I spindle aligns the bivalents at the equatorial plane (Fig. 2A) and subsequently segregates the homologues towards the opposite poles at anaphase I (Fig. 2B). Chromosomes decondensed at telophase I at the end of meiosis I (Fig. 2C). Perpendicular spindles are usually observed at metaphase II (Fig. 2D), leading to the formation of four well-separated poles at anaphase II and telophase II (Fig. 2E, F) and tetrads at the end of meiosis. In abnormal meiosis, abnormalities were observed in meiosis I and II, including an absence of the second meiotic division (Fig. 2G-I), tripolar spindles (Fig. 2J), fused spindles (Fig. 2L) and lagged chromosomes at metaphase II (Fig. 2N, O) and telophase II (Fig. 2P). These defects might be the main reason for dyad and triad formation in the DcPS1 RNAi line, which occurred at a higher frequency than in the wild type.
The OSDLa RNAi line produces diploid male gametes in carnation
PCR amplification was used to screen positive transgenic lines (Supplementary Figs. 4A, 5). Four transgenic plants were obtained in which OSDLa expression was suppressed to varying degrees in different RNAi lines. The results showed that the expression of OSDLa was significantly decreased (t test, P < 0.05) in transgenic plants compared with the wild type (Supplementary Fig. 4B).
The meiotic division of some microspores in transgenic carnation plants proceeded normally (Fig. 3). At the leptotene stage, the chromatin condensed into a filamentous structure due to spiral curling (Fig. 3A, B). At the zygotene stage, homologous chromosomes gradually moved closer and started pairing, and the synaptonemal complex (SC) also began to form at this stage (Fig. 3C, D). The chromosomes were further condensed and fully synapsed along the SC at pachytene, and the nonsister chromatids of homologous chromosomes recombined locally (Fig. 3E). Homologous chromosomes were further condensed, and the SC began to disintegrate and fall off chromosomes, but homologous chromosomes were still cross-linked (Fig. 3F, G). Chromosomes were neatly arranged on the metaphase I equatorial plate (Fig. 3H), and homologous chromosomes were separated and moved to the two poles at anaphase I (Fig. 3I). The two sets of fifteen homologues aligned neatly on the two metaphase II plates (Fig. 3J). The second round of segregation at anaphase II (Fig. 3K) led to the formation of four sets of chromosomes that decondensed to form the spore nuclei and finally formed a tetrad (Fig. 3L). In OSDLa RNAi line meiocytes, abnormal chromosome bridges appeared at meiosis I (Fig. 4A). Homologous chromosomes separated at telophase I (Fig. 4B) and formed a dyad (Fig. 4C, D), suggesting that dyad production is due to an absence of the second meiotic division. Lagging chromosomes at anaphase II (Fig. 4E) maybe lead to unbalanced polyads (Fig. 4F). Oriented metaphase II/anaphase II spindles were aberrant and formed tripolar spindles (Fig. 4G, H). This defect in spindle orientation explained the appearance of triads (Fig. 4I, J).
At the tetrad stage of pollen mother cells (PMCs) in WT, OSDLa RNAi plants C1, C2, G3 and B1, a total of 7360, 1326, 4254, 363 and 509 PMCs were observed, respectively. In wild-type diploid plants, 97.53% tetrads were observed (Fig. 4K). In addition, we found 0.56% triads, 1.79% dyads, 0.11% monads and 0.01% polyads (Fig. 4K). In the OSDLa RNAi plant C1, 93.21% tetrads, 1.28% triads, 5.51% dyads were observed. In the OSDLa RNAi plant C2, we found 95.79% tetrads, 0.68% triads, 3.31% dyads, 0.19% monads and 0.02% polyads. In the OSDLa RNAi plant G3, 75.48% tetrads, 1.93% triads, 22.59% dyads were found. In the OSDLa RNAi plant B1, 89.98% tetrads, 1.96% triads, 7.86% dyads and 0.2% polyads were observed (Fig. 4K). Compared to the WT, in the OSDLa RNAi lines, a high frequency of dyads and triads was observed, confirming that the diploid microspores in OSDLa RNAi lines are produced by the defects in the meiotic process.
Carnation OSDLa interacts with DcPS1 and DcRAD51D
We chose full-length DcPS1, DcRAD51D and OSDLa to further analyse protein-protein interactions. The vectors pGBKT7-DcPS1, pGBKT7-DcRAD51D and pGBKT7-OSDLa were transformed into the Y2HGold yeast strain, and the resulting transformants were plated on SD medium lacking tryptophan but containing a chromogenic substrate for yeast galactosidase and aureobasidin A (SD/−Trp/X-α-gal/AbA). The pGBKT7-DcPS1 colonies turned blue on SD/ -Trp/X-α-gal/AbA plates (Supplementary Fig. 6A), while the pGBKT7-DcRAD51D and pGBKT7-OSDLa colonies did not turn blue (Supplementary Fig. 6B, C). This result indicated that the bait pGBKT7-DcPS1 could autonomously activate the reporter genes,whereas pGBKT7-DcRAD51D and pGBKT7-OSDLa could not autonomously activate the reporter genes in the absence of prey protein and were, therefore, suitable for screening in the Y2H assay. Y2HGold cells were cotransformed with pGBKT7-DcRAD51D and pGADT7-OSDLa and plated on QDO/X/A plates, and a positive interaction was indicated by the presence of blue colonies (Fig. 5A). The results showed that carnation DcRAD51D interacted with OSDLa.
The DcRAD51D-OSDLa interaction was also confirmed via bimolecular fluorescence complementation (BiFC) assays in Nicotiana protoplasts in vivo, in which DcRAD51D was fused to the amino-terminal half of yellow fluorescent protein (YFP) (DcRAD51D-YFPN) and OSDLa was fused to the carboxy-terminal half of YFP (OSDLa-YFPC). For two combinations (DcRAD51D-YFPN + OSDLa-YFPC and DcPS1-YFPN + OSDLa-YFPC), YFP signals were observed in the nucleus and cytoplasm (Fig. 5B, C). There was no YFP signal when DcPS1-YFPN was cotransformed with DcRAD51D -YFPC (Fig. 5D). These results showed that carnation OSDLa interacted with DcPS1 in vivo, and OSDLa interacted with DcRAD51D both in vivo and in vitro, whereas DcPS1 was unable to associate with DcRAD51D.
DcPS1 positively regulates the expression levels of OSDLa and DcRAD51D
Similar to the DcPS1 RNAi lines, RNAi lines for OSDLa also produced high numbers of unreduced microspores through the absence of the second meiotic division and a tripolar spindle mechanism at meiosis II (MII) (Fig. 4). Moreover, the frequencies of dyad/triad formation in DcPS1 and OSDLa RNAi lines were higher than those in the wild type (Figs. 1 and 4). BiFC analysis also confirmed that OSDLa interacted with DcPS1 (Fig. 5), suggesting that both genes are involved in the same molecular pathway.
To analyse the potential regulatory effects of the DcPS1 protein on OSDLa and DcRAD51D gene expression and the OSDLa protein on DcPS1 and DcRAD51D gene expression, RT–qPCR analyses using gene-specific primers were performed on RNA harvested from both DcPS1 and OSDLa RNAi lines. In the DcPS1 RNAi lines, a decrease in OSDLa and DcRAD51D expression was observed (Fig. 6A), indicating that DcPS1 either positively regulates the expression of OSDLa and DcRAD51D or prevents the degradation of the corresponding transcript. Compared with the wild type, DcPS1 RNAi lines showed a significant decrease in DcPS1 gene expression (Fig. 6A), which led to a larger increase in the frequency of dyad and triad formation (Fig. 1K), similar to that in OSDLa-RNAi lines (Fig. 6B-E, Fig. 4K). In the OSDLa RNAi lines C1 and C2,there was no significant decrease in DcPS1 gene expression (Fig. 6C, D) with a smaller increase in the frequency of dyad and triad formation than in the WT (Fig. 4K). However, in the OSDLa RNAi lines B1 and G3, there was a significant decrease in DcPS1 gene expression (Fig. 6B, E) and a greater increase in the frequency of dyad and triad formation than in the WT (Fig. 4K). Because the penetrance of the 2n pollen phenotype (dyad and triad formation) for all lines was reflected by the decrease in DcPS1 expression, there appeared to be a close correlation between the DcPS1 expression decrease and the frequency of 2n formation. In OSDLa RNAi lines, DcRAD51D gene expression was decreased compared to the wild type (Fig. 6), and the Y2H assay and BiFC analysis also confirmed that OSDLa interacted with DcRAD51D, which suggested a meiosis-specific function of OSDLa in regulating DcRAD51D expression.