ddp1 mutants exhibit partial pollen abortion
The rice genome contains 22 KCS genes (Yang et al. 2023). To investigate the functions of these KCSs, we employed CRISPR/Cas9 technology to generate a series of knockout (KO) mutants targeting these genes (Table S1). Among all developed mutants, the knockout mutants of LOC_Os02g49920 exhibited a partial male sterility attributed to defective anther dehiscence and pollen fertility. Thus, the mutant of LOC_Os02g49920 was designated as ddp1 (defective in dehiscence and pollen1). two alleles ddp1 namely ddp1-1 and ddp1-2 were identified in T1 progeny. ddp1-1 and ddp1-2 carried 4-bp and 2-bp deletion in the sole exon of LOC_Os02g49920, respectively (Fig. 1b). These deletions caused reading-frame shifts, generating premature and delayed stop codon, respectively.
Compared to WT, ddp1-1/2 mutants exhibited no significant differences in whole-plant morphology and the structure of flower organs such as spikelets and anthers (Fig. 1c-f). However, at the mature stage, the seed-setting rates of ddp1-1 (41.533 ± 1.050%) and ddp1-2 (32.100 ± 1.400%) were significantly lower than that of WT (95.000 ± 0.816%) (Table 1)(Fig. 1g, k). To determine if the lower seed-setting rate in ddp1-1/2 mutants resulted from a male or female gametophyte impairment, we evaluated the pollen and female fertility. Notably, the fertile pollen rate was significantly lower in ddp1-1 (54.67% ± 3.68%) and ddp1-2 (49.17% ± 1.65%) compared to WT, which exhibited a fertile pollen rate of 95.13% ± 0.84% (Fig. 1h-j, l). When ddp1 mutants were self-pollinated with its pollen, their seed setting rate was only 24.45% ± 2.53%. However, after pollination with WT pollen, the seed setting rate of ddp1 mutants reached 58.60% ± 5.88% (Fig. S1). The findings suggest that the reduced seed-setting rate observed in ddp1 mutants is due to abnormalities in pollen development, rather than any deficiency in the female gametophyte.
Table 1
The agronomic traits of WT and ddp1 mutant.
Agronomic traits | WT Mean ± SD | ddp1-1 Mean ± SD | ddp1-2 Mean ± SD |
Plant height (cm) | 86.333 ± 2.055 | 83.667 ± 1.247 | 88.000 ± 0.816 |
Panicle length (cm) | 21.500 ± 0.408 | 23.533 ± 0.500 | 21.400 ± 0.638 |
1000-grain weight (g) | 26.647 ± 0.509 | 27.406 ± 0.561 | 26.810 ± 0.752 |
Grain length (cm) | 0.811 ± 0.009 | 0.797 ± 0.005 | 0.811 ± 0.009 |
Grain width (cm) | 0.351 ± 0.002 | 0.354 ± 0.003 | 0.352 ± 0.002 |
Grain thickness (mm) | 2.32 ± 0.037 | 2.303 ± 0.054 | 2.317 ± 0.049 |
Seed-setting Rate (%) | 95.000 ± 0.816 | 41.533 ± 1.050** | 32.100 ± 1.400* |
Number of tillers | 10.333 ± 1.247 | 9.667 ± 1.247 | 9.667 ± 0923 |
Number of primary branches | 10.000 ± 0.816 | 10.333 ± 1.247 | 10.000 ± 1.633 |
Number of secondary branches | 27.667 ± 1.700 | 23.000 ± 2.60* | 17.333 ± 1.700* |
Yield per plant (g) | 16.000 ± 1.633 | 6.233 ± 1.481** | 4.433 ± 0.419** |
*, ** significant at 5% and 1% levels of probability, respectively. |
In addition to partial pollen abortion, ddp1-1/2 mutants also exhibited significant differences in certain yield-related traits compared to WT. For instance, the number of secondary branches and yield per plant were significantly lower in ddp1-1 mutant (23.00 ± 2.16% and 6.23 ± 1.48%) compared to WT (27.67 ± 1.70% and 16.00 ± 1.62%) (Table 1). The ddp1-2 mutant showed similar trends to those observed in ddp1-1 (Table 1). Nevertheless, the comparison revealed no notable differences in plant height, number of tillers and primary branches, panicle length, 1000-grain weight, and grain dimensions (length, width, thickness) between ddp1-1/2 and WT (Table 1).
To confirm whether LOC_Os02g49920 was the target gene responsible for the ddp1 mutants phenotype, a complementation test was performed by transforming a construct DPP1pro:DPP1CDS-mCitrine containing a 2.5-kb DDP1 promoter region, 1.458-kb DDP1 coding region and 0.72-kb mCitrine coding region into Cas9-free ddp1-1 mutants (Table S1). As expected, the complementary transgenic lines exhibited a restoration of seed-setting rates and pollen fertility to WT levels (Fig. S2), thereby confirming that LOC_Os02g49920 is indeed the target gene whose knockout causes partial male sterility in ddp1 mutants.
Anther dehiscence is defective in ddp1 mutants
To determine whether partial male sterility of ddp1-1/2 mutants was associated with the defects in anther dehiscence, pollen germination, tube growth, and dehydration, we compared these critical reproductive processes in ddp1-1/2 mutants to those in WT. For anther dehiscence analysis, we examined the anther dehiscence phenotype in both ddp1-1/2 and WT. WT anthers typically undergo dehiscence, facilitating the dispersal of mature pollen (Fig. 2a, f), while ddp1-1/2 mutants exhibited various defects in anther dehiscence and pollen grain release (Fig. 2b-e, g-j). The anther dehiscence phenotype was quantified using the dehiscence index as reported previously (Zhang et al. 2004). A higher dehiscence index indicates that anthers have dehisced more completely. Our findings indicated that the anther dehiscence index was significantly lower in ddp1-1 (2.89) and ddp1-2 (2.73) compared to WT (4.09) (Table 2). This impaired anther dehiscence in ddp1-1/2 mutants may result in reduced adhesion and germination of fertile pollen grains on the stigmas. To validate this finding, we performed an in vivo pollen germination experiment to evaluate the pollen germination rates on the stigmas of ddp1-1/2 and WT. The assay demonstrated a markedly lower pollen germination rate on the stigmas in ddp1-1/2 mutants relative to WT (Fig. 2n-r, u). Collectively, these results indicate that defective anther dehiscence in ddp1-1/2 mutants hinders the adhesion and subsequent germination of fertile pollen grains on the stigma, thereby reducing the seed-setting rate.
Table 2
The anther dehiscence index of WT and ddp1.
| Level of anther dehiscence | Total number of anther | Dehiscence index |
| 5 | 4 | 3 | 2 | 1 |
WT | 8 | 84 | 0 | 0 | 0 | 92 | 4.09 |
ddp1-1 | 10 | 34 | 26 | 18 | 24 | 112 | 2.89** |
ddp1-2 | 14 | 41 | 36 | 30 | 41 | 162 | 2.73** |
** significant at 1% levels of probability. |
For pollen germination analysis, we compared pollen germination rate between ddp1-1/2 and WT using an in vitro pollen germination assay. Because ddp1-1/2 mutants contained sterile pollen grains which could not normally germinate (Fig. 2k-m), the overall in vitro pollen germination rates were significantly lower in ddp1-1 (50.65 ± 1.86%) and ddp1-2 (42.23 ± 2.36%) compared to WT (89.4 ± 1.11%) (Fig. 2s). However, when only considering the germination rate of fertile pollen grains, the comparison revealed no notable differences between ddp1-1/2 and WT (with 82.03 ± 0.82% of ddp1-1 and 81.37 ± 0.45% of ddp1-2 versus 84.90 ± 1.35% of WT) (Fig. 2t). These findings suggest that the germination of fertile pollen in ddp1-1/2 mutants remains unaffected.
For pollen tube growth analysis, we compared pollen tube growth within the pistil of ddp1-1/2 and WT using an in vivo pollen germination assay. In self-pollinated WT pistils, at 5 hours post-pollination, pollen tubes had successfully reached the ovule micropyle (Fig. 2n). However, self-pollinated pistils of the ddp1-1/2 mutants showed a reduced germination rate of pollen grains, yet those that did germinate were able to extend their pollen tubes to the ovule micropyle (Fig. 2o, q), suggesting normal pollen tube growth within ddp1-1/2 pistils.
For pollen dehydration analysis, we compared the in vitro dehydration rate of fertile pollen between ddp1-1 and WT. However, no significant differences in the dehydration rate of fertile pollen were observed between ddp1-1 and WT from 0 second to 30 minutes (Fig. S3). Altogether, these findings suggest that while pollen germination, pollen tube growth and pollen dehydration are not compromised in ddp1-1/2 mutants, anther dehiscence is abnormal. Consequently, we conclude that the decrease in seed-setting rate of ddp1-1/2 mutants was attributed to the aborted pollen and aberrant anther dehiscence.
Cytological comparison of anther development between ddp1 and WT
To determine the morphological defect of ddp1-1/2 anthers, semi-thin sections of both ddp1-1/2 and WT anthers at different development stages were examined using light microscopy. Initial observations revealed that no significant differences were detected between ddp1-1/2 and WT anthers from stage 6 to 9 (Fig. 3a-c, h-j, p-r). However, at stage 10, the epidermis and endothecium of ddp1-1/2 anthers appeared to exhibit a greater thickness compared to their counterparts in WT anthers, and microspores of ddp1-1/2 exhibited collapse and a defect in vacuolation (Fig. 3k, s). At stage 11, the epidermis and endothecium of ddp1-1/2 anthers continued to show swelling, and tapetum degradation seemed to be delayed. Only a fraction of microspores formed falcate-shaped pollen grains filled with starch granules, while the rest developed into irregularly shaped pollen grains lacking starch accumulation (Fig. 3l, t). At stage 12, the WT tapetum had completely degraded, allowing microspores to mature into pollen grains (Fig. 3f). In contrast, in ddp1-1/2 anthers, the incompletely degraded tapetum residues were observed, with only some microspores maturing into fertile pollen grains, while others aborted (Fig. 3m, u). At stage 14, WT anthers were able to dehisce normally as the connective tissue between anthers cracked (Fig. 3g). However, in ddp1-1/2, only a subset of anthers dehisced, while others remained closed (Fig. 3n-o, v-w). These findings show that ddp1-1/2 mutations lead to developmental defects in anther and pollen from stage 10 to 14.
To precisely characterize the differences between ddp1-1/2 and WT, SEM was employed to compare the anthers and pollen grains. The anthers of ddp1-1/2 and WT were comparable in size (Fig. 4a, e, i). WT anthers exhibited an outer epidermis covered by spaghetti-like cuticle layers (Fig. 4b) and an inner epidermis with Ubisch bodies (Fig. 4c, d). In contrast, the cuticle layers of the outer epidermis and the Ubisch bodies of the inner epidermis in ddp1-1/2 anthers seemed to be denser (Fig. 4f-h, g-l). In addition, WT pollen grains were spherical and smooth, uniformly coated by a layer of sporopollenin (Fig. 5a, d). ddp1-1/2 mutants had two types of pollen grains, namely spherical pollen grains and shrunken pollen grains (Fig. 5b-c, e-h). The spherical pollen grains in ddp1-1/2 mutants exhibited a slightly sparse sporopollenin layer (Fig. 5j, l), while the sporopollenin on the shrunken pollen grains in ddp1-1/2 mutants appeared to be denser (Fig. 5k, m). The increased accumulation of sporopollenin on the surfaces of shrunken ddp1-1/2 pollen may be attributed to the shrinkage of the pollen grain surface area. Given the similar phenotypes exhibited by ddp1-1 and ddp1-2 under identical growth conditions, the following analyses within this study were conducted solely on ddp1-1.
To further elucidate the cytological differences between ddp1 and WT, TEM was also employed to observe the ultrastructures of anthers and pollen grains. At stage 11, in accordance with our light microscopy results, the anther wall, particularly the tapetum of ddp1-1 mutant was thicker compared to that of WT (Fig. 6a, e). At stage 13, the WT tapetum fully degraded while the incompletely degraded tapetum residue was still observed in ddp1-1 mutant (Fig. 6i, m). At stage 11 and 13, the cuticle layer and the Ubisch body in ddp1-1 anthers were denser than those of WT (Fig. 6f-g, n-o), aligning with SEM observations. Furthermore, the aborted ddp1-1 pollen grains exhibited a thicker exine and smaller bacula compared to WT (Fig. 6h, p), indicating the alterations in the pollen wall ultrastructure. At stage 10, neither ddp1-1 nor WT microspores showed signs of starch accumulation and intine formation (Fig. 6q, r). By stage 12, WT microspores had developed an intine and were filled with a substantial amount of starch granules (Fig. 6s), whereas aborted pollen grains of ddp1-1 mutant lacked both features (Fig. 6t). In summary, ddp1-1 mutant exhibited abnormalities in their anther cuticle layer, Ubisch bodies, tapetum degradation, anther dehiscence, pollen wall structure, and starch accumulation.
Gametogenesis at stage 11 and 12 is defective in ddp1-1 anthers
To determine whether gametogenesis was affected in ddp1 anthers, we performed the 4’,6-diamidino-2-phenylindole (DAPI) staining on microspores from both ddp1-1 and WT at different stages. At stage 9, WT microspores typically developed into uninucleate microspores. At stage 11, they underwent the first mitotic division, resulting in binucleate microspores. At stage 12, they underwent the second mitotic division, yielding trinucleate microspores with one vegetative nucleus and two generative nuclei. Accordingly, a portion of ddp1-1 microspores successfully developed into uninucleate, binucleate and trinucleate microspores at stage 9, 11 and 12, respectively, similar to WT (Fig. 7a). However, other ddp1-1 microspores exhibited a brightly stained nucleus at stage 9 and 11 but failed to complete the first and second mitotic divisions to form binucleate and trinucleate microspores at stage 11 and 12 (Fig. 7a). The DAPI staining analysis indicates that defective gametogenesis in the aborted ddp1-1 microspores begins at stage 11, which is roughly consistent with the semi-thin section observations that anther and pollen abnormalities occur at stage 10.
Defective gametogenesis was often closely associated with abnormalities in callose biosynthesis and degradation. Callose biosynthesis initiated at stage 7, and its degradation at stage 9 facilitated the release of microspores from the tetrads. In order to investigate whether these processes were affected in ddp1-1 anthers, we performed aniline blue staining on anthers from both ddp1-1 and WT from stage 7 to 9. However, no obvious differences in the fluorescence signal surrounding microspores were observed between ddp1-1 and WT anthers (Fig. 7b). This result suggests that the callose formation and degradation steps in ddp1-1 anthers from stage 7 to 9 are not impaired, consistent with the semi-thin section observations that anther and microspore development in ddp1-1 mutant appeared normal from stage 6 to 9.
Altered wax content and composition of the anther epidermis and pollen wall in the ddp1 mutant
DDP1 encodes a β-ketoacyl-CoA synthase with 485 amino acids residues, catalyzing cuticular wax synthesis. To elucidate the function of DDP1, a phylogenetic analysis was performed using DDP1 and other plant KCS proteins with both known and unknown biochemical and biological functions. Phylogenic analysis revealed that DDP1 shares high sequence similarity with its homologs (Fig. S4a). Domain homology analysis using SMART demonstrated that the DDP1 protein contains three conserved domains: an ACP_syn_III_C domain, a transmembrane region, and a FAE1_CUT1_RppA domain, similar to its counterparts in other plant species (Fig. S4b). The analysis suggests that DDP1 may play a conserved function in cuticular wax synthesis. Given that ddp1-1 mutant also showed structural abnormalities in the lipidic anther epidermis and pollen wall, we speculated that the wax content and composition of these structures in the ddp1-1 mutant might be altered. To test this speculation, we initially examined the lipidic compounds by staining anthers and pollen from ddp1-1 and WT at stage 12 with the lipophilic dye Sudan Red 7B. While the WT anther epidermis and pollen wall showed intense staining (Fig. 8a1-a3,), the ddp1-1 mutant displayed strong staining in anther epidermis and fertile pollen wall, and weak staining in the wall of aborted pollen grains (Fig. 8b1-b3). This suggests a reduction in lipidic compounds within the wall of ddp1-1 aborted pollen. Subsequently, we analyzed the phenolic compounds by comparing the fluorescence intensity emitted by phenolic compounds in pollen walls under ultraviolet radiation between ddp1-1 and WT. The fertile pollen grains of ddp1-1 emitted fluorescence similar in intensity to that of WT, whereas the aborted pollen grains exhibited significantly weaker fluorescence (Fig. 8b4), indicating a substantial decrease in phenolic compounds in the wall of ddp1-1 aborted pollen. In conclusion, these findings indicate that precursors to sporopollenin, including lipidic and phenolic compounds, are significantly reduced in the walls of ddp1-1 aborted pollen.
To ascertain the chemical profiles of anther epidermis and pollen wall, we quantified wax constituents at stage 12 using gas GC-MS. This result indicated that the total wax content in ddp1-1 anthers increased by 34%±0.08% compared to WT (Fig. 8c), primarily due to significant elevations in free fatty acids (C16:0, C18:2, C18:3, C18:0, C20:0, C22:0, C26:0,C28:0 and C30:0), alkenes (C25:0, C27:0, C33:1, C33:0, C35:0 and C35:1), and alcohols (C26:0 and C28:0), as well as sterols such as campersterol, stigmasterol and β-sitosterol (Fig. 8d, e). In contrast, the levels of certain free fatty acid (C24:0), alkenes (C23:0 and C27:1) and inositol were dramatically decreased in ddp1-1 mutants, and other wax constituents remained largely unchanged (Fig. 8e). The altered wax profile in ddp1-1 mutant indicated that DDP1 is essential for maintaining lipid metabolic homeostasis. The alteration in the content and composition of wax in the anther epidermis and pollen wall may underlies the structural abnormalities observed in the ddp1-1 anther epidermis and pollen wall.
Expression pattern and subcellular location of DDP1
To determine DDP1 expression pattern, RT-qPCR was conducted to quantify its expression in various tissues of WT. The result suggested that DDP1 transcription was uniformly detected across all assessed tissues, with peak transcription observed in the leaves (Fig. 9a). DDP1 expression was observed throughout the stages of anther development, peaking at stage 8 (Fig. 9a). To further analyze the expression of DDP1 in anthers, RNA in situ hybridization was conducted on sections of WT anthers. The hybridization signal was first observed in tapetum and pollen mother cells (PMCs) at stage 7, then in tapetum and tetrads at stage 8, and in tapetum and microspores from stage 9 to 11 (Fig. 9b). In control experiments, the sense probe failed to produce any detectable signal, confirming the specificity of the hybridization. These findings indicate that DDP1 is predominantly expressed in the tapetum, PMCs and microspores, which correlates with the phenotypes of delayed tapetum degradation and impaired microspore development in ddp1 mutants.
To further elucidate the spatial distribution of DDP1 protein, we utilized confocal microscopy to observe the mCitrine signal in ddp1-complementary transgenic plants expressing DDP1pro:DDP1cDNA-mCitrine. At stage 10, yellow mCitrine fluorescence was detected in the tapetum and microspores (Fig. 10a-c). By stage 11, yellow mCitrine fluorescence was noted in the partially degraded tapetum residues and microspores (Fig. 10d-f). At stage 12, with the tapetum fully degraded, mCitrine signal was exclusively detected in pollen grains (Fig. 10g-i). These observations indicate that DDP1 expression aligns with DDP1 protein distribution in the tapetum, microspores, and pollen grains.
Protein domain analysis using SMART revealed that the DDP1 protein featured a transmembrane domain at its N-terminus (Fig. S4b). Predictions of secondary structure and three-dimensional models using PSIPRED and SWISS-MODEL indicated that DDP1 possesses coil, helix, and strand domains (Fig. S5a). Mutations in the first helix domain of the ddp1-1/2 mutants led to noticeable alterations in the protein's three-dimensional structure (Fig. S5b, c).
To confirm the localization of DDP1 and ddp1-1 in rice protoplasts, we produced DDP1-GFP and ddp1-1-GFP fusion constructs and introduced them into rice protoplasts alongside the ER marker mCherry-HDEL. DDP1-GFP fluorescence was found to almost completely overlap with the mCherry-HDEL fluorescence (Fig. 9c), suggesting that DDP1 is specifically localized to the ER, as was expected. In contrast, the ddp1-1-GFP fluorescence signal only partially co-localized with the ER marker, with a significant portion of the green fluorescence observed in the cytoplasm (Fig. 9c).
The expression levels of known genes complicated in lipid metabolism, anther and pollen development, and anther dehiscence are significantly changed in ddp1 anthers
Observations of delayed tapetum degradation, defective anther and pollen development, impaired lipid metabolism, and abnormal anther dehiscence in the ddp1-1 mutant have led us to hypothesize that DDP1 mutations may disrupt the expression of genes associated with these processes. To test this speculation, we performed a comparative analysis to evaluate the expression profiles of known genes associated with these processes between ddp1-1 and WT from stage 7 to 10. Compared to WT, the expression levels of six lipid metabolism genes (OsABCG26, OsPKS2, OsNP1, RMS2, CYP703A3 and OsCER2), five genes complicated in the anther, pollen and tapetum development (OsCP1, GAMYB, OsSTRL2, STS1 and DPW3), and two genes involved in anther dehiscence (DAO and OsH1) were significantly decreased in the ddp1-1 mutant (Fig. 11). In contrast, the ddp1-1 mutant displayed a pronounced upregulation in the expression levels of two lipid metabolism genes (OsABCG3 and OsCER1), along with one tapetum development gene (OsAP25), specifically at stages 8 and 9 (Fig. 11). In addition, considering the close association of anther dehiscence with auxin and jasmonic acid (JA) signaling, we also conducted a comparative analysis of gene expression levels related to these signaling pathways between ddp1-1 and WT. Six genes (OsETTIN1, OsETTIN2, OsETTIN3, OsYUCCA1, OsYUCCA4 and OsMP) required for auxin synthesis and response and three JA signaling suppressor genes (JAZ1, JAZ6 and JAZ8) were prominently up-regulated expressed in the ddp1-1 mutant, while a significant downregulation of two JA synthesis genes (LOX2 and AOS3) was detected in the ddp1-1 mutant at stage 9 and 10, respectively (Fig. 11). Together, these findings suggest that mutations in DDP1 alter the expression of known genes complicated in lipid metabolism, anther development and dehiscence, thus causing defective phenotypes associated with these processes.