Two allelic recessive mutants of DH2, designated dh2-1, and dh2-2, were identified in this study. Given that dh2-1 showed typical defects, the remainder of this study focuses on the dh2-1 mutants.
3.1. dh2 mutantsexhibited defective lemmas
In rice, the wild-type (WT) spikelet was composed of a pair of rudimentary glumes, a pair of sterile lemmas and a fertile terminal floret, which contained four whorls of floral organs from outside to inside. Lemma and palea were shell-shaped and have 5 and 3 vascular bundles, respectively. The lemma had two inward hook-like edges, when the palea had two inward hook-like structures between the body of palea and margin region of palea. Then they hooked together perfectly to provided a protective environment for the internal floral organs (Fig. 1-A1, A2, A3).
In all the dh2-1 spikelets, it was observed that the lemma and palea could not hooked together mainly because of defective lemma. Based on the lemma phenotype, then dh2-1 spikelets were classified into two types. In type Ⅰ spikelets, although it still showed shell-shaped, the lemma becomed significantly narrower than that of the WT. (Fig. 1-B1, B2, B3). In type Ⅱ spikelets, the whole lemma almost degraded into a rod-shaped awn structure, losing its shell-shaped structural characteristics and exposing the inner floral organs (Fig. 1-C1, C2, C3).
For further clarifying the defects in dh2-1 lemma, the histology was analyzed between WT and dh2-1 mutants. Firstly, the lemma of WT had 5 vascular bundles and wider than the palea. However, the dh2-1 lemma became narrower, because the numbers of cells and vascular bundles were significantly reduced (Fig. 1-A8, A9). Secondly, the WT lemma had different epidermis cell types between the abaxial and the adaxial sides. The abaxial epidermis covered a silicified cell layer bearing trichomes and protrusions, when the adaxial epidermis consist of a layer of vacuolated cells (Fig. 1-A7, A8, A9). In lemma of type I spikelets, the abaxial and the adaxial sides were similar to the wild type. However, in lemma of type II spikelets, the ectopic silicified cells were observed in the adaxial epidermis, even some lemmas completely lost epidermis and showed rod-shaped with a silicified abaxial epidermis (Fig. 1-B7-B9, C7-C9). In general, the dh2-1 mutant showed defectivemedial-lateral and adaxial-abaxial polarities of lemma to different degrees, which suggests that DH2 plays a vital role in the establish of lemma polarity.
3.2 The defects of lemma in dh2-1 occured in the early stage of spikelet development
The development of WT and dh2-1 mutants during the early spikelet development stage was observed by SEM. In the Sp4 stage, the lemma primordia had already formed and the palea primordia appears. In the stage of Sp5-6, while the lodicules and stamens primordia of the WT spikelet began to formed, the lemma and palea primordia continue to differentiate, and the bottom edges of the palea and lemma began to hook together. Then, the lemma and palea primordia continued to differentiate and finally closed completely at the Sp8 (Fig. 2-A, B, C, D). It should be pointed out that the lemma and palea showed shell-shaped from the beginning of their primordia formation.
At the stage of Sp4 when lemma just formed, compared with WT, there were no obvious abnormalities in the dh2-1 mutants (Fig. 2-E, I). However, from the stage of Sp5, the dh2-1 spikelet displayed obvious defects. In type I spikelet, the lemma began to become narrower in the medial-lateral direction than that of WT at the stage of Sp5, although it still showed the shell-shaped (Fig. 2-F). Next, from Sp6 to Sp8, it was observed that the lemma primordia could not differentiate enough in medial-lateral direction and showed narrower than that of WT, resulting the inner organ primordia partially exposed (Fig. 2-G). In type II spikelet, the lemma almost fail to differentiate in medial-lateral direction at the Sp5 and showed rod-like in its upper part (Fig. 2-F). Then, the lemma primordia of the dh2-1 mutants further developed into a rod-like structure, and the stamens and pistil primordia were completely exposed from Sp6 to Sp8 (Fig. 2-K, L). Therefore, the entire data above indicated that the establish of medial-lateral and adaxial-abaxial polarities of lemma were seriously affected at the early stage of spikelet development, which was consistent with the observations at the mature stage.
3.3 dh2-1 mutants exhibited defective leaves
Under normal physiological conditions, the WT leaf blade was flat (Fig. 3- A, B, C), whereas the dh2-1 mutants showed narrow and abaxially-rolled leaves throughout the whole growth period (Fig. 3-A, I, J). At the mature period, the LRIs of the top three leaves was investigated between the WT and dh2-1 mutants. The results showed that the abaxial LRIs of the top three leaves of dh2-1 mutants were significantly increased compared with that of the WT (Fig. 3-P). Meanwhile, it was observed that the width of all top three leaves of dh2-1 mutants were significantly decreased compared with that of the WT (Fig. 3-Q, R). In addition, it was also found that the length of flag leaves in dh2-1 mutants was significantly decreased while that of other leaves showed no obvious difference compared with that of the WT (Fig. 3-Q, R).
For further clarifying the defects of dh2-1 leaves at the histological level, the SEM and freezing microtome section were conducted. Firstly, by using SEM, there were no obvious differences in both adaxial and abaxial surface between the WT and dh2-1 mutants leaves (Fig. 3-D, E, K, L). Next, by analyzing the cross sections, it was found that both the number and area of bulliform cells on the adaxial surface in dh2-1 mutants were increased significantly, between large and small veins, and among small veins, compared with those of the WT (Fig. 3-F, G, H, M, N, O, T, U). As a result, the leaves in dh2-1 were rolled abaxially. We also observed that the numbers of both large and small veins in dh2-1 were significantly decreased compared with that of the WT, which accounting for the narrow leaves. In general, in dh2-1 leaves, the adaxial cell feature was disordered and cell differentiation along medial-lateral direction was repressed, which means DH2 also played crucially important role in the establish of leaf polarity in rice.
Besides of exhibiting abaxially rolled leaves and rod-like lemmas, dh2-1 showed a series of growth and development defects. Compared with WT, the plant height of dh2-1 was lightly reduced (Fig. S3-A), the panicle length was significantly reduced (Fig. S3-B, G), the number of primary and secondary branches, the number of spikelets per panicle , the number of grains per panicle, and the seed setting rate were both decreased (Fig. S3-H-L). Meanwhile, we noted that mature grains of dh2-1 were also affected. Both the length and width of grain and brown grain of dh2-1 were significantly reduced, which resulted the decrease of the 1000-grain weight of dh2-1. (Fig. S3-C, D, E, F, M, N, O, P, Q, R). Taken together, the abnormal development of lateral organs has a significant effect on yield in rice.
3.4 Similarity of Defects in the dh2-2 Mutants to Those of the dh2-1 Mutant
The dh2-2 mutants showed very similar defects in morphology, histology, and agronomy, compared with those of the dh2-1 mutant (Fig. S1, Fig. S2, Fig. S4, Fig. S5). However, some differences were observed. dh2-2 spikelets were also classified into two types. In type Ⅰ spikelets, the degree of narrowing of the lemma in dh2-2 was stronger than that of dh2-1, and dh2-2 produced the elongated rod-like awn at the top of lemma (Fig. S1-B1, Fig. 1-B1). In the same planting environment, the changes in dh2-2 mutants are more severe compared to that of dh2-1 mutants. The spikelet of dh2-2 mutants usually presented a completely rod-shaped lemma, while the lemma of dh2-1 mutants was only slightly narrowed (Fig. S1-C1, Fig. 1-C1). The LRIs of the top three leaves of dh2-2 mutants leaves were higher than that of dh2-1 mutants leaves (Fig. S4-P, Fig. 3-P1). At the same time, the mature brown grains of type I dh2-1 spikelets had a complete grain morphology, while the morphology of the mature brown grains of dh2-2 mutants had changed. The shape of the brown grains of dh2-2 mutants were similar to that of water droplets (Fig. S5-D, E; Fig. S3-D, F). Taken together, these observations demonstrate that the dh2-2 mutant showed more severe defects than the dh2-1 mutant.
3.5. DH2 encode a Rice Argonaute7 (AGO7) protein
Previous studies showed that DH2 was located between the INDEL markers of IND7 and IND12 on chromosome 3, and there are total of 13 annotation genes related to flower development within this interval (Guo et al. 2013). In the present study, we performed gene cloning and sequence alignment analysis among those 13 annotation genes between WT and dh2 mutants. Sequencing analysis showed that a C to A transversion in the 1661th base of second exon of LOC_Os03g33650 occurred in dh2-1, which causes the replacement of amino acids Pro to Gln. Meanwhile, in dh2-2, there was a C-T transition at the 2861th base of the third exon of LOC_Os03g33650, leading to Ser-Phe substitution of the 221th amino acid (Fig. 4-A).
To further determine whether the phenotype of dh2 were due to the mutation of LOC_Os03g33650, we constructed the complementary vector containing 7623bp genomic sequence of wild-type LOC_Os03g33650 and introduced it into dh2-1 mutants. It was found that in positive transgenic lines, the mutation phenotype of dh2-1 was restored completely. Therefore, the results verified that the DH2 gene was LOC_Os03g33650 (Fig. 4-B).
We also utilize CRISP-CAS9 technology to perform targeted mutagenesis in coding frame of the LOC_Os03g33650 in the ZH11 (a geng-type) background. Compared with the ZH11, the leaves of DH2-cas9-1 lines were slightly abaxial curled (Fig. S4-A-E2). Meanwhile, the width and length of DH2-cas9-1 leaves were significantly decreased (Fig. S4-O, P, Q). Furthermore, the number and area of bulliform cells on the adaxial surface of DH2-cas9-1 were also increased (Fig. S4-F1, G1, F2, G2). Compared with the spikelet of ZH11 (Fig. S4-H1, I1, J1), it was obvious that the development of lemma in DH2-cas9-1 was affected, which resulted in lemma and palea primordia could not closed completely (Fig. S4-H2, I2, J2). It was worth noting that we found the plant height, and the length of each internode of DH2-cas9-1 were significantly decreased. Meanwhile, the panicle length, the number of primary and secondary branches, the number of spikelets per panicle, the number of grains per panicle, and the setting rate of DH2-cas9-1 were both decreased. Both the length and width of grain and brown grain of DH2-cas9-1 were significantly reduced, which resulted the decrease of the 1000-grain weight of DH2-cas9-1 (Fig. S7-A-P). Those phenotypes of DH2-cas9-1 were similar to the phenotypes in type I of dh2-1, and dh2-2, further verifying that LOC_Os03g33650 was the DH2 gene (Fig. S4-K1, L1, M1, K2, L2, M2).
DH2 encodes an Argonaute-like protein, which has four conserved domains, Piwi_ago-like domain, PAZ_argonaute_like domain, ArgoL1 domain, and ArgoN domain (Fig.4-C). Piwi_ago-like domain is the C-terminal portion of Argonaute7, which provides the 5' anchoring of the guide RNA and the catalytic site for slicing. PAZ_argonaute_like domain functions as a nucleic acid binding domain, with a strong preference for single-stranded nucleic acids (RNA or DNA) or RNA duplexes with single-stranded 3' overhangs. Phylogenetic tree analysis showed that DH2 and its homologous genes were widely present, and highly conserved in monocotyledonous and dicotyledonous plants. Those homologous genes could be classified into four Types, and DH2 protein belongs to the Type Ⅰ. Therefore, it indicated that DH2 encodes an OsAGO7 protein and may be act as the similar role with AtAGO7 (Fig. 4-D).
3.6. Overexpression of DH2 Causes Adaxially Rolled Leaves
As analyzing the phenotype of complementary transgenic plants, it was interesting that some of these lines showed adaxially curled leaves, which was just the opposite of abaxially curled leaves in loss-of-function mutants of DH2, in addition to being able to restore the phenotype of the dh2-1 mutant (Fig. 5-A, B, C, D). We further detected the expression level of DH2 in these complementary transgenic lines, OE1 to OE6 and WT. The result displayed that in OE1, OE3, OE5, and OE6 lines, the expression of DH2 increased by at least 100 times, while OE2 and OE4 increased about 50 times (Fig. 5-J). Then, it was suggested that adaxially curled leaves in these transgenic plants might be caused by overexpression of DH2.
For further clarifying the function of DH2, we analyzed the phenotypes of overexpression lines OE-2 and OE-6 in detail. During the mature period, the LRIs of the top three leaves of OE lines were significantly increased. In OE-2 , OE-4 and OE-6 lines, the average values of the LRIs of the top three leaves are 50%, 70% and 75%, respectively (Fig. 5-K), but the LEIs of the top three leaves of OE lines have no obvious changes (Fig. 5-L). In addition, the width and number of large or small vein of leaves in OE lines were extremely significantly decreased while the length of leaves were slightly decreased or unchanged, compared with those in the WT (Fig. 5-M, N, O). Next, by using SEM, it was found that the adaxial surface of OE-2, and OE-6 was similar to that of WT. However, in the abaxial surface of OE-2 and OE-6, prickly hair cells massively and regularly produced on the surface of dumbbell cells, while it was rare in the abaxial surface of WT but generally occurred in the adaxial surface of WT(Fig. 5-E1, E2, E3, F1, F2, F3). Finally, by the freezing microtome sections of the leaves in WT and OE-2, OE-6, it was observed that the number and area of bulliform cells on the adaxial surface between large and small vein, and among small veins, decreased significantly compared with that in WT (Fig. 5-H1, H2, H3, I1, I2, I3, P, Q).
Therefore, the results mentioned above showed that overexpression of DH2 could cause defects on development of abaxial-adaxial and medial-lateral axiles, further verifying that DH2 played crucially important role in the establish of leaf polarity in rice.
3.7. Overexpression of DH2 affects the panicle and grain development
Moderate leaf rolling, which is helpful to establish ideal plant architecture, plays a key role contributing to the high yields in rice. For evaluating the leaf rolling of overexpression of DH2 in breeding, we investigated the agronomy traits of the WT and overexpression lines.
It was found that in OE lines the panicle morphology and grains development were changed. Compared with the WT, The panicle length of OE-2, OE-4, OE-6 significantly decreased (Fig. 6-A, F). The number of primary branches, secondary branches, spikelets per panicle, grain per panicle of OE-2, OE-4 and OE-6 increased significantly. The setting rates of OE-2, OE-4, OE-6 were gently decreased (Fig. 6-G, H, I, J, K). The length and width of grain and brown grain of OE-2, OE-4, OE-6 were delightedly decreased, while the weight of 1000 grains and brown grains were decreased in the OE lines (Fig. 6- B, C, D, E, L, M, N, O, P, Q). The results showed that overexpression of DH2 may affect the development of rice panicles and grain.
3.8. DH2 affects the expressions of tasiR-ARF and OsARFs on the lemma
In order to further understand the functions of DH2 gene, we examined the spatial expression patterns of DH2 by RT-qPCR. The results indicated that DH2 were expressed mainly in leave blades and young panicle, while it showed very low expression level in root, shoot and leave sheath. Further, in spikelets, it was found that DH2 expressed mainly in lemma but not palea, lodicule, pistil (Fig. S8- A, B). Then, the expression pattern was constant with the role of DH2 in regulating the development of leave blades and lemma.
In Arabidopsis, AtAGO7, the DH2 homologous, had been verified to function on the production and silencing complex formation of tasiR-ARFs, then mediate the post-transcriptional gene silencing of AtARF3/4 by specially degrading the mRNA of ARF3/4 on adaxial side of leaves to function on the adaxial identity establishment (Allen et al. 2005; Axtell et al. 2006; Montgomery et al. 2008; Peragine et al. 2004; Iwakawa et al. 2021; Sakurai et al. 2021).
In this study, the expression of tasiR-ARFs was detected firstly by qPCR. The results showed that it was decreased significantly in the leaf, lemma and panicle of dh2-1 compared with that in WT (Fig. 8-A). Then, we detected the expression of OsARF2, OsARF3, OsARF14, and OsARF15, the homologous of the AtARF3/4, between dh2-1 mutants and the WT. The results showed that their expressions were all increased significantly in the leaf and lemma in dh2-1 compared with those in WT (Fig. 8-B, C, D), whereas their expression in the leaf of OE-2, OE-4 were decreased significantly compared with those in the WT (Fig. 8-D). Meanwhile, expression levels of these genes in young panicles were detected. The results showed that their expression were all increased significantly in panicles with different length, compared with those in the WT (Fig. 8-E, F, G, H). Therefore, these data indicated that DH2/OsAGO7 also mediated the formation of tasiR-ARFs, and then repress the expression of ARFs in rice, which may be quite conserved in monocots and dicots.
For further clarify the expression pattern of tasiR-ARFs and its targets OsARFs, we conducted in situ hybridization. The results showed that tasiR-ARFs is widely expressed in the floral organ primordia of the WT (Fig. 8-I1). However, there was no obvious signal of tasiR-ARFs in the rod-like lemma primordia in dh2-1 mutants, while it was still expressed in the other floral organ primordia in dh2-1 mutants (Fig. 8-I2). In the WT spikelets, OsARF2, OsARF3, OsARF14, and OsARF15 shared a similar expression pattern, mainly expressing on the abaxial surface of the margin region of palea, and weak signals of these genes could be detected in the stamens and pistils, whereas they were not almost expressed in the lemma primordia (Fig. 8-K1, L1, M1, N1). In dh2-1 spikelets, the signal of OsARF2, OsARF3, OsARF14, or OsARF15 was slightly enhanced in the abaxial surface of the margin region of palea, pistil and stamen. Most importantly, strong ectopic signals of these genes were detected in the rod-shaped lemma primordia (Fig. 8-K2, L2, M2, N2). In addition, we also detect the expression of OsARF4. In the WT flowers, OsARF4 specifically expresses in the lemma and palea primordia. However, it could not be detected in the rod-like lemma primordia in the dh2-1 spikelets (Fig. 8-J1, J2). Together, these results indicated DH2 was mainly responsible for the repressing the expression of OsARF2, OsARF3, OsARF14 and OsARF15 in lemma in a way of mediating the production of ta-siRNA.