miR167d-ARFs module regulates flower opening and stigma size in rice

doi:10.21203/rs.3.rs-1631437/v1

Download PDF

Research Article

miR167d-ARFs module regulates flower opening and stigma size in rice

https://doi.org/10.21203/rs.3.rs-1631437/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Flower opening and stigma exertion are two critical traits for cross-pollination during seed production of hybrid rice (Oryza sativa L.). In this study, we demonstrate that the miR167d-ARFs module regulates stigma size and flower opening that is associated with the elongation of stamen filaments and the cell arrangement of lodicules, resulting in cleistogamy. The overexpression of miR167d (OX167d) resulted in failed elongation of stamen filaments, increased stigma size, and morphological alteration of lodicule. Blocking miR167 by target mimicry also led to a morphological alteration of the individual floral organs, including a reduction in stigma size and alteration of lodicule cell morphology, but did not show the cleistogamous phenotype. In addition, the four ARFs target genes of miR167d, namely ARF6, ARF12, ARF17, and ARF25, have overlapping functions in flower opening and stigma size. The loss-of-function of a single ARF gene did not influence the flower opening and stigma size, but arf12 single mutant showed a reduced plant height and aborted apical spikelets. However, mutation in ARF12 together with mutation in either ARF6, ARF17, or ARF25 led to the same defective phenotypes that were observed in OX167d, including the failed elongation of stamen filaments, increased stigma size, and morphological alteration of lodicule. These findings indicate that the appropriate expression of miR167d is crucial and the miR167d-ARFs module plays important roles in the regulation of flower opening and stigma size in rice.

Oryza sativa

miR167d

auxin response factor

flower opening

stigma size

Rice (Oryza sativa L.) is one of the most important food crops in the world. The improvement in grain yield has been a major focus of crop breeding programs. The yield of rice grains is determined by four major components, including panicle numbers per plant, spikelet number per panicle, filling rate of the grains, and the grain weight (Sakamoto & Matsuoka, 2008). Spikelets and the inner floral organs are crucial to the filling rate of grains and grain weight (Chen et al., 2022). In addition, successful hybrid rice production needs a sufficient number of pollen grains deposited on the stigma lobes of the seed parent (Virmani et al., 2002). The floral traits and flowering behavior traits influence outcrossing in rice (Virmani et al., 2002), including sigma size, style length, stigma exertion, stamen filament length, and flower opening. (Virmani et al., 2002). Thus, elucidating the regulatory mechanisms of flower opening and stigma size are important to outcrossing in rice and genetically improve seed production of hybrid rice.

microRNAs (miRNAs), a class of short non-coding RNAs, play essential roles in various biological processes by regulating the expression of their target genes (Zheng & Qu, 2015). A variety of miRNAs have been reported to determine various agronomic traits of rice. For example, miR156 targets SQUAMOSA promoter-binding protein-like transcription factor (SPL) genes that are involved in the regulation of various agronomic traits (Zheng & Qu, 2015). Overexpression of miR156 results in rapid leaf/tiller initiation and precocious leaf maturation (Xie et al., 2012). A point mutation of OsSPL14 at the miR156-target site, which perturbs the cleavage of the OsSPL14 transcripts by miR156, showed enhanced grain yields (Jiao et al., 2010; Miura et al., 2010). OsmiR397 promotes panicle branching and enlarges grain size, leading to an increase in the grain yield of up to 25% in the field by targeting OsLAC, which encodes a laccase-like protein involved in the sensitivity of plants to brassinosteroids (Zhang et al., 2013). OsmiR444a suppresses tillering by negatively regulating OsMADS57, which encodes a MADS-domain family transcription factor (Guo et al., 2013). Overexpression of miR172 leads to a decrease in seed weight and floral defects by delaying the transition from spikelet meristem to floral meristem (Zhu et al., 2009). In addition, a coordinated gene network formed by miR172-APETALA2 (AP2) and miR156/miR529-SPL2 is involved in the regulation of panicle branching and rice tillering (Wang et al., 2015). The OsmiR396-growth regulating factor (GRF) module regulates grain size by directly regulating miR408 (Yang et al., 2021). A novel miR167a-OsARF6-OsAUX3 module regulates grain length and width (Qiao et al., 2021). OsmiR167a-targeted auxin response factors (OsARF12, -17, and − 25) control the tiller angle (Li et al., 2020). However, the full roles of miR167 in the regulation of agronomic traits are largely unclear.

The phytohormone auxin is known to play an essential role in regulating agricultural traits in rice, including root system, tillers, panicle architecture, flower development, and seed development (Wang et al., 2017). For example, the overexpression of OsAUX1, the putative auxin influx carrier gene, led to an increase in the number of lateral roots (Zhao et al., 2015). Tillering in rice is affected by an alteration in auxin signaling. More tillers were produced following the accumulation of miR393, which targets the OsTIR1/OsAFB auxin receptor genes (Xia et al., 2012). Loss of function in the auxin biosynthesis gene FIB led to small panicles and a reduced number of spikelets (Yoshikawa et al., 2014).

The auxin response factors (ARFs) are a transcription factor family that play crucial roles in auxin-signaling pathway through their specifical binding to auxin response elements (AuxRE) in the promoter of auxin response genes (Li et al., 2016). Based on genome-wide analysis studies, ARFs from 15 plant species have been identified (Li et al., 2016). For example, 22 ARF genes and a pseudogene have been identified in Arabidopsis, and 25 genes have been identified in rice (Li et al., 2016). It has been reported that ARFs play an essential role in several aspects of plant growth and development. In Arabidopsis, ARF1 and ARF2 are involved in controlling leaf senescence and the abscission of floral organs (Ellis et al., 2005). AFR8 negatively regulates fertilization and fruit initiation (Goetz et al., 2006). In rice, OsARF6 and OsARF17 are involved in controlling flag leaf angle by regulating secondary cell wall biosynthesis of the lamina joints (Huang et al., 2021). OsARF19 regulates rice leaf angles by positively regulating brassinosteroid insenstive1 (OsBRI1) and OsGH3-5 (Zhang et al., 2015a). In tomato (Solanum lycopersicum L.), SlARF3 is reported to play multiple roles in development and is associated with the formation of trichomes and epidermal cells (Zhang et al., 2015b).

A rice spikelet contains a flower subtended by a pair of glumes that are aborted florets. Spikelets are the primary reproductive organs that determine grain yield. A rice flower consists of a lemma and palea, two lodicules, six stamens, and one pistil with two stigmas. Several genes involved in spikelet development have been characterized. For example, OsMADS34, a SEPALLATA-LIKE gene, controls rice inflorescence and spikelet development (Gao et al., 2010). The SEPALLATA MADS box gene PAP2, together with three APETALA1 (AP1)-like genes, plays a role in the specification of inflorescence meristems (Kobayashi et al., 2012). OsMADS1 functions in the development of rice flowers (Jeon et al., 2000). In addition, mutants showing cleistogamy have been characterized. A mutant lacking of lodicules inside the spikelets showed cleistogamy (Maeng et al., 2006). A rice mutant cl7(t) developed by ethyl methane sulfonate mutagenesis (EMS) exhibited cleistogamy, contained thicker stems and more vascular bundles, had closed spikelets and a reduction in plant height, and the weak swelling ability of the lodicules led to cleistogamy (Ni et al., 2014). However, our knowledge of cleistogamy is limited.

In a previous study, we demonstrated that miR167d functions in rice immunity against Magnaporthe oryzae (Zhao et al., 2020). Here, we investigated the roles of miR167d in regulating flower opening and stigma size. First, we examined the flower opening and morphology of individual flower organs in overexpressing and blocking miR167d transgenic lines, including stigma size, elongation of stamen filaments, and morphological features of lodicules. Further, we obtained single or double mutants for the four ARFs that are the target genes of miR167d using CRISPR/Cas9 system or genetic crosses. Then, we investigated the flower opening and morphology of individual flower organs in these mutants. Our data demonstrate that miR167, together with four ARFs, had crucial roles in regulating stigma size, elongation of stamen filaments, and morphological features of lodicules. These findings establish a novel function of the miR167-ARFs module in regulating flower opening and stigma size.

Plant Materials and Growth Conditions

The rice lines used in this study included OX167d, MIM167d (transgenic line expressing a target mimic of miR167d), arf6, arf12, arf17, arf25, arf6 arf12, arf12 arf17, and arf12 arf25. All these lines were in the Kasalath (KA) background. Besides, we also included OX167d in Zhonghua 11 backgrounds (ZH11, a japonica cultivar), which was generated by using the same construct and method as previously described (Zhao et al., 2020). OX167d/KA, MIM167d, arf12, and arf25 were from a previous study (Zhao et al., 2020). The arf6, arf17, and arf6 arf12 mutants were generated using CRISPR/Cas9 technology in the KA background. The arf12 arf17 and arf12 arf25 double mutants were obtained by genetic crosses and confirmed by PCR-based genotyping.

CRISPR/Cas9 Plasmids Construction and Mutant Screening

To generate the ARF6, ARF17, and ARF6 ARF12 knockout mutants, we constructed the plasmid as previously described (Xie et al., 2015) with slight modifications. Briefly, the target sites for each gene were screened and designed using the CRISPR-PLANT system (https://www.genome.arizona.edu/crispr/CRISPRsearch.html). The primers that contained the gene-specific sequences were designed using standards that were previously described (Xie et al., 2015). The DNA fragments to generate the ARF6 knockout construction were amplified using the following primer pairs: PRGEB32-S5AD5-F/ARF6-gR1-R, ARF6-gR1-F/ARF6-gR2-R, and ARF6-gR2-F/PRGEB32-S3AD5-R. Together with the pRGEB32 plasmid, they were assembled following the manufacturer’s instructions for the NEB® Golden Gate Assembly Kit (New England Biolabs, Ipswich, MA, USA). The ARF17 and ARF6 ARF12 knockout constructs were generated using the same methods as those for ARF6 by just changing the gene-specific primers. All the positive constructs were confirmed using PCR and sequencing and were transformed into KA using Agrobacterium-mediated transformation with strain GV3101. To confirm the genotypes of knockout lines, all the transgenic T0 lines were verified using hygromycin (Hyg). DNA was then extracted from the Hyg-positive lines using the CTAB method and used for PCR amplification with gene-specific primers. The PCR products were sequenced, and the genotype of each mutation line was confirmed. All the primers for constructions are shown in Additional file 1: Table S1.

Morphological Analysis and Microscopic Observations

The phenotypic characteristics of panicles, spikelets, and the floral organs were investigated using a DSLR camera (Nikon, Tokyo, Japan), an Olympus (Tokyo, Japan) stereoscope, and a microscope (Carl Zeiss AG, Jena, Germany), respectively. During the mature stage, the phenotypes of all the lines were investigated using a DSLR camera. The traits were measured as previously described (Li et al., 2022).

For paraffin sectioning, the spikelets from the different lines at the heading stage were fixed in FAA (50% [v/v] ethanol, 5% [v/v] acetic acid, and 3.7% [v/v] formaldehyde), vacuumed for 15 min, and incubated overnight at 4 ℃. The fixed spikelets were dehydrated by gradient ethanol, infiltrated with xylene, embedded into paraffin (Sigma-Aldrich, St. Louis, MO, USA), cut into thick slices, and dyed with 1% (w/v) safranin after they were pasted on microscope slides. Next, the slices were dewaxed, rehydrated, and dyed with 0.5% (w/v) toluidine blue. The microscopic observations were performed using a Zeiss microscope.

RNA Isolation and Reverse Transcription-Quantitative Real-Time PCR (RT-qPCR) Assay

The plant tissues were collected and used for total RNA extraction as previously described (Zhao et al., 2021; Hassan et al., 2022). cDNA was synthesized using NovoScript® Plus All-in-one 1st Strand cDNA Synthesis SuperMix (gDNA Purge) (Novoprotein, Shanghai, China). RT-qPCR was performed using Taq Pro Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) and the indicated primers to determine the abundance of mRNA on a Bio-Rad CFX96^TM Real-Time system (Hercules, CA, USA). The abundance of gene expression was normalized using the ubiquitin (UBI) gene as an internal standard. The same method was used to determine the abundance of miR167d in ZH11 as previously described (Zhao et al., 2020). All the primers for the RT-qPCR assay are shown in Additional file 1: Table S1.

Overexpression of miR167d Results in Cleistogamy in KA

The miR167 family is conserved in planta based on annotations of the miRBase (https://www.mirbase.org/) (Kozomara et al., 2018), which is a searchable database of miRNA. There are 10 MIR167 genes in rice that generate two types of mature miR167 sequences with only a single nucleotide difference at the 3’ end, i.e., miR167a-c and miR167d-j. In a previous study, we reported that miR167d functions in rice immunity against rice blast (Zhao et al., 2020). In this study, we describe our observation of the transgenic plants overexpressing miR167d (OX167d).

One of the most striking phenotypes of the OX167d plants was cleistogamy. During flowering time, the anthers were not observed outside of the lemma and palea in OX167d compared with those in KA. Continual observation of spikelets in OX167d from the heading to the filling stages indicated that the anthers were pushed out of the lemma and palea in KA flowers by the elongation of filaments and remained outside of them for a few days (Fig. 1a). In contrast, no anthers were observed outside of the lemma and palea in OX167d (Fig. 1a) except for the black spots inside the spikelets (Fig. 1a). When we peeled off parts of the lemma and palea in OX167d, the dead anthers were clearly displayed (Fig. 1b). The anthers always remained inside the spikelets even during the early, middle, and late filling stages of OX167d (Fig. 1c). Thus, these results suggest that the overexpression of miR167d results in cleistogamy.

In addition, the seedling growth and plant height were impacted by the overexpression of miR167. The seedling height of the OX167d lines was significantly reduced compared with that of KA (Additional file 2: Fig. S1a). At the heading stage, the plant height of OX167d lines was reduced to approximately 50% that of KA (Additional file 2: Fig. S1b, c). Each of the four internodes was obviously shorter in OX167d than those in KA (Additional file 2: Fig. S1d). The OX167d lines showed larger tiller angles than those of KA at both the seedling and heading stages (Additional file 2: Fig. S1a, b, e, g), which were associated with the curve at the base of each tiller (Additional file 2: Fig. S1f). The number of tillers was significantly reduced in OX167d compared with KA (Additional file 2: Fig. S1h). Thus, our results suggest that overexpression of miR167d results in a reduction in plant height, larger tiller angle, and fewer tillers in rice.

Furthermore, the OX167d displayed aborted apical spikelets, panicle enclosure, and delayed maturity (Additional file 2: Fig. S2a). Moreover, the size of panicles was significantly reduced in OX167d compared with those in KA (Additional file 2: Fig. S2b, c). The seed setting rate was significantly reduced in OX167d compared with KA (Additional file 2: Fig. S2d). The 1,000-grain weight decreased significantly in OX167d compared with that in KA (Additional file 2: Fig. S2e). The grain width and length in OX167d were comparable with those in KA (Additional file 2: Fig. S2f, h, j, l). However, the width and length of brown rice grain were significantly reduced compared with those in KA (Additional file 2: Fig. S2g, i, k, m). These results suggest that the overexpression of miR167d affects plant morphology and yield components.

Overexpression of miR167d Alters the Elongation of Stamen Filaments, Stigma Size, and Lodicule

To elucidate the reason for the cleistogamy of OX167d flowers, we examined the floral organs in detail. The stamens of OX167d had no apparent significant difference from that of KA at the heading stage (Fig. 2a), except that the lodicules in OX167d seemed white in contrast to the watery and transparent lodicules in KA (indicated by the red arrows) (Fig. 2b). However, the stamen filament failed in elongation at the flowering stage in OX167d (Fig. 2c). The size of stigma significantly increased in OX167d compared with that in KA (Fig. 2d, e). Since the lodicule enables flowers to open by swelling to push the lemma and the palea (Honda et al., 2005), its morphology is associated with cleistogamy (Yoshida et al., 2007). The microscopic observations showed that the size of the lodicule in OX167d was comparable with that in KA (Fig. 2f). However, the cell arrangement of lodicules in OX167d was significantly crowded compared with that in KA (Fig. 2g), and the cell width significantly decreased in OX167d, but not the cell length (Fig. 2h). The cross-sections revealed that OX167d showed a narrow shape, but the KA showed the plump shape of lodicule (Fig. 2i). Therefore, the cleistogamy of OX167d resulted from the defective lodicule. Collectively, overexpression of miR167d leads to defects in stamen filament elongation and lodicules narrowing, but increased stigma size.

Overexpression of miR167d Results in Cleistogamy in ZH11

To exclude that the effects of miR167d on flower opening and stigma size are dependent on genetic background, we constructed the overexpressing miR167d lines in ZH11 (a japonica cultivar) (hereafter, named OX167d/ZH11) and obtained 17 transgenic lines that showed the same phenotypes. The plant height was significantly reduced compared with that of ZH11 (Additional file 2: Fig. S3a). The lengths of four internodes from the top were shorter in OX167d/ZH11 than in ZH11 (Additional file 2: Fig. S3b). Thus, we selected three lines in which there was a significant increase in the accumulation of mature miR167d for further study (Additional file 2: Fig. S3c). As previously described (Zhao et al., 2020), miR167d suppresses its target genes at the transcriptional level (Zhao et al., 2020). Thus, we subsequently examined the expression of miR167d target genes and found that the abundance of mRNA for each gene was significantly reduced compared with that in ZH11 (Additional file 2: Fig. S3d), indicating that miR167d was successfully overexpressed and functioned.

Next, we examined the inner floral organs in detail. The stamens in OX167d/ZH11 had no apparent difference from those of ZH11 (Additional file 2: Fig. S4a). However, the stigma size of OX167d/ZH11 significantly increased (Additional file 2: Fig. S4b), and the stamen filament failed in elongation compared with that in ZH11 at the flowering stage (Additional file 2: Fig. S4c, d). In addition, the anthers remained inside the spikelets (Additional file 2: Fig. S4e). These results indicate that the influence of miR167d on flower opening and stigma size is independent on genetic background. Therefore, we used KA to conduct in-depth research.

Blocking miR167 by Target Mimicry Results in Morphological Alteration of Stigma and Lodicule

To further clarify the functions of miR167d in flower opening and stigma size, we used two independent transgenic lines designated MIM167d that overexpressed a target mimic of miR167 from a previous study (Zhao et al., 2020), which led to a significant reduction in the accumulation of miR167d. In MIM167d and KA flowers, the anthers were pushed out of the lemma and palea by the elongation of filaments (Fig. 3a), indicating normal flowering. The anthers of MIM167d exhibited no significant difference from those of KA (Fig. 3b). The filament elongation of MIM167d was comparable to that of KA (Fig. 3c). However, the size of stigma in MIM167d was significantly reduced compared with that in KA (Fig. 3d, e). In addition, microscopic observation showed that the size of the lodicule in MIM167d was comparable with that in KA (Fig. 3f). However, the cell width of the lodicule, but not the cell length, was significantly increased in MIM167d compared with that in KA (Fig. 3g, h). The cross-sections revealed that the shape of lodicule in MIM167d was similar to that in KA (Fig. 3i). These results suggest that blocking miR167d results in morphological alteration of stigma and lodicule but cannot result in cleistogamy in rice.

In addition, the survey of agronomic traits demonstrated that the plant height of MIM167d was significantly reduced compared with that in KA (Additional file 2: Fig. S5a, b). The tiller angles were comparable with those of KA (Additional file 2: Fig. S5c, d). However, the tiller number of MIM167d was significantly higher than that of KA (Additional file 2: Fig. S5e). The size of the panicles was significantly smaller in MIM167d than in KA (Additional file 2: Fig. S6a, b). The seed setting rate and 1,000-grain weight were significantly decreased in MIM167d compared with those in KA (Additional file 2: Fig. S6c, d), leading to straight panicles at the mature stage (Additional file 2: Fig. S5a). In addition, the width and length of grain and brown rice grain were significantly reduced in MIM167d compared with those in KA (Additional file 2: Fig. S6e-h, i-l). These results suggest that blocking miR167d causes defects in agronomic traits.

Four AFR Genes Have Overlapping Functions in Flower opening and Stigma Size

Previous studies have shown that miR167d has four target genes, namely ARF6, ARF12, ARF17, and ARF25, which encode auxin response factors (Zhao et al., 2020). This prompted us to identify which target genes function in flower opening and stigma size, particularly those that result in cleistogamy. Therefore, we generated single mutants for each of these four ARF genes. Among them, arf12-1, arf12-2, arf25-1, and arf25-2 were created in a previous study (Zhao et al., 2020). The mutants of ARF6 and ARF17 were constructed using CRISPR/Cas9 technology (Additional file 2: Fig. S7). Fortunately, we obtained two independent homozygous mutants for each gene, including arf6-1 and arf6-2 for ARF6, and arf17-1 and arf17-2 for ARF17 (Additional file 2: Fig. S7). All the mutants harbored deletions that resulted in a frameshift and caused protein truncation (Additional file 2: Fig. S7). Next, we assessed the phenotype of spikelet organs in these mutants. In KA and all the ARF single mutant flowers, the anthers were out of the lemma and palea and remained outside (Fig. 4a-c), and no anthers remained inside the spikelets (Fig. 4d). The elongation of filaments was normal in the single mutants compared with that in KA (Fig. 5a). Microscopic observation showed that the size of stigma in the ARF single mutants was comparable with that in KA (Fig. 5b, c). The cross-sections revealed that the shape of the lodicule in the ARF single mutants was similar to that in KA (Fig. 5d). Thus, these results suggest that ARF6, ARF12, ARF17, and ARF25 may have functional redundancy in regulating flower opening and stigma size.

In addition, the plant height and the number of tillers per plant were similar in the arf6, arf17, and arf25 mutants compared with those in KA (Additional file 2: Fig. S8). However, the plant height of the arf12 mutants was significantly reduced compared with that of KA (Additional file 2: Fig. S8a, b), but there were no significant differences in the tiller numbers (Additional file 2: Fig. S8c). Meanwhile, the arf12 mutants displayed aborted apical spikelets (Fig. 4b). These results indicate that ARF12 may play a major role in plant height and the development of spikelets among the ARFs targeted by miR167.

To confirm whether the four ARFs were functionally redundant, we constructed double mutants of the ARF genes, i.e., arf6 arf12, arf12 arf17, and arf12 arf25, because only the ARF12 single mutant lines resulted in defects in development of spikelets and the plant height. Among the double mutants, arf12-1 arf17-1 and arf12-1 arf25-1 were obtained by genetic crosses using arf12-1, arf17-1, and arf25-1 single mutants. arf6 arf12 was constructed using CRISPR/Cas9 technology. We obtained two independent double homozygous mutants, i.e., arf6 arf12-1 and arf6 arf12-2 (Additional file 2: Fig. S7). Both of the mutant lines harbored deletion or insertion that resulted in a frameshift and caused the proteins to be truncated (Additional file 2: Fig. S7). Consistently, the aborted apical spikelets were observed in all the double homozygous mutants compared with that in KA (Fig. 6a, b). Moreover, no anthers were observed outside of the lemma and palea (Fig. 6c, d) but remained inside the spikelets in the double mutant (Fig. 6e). Conversely, the stamen filaments elongated out of the lemma and palea in KA (Fig. 6c, d). In addition, the filaments of all three double mutants failed in elongation at the flowering stage (Fig. 7a). Microscopic observations showed that the size of the stigma had increased significantly compared with that of KA (Fig. 7b, c). The cross-sections revealed that the double mutants showed a narrow, but the KA showed the plump lodicule (Fig. 7d). These results suggest that ARF6, ARF12, ARF17, and ARF25 have an overlapping function in regulating flower opening, stigma size, and cleistogamy.

Furthermore, we assessed the plant phenotypes of these double mutants. It showed that plant height and the number of tillers per plant were significantly reduced in the double mutants compared with those in KA (Additional file 2: Fig. S9). These results suggest that ARF6, ARF12, ARF17, and ARF25 have an overlapping function in regulating plant height and tiller number in rice.

The miR167d-ARFs Module is Crucial to Regulate Flower opening and Stigma Size

miR167 is ubiquitous across all the terrestrial plant species from ferns to monocotyledons and eudicotyledons, and it is generally highly expressed next only to miR156 and miR166 (Chávez Montes et al., 2014). This suggests that miR167 may have conservative and crucial functions. Flower opening is essential for cross-pollination in seed production of hybrid rice. A previous study suggested that miR167 regulates pistil development by targeting InARF8 in Japanese morning glory (Ipomoea nil) (Glazińska et al., 2014). In this study, we identified the miR167d-ARFs as key regulatory elements in flower opening and stigma size in rice. Overexpression of miR167d led to increased stigma size and cleistogamy that is associated with failure of filament elongation and morphological alteration of lodicule (Figs. 1 and 2). In contrast, blocking miR167d led to decreased stigma size and alteration of lodicule cell morphology (Fig. 3). In addition, double mutants of ARFs showed morphological alteration similar to those of OX167d compared with wild type (Figs. 6 and 7). Thus, the miR167d-ARFs module is crucial to regulating the flower opening and stigma size.

In flowering plants, anther dehiscence following the release of pollen is essential for sexual reproduction. In Arabidopsis, miR167 and its target genes, ARF6 and ARF8, function in anther growth and dehiscence (Zheng et al., 2019). Thus, future studies should focus on examining the function of the miR167d-ARFs module in anther growth and dehiscence in rice.

ARF6, ARF12, ARF17, and ARF25 Have Overlapping Functions

There are 25 ARFs genes in the rice genome (Wang et al., 2007). These genes encode the auxin response factors that specifically bind to the AuxREs found in the promoters of auxin response genes that contain TGTCTC. These ARFs contain a conserved DNA binding domain at the N-terminal, a central region, and a protein-protein interaction domain at the C-terminal (Wang et al., 2007; Li et al., 2020). The function of most remains unclear (Li et al., 2020). In this study, we showed that ARF6/12/17/25 have overlapping functions in regulating flower opening, stigma size, and agronomic traits in rice.

Single knockout mutants for these four ARFs showed the normal flower opening and filament elongation, and comparable stigma size with KA (Figs. 4 and 5). In addition, the agronomic trait survey showed that the plant height and tiller number in the mutants were similar to those in KA (Fig. S8), except for arf12, which showed a significant reduction in plant height and aborted apical spikelets (Fig. 4b and Additional file 2: Fig. S8), indicating that ARF12 has a major role in regulating plant height and apical spikelet development among these ARFs. However, a previous study showed that the plant height of arf12 in the ZH11 background did not differ significantly compared with that of the WT (Li et al., 2020). We hypothesized that this could result from the difference in plant heights between KA and ZH11. Thus, the arf12 mutant in ZH11 did not show a significant difference in plant height compared with that of ZH11.

These four ARFs belong to the same subgroup in rice (Wang et al., 2007), also suggesting they are redundantly functional. In this study, mutation in ARF12 together with the mutation in either ARF6, ARF17, or ARF25 exhibited aborted apical spikelets, failed elongation of stamen filaments, increased stigma size, morphological alteration of lodicule, and cleistogamy (Figs. 6 and 7). In addition, the plant height and tillers per plant were significantly reduced in the double homozygous mutants compared with those in KA (Additional file 2: Fig. S9). Therefore, ARF6, ARF12, ARF17, and ARF25 have overlapping functions in flower opening, stigma size, plant height, and tiller number in rice.

Cleistogamy Maybe due to the Modulation of the Lodicule Development Genes

Cleistogamy is defined as fertilization without flower opening (Ohmori et al., 2018). It is an effective strategy to prevent “gene flow” from the cultivated species to their wild relatives and helps the plants to survive under unfavorable conditions (Yoshida et al., 2007; Lombardo et al., 2017). It has been reported that altering lodicule morphology can result in cleistogamy because the defective organs are unable to exert sufficient outward pressure to trigger the flower opening (Lombardo et al., 2017). In this study, we demonstrated that overexpressing miR167d or knocking out ARF6 ARF12, ARF12 ARF17, and ARF12 ARF25 results in morphological alteration of the lodicule and leads to cleistogamy (Figs. 1, 2, 6 and 7). To our knowledge, this study represents the first characterization of a miR167d-ARFs module involved in cleistogamy in rice.

To date, more than 10 genes have been identified that are involved in cleistogamy. However, only two of them had been cloned, including SUPERWOMAN 1 (SPW1) (Yoshida et al., 2007) and CL7(t)/DEP2 (Ni et al., 2014). SPW1 is a class-B MADS-box gene, which specifies the identities of lodicules and stamens (Yoshida et al., 2007), and CL7(t) could participate in the development of lodicules (Ni et al., 2014). Based on genetic and molecular studies, several MADS-box transcription factors were involved in the development of lodicules, including OsMADS1, OsMADS2, and OsMADS3 (Yoshida et al., 2007). Suppression of the expression OsMADS2 results in flowers that do not open (Yadav et al., 2007), suggesting that manipulating the expression of lodicule development genes could result in cleistogamy. Thus, future studies should focus on examining the function of ARFs in regulating lodicule development genes.

The Appropriate Expression of miR167d is Crucial for Agronomic Traits

It has been reported that miR167 is involved in the regulation of plant vegetation, flowering time, reproductive organ development, and stress response by the regulation of its target genes, which encode the ARFs (Liu et al., 2021). For example, in Arabidopsis, miR167 targets ARF6 and ARF8 and controls somatic embryogenesis, seed development, and adventitious rooting (Gutierrez et al., 2009; Su et al., 2016; Yao et al., 2019). Overexpression of target mimics of miR167 results in a late flowering phenotype (Todesco et al., 2010). In rice, overexpression of miR167b results in a decreased abundance of mRNA of the four ARF genes, and the transgenic lines were small in stature and had remarkably reduced tiller numbers (Liu et al., 2012). In tomatoes, the downregulation of ARF6 and ARF8 by miR167 leads to female sterility and defects in floral development (Liu et al., 2014). In Japanese morning glory (Ipomoea nil), miR167 and InARF8 participate in vegetative and generative development (Glazińska et al., 2014). In this study, we identified miR167d as a key regulator of agronomic traits in rice.

Previous studies suggest that the overexpression of miR167a in the ZH11 background results in reduced plant height, fewer tiller numbers, and larger tiller angles (Li et al., 2020). Overexpression of miR167b in the cultivar Nipponbare resulted in a shorter stature and fewer tillers (Liu et al., 2012). In this study, similar phenotypes were observed in the miR167d lines that were overexpressed in the KA background (Additional file 2: Fig. S1). In addition, we observed that the overexpression or inhibition of the abundance of miR167d affected rice agronomic traits, including plant height, the tiller number, panicle and spikelet development, grain width and length, seed setting rate, and 1,000-grain weight (Additional file 2: Figs. S1, S2, S5 and S6). Among them, some agronomic traits were defective in both OX167d and MIM167d. For example, the plant height was significantly reduced in both OX167d and MIM167d (Additional file 2: Figs. S1 and S5). Yield-related traits were defective in both lines (Additional file 2: Figs. S2 and S6). Thus, our research demonstrates that the appropriate expression of miR167d is crucial for agronomic traits.

miR167d May Function Differently from miR167a-c and/or in Different Genetic Backgrounds

Though miR167 is a highly conserved miRNA family in plants, the number of MIR167 genes are various among different plant species. For example, there are four MIR167 genes in Arabidopsis, designated MIR167a to MIR167d, and miR167a is highly abundant and acts as a key regulator in the development of female and male organs (Yao et al., 2019). There are 10 MIR167 genes in rice, each of which has different abundance. MIR167a-c probably plays a primary role because the abundance of miR167a-c was higher than that of miR167d-j when the genomic sequences covering the stem-loop region of each MIR167 gene were transiently expressed in tobacco leaves (Liu et al., 2012). However, miR167d-j is induced higher than miR167a-c in Lijiangxin Tuan Heigu upon biotic stress (Li et al., 2014). Previous studies suggested that miR167a overexpression lines produced longer grains and had increased 1,000-grain weight in the WT/DJ background (Qiao et al., 2021). The overexpression of MIM167a has no effect on the plant height in the ZH11 background (Li et al., 2020). However, in this study, overexpressing miR167d lines significantly reduced the 1,000-grain weight in the KA background (Additional file 2: Fig. S2). Moreover, the brown rice grain width and length were significantly reduced compared with those of KA (Additional file 2: Fig. S2), and the MIM167d lines significantly reduced the plant height (Additional file 2: Fig. S5a, b). Thus, miR167d may function differently with miR167a-c and/or in different genetic backgrounds.

In the present study, we demonstrated that overexpressing miR167d, or knockout of ARF12, together with ARF6, ARF17, or ARF25, resulted in failed elongation of stamen filaments, increased stigma size, morphological alteration of lodicule. In contrast, blocking miR167 expression led to a reduction in stigma size and alteration of lodicule cell morphology. Furthermore, ARF6, ARF12, ARF17, and ARF25 have overlapping functions in flower opening and stigma size. Thus, appropriate expression of miR167d and miR167d-ARFs module are crucial for agronomic traits in rice.

miRNA: microRNA; SPL: SQUAMOSA promoter-binding protein-like transcription factor; GRF: growth-regulating factor; ARF: auxin response factor; AuxRE: auxin response elements; AP1: APETALA1; EMS: methane sulfonate mutagenesis; KA: Kasalath; ZH11: Zhonghua 11; Hyg: hygromycin; RT-qPCR: Reverse transcription-quantitative real-time PCR.

Acknowledgments

We thank Xin-Hui Zhou (Sichuan Agricultural University) for technical support.

Authors’ contributions

W.-M. W. and Z.-X. Z. conceived this project, designed the experiments, and wrote the manuscript. Z.-X. Z. and X.-X. Y. performed most of the experiments with the support of S. L., Y.-T. P., X.-L. Y., C. C., B. H., S.-X. Z., M. P., J.-H. Z., X.-H. H., G.-B. L., H. W., J.-W. Z., Y.-Y. H., J. F., and Y. L.. All the authors read and approved the final manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (U19A2033 to W.-M. W., 32101728 to Z.-X. Z.), the China Postdoctoral Science Foundation (2020M683346 to Z.-X. Z.), and the Sichuan Youth Science and Technology Innovation Research Team Foundation (2022JDTD0023).

Availability of Data and Materials

Not applicable.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

All the authors declare no conflict of interests.

Chávez Montes RA, Rosas-Cárdenas dFF, De Paoli E, Accerbi M, Rymarquis LA, Mahalingam G, Marsch-Martínez N, Meyers BC, Green PJ, de Folter S (2014) Sample sequencing of vascular plants demonstrates widespread conservation and divergence of microRNAs. Nat Commun 5: 3722
Chen S, Lou S, Zhao X, Zhang S, Chen L, Huang P, Li G, Li Y, Liu Y-G, Chen Y (2022) Ectopic expression of a male fertility gene, LOGL8, represses LOG and hinders panicle and ovule development. Crop J
Ellis CM, Nagpal P, Young JC, Hagen G, Guilfoyle TJ, Reed JW (2005) AUXIN RESPONSE FACTOR1 and AUXIN RESPONSE FACTOR2 regulate senescence and floral organ abscission in Arabidopsis thaliana. Development 132: 4563-4574
Gao X, Liang W, Yin C, Ji S, Wang H, Su X, Guo C, Kong H, Xue H, Zhang D (2010) The SEPALLATA-Like Gene OsMADS34 Is Required for Rice Inflorescence and Spikelet Development. Plant Physiol 153: 728-740
Glazińska P, Wojciechowski W, Wilmowicz E, Zienkiewicz A, Frankowski K, Kopcewicz J (2014) The involvement of InMIR167 in the regulation of expression of its target gene InARF8, and their participation in the vegetative and generative development of Ipomoea nil plants. J Plant Physiol 171: 225-234
Goetz M, Vivian-Smith A, Johnson SD, Koltunow AM (2006) AUXIN RESPONSE FACTOR8 Is a Negative Regulator of Fruit Initiation in Arabidopsis. Plant Cell 18: 1873-1886
Guo S, Xu Y, Liu H, Mao Z, Zhang C, Ma Y, Zhang Q, Meng Z, Chong K (2013) The interaction between OsMADS57 and OsTB1 modulates rice tillering via DWARF14. Nat Commun 4: 1566
Gutierrez L, Bussell JD, Păcurar DI, Schwambach Jl, Păcurar M, Bellini C (2009) Phenotypic Plasticity of Adventitious Rooting in Arabidopsis Is Controlled by Complex Regulation of AUXIN RESPONSE FACTOR Transcripts and MicroRNA Abundance. Plant Cell 21: 3119-3132
Hassan B, Peng Y-T, Li S, Yin X-X, Chen C, Gulzar F, Zhou S-X, Pu M, Ji Y-P, Wang Y-P, Zhao W, Huang F, Peng Y-L, Zhao Z-X, Wang W-M (2022) Identification of the blast resistance genes in three elite restorer lines of hybrid rice. Phytopathol Res 4: 15
Honda I, Turuspekov Y, Komatsuda T, Watanabe Y (2005) Morphological and physiological analysis of cleistogamy in barley (Hordeum vulgare). Physiol Plant 124: 524-531
Huang G, Hu H, van de Meene A, Zhang J, Dong L, Zheng S, Zhang F, Betts NS, Liang W, Bennett MJ, Persson S, Zhang D (2021) AUXIN RESPONSE FACTORS 6 and 17 control the flag leaf angle in rice by regulating secondary cell wall biosynthesis of lamina joints. Plant Cell 33: 3120-3133
Jeon J-S, Jang S, Lee S, Nam J, Kim C, Lee S-H, Chung Y-Y, Kim S-R, Lee YH, Cho Y-G, An G (2000) leafy hull sterile1 Is a Homeotic Mutation in a Rice MADS Box Gene Affecting Rice Flower Development. Plant Cell 12: 871-884
Jiao Y, Wang Y, Xue D, Wang J, Yan M, Liu G, Dong G, Zeng D, Lu Z, Zhu X, Qian Q, Li J (2010) Regulation of OsSPL14 by OsmiR156 defines ideal plant architecture in rice. Nat Genet 42: 541-544
Kobayashi K, Yasuno N, Sato Y, Yoda M, Yamazaki R, Kimizu M, Yoshida H, Nagamura Y, Kyozuka J (2012) Inflorescence Meristem Identity in Rice Is Specified by Overlapping Functions of Three AP1/FUL-Like MADS Box Genes and PAP2, a SEPALLATA MADS Box Gene. Plant Cell 24: 1848-1859
Kozomara A, Birgaoanu M, Griffiths-Jones S (2018) miRBase: from microRNA sequences to function. Nucleic Acids Res 47: D155-D162
Li S-B, Xie Z-Z, Hu C-G, Zhang J-Z (2016) A Review of Auxin Response Factors (ARFs) in Plants. Front Plant Sci 7
Li Y, Li J, Chen Z, Wei Y, Qi Y, Wu C (2020) OsmiR167a-targeted auxin response factors modulate tiller angle via fine-tuning auxin distribution in rice. Plant Biotechnol J 18: 2015-2026
Li Y, Li T-T, He X-R, Zhu Y, Feng Q, Yang X-M, Zhou X-H, Li G-B, Ji Y-P, Zhao J-H, Zhao Z-X, Pu M, Zhou S-X, Zhang J-W, Huang Y-Y, Fan J, Wang W-M (2022) Blocking Osa-miR1871 enhances rice resistance against Magnaporthe oryzae and yield. Plant Biotechnol J n/a
Li Y, Lu Y-G, Shi Y, Wu L, Xu Y-J, Huang F, Guo X-Y, Zhang Y, Fan J, Zhao J-Q, Zhang H-Y, Xu P-Z, Zhou J-M, Wu X-J, Wang P-R, Wang W-M (2014) Multiple Rice MicroRNAs Are Involved in Immunity against the Blast Fungus Magnaporthe oryzae. Plant Physiol 164: 1077-1092
Liu H, Jia S, Shen D, Liu J, Li J, Zhao H, Han S, Wang Y (2012) Four AUXIN RESPONSE FACTOR genes downregulated by microRNA167 are associated with growth and development in Oryza sativa. Funct Plant Biol 39: 736-744
Liu N, Wu S, Van Houten J, Wang Y, Ding B, Fei Z, Clarke TH, Reed JW, van der Knaap E (2014) Down-regulation of AUXIN RESPONSE FACTORS 6 and 8 by microRNA 167 leads to floral development defects and female sterility in tomato. J Exp Bot 65: 2507-2520
Liu X, Huang S, Xie H (2021) Advances in the regulation of plant development and stress response by miR167. Front Biosci 26: 655-665
Lombardo F, Kuroki M, Yao S-G, Shimizu H, Ikegaya T, Kimizu M, Ohmori S, Akiyama T, Hayashi T, Yamaguchi T, Koike S, Yatou O, Yoshida H (2017) The superwoman1-cleistogamy2 mutant is a novel resource for gene containment in rice. Plant Biotechnol J 15: 97-106
Maeng J-Y, Won Y-J, Piao R, Cho Y-I, Jiang W, Chin J-H, Koh H-J (2006) Molecular mapping of a gene ‘ld(t)’ controlling cleistogamy in rice. Theor Appl Genet 112: 1429-1433
Miura K, Ikeda M, Matsubara A, Song X-J, Ito M, Asano K, Matsuoka M, Kitano H, Ashikari M (2010) OsSPL14 promotes panicle branching and higher grain productivity in rice. Nat Genet 42: 545-549
Ni D-H, Li J, Duan Y-B, Yang Y-C, Wei P-C, Xu R-F, Li C-R, Liang D-D, Li H, Song F-S, Ni J-L, Li L, Yang J-B (2014) Identification and utilization of cleistogamy gene cl7(t) in rice (Oryza sativa L.). J Exp Bot 65: 2107-2117
Ohmori S, Koike S, Hayashi T, Yamaguchi T, Kuroki M, Yoshida H (2018) The cleistogamy of the superwoman1-cleistogamy1 mutation is sensitive to low temperatures during the lodicule-forming stage. Breeding Science 68: 432-441
Qiao J, Jiang H, Lin Y, Shang L, Wang M, Li D, Fu X, Geisler M, Qi Y, Gao Z, Qian Q (2021) A novel miR167a-OsARF6-OsAUX3 module regulates grain length and weight in rice. Mol Plant 14: 1683-1698
Sakamoto T, Matsuoka M (2008) Identifying and exploiting grain yield genes in rice. Current Opinion in Plant Biology 11: 209-214
Su YH, Liu YB, Zhou C, Li XM, Zhang XS (2016) The microRNA167 controls somatic embryogenesis in Arabidopsis through regulating its target genes ARF6 and ARF8. Plant Cell, Tissue Organ Cult 124: 405-417
Todesco M, Rubio-Somoza I, Paz-Ares J, Weigel D (2010) A Collection of Target Mimics for Comprehensive Analysis of MicroRNA Function in Arabidopsis thaliana. PLoS Genet 6: e1001031
Virmani SS, Mao C, Toledo R, Hossain M, Janaiah A (2002). Hybrid rice seed production technology and its impact on seed industries and rural employment opportunities in Asia: FFTC
Wang D, Pei K, Fu Y, Sun Z, Li S, Liu H, Tang K, Han B, Tao Y (2007) Genome-wide analysis of the auxin response factors (ARF) gene family in rice (Oryza sativa). Gene 394: 13-24
Wang L, Sun S, Jin J, Fu D, Yang X, Weng X, Xu C, Li X, Xiao J, Zhang Q (2015) Coordinated regulation of vegetative and reproductive branching in rice. Proc Natl Acad Sci USA 112: 15504-15509
Wang Y, Zhang T, Wang R, Zhao Y (2017) Recent advances in auxin research in rice and their implications for crop improvement. J Exp Bot 69: 255-263
Xia K, Wang R, Ou X, Fang Z, Tian C, Duan J, Wang Y, Zhang M (2012) OsTIR1 and OsAFB2 Downregulation via OsmiR393 Overexpression Leads to More Tillers, Early Flowering and Less Tolerance to Salt and Drought in Rice. PloS One 7: e30039
Xie K, Minkenberg B, Yang Y (2015) Boosting CRISPR/Cas9 multiplex editing capability with the endogenous tRNA-processing system. Proc Natl Acad Sci USA 112: 3570-3575
Xie K, Shen J, Hou X, Yao J, Li X, Xiao J, Xiong L (2012) Gradual Increase of miR156 Regulates Temporal Expression Changes of Numerous Genes during Leaf Development in Rice. Plant Physiol 158: 1382-1394
Yadav SR, Prasad K, Vijayraghavan U (2007) Divergent Regulatory OsMADS2 Functions Control Size, Shape and Differentiation of the Highly Derived Rice Floret Second-Whorl Organ. Genetics 176: 283-294
Yang X, Zhao X, Dai Z, Ma F, Miao X, Shi Z (2021) OsmiR396/growth regulating factor modulate rice grain size through direct regulation of embryo-specific miR408. Plant Physiol 186: 519-533
Yao X, Chen J, Zhou J, Yu H, Ge C, Zhang M, Gao X, Dai X, Yang Z-N, Zhao Y (2019) An Essential Role for miRNA167 in Maternal Control of Embryonic and Seed Development. Plant Physiol 180: 453-464
Yoshida H, Itoh J-I, Ohmori S, Miyoshi K, Horigome A, Uchida E, Kimizu M, Matsumura Y, Kusaba M, Satoh H, Nagato Y (2007) superwoman1-cleistogamy, a hopeful allele for gene containment in GM rice. Plant Biotechnol J 5: 835-846
Yoshikawa T, Ito M, Sumikura T, Nakayama A, Nishimura T, Kitano H, Yamaguchi I, Koshiba T, Hibara K-I, Nagato Y, Itoh J-I (2014) The rice FISH BONE gene encodes a tryptophan aminotransferase, which affects pleiotropic auxin-related processes. Plant J 78: 927-936
Zhang S, Wang S, Xu Y, Yu C, Shen C, Qian Q, Geisler M, Jiang Da, Qi Y (2015) The auxin response factor, OsARF19, controls rice leaf angles through positively regulating OsGH3-5 and OsBRI1. Plant Cell Environ 38: 638-654
Zhang X, Yan F, Tang Y, Yuan Y, Deng W, Li Z (2015) Auxin Response Gene SlARF3 Plays Multiple Roles in Tomato Development and is Involved in the Formation of Epidermal Cells and Trichomes. Plant Cell Physiol 56: 2110-2124
Zhang Y-C, Yu Y, Wang C-Y, Li Z-Y, Liu Q, Xu J, Liao J-Y, Wang X-J, Qu L-H, Chen F, Xin P, Yan C, Chu J, Li H-Q, Chen Y-Q (2013) Overexpression of microRNA OsmiR397 improves rice yield by increasing grain size and promoting panicle branching. Nat Biotechnol 31: 848-852
Zhao H, Ma T, Wang X, Deng Y, Ma H, Zhang R, Zhao J (2015) OsAUX1 controls lateral root initiation in rice (Oryza sativa L.). Plant Cell Environ 38: 2208-2222
Zhao Z-X, Feng Q, Cao X-L, Zhu Y, Wang H, Chandran V, Fan J, Zhao J-Q, Pu M, Li Y, Wang W-M (2020) Osa-miR167d facilitates infection of Magnaporthe oryzae in rice. J Integr Plant Biol 62: 702-715
Zhao Z-X, Xu Y-J, Lei Y, Li Q, Zhao J-Q, Li Y, Fan J, Xiao S, Wang W-M (2021) ANNEXIN 8 negatively regulates RPW8.1-mediated cell death and disease resistance in Arabidopsis. J Integr Plant Biol 63: 378-392
Zheng L-L, Qu L-H (2015) Application of microRNA gene resources in the improvement of agronomic traits in rice. Plant Biotechnol J 13: 329-336
Zheng L, Nagpal P, Villarino G, Trinidad B, Bird L, Huang Y, Reed JW (2019) miR167 limits anther growth to potentiate anther dehiscence. Development 146
Zhu Q-H, Upadhyaya NM, Gubler F, Helliwell CA (2009) Over-expression of miR172 causes loss of spikelet determinacy and floral organ abnormalities in rice (Oryza sativa). BMC Plant Biol 9: 149

No competing interests reported.

TableS1.Theprimersusedinthisstudy.xlsx
Additional file 1: Table S1. The Primers used in this study.
2022.5.5Fig.S.pptx
Additional file 2: Fig. S1. Phenotypes of OX167d in the KA background. a, b Plant morphology of KA and OX167d lines at the seedling stage (a) and heading stage (b). c Comparison of the plant height between KA and OX167d lines. d Comparison of the internode length. e Vertical observation. f The curve observation. g, h Comparison of the tillers angle (g) and tiller number (h). In (c), (g), (h), error bars indicate the standard deviation (SD) (n=10). **P<0.01 (Student’s t-test). Fig. S2. Panicle morphology and yield components of OX167d in the KA background. a, b Panicle morphology of the indicated lines. Pictures were taken before (a) and at the mature stage (b). In (a), red arrows indicate panicle enclosure, and white arrows indicate aborted apical spikelets. Bar, 8 cm in (b). c-e Comparison of the panicle length (c), seed setting rate (d), and 1,000-grain weight (e). f-i Analysis of grain and brown rice grain phenotypes in the indicated lines. Grain width (f), brown rice grain width (g), grain length (h), and brown rice grain length (i). j-m Comparison of the grain width (j), brown rice grain width (k), grain length (l), and brown rice grain length (m) in the indicated lines. In (c-e, j-m), error bars indicate the standard deviation (SD) (n=10). **P<0.01 (Student’s t-test). Fig. S3. Phenotypes of OX167d in the ZH11 background. a Plant morphology of ZH11 and OX167d lines at the heading stage. b The internode length comparison of ZH11 and OX167d lines. c Analysis of the expression of miR167d in the indicated lines by RT-qPCR. U6 was used as an internal reference. d Relative mRNA abundance of the genes targeted by miR167d in the indicated lines. ARF6: Os02G0164900; ARF12: Os04g0671900; ARF17: Os06g0677800; ARF25: Os12g0613700. In (c, d), error bars indicate the standard deviation (SD) (n=3). **P<0.01 (Student’s t-test). Fig. S4. Spikelet characteristics of ZH11 and OX167d lines. a, b Spikelets of the ZH11 and OX167d lines. The stamen was removed in (b) from (a). Red arrows indicate stigmas. c, d Comparison of filament elongation in the indicated lines. e Comparison of the inner organs at the grain filling stage. Parts of the palea and lemma have been removed to allow the observation of inner organs. White arrows indicate the stigmas, and green arrows indicate the stamens. Fig. S5. Phenotypes of MIM167d in the KA background. a Plant morphology of KA and MIM167d lines at the heading stage. b Comparison of the plant height. c The curve observation. d, e Comparison of the tillers angle (d) and tiller number (e) in the indicated lines. In (b, d, e), error bars indicate the standard deviation (SD) (n=10). **P<0.01 (Student’s t-test). Fig. S6. Panicle morphology and yield components in the MIM167d lines. a Panicle morphology of KA and MIM167d lines. Pictures were taken at the mature stage. Bar, 8 cm. b-d Comparison of the panicle length (b), seed setting rate (c), and 1,000-grain weight (d) of the indicated lines. e-h Analysis of grain and brown rice grain phenotypes in the indicated lines. Grain width (e), brown rice grain width (g), grain length (f), and brown rice grain length (h). i-l Comparison of the grain width (i), grain length (j), brown rice grain length (k), and brown rice grain width (l) in the indicated lines. In (b-d, i-l), the error bars indicate the standard deviation (SD) (n=10). **P<0.01 (Student’s t-test). Fig. S7. The genotypes of the indicated knockout mutants. Alignment of the gDNA sequences around the gRNA sites in the indicated mutants. Two typical lines were obtained for each single or double mutant. The target sites are highlighted in yellow. The protospacer adjacent motifs are highlighted in red. Numbers on the right side indicate the deleted or inserted base pairs of the mutant lines indicated. Fig. S8. Phenotypes of ARF single mutants. a Plant morphology of the indicated lines at the heading stage, including arf6, arf12, arf17, and arf25. b, c Comparison of the plant height (b) and tiller number (c) in the indicated lines. Error bars indicate the standard deviation (SD) (n=10). **P<0.01 (Student’s t-test). Fig. S9. Phenotypes of ARFs double mutants. a Plant morphology of the indicated lines at the heading stage, including arf6 arf12, arf12 arf17, and arf12 arf25. b, c Comparison of the plant height (b) and tiller number (c) in the indicated lines. Error bars indicate the standard deviation (SD) (n=10). **P<0.01 (Student’s t-test).

Download PDF

Editorial decision: Major revision
22 May, 2022
Reviews received at journal
17 May, 2022
Reviewers agreed at journal
11 May, 2022
Reviewers invited by journal
10 May, 2022
Editor assigned by journal
09 May, 2022
Submission checks completed at journal
09 May, 2022
First submitted to journal
06 May, 2022

You are reading this latest preprint version

miR167d-ARFs module regulates flower opening and stigma size in rice

Status:

Version 1

Abstract

Figures

Background

Materials And Methods

Results

Discussion

Conclusions

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1