#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Auxin apical dominance governed by the OsAsp1-OsTIF1 complex determines distinctive rice caryopses development on different branches


Authors: Shu Chang aff001;  Yixing Chen aff001;  Shenghua Jia aff001;  Yihao Li aff001;  Kun Liu aff001;  Zhouhua Lin aff001;  Hanmeng Wang aff001;  Zhilin Chu aff001;  Jin Liu aff001;  Chao Xi aff001;  Heping Zhao aff001;  Shengcheng Han aff001;  Yingdian Wang aff001
Authors place of work: Beijing Key Laboratory of Gene Resource and Molecular Development, College of Life Sciences, Beijing Normal University, China aff001;  Academy of Plateau Science and Sustainability of the People's Government of Qinghai Province & Beijing Normal University, Qinghai Normal University, Qinghai, China aff002
Published in the journal: Auxin apical dominance governed by the OsAsp1-OsTIF1 complex determines distinctive rice caryopses development on different branches. PLoS Genet 16(10): e1009157. doi:10.1371/journal.pgen.1009157
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1009157

Summary

In rice (Oryza sativa), caryopses located on proximal secondary branches (CSBs) have smaller grain size and poorer grain filling than those located on apical primary branches (CPBs), greatly limiting grain yield. However, the molecular mechanism responsible for developmental differences between CPBs and CSBs remains elusive. In this transcriptome-wide expression study, we identified the gene Aspartic Protease 1 (OsAsp1), which reaches an earlier and higher transcriptional peak in CPBs than in CSBs after pollination. Disruption of OsAsp1 expression in the heterozygous T-DNA line asp1-1+/–eliminated developmental differences between CPBs and CSBs. OsAsp1 negatively regulated the transcriptional inhibitor of auxin biosynthesis, OsTAA1 transcriptional inhibition factor 1 (OsTIF1), to preserve indole-3-acetic acid (IAA) apical dominance in CPBs and CSBs. IAA also facilitated OsTIF1 translocation from the endoplasmic reticulum (ER) to the nucleus by releasing the interaction of OsTIF1 with OsAsp1 to regulate caryopses IAA levels via a feedback loop. IAA promoted transcription of OsAsp1 through MADS29 to maintain an OsAsp1 differential between CPBs and CSBs during pollination. Together, these findings provide a mechanistic explanation for the distributed auxin differential between CPBs and CSBs to regulate distinct caryopses development in different rice branches and potential targets for engineering yield improvement in crops.

Keywords:

Arabidopsis thaliana – Auxins – Endoplasmic reticulum – Endosperm – Gene expression – Mesophyll – Rice – Transcriptional control

Introduction

Rice (Oryza sativa) is a staple food for approximately two-thirds of the global population. Due to the application of dwarf breeding and hybrid breeding in the past 50 years, rice grain yield has more than doubled in most parts of the world, and even tripled in certain countries and regions [1]. However, these newer rice cultivars have not entirely fulfilled their high-yield potential due to smaller grain size and poorer grain filling in caryopses located on proximal secondary branches (CSBs) than in those of apical primary branches (CPBs). This problem is exacerbated in recently bred ‘super’ rice cultivars [2].

Rice caryopsis development begins with a double-fertilisation event, following which endosperm mainly contributes to grain size and weight in two continuous developmental progresses: early-stage endosperm cell proliferation and differentiation, and late-stage grain filling. Extensive studies have shown that grain size and weight on a single spikelet are regulated by the histone acetylation and deacetylation pathway [3], the ubiquitin–proteasome pathway [46], G-protein signalling [7, 8], MAPK signalling [9], and phytohormone signalling, such as that involving brassinosteroids [1012] or auxins [13, 14], which influence spikelet hull and endosperm growth. Grain filling is associated with carbohydrate supply, activity of the enzymes involved in starch metabolism, levels of various endogenous hormones, and environmental conditions [1518]. Studies using omics approaches have revealed differences in the expression of various proteins [19], in the expression of genes related to starch synthesis and hormone signalling [20], and in miRNA levels between CPBs and CSBs during the middle and late stages of grain filling after fertilisation [21, 22]. However, the molecular mechanisms determining the developmental differences between CPBs and CSBs remain elusive. In this transcriptome-wide expression study, we sought to identify the key genes that are differentially expressed in post-anthesis CPBs and CSBs, and further characterised the role of one of these genes, Aspartic Protease 1 (OsAsp1), in determining the distinctive development of CPBs and CSBs in rice.

Results

Transcriptome-wide expression analysis of CPBs and CSBs to identify OsAsp1

To identify the key factors that determine the developmental differences between CPBs and CSBs in rice, we performed a transcriptome-wide expression analysis of CPB and CSB samples obtained at four different dates. Samples were chosen separately for each branch type to ensure morphological similarity of CPBs and CSBs, given that primary branches develop earlier than secondary branches. Specifically, samples were collected at 0, 5, 12, and 20 days after heading (DAH) for CPBs (CPB-0, -5, -12, and -20) and 5, 12, 25, and 35 DAH for CSBs (CSB-5, -12, -25 and -35) (S1A Fig).

The gene expression profiles formed two separate clusters, one consisting of CPB-0, CSB-5, CPB-5, and CSB-12 and the other of CPB-12, CSB-25, CPB-20, and CSB-35, among which two pairs, CPB-5/CSB-12 and CPB-20/CSB-35, were closer than the others (Fig 1A). This pattern is consistent with the existence of two continuous phases of caryopsis development: endosperm cell proliferation, which mainly contributes to grain enlargement, followed by grain filling. We also found that CPB endosperm had a higher cell reproduction rate and faster grain filling rate than that of CSBs (S1B and S1C Fig), which is consistent with the findings of previous studies [17, 23]. We propose that the difference in endosperm cell proliferation rates between CPBs and CSBs is the critical factor for grain production because it determines the sink capacity for caryopsis. Therefore, we further compared the genes that were more highly expressed in CPB-0 and CSB-12 than in CSB-5 and CPB-5 (Fig 1A) with the results of the gene expression ratio distribution analysis, based on reads per kb per million mapped reads (RPKM) values. Then, we identified candidate genes with expression ratios higher than 4.8 in both comparisons (Fig 1B and S1 Table) and used quantitative reverse-transcription polymerase chain reaction (qRT–PCR) to confirm the relative transcript levels of the top three genes (S2 Fig): Os07g39020, encoding the SUBTILISIN-LIKE PROTEASE OsSUB53; Os06g41030, encoding an unknown protein containing a DUF1680 domain; and OsAsp1, a previously characterised nucellin gene (aspartic protease1, Os11g08200) [24]. Next, we characterised the role of OsAsp1 in determining differential development between CPBs and CSBs.

Fig. 1. Disruption of differential OsAsp1 expression eliminated the distinct development of caryopses on different branches.
Disruption of differential <i>OsAsp1</i> expression eliminated the distinct development of caryopses on different branches.
A, Heat map showing differentially expressed genes in CPBs at 0, 5, 12, and 20 days after heading (DAH) and in CSBs at 5, 12, 25, and 35 DAH. The gene expression profiles were deposited in the Sequence Read Archive (accession no. SRP083496). Blue rectangles on heat map indicate genes with higher transcriptional levels in CPB-0 and CSB-12 than in CSB-5 and CPBs5. Colour scale shows log ratio of fold change. B, Distributions of ratios of gene expression levels (blue rectangles, Fig 1A) based on the reads per kb per million mapped reads (RPKM) values for CPB-0 and CSB-5 (red line) and CSB-12 and CPB-5 (green line). *Ratio distribution of a significant coincidence for genes of CPB-0 vs. CSB-5 and CSB-12 vs. CPB-5, centred on a value of 4.8 (≥ 4.8). C, Schematic representation of the two OsAsp1 T-DNA insertions found in the asp1-1 and asp1-2 mutant rice lines. D, Identification of the genotypes of asp1-1 and asp1-2 by polymerase chain reaction (PCR). Wild-type (WT) rice was used as a control. LP, LB, and RP primers used for PCR are listed in S4 Table. E, Fertility of the WT, asp1-1+/–, and asp1-2+/–. Unfertilised and fertilised seeds were counted for each plant. F, Transcription levels of OsAsp1 detected by quantitative reverse-transcription (qRT)-PCR of CPBs and CSBs in the WT, asp1-1, and asp1-2 at 0, 3, 5, and 10 DAF. Solid and dotted lines indicate the transcript levels of OsAsp1 in CPBs and CSBs, respectively. Red, blue, and green lines indicate the transcript levels of OsAsp1 in the WT, asp1-1, and asp1-2, respectively. The expression of OsAsp1 was normalised to that of OsActin1, and relative expression in CPBs in the WT at 0 DAF was set at 1.0. qRT-PCR data are means ± standard deviation (SD) of three independent experiments. G, Dry weights of developed grain from CPBs and CSBs in one whole plant of the WT, asp1-1+/–, and asp1-2+/–. Values are means ± SD (n = 30). Different letters indicate significant differences (P < 0.05, Student’s t-test).

Disrupting OsAsp1 eliminated CPB–CSB developmental differences

To explore the role of OsAsp1 in caryopsis development, we obtained two OsAsp1 T-DNA mutation lines from the Rice Functional Genomic Express Database (RiceGE, http://signal.salk.edu): asp1-1 (PFG_3A-60042) and asp1-2 (PFG_3A-04889) (Fig 1C). Through self-crossing and back-crossing with wild-type (WT) plants, we confirmed that both asp1-1 and asp1-2 are single-insertion and homozygous-lethal mutants, because about 75% of spikelets developed into fertilised seeds and none seeds of these seeds was homozygous (S2 Table, Fig 1D and 1E). We further detected OsAsp1 transcript levels in CPBs and CSBs of the WT, asp1-1+/–, and asp1-2+/–on four different dates, and found that OsAsp1 expression peaked in CPBs at 3 days after flowering (DAF) and in CSBs at 5 DAF in the WT, and that the peak OsAsp1 transcript level of WT CPBs at 3 DAF was higher than that of WT CSBs at 5 DAF (Fig 1F). However, in both CPBs and CSBs of asp1-1+/–and asp1-2+/–, OsAsp1 expression was stable and expression values at different DAF were similar to those at 0 DAF in both CPBs and CSBs of the WT (Fig 1F). These results indicate that T-DNA insertion decreased expression of OsAsp1 in CPBs and CSBs in the independent lines asp1-1+/–and asp1-2+/–. Then, we assessed the developmental parameters of both lines, and found that plant height was lower in both asp1-1+/–and asp1-2+/–than in the WT, but that spikelet number per panicle in asp1-1+/–was similar to that in the WT and greater than that in asp1-2+/–(S3 Fig). Of note, compared to the WT, the difference in the grain weight of the developed seeds between CPBs and CSBs disappeared in both asp1-1+/–and asp1-2+/–, whereas seed grain weights in asp1-1+/–and asp1-2+/–were the same or closer to those of CPBs in the WT (S4 Fig and Fig 1G). These results indicate that OsAsp1 is essential for rice growth and caryopsis development; therefore, we focused on asp1-1+/- in the subsequent analysis.

Next, we observed the ovary development of CPBs and CSBs in the WT and asp1-1+/–at 3 and 5 DAF, respectively, and found that at both time points, asp1-1+/–displayed faster ovary development in CSBs than in the WT, and that the ovary size of CSBs in asp1-1+/–was similar to that of CPBs in both the WT and asp1-1+/–(Fig 2A and S5 Fig). We further used Eosin-B staining to assess the ovary development of CPBs and CSBs in the WT and asp1-1+/–at 3 and 5 DAF, and found that the CSBs in asp1-1+/–showed a higher level of endosperm cellularisation than either CSBs in the WT or CPBs in the WT and asp1-1+/–at both time points (Fig 2B). In addition, CSBs in asp1-1+/–displayed more developed embryos at 5 DAF than either CSBs in the WT or CPBs in the WT and asp1-1+/–(Fig 2B). Moreover, DNA ploidy analysis of CPBs and CSBs in the WT and asp1-1+/–at 3 and 5 DAF by flow cytometry showed that at both time points, the relative 3C plus 6C DNA contents, which represent the numbers of endosperm cells, were significantly higher in CSBs of asp1-1+/–than in CPBs of asp1-1+/–or either CPBs or CSBs of the WT, indicating that CSBs in asp1-1+/–developed better than did CPBs of asp1-1+/–or CPBs and CSBs of the WT (Fig 2C). Previous studies have identified Orysa;CycB1;1 (OsCycB1), Orysa;KRP1 (OsKRP1) and the rice cell cycle switch 52A gene (OsCCS52A) as three marker genes associated with cell cycle and division in rice endosperm [2527]. qRT-PCR analysis showed that these three genes displayed differential expression patterns between CPBs and CSBs in the WT and asp1-1+/–at both 3 and 5 DAF (Fig 2D). Together, these data confirm that disruption of OsAsp1 transcription accelerated CSBs’ developmental progress, thereby relieving the developmental differences between CPBs and CSBs in asp1-1+/–.

Fig. 2. OsAsp1 is required for differential ovary development in CPBs and CSBs.
OsAsp1 is required for differential ovary development in CPBs and CSBs.
A, Ovary phenotypes of CPBs and CSBs in the WT and asp1-1+/–at 3 and 5 days after flowering (DAF). Scale bars, 0.2 cm. B, Ovary developmental status of CPBs and CSBs in the WT and asp1-1+/–at 3 and 5 DAF detected by Eosin B staining. Em, embryo; Es, endosperm. Scale bars, 100 μm. C, Nuclear DNA ploidy distributions of cells from CPBs and CSBs (n = 6) of the WT and asp1-1+/–at 3 and 5 DAF detected by flow cytometry. MC indicates nuclear DNA ploidy exceeding 6C. D, Relative expression levels of OsCycB1, OsKRP1, and OsCCS25A in CPBs and CSBs of the WT and asp1-1+/–at 3 and 5 DAF detected by RT-qPCR. Target gene expression was normalised to that of OsActin1, and relative expression of the target gene in CPBs of the WT was set at 1.0. Values are means ± SD of three independent experiments.

OsAsp1 activated TAA1 transcription to upregulate IAA levels

A prior study showed that auxin biosynthesis couples central cell division for endosperm development to fertilisation in Arabidopsis [28]. Therefore, we monitored the levels of indole-3-acetic acid (IAA) in CPBs and CSBs of the WT and asp1-1+/–using high-performance liquid chromatography and found that in the WT, free IAA content was significantly higher in CPBs than in CSBs at 0, 3, and 5 DAF (Fig 3A), which is consistent with a previous report [16]. However, in asp1-1+/–, there was no difference in IAA levels between CPBs and CSBs; IAA content was very low and almost identical in CPBs and CSBs at 0, 3, and 5 DAF (Fig 3A). In rice, IAA is mainly synthesised via the indole-3-pyruvic acid (IPA) pathway, catalysed by tryptophan aminotransferase (OsTAA) and flavin monooxygenase-like enzymes (OsYUCCAs) [29, 30]. As expected, in the WT, expression of OsTAA1 gradually increased after fertilisation, reaching a peak at 5 DAF in both CPBs and CSBs, with a higher peak in CPBs. However, in asp1-1+/–, the expression peaks of OsTAA1 in both CPBs and CSBs disappeared at 5 DAF, but occurred at 3 DAF. Of note, OsTAA1 expression was lower in CPBs than in CSBs of asp1-1+/–(Fig 3B). By contrast, the expression patterns of OsYUCCA1, which was more highly expressed in CPBs than in CSBs at 5 and 10 DAF, respectively, were similar in CPBs and CSBs in the WT and asp1-1+/–, except that expression levels of OsYUCCA1 were lower in CPBs of asp1-1+/–than in the WT at 10 DAF (Fig 3C). In addition, we heterologously expressed OsAsp1 in DR5::GUS transgenic Arabidopsis thaliana plants, and found upregulated transcription of AtTAA1 and AtYUCCA8, but not AtYUCCA9 (Fig 3D); GUS activities were also strongly increased in root apical meristem of the main and lateral roots of 3-day-old seedlings (Fig 3E and 3F). These results imply that OsTAA1 is upregulated by OsAsp1 to preserve IAA levels and ensure normal development of CPBs and CSBs.

Fig. 3. OsAsp1 activates TAA1 transcription to upregulate IAA levels.
OsAsp1 activates <i>TAA1</i> transcription to upregulate IAA levels.
A, IAA content of CPBs and CSBs of the WT and asp1-1+/–at 0, 3, and 5 DAF. Values are means ± SD of three independent experiments. B, Transcript levels of OsTAA1 detected by qRT-PCR in CPBs and CSBs of the WT and asp1-1+/–at 0, 3, 5, and 10 DAF. OsTAA1 expression was normalised to that of OsActin1, and relative expression in CPBs of the WT at 0 DAF was set at 1.0. qRT–PCR data are means ± SD of three independent experiments. C, Transcription levels of OsYUCCA1 in CPBs and CSBs of the WT and asp1-1+/–on 5 and 10 DAF. OsYUCCA1 expression was normalised to that of OsActin1, and relative expression of OsYUCCA1 in CPBs of the WT at 5 DAF was set at 1.0. Values are means ± SD of three independent experiments. NS, no significant difference; *0.01 < P < 0.5; **0.001 < P < 0.01 (Student’s t-test). D, Expression of OsAsp1, AtTAA1, AtYUCCA8, and AtYUCCA9 in Arabidopsis WT and OsAsp1-overexpressed lines OE-5, OE-15, and OE-19. Target gene expression was normalised to that of AtUBQ10, and relative expression in the WT was set at 1.0. Values are means ± SD of three independent experiments. E, GUS staining of root tips from 3-day-old seedlings of DR5:GUS Arabidopsis, and DR5:GUS and OsAsp1-overexpressed crossed lines OE-5, OE-15, and OE-19. Scale bars, 20 μm. F, GUS staining of lateral roots of 3-day-old seedlings of DR5:GUS, and DR5:GUS and OsAsp1-overexpressed crossed lines OE-5, OE-15, and OE-19. Scale bars, 50 μm.

OsAsp1 dismissed transcriptional inhibition of OsTIF1 to OsTAA1 to sustain IAA levels

OsAsp1 encodes an aspartic protease [24], implying that it may regulate the transcription of OsTAA1 via an uncharacterised transcriptional regulator. As shown in our flowchart of OsTAA1 transcription regulator prediction in rice (Fig 4A), we first used the Arabidopsis co-expression networks of AtTAA1 (http://atted.jp/) [31] to identify a C2H2-type zinc-finger protein (At1G75710) as the candidate (S6A Fig), and then performed a BLAST search of the Rice Annotation Project (RAP) rice protein database (http://rapdb.dna.affrc.go.jp/) to obtain five homologs of At1G75710 that were presumed to regulate OsTAA1 expression (S6B Fig). Based on their sequence homology and high expression in pistils, seeds, and endosperm (S6C Fig), we identified two candidate genes, Os04g59380 and Os04g02510, which we named OsTAA1 Transcription Inhibition Factor 1 (OsTIF1) and OsTIF2, respectively. We further detected the transcription levels of OsTIF1 and OsTIF2 in CPBs and CSBs of the WT and asp1-1+/–, and found that expression of both genes was higher in both CPBs and CSBs of asp1-1+/–than in the WT (Fig 4B). We transiently co-transformed POsTAA1:GUS with either OsAsp1 and OsTIF1 or OsTIF2 (S7A Fig) in rice mesophyll protoplasts, and found that both OsTIF1 and OsTIF2 inhibited the transcription activity of POsTAA1, and OsTIF1 more strongly (Fig 4C). In addition, OsAsp1 relieved the inhibitory effects of OsTIF1 or OsTIF2 on the activity of POsTAA1 in rice mesophyll protoplasts (Fig 4C).

Fig. 4. OsAsp1 dismisses OsTIF1 binding to the OsTAA1 promoter to regulate the transcription of OsTAA1.
OsAsp1 dismisses OsTIF1 binding to the <i>OsTAA1</i> promoter to regulate the transcription of <i>OsTAA1</i>.
A, Bioinformatics flowchart for the prediction of OsTAA1 transcription inhibition factors (OsTIF1 and OsTIF2) in rice. B, Transcript levels of OsTIF1 and OsTIF2 in CPBs and CSBs of the WT and asp1-1+/–at 5 DAF. Target gene expression was normalised to that of OsActin1, and relative expression of the target gene in CPBs of the WT was set at 1.0. Values are means ± SD of three independent experiments. *0.01 < P < 0.5; ***P < 0.001 (Student’s t-test). C, Promoter activity of POsTAA1 under the regulation of OsAsp1 and/or OsTIFs in the protoplast of rice leaves. The POsTAA1:GUS construct was transfected along with PUBI:OsAsp1, PUBI:OsTIFs, or PUBI:OsAsp1 and PUBI:OsTIFs, or with empty vector as control. We included 35S:LUC as an internal control. Relative GUS activity normalised to luciferase activity was adopted as the promoter activity of POsTAA1. Values are means ± SD of three independent experiments. D, Chromatin immunoprecipitation (ChIP)–qPCR was used to detect the enrichment of OsTIF1-binding motifs at positions –2349, –1748, or –377 nt in the promoter region of OsTAA1. OsTIF1-HA was overexpressed in rice calli and purified with HA antibody. NC, negative control. Values are means ± SD of three independent experiments. E, EMSA showing recombinant OsTIF1-His directly binding the promoter region of OsTAA1 at its potential motifs. The predicted OsTIF1 binding motifs at positions –2349, –1748, or –377 nt were labelled with biotin as the probe, and the unlabelled DNA fragments were used as the competitor. Red and black arrows indicate the shifted band and free probe, respectively. F, and G, Transcript levels of OsAsp1, OsTIF1, and OsTAA1 in rice calli overexpressing OsAsp1-GFP (F) or OsTIF1-HA (G). The expression of each target gene was normalised to that of OsActin1, and calli transformed with empty vector were used as a control. qRT-PCR data are means ± SD of three independent experiments.

Next, according to the target sequences of C2H2-type transcription factors reported previously [32], three potential OsTIF1 binding sites at –2349, –1748, and –377 nt were predicted in the 5’ region upstream from the ATG site of OsTAA1 (S7B Fig). Chromatin immunoprecipitation (ChIP)–qPCR showed that sites at positions –2349 and –1748 nt, but not at –377 nt, were enriched by OsTIF1 in transgenic rice calli (Fig 4D). An electrophoretic mobility shift assay (EMSA) showed that recombination OsTIF1-His could separately bind these three sites to produce the band shift visible in Fig 4E. We also found that overexpression of OsAsp1 effectively upregulated transcription of OsTAA1 (Fig 4F), whereas overexpression of OsTIF1 strongly downregulated expression of OsTAA1 in rice calli (Fig 4G). However, overexpression of OsAsp1 and OsTIF1 in rice calli slightly activated transcription of OsTIF1 and OsAsp1, respectively (Fig 4F and 4G). These results demonstrate that OsAsp1 activates OsTAA1 transcription by relieving the inhibition of OsTIF1 to ensure IAA biosynthesis.

IAA induced OsTIF1 ER-nucleus transfer by reducing its interaction with OsAsp1 to form a feedback loop

We further investigated the possibility of an interaction between OsAsp1 and OsTIF1 to characterise the regulation mechanism of OsTAA1 transcription. First, after separately co-transforming OsAsp1-GFP with an ER marker, OsTIF1-GFP with an ER marker, and OsAsp1-GFP with OsTIF1-mCherry, we showed that OsAsp1 and OsTIF1 were co-localised in the ER of rice mesophyll protoplasts (Fig 5A), which raised the question of how OsTIF1 enters the nucleus to inhibit OsTAA1 expression. Due to the involvement of OsTAA1 in IAA formation in rice caryopses, we treated rice mesophyll protoplasts transformed with OsAsp1 and/or OsTIF1 with 50 μm IAA, and found that OsTIF1 was transferred into the nucleus, but that OsAsp1 remained localised in the ER after IAA treatment (Fig 5B). Then, using a yeast two-hybrid assay, we showed that OsAsp1 and OsTIF1 separately formed homodimers and interacted with each other (Fig 5C). Furthermore, IAA inhibited the interaction between OsAsp1 and OsTIF1 in the yeast two-hybrid system (Fig 5C). A co-immunoprecipitation (Co-IP) assay also demonstrated that OsAsp1-GFP interacted with OsTIF1-HA in rice mesophyll protoplasts and that their interaction was significantly depressed upon the addition of IAA (Fig 5D). These results indicate that IAA mediates a regulatory feedback loop involving the OsAsp1–OsTIF1 complex to fine-tune IAA levels between CPBs and CSBs during caryopsis development.

Fig. 5. OsTIF1 transfers from the endoplasmic reticulum (ER) into the nucleus after IAA releases its interaction with OsAsp1.
OsTIF1 transfers from the endoplasmic reticulum (ER) into the nucleus after IAA releases its interaction with OsAsp1.
A, OsAsp1 co-localised with OsTIF1 in the ER of rice mesophyll protoplast. The ER marker was ER-rk CD3-959, purchased from the Arabidopsis Stock Center (http://www.arabidopsis.org/). Scale bars, 5 μm. B, OsTIF1, but not OsAsp1, was transferred from the ER into the nucleus in rice mesophyll protoplasts in response to IAA treatment (50 μM). The nuclear marker was AtBZR2(1–158 aa)-mCherry. Scale bars, 5 μm. C, Interaction between OsAsp1 and OsTIF1 detected by yeast two-hybrid assay.–LH indicates SD medium without LEU and TRP, and–LTH indicates SD medium without LEU, TRP, and HIS. D, Interaction between OsAsp1 and OsTIF1 detected by a Co-IP assay. PUBI:OsTIF1-HA was co-transformed with PUBI:OsAsp1-GFP or PUBI:GFP into rice mesophyll protoplast, and then IAA (50 μM) was added, along with OsTIF1-HA as bait to pull down OsAsp1-GFP. Anti-GFP (a-GFP) was used to detect OsAsp1-GFP or green fluorescent protein (GFP), and anti- hemagglutinin (a-HA) was used to detect OsTIF1-HA.

OsAsp1 transcription was regulated by auxin-inducible MADS29 after pollination

Our analyses showed that IAA levels were always higher in CPBs than in CSBs from the day of pollination (Fig 3A), which is consistent with the phenomenon of auxin apical dominance in plants. Therefore, we hypothesised that IAA regulates OsAsp1 expression to maintain its differential level between CPBs and CSBs in rice. To test this hypothesis, we first cultured unpollinated spikelets (-1 DAF) in vitro with IAA, and found that transcription of OsAsp1 was induced in a concentration-dependent manner (Fig 6A). In addition, in vitro culture of unpollinated spikelets with 2,4-dichlorophenoxyacetic acid (2,4-D) also induced OsAsp1 expression (S8 Fig). Moreover, OsAsp1 transcription in spikelets of CSBs strongly increased after treatment with IAA at 3 days before fertilisation (-3 DAF) compared with OsAsp1 transcription in spikelets given a mock treatment (methanol) (Fig 6B). The grain weight of mature CSB seeds after IAA treatment was also greater than that of CSB seeds treated with mock solution and similar to that of CPB seeds (Fig 6C and S9 Fig). We further treated CPBs of asp1-1+/–with IAA to restore the gradient of IAA levels between CPBs and CSBs, and found that the grain weight of mature CPB seeds of asp1-1+/–after IAA treatment was greater than that of CSB seeds (S10A Fig), and that this difference was similar to the grain weight difference between CPB and CSB seeds in the WT. In addition, the difference in OsAsp1 expression between CPBs and CSBs in asp1-1+/- was also complemented by IAA treatment of CPBs in asp1-1+/- (S10B Fig).

Fig. 6. IAA promotes the transcription of OsAsp1 through MADS29 to maintain differential levels of OsAsp1 in CPBs and CSBs after pollination.
IAA promotes the transcription of <i>OsAsp1</i> through MADS29 to maintain differential levels of OsAsp1 in CPBs and CSBs after pollination.
A, OsAsp1 levels detected by semi-quantitative RT-PCR in spikelets at –1, 0, and 1 DAF, and in detached spikelets at –1 DAF after culture in Murashige and Skoog (MS) medium plus IAA for 24 h. OsActin1 was used as a control. B, Quantification of OsAsp1 levels in CSBs of the WT treated at –3 and 0 DAF with IAA, or with methanol as a mock treatment. OsAsp1 expression was normalised to that of OsActin1, and relative expression of OsAsp1 at –3 DAF in mock-treated samples was set at 1.0. Values are means ± SD of three independent experiments. C, Mature grain dry weight of CPBs and CSBs after the IAA treatment described in B. Values are means ± SD (n = 30). NS, no significant difference; ***P < 0.001 (Student’s t-test). D, Transcript levels of MADS29 and OsAsp1 detected by RT-PCR in CPBs and CSBs of the WT at 0, 3, 5, 12, and 20 DAF. OsActin1 was used as a control. E, Transcript level of MADS29 in CSBs of the WT in response to IAA treatment. CSB spikelets on panicles were treated at –3 and 0 DAF with IAA, or with methanol as mock treatment. MADS29 expresion was normalised to that of OsActin1, and relative expression of MADS29 at –3 DAF in mock-treated samples was set at 1.0. Values are means ± SD of three independent experiments. F, Transcript levels of OsAsp1 in CPBs and CSBs of the WT and Antisense-MADS29 transgenic rice A14 (A-MADS29 A-14) at 5 DAF. OsAsp1 expression was normalised to that of OsActin1, and relative expression in CPBs of the WT was set at 1.0. Values are means ± SD of three independent experiments. NS, no significant difference; ***P < 0.001 (Student’s t-test). G, Mature grain dry weight of CPBs and CSBs of the WT and A-MADS29 A-14. Values are means ± SD (n = 30). ***P < 0.001 (Student’s t-test); NS, no significant difference. H, Constructs used for GUS assay of OsAsp1 promoter activity. The three predicted MADS29 binding sites (CArG box) were labelled in the 1500-bp promoter region of OsAsp1. LUC, firefly luciferase; GUS, beta-glucuronidase. Scale bars, 300 bp. I, Promoter activity of POsAsp1 activated by MADS29 in the protoplast of Arabidopsis leaves. POsAsp1:GUS was co-transfected with MADS29 or an empty vector as a control. 35S:LUC was used as an internal control. Relative GUS activity normalised to luciferase activity was adopted as the promoter activity of POsAsp1. Values are means ± SD of three independent experiments. J, ChIP-qPCR was used to detect the enrichment of CArG boxes by MADS29 at sites –1421, –1217, and –770 of the OsAsp1 promoter region. 35S:MADS29-GFP and POsAsp1:GUS were co-transformed into the protoplast of Arabidopsis leaves, and MADS29-GFP was purified with GFP antibody. NC, negative control. Values are means ± SD of three independent experiments. K, EMSA showing the recombinant MADS29-His directly binding the promoter region of OsAsp1 on its CArG boxes. The CArG boxes at positions –1421, –1217, and –770 were labelled with biotin as the probe, and the unlabelled DNA fragments were used as the competitor. Red and black arrows indicate the shifted band and free probe, respectively.

To further investigate the molecular mechanism regulating OsAsp1 expression, we identified genes that were co-expressed with OsAsp1 using the Rice Oligo Array Database [33], and found five MADS family members with enriched expression (S3 Table). A prior study showed that MADS29 is induced by IAA and is involved in rice seed development [34]. We found that the expression pattern of MADS29 was similar to that of OsAsp1 during the development of CPBs and CSBs after fertilisation (Fig 6D), and that IAA treatment of CSB spikelets at –3 DAF increased MADS29 transcription compared to that in mock-treated spikelets (Fig 6E). Furthermore, OsAsp1 expression was almost the same and downregulated in the CPBs and CSBs of Antisense-MADS29 transgenic rice lines A2 and A14, respectively, but not in A33. (Fig 6F and S11 Fig). The difference in mature grain dry weight between CPBs and CSBs was absent in Antisense-MADS29 lines A2 and A14, but not in A33, which was similar to that of asp1-1+/–(Fig 6G and S12 Fig). Three MADS protein binding sites, CArG boxes at positions –1412, –1217, and –770 nt were predicted in the 5’ region upstream from the ATG site of OsAsp1 (Fig 6H). We found that MADS29 effectively promoted the activity of POsAsp1:GUS in Arabidopsis mesophyll protoplasts (Fig 6I). ChIP–qPCR analysis showed that the CArG box at position –770 nt, but not those at positions –1412 or –1217 nt, was enriched by MADS29-His in transgenic rice calli (Fig 6J). An EMSA further showed that recombinant MADS29-His bound all three of these CArG boxes to produce a specific shifted band (Fig 6K). These results indicate that IAA promotes the expression of OsAsp1 via MADS29 to sustain higher OsAsp1 levels in CPBs than in CSBs during pollination.

Based on these results, we propose a model to explain how OsAsp1 maintains auxin apical dominance to regulate the developmental differential between CPBs and CSBs (Fig 7). In this model, high IAA levels upregulate OsAsp1 transcription via MADS29 at 0 DAF to produce higher OsAsp1 levels in CPBs than in CSBs, which explains why OsAsp1 was identified as differentially expressed in our transcriptome-wide expression analysis. Then, OsAsp1 interacts with OsTIF1 in the ER to dispel OsTIF1 transcriptional inhibition of OsTAA1, causing IAA content to peak at 5 DAF. At the same time, IAA induces transfer of OsTIF1 into the nucleus by reducing its interaction with OsAsp1, causing IAA biosynthesis to decrease, and thereby forming a feedback regulation loop controlling IAA levels. From this model, we conclude that differential IAA levels between CPBs and CSBs are determined by the OsAsp1–OsTIF1 complex, which is essential for the distinct development of CPBs and CSBs.

Fig. 7. Schematic model of the OsAsp1–OsTIF1 complex maintaining IAA apical dominance to establish differential development between CPBs and CSBs in rice.
Schematic model of the OsAsp1–OsTIF1 complex maintaining IAA apical dominance to establish differential development between CPBs and CSBs in rice.

Discussion

Auxin apical dominance has evolutionary advantages for plant growth and reproduction because it plays an essential role in regulating plant shoot architecture and ensuring the production of at least some high-quality seed within a short period [35, 36]. In this study, we demonstrated that differential distribution of auxin determines developmental differences between CPBs and CSBs in rice, causing the grain-filling problem observed in CSBs and thereby reducing total grain production and quality. We found that auxin-dependent transcription of OsAsp1 and interaction of OsAsp1 with OsTIF1 regulate the transcription of OsTAA1 to maintain auxin apical dominance between CPBs and CSBs after fertilisation, providing a potential strategy for enhancing grain size and grain filling and thereby increasing rice yield.

Asps comprise a large endopeptidase family that is widely distributed in all three domains of life [37, 38]. The rice genome contains 96 putative OsAsp genes [39], among which atypical and nucellin-like Asps are characterised by the absence of a plant-specific insert and an unusually high number of cysteine residues, which play regulatory roles in biotic and abiotic stress responses, chloroplast metabolism, and reproductive development [40]. A prior study showed that EAT1 (ETERNAL TAPETUM 1), a basic helix-loop-helix transcription factor, promotes tapetal programmed cell death (PCD) by directly regulating the expression of OsAsp25 and OsAsp37 in rice anthers [41]. However, caspase activity of OsAsp25 or OsAsp37 was not detected, indicating that OsAsp25 and OsAsp37 regulate downstream caspase-like protease(s) involved in PCD rather than directly acting on caspase-specific substrates [41]. Rice nucellin gene OsAsp1 is strongly expressed in zygotic embryos at 1–2 DAF [24]. However, the biological function of OsAsp1 remains unclear. Auxin has been found to induce OsMADS29 transcription to promote PCD of the nucellus and nucellar projection during rice seed development [34]. The present study was the first to determine that OsAsp1 is essential for rice caryopse development because both the asp1-1 and asp1-2 T-DNA insertion mutant lines were homozygous-lethal and, compared to the WT, the difference in mature seed grain weight between CPBs and CSBs disappeared in both lines. Furthermore, we showed that auxin regulates expression of OsAsp1 through OsMADS29 to maintain high levels of OsAsp1 in CPBs, but not in CSBs, after fertilisation, and that OsAsp1 then interacts with OsTIF1 and dismisses its transcriptional inhibition for OsTAA1 to promote IAA biosynthesis. In addition, IAA facilitates OsTIF1 translocation from the ER to the nucleus by releasing the interaction of OsTIF1 with OsAsp1 to form a feedback regulation of IAA levels in caryopses. The disruption of OsAsp1 levels in asp1-1+/- and overexpression of OsAsp1 in rice calli activates transcription of OsTIF1, indicating that OsTIF1 expression is regulated by an uncharacterised and sophisticated mechanism. Therefore, our results demonstrate that OsAsp1 functions as a transcriptional regulator to form a complex with OsTIF1 to regulate the differential development of CPBs and CSBs after fertilisation.

As a critical phytohormone, auxin regulates nearly all aspects of plant development including morphogenesis and adaptive responses [42]. Auxin distribution depends on both auxin metabolism (biosynthesis, conjugation, and degradation) and cellular auxin transport [43, 44]. Auxin is mainly synthesised in shoot apical meristem to maintain increased shoot apical dominance, which transfers to inflorescence meristem when plants flower. In higher plants such as Arabidopsis, maize, and rice, the IPA pathway, catalysed by TAA and YUCCA, has been shown to be the major auxin biosynthetic pathway [29, 30, 4547]. The rice genome contains two OsTAA (synonym of OsTAR) and 14 OsYUCCA genes; increased IAA content during development of rice grains was strongly correlated with the expression of OsTAR1, OsYUCCA9, and OsYUCCA11 [48]. IAA content has been shown to be higher in CPBs than in CSBs at the early grain-filling stage [16], which is consistent with our results. However, few studies have focused on the differential expression of OsTAAs and OsYUCCAs in CPBs and CSBs. In this study, we showed that OsTAA1/OsTAR2 expression gradually increased after fertilisation, reaching a peak at 5 DAF in both CPBs and CSBs in the WT, and with a higher peak in CPBs, which was correlated with the IAA content of CPBs and CSBs in the WT. However, in asp1-1+/–, expression peaks of OsTAA1 disappeared at 5 DAF, but occurred at 3 DAF in both CPBs and CSBs, indicating that OsTAA1 transcription is regulated by OsAsp1 and that OsTAA1 is a key enzyme for IAA biosynthesis, maintaining differential auxin distribution between CPBs and CSBs. By contrast, the expression pattern of OsYUCCA1, which was more highly expressed in CPBs than in CSBs at 5 and 10 DAF, respectively, was similar in the WT and asp1-1+/–, except that OsYUCCA1 expression was lower in CPBs in asp1-1+/–than in the WT at 10 DAF. These results imply that OsYUCCA1 is not a main enzyme for IAA biosynthesis in CPBs and CSBs. Therefore, we will aim to detect the expression levels of 13 other OsYUCCAs of CPBs and CSBs in asp1-1+/–in a future study to determine which OsYUCCA is regulated by OsAsp1 and associated with IAA biosynthesis in CPBs and CSBs. The heterologous expression of OsAsp1 in Arabidopsis improved IAA levels in root apical meristem in 3-day-old seedlings, as well as transcription levels of AtTAA1 and AtYUCCA8, but not AtYUCCA9, providing additional support for the detection of expression of other OsYUCCAs in asp1-1+/-.

In this study, we characterised the effect of differential auxin distribution mediation by OsAsp1 on the developmental difference between CPBs and CSBs. First, we found that two independent T-DNA insertion mutant lines, asp1-1 and asp1-2, present a homozygous-lethal phenotype and that differences in ovary development and size, endosperm cellularisation and development, and mature seed grain weight between CPBs and CSBs disappeared, and were closer to CPBs in asp1-1+/–. These results indicate that OsAsp1 is critical for grain development and maintenance of the developmental difference between CPBs and CSBs. We further showed that IAA content was very low and almost identical in CPBs and CSBs in asp1-1+/–during the early stages of caryopse development, indicating that low auxin levels are sufficient to trigger embryo and endosperm development in CPBs and CSBs, and promote similar development of CSBs and CPBs in asp1-1+/–. In vitro treatment of CPBs in asp1-1+/–with IAA phenocopied the expression difference of OsAsp1 and the development difference between CPBs and CSBs in the WT. These results imply that the higher auxin levels in CPBs than in CSBs, but not absolute IAA concentration, are critical for maintaining differential development of CPBs and CSBs. Treatment of CSBs in the WT with extra IAA before pollination promoted the development of CSBs, allowing the grain weight of CSBs to increase to that of CPBs, which provides strong evidence for differential IAA distribution, but not concentration, determining the developmental difference between CPBs and CSBs in rice. Auxin is produced post-fertilisation in seeds, which drives central cell division and is necessary for MADS-box transcription factor AGL62 (AGAMOUS-LIKE 62)-dependent endosperm development in Arabidopsis [28]. In the present study, we showed that auxin upregulates OsAsp1 transcription via MADS29 after fertilisation to produce higher levels of OsAsp1 in CPBs than in CSBs. The interaction of OsAsp1–OsTIF1 regulated by auxin formed a feedback regulation loop to maintain auxin apical dominance during rice caryopsis development. Therefore, our results present a mechanistic explanation for differential distribution of auxin between CPBs and CSBs to regulate distinct caryopse development patterns in rice.

Materials and methods

Plant materials and growth conditions

The WT rice (Oryza sativa L. subspecies japonica) varieties used in this study were Zhong Hua 15, for RNA sequencing (RNA-seq) and auxin treatment; Zhong Hua 11 for A-MADS29 lines; and Dong Jin for asp1-1 and asp1-2. Two T-DNA mutant lines, PFG_3A-60042 and PFG_3A-04889, were obtained from http://signal.salk.edu, named asp1-1 and asp1-2, respectively, and confirmed by genomic PCR and RT-PCR. Both rice lines were cultivated in a field at Beijing Normal University during the natural growing season for observation of their growth phenotypes during the vegetative and reproductive periods.

RNA-seq and data analysis

Caryopses were collected from growing Zhong Hua 15 rice at different stages after heading for total RNA extraction. Libraries for each sample were constructed using a Next Multiplex RNA Library Prep Set (NEB) and sequenced on an Illumina HiSeq 2000 Platform. Clean reads were mapped to the rice genome (version 7.0) [49] using Cufflink software (http://cufflinks.cbcb.umd.edu/). For gene expression analysis, the numbers of reads per kilobase per million mapped reads (RPKM) were calculated. Differentially expressed genes were defined as having a fold change of ≥ 2 or ≤ 0.5 and a false discovery rate (FDR) ≤ 0.05.

qRT-PCR and semi-quantitative RT-PCR

Total RNA was separately extracted from caryopses at different developmental stages, transgenic rice calli, and 5-day-old Arabidopsis seedlings using TRIzol reagent (Invitrogen, USA) and purified using a PureLink RNA Mini Kit (Invitrogen) and a PureLink DNase Kit (Invitrogen), according to the manufacturer’s protocol. Approximately 2 μg RNA was reverse transcribed using a Reverse-Aid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, USA) and used as template for RT-PCR and RT-qPCR. PCR reactions were carried out in a MyCycler thermal cycler (Bio-Rad, USA) using extension cycling conditions of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s. The amplifications were performed for 28–33 cycles. qRT-PCR was performed using a 7500 Fast Real-Time PCR System (Applied Biosystems, USA) with Power SYBR Green PCR Master Mix (Applied Biosystems). The thermal program was 2 min at 50°C, 10 min at 95°C, followed by 40 cycles of 15 s at 95°C and 60 s at 60°C. The dissociation curve program was used to confirm the specificity of the target amplification product. Three independent biological replicates were performed for RT-qPCR analyses. OsActin1 or AtUBQ10 was used as an internal control in both qRT-PCR and RT-PCR analyses. The primers for the target genes are listed in S4 Table.

Eosin staining

Whole-mount stain-clearing laser-scanning confocal microscopy (WCLSM) was used to evaluate the structure of mature embryo sacs of WT and asp1+/–mutant rice. The spikelets of WT and asp1+/–mutants were collected, fixed, hydrated, stained in Eosin Y water solution, dehydrated, and cleared as previously described [50]. The developmental phenotype of the embryo sacs was examined using confocal laser-scanning microscopy (LSM 700, Zeiss, Germany).

Flow cytometry analysis

Ovaries of WT and asp1-1+/–rice sampled at 3 or 5 DAF were chopped with a razor blade (six ovaries per group). We added 500 μL GS buffer (45 mM MgCl2, 30 mM sodium citrate, 20 mM MOPS, 0.1 v/v Triton X-100) to each group and the mixture was filtered by passing it through a sieve (48 μm2 mesh). DNA was stained by adding 10 mM propidium iodide. The DNA content of the nuclei was measured using a flow cytometer (BD FACSAria II; BD Biosciences) as previously reported [51].

Liquid chromatography–mass spectrometry (LC-MS) analysis

IAA content was quantified using LC–MS as previously reported [52]. Fresh caryopses were collected and ground into powder using liquid nitrogen. Approximately 100 ± 5 mg powder was extracted using 1.5 mL methanol. The supernatant was lyophilised in a freeze-drier (Thermo Fisher Scientific) and then resuspended in 85% (v/v) methanol in water for LC–MS analysis. IAA (Sigma-Aldrich) was used as a standard to prepare the standard curve.

EMSA analysis

The coding sequences (CDSs) of MADS29 and OsTIF1 were separately amplified and cloned into the vector pET32a. MADS29-His and OsTIF1-His were expressed in Escherichia coli (BL21) and affinity-purified using Ni-NTA Resin. Synthesised fragments of the DNA probes (S4 Table) were labelled using a Biotin 3'-End DNA Labeling Kit according to the manufacturer’s instructions (Thermo Fisher Scientific). The protein–DNA binding reactions were performed using a LightShift Chemiluminescent EMSA Kit according to the manufacturer’s instructions (Thermo Fisher Scientific).

ChIP-qPCR

Ubi1:OsTIF1-HA and Ubi1:MADS29-GFP were transformed into rice callus by particle bombardment using a PDS-1000/He biolistic particle delivery system (Bio-Rad) as described previously [53]. ChIP of OsTIF1-HA and MADS29-GFP was performed as previously described, with minor modifications [54]. Briefly, 1 g transgenic callus was ground to a fine powder with liquid nitrogen and combined with 10 mL ChIP extraction buffer (0.4 M sucrose, 10 mM Tris-HCl, 10 mM MgCl2, 1 mM DTT, 0.1 mM PMSF and protease inhibitor cocktail, pH 8.0). Then, 270 μL 37% formaldehyde solution (to a final concentration of 1%) was added and the mixture was incubated at 4°C for 10 min to cross-link DNA to protein. The cross-linking reaction was quenched by adding 0.63 mL 2 M glycine, and the mixture was incubated at 4°C for 5 min. Finally, the nuclear pellet was isolated by centrifugation and resuspended with ChIP lysis buffer (50 mM Tris-HCl, 10 mM EDTA, and 1% sodium dodecyl sulfate [SDS], pH 8.0) and kept on ice for 30 min. One volume of ChIP dilution buffer (16.7 mM Tris-HCl, 167 mM NaCl and 1.1% Triton X-100, pH 8.0) was added to the samples before sonication for 12 min (30 s on, 30 s off, high level) in a Bioruptor (Diagenode) to yield DNA fragments of 0.2–1.0 kb in length. The lysates were diluted fivefold in ChIP dilution buffer to decrease the concentration of SDS to 0.1% and cleared by centrifugation (16,000 × g for 5 min at 4°C). After keeping 5% of the sample as an input control, the remaining supernatant was incubated with antibody-bound Dynabeads Protein A or G (Invitrogen) overnight at 4°C. Washing, elution, reverse cross-linking, and DNA purification were performed according to the methods of a previous study [55]. The antibodies used for ChIP were hemagglutinin (HA, Sigma) and green fluorescent protein (GFP, Sigma). DNA isolated by ChIP was used for qPCR analysis. The 100-bp 5’-region of OsTAA1 was used as a negative control for OsTIF1, and the 100-bp N-terminal region of GUS was used as a negative control for MADS29. The primers used for qPCR are listed in S4 Table.

Transient transactivation assay

The promoters of OsTAA1 or OsAsp1 were amplified by PCR from rice genomic DNA and cloned into the pGPTV vector to generate the POsTAA1:GUS and POsAsp1:GUS constructs, respectively. POsTAA1:GUS was co-transformed with either 35S:Luc and 35S:OsAsp1-GFP or 35S:OsTIF1-HA into rice mesophyll protoplasts [56, 57]. POsAsp1:GUS was co-transformed with 35S:Luc and 35S:MADS29-GFP into Arabidopsis mesophyll protoplasts [58]. 35S:Luc was used as an internal control. GUS and luciferase (Luc) activities were measured as previously described [59], and relative GUS activity (GUS/Luc) was calculated to determine the activity of POsTAA1 or POsAsp1.

Yeast two-hybrid assay

The CDS of OsAsp1 or OsTIF1 was subcloned into pGAD7 or pGBK7 (Clontech) as prey and bait constructs, respectively. Then, the pairs of bait and prey constructs were co-transformed into yeast strain AH109 and screened by growth on selective medium (SD–Leu–Trp–His). The same strains were grown on a control medium (SD–Leu–Trp) for use as controls.

Co-IP assays

Ubi1:OsTIF1-HA was co-transformed with Ubi1:OsAsp1-GFP or Ubi1:GFP into rice mesophyll protoplasts, which then were incubated with or without 50 μM IAA for 16 h. Total protein was extracted and immunoprecipitated with GPF or HA antibody incubated with antibody-bound Dynabeads Protein A or G. The immunoprecipitates were separated by 10% SDS–polyacrylamide gel electrophoresis (PAGE) and detected with GFP antibody. Ubi1:HA and Ubi1:OsAsp1-GFP transgenic protoplasts were used as a negative control.

IAA and 2,4-D treatments

Spikelets were collected from panicles at -1 DAF and cultured in Murashige and Skoog (MS) medium plus IAA at different concentrations (100, 1.0, or 0.1 μM) or 2,4-dichlorophenoxyacetic acid (2,4-D, Sigma) at different concentrations (10, 1.0, or 0.1 μM) in 0.1 v/v Triton X-100. Some of the collected spikelets were cultured in the MS medium with methanol in 0.1 v/v Triton X-100 as a control, as previously described [34, 60]. Spikelets located on CSBs in the WT and on CPBs in asp1-1+/- were separately treated with 100 μM IAA (Sigma) in 0.1 v/v Triton X-100 using a fine brush once per day from –3 to 0 DAF. Methanol treatment used as a negative control (mock), following a method modified from a previous study [13].

Generation of transgenic Arabidopsis plants and GUS staining

The constructs 35S:OsAsp1-GFP and DR5:GUS [61] were separately introduced into Agrobacterium GV3101, and transformed into Arabidopsis WT (Col-0) using the floral dip transformation method [62]. The transformants were screened with 1/2 MS medium containing 50 mg/L hygromycin B (Amresco) and 50 mg/L kanamycin (INALCO) to obtain T3-generation homozygous transgenic lines. DR5:GUS transgenic lines were crossed with OsAsp1-overexpressed lines to generate DR5:GUS x OEOsAsp1-GFP plants for GUS staining The 3-day-old seedlings were stained in the GUS stain buffer (0.5 mg/mL X-Gluc, 0.01 v/v TritonX-100, 0.01 v/v DMSO, 10 mmol/L EDTA in sodium phosphate buffer, pH 7.0) for 2 h at 37°C, and then destained in ethanol solution for photographing with confocal laser-scanning microscopy (LSM 700, Zeiss).

Supporting information

S1 Table [docx]
Candidate genes were more highly expressed in apical primary branches (CPBs) at 0 days after heading (DAH) than in 5-DAH proximal secondary branches (CSBs), and expression was higher in 12-DAH CSBs than in 5-DAH CPBs.

S2 Table [docx]
Genotype of the F1 progeny of two independent lines self-crossed or back-crossed with wild-type (WT) plants.

S3 Table [docx]
Co-expression genes of from the Rice Oligo Array Database (ROAD).

S4 Table [docx]
Primers and probe sequences used in the study.

S1 Fig [jpg]
Development status of apical primary branches (CPBs) and proximal secondary branches (CSBs) at different days after heading (DAH).

S2 Fig [sd]
Transcription levels of (), and in CPBs and CSBs detected by quantitative reverse-transcription polymerase chain reaction (qRT-PCR) at different DAH.

S3 Fig [wt]
Growth phenotypes of the rice T-DNA insertion mutant lines and .

S4 Fig [jpg]
Mature grain phenotypes of CPBs and CSBs in the WT and mutant lines and .

S5 Fig [jpg]
Ovary length and width of CPBs and CSBs in the WT and at 3 and 5 days after flowering (DAF).

S6 Fig [jpg]
Potential regulation factors of predicted by the co-expression network.

S7 Fig [emsa]
Constructs for the GUS assay of promoter activity of and predicted OsTIF1 binding sites located in the promoter region of .

S8 Fig [tif]
Transcription levels of OsAsp1 detected by semi-quantitate RT-PCR in spikelets at -1, 0 and 1 DAF, and in detached spikelets at -1 DAF, after culture in Murashige and Skoog (MS) medium plus 2,4-D for 24 h. OsActin1 was used as a control.

S9 Fig [jpg]
Mature grain phenotypes of CPBs and CSBs in the WT in response to IAA treatment. CSBs were treated with IAA at each time point from -3 to 0 DAF, or with methanol as a control (mock).

S10 Fig [mock]
Mature grain weight and expression of in CPBs and CSBs in one whole plant of following IAA treatment.

S11 Fig [jpg]
Transcription levels of of CPBs and CSBs in the WT and () A-2 and A-33 at 5 DAF were detected by qRT-PCR.

S12 Fig [b]
Mature grain phenotype of CPBs and CSBs in the WT and () lines.


Zdroje

1. Xing Y, Zhang Q. Genetic and molecular bases of rice yield. Annual Review of Plant Biology. 2010;61(1):421–42. doi: 10.1146/annurev-arplant-042809-112209 20192739

2. Yang J, Zhang J. Grain-filling problem in ‘super’ rice. Journal of Experimental Botany. 2010;61(1):1–5. doi: 10.1093/jxb/erp348 19959608

3. Song XJ, Kuroha T, Ayano M, Furuta T, Nagai K, Komeda N, et al. Rare allele of a previously unidentified histone H4 acetyltransferase enhances grain weight, yield, and plant biomass in rice. Proceedings of the National Academy of Sciences. 2015;112(1):76–81. doi: 10.1073/pnas.1421127112 25535376

4. Song X-J, Huang W, Shi M, Zhu M-Z, Lin H-X. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nature Genetics. 2007;39:623. doi: 10.1038/ng2014 17417637

5. Shomura A, Izawa T, Ebana K, Ebitani T, Kanegae H, Konishi S, et al. Deletion in a gene associated with grain size increased yields during rice domestication. Nature Genetics. 2008;40:1023. doi: 10.1038/ng.169 18604208

6. Weng J, Gu S, Wan X, Gao H, Guo T, Su N, et al. Isolation and initial characterization of GW5, a major QTL associated with rice grain width and weight. Cell Research. 2008;18:1199. doi: 10.1038/cr.2008.307 19015668

7. Miao J, Yang Z, Zhang D, Wang Y, Xu M, Zhou L, et al. Mutation of RGG2, which encodes a type B heterotrimeric G protein γ subunit, increases grain size and yield production in rice. Plant Biotechnology Journal. 2018. doi: 10.1111/pbi.13005 30160362

8. Sun S, Wang L, Mao H, Shao L, Li X, Xiao J, et al. A G-protein pathway determines grain size in rice. Nature Communications. 2018;9(1):851. doi: 10.1038/s41467-018-03141-y 29487318

9. Xu R, Duan P, Yu H, Zhou Z, Zhang B, Wang R, et al. Control of grain size and weight by the OsMKKK10-OsMKK4-OsMAPK6 signaling pathway in rice. Molecular Plant. 2018;11(6):860–73. doi: 10.1016/j.molp.2018.04.004 29702261

10. Li Y, Fan C, Xing Y, Jiang Y, Luo L, Sun L, et al. Natural variation in GS5 plays an important role in regulating grain size and yield in rice. Nature Genetics. 2011;43:1266. doi: 10.1038/ng.977 22019783

11. Qi P, Lin Y-S, Song X-J, Shen J-B, Huang W, Shan J-X, et al. The novel quantitative trait locus GL3.1 controls rice grain size and yield by regulating Cyclin-T1;3. Cell Research. 2012;22:1666. doi: 10.1038/cr.2012.151 23147796

12. Zhang X, Wang J, Huang J, Lan H, Wang C, Yin C, et al. Rare allele of OsPPKL1 associated with grain length causes extra-large grain and a significant yield increase in rice. Proceedings of the National Academy of Sciences. 2012;109(52):21534–9. doi: 10.1073/pnas.1219776110 23236132

13. Ishimaru K, Hirotsu N, Madoka Y, Murakami N, Hara N, Onodera H, et al. Loss of function of the IAA-glucose hydrolase gene TGW6 enhances rice grain weight and increases yield. Nature Genetics. 2013;45(6):707–11. doi: 10.1038/ng.2612 23583977

14. Liu L, Tong H, Xiao Y, Che R, Xu F, Hu B, et al. Activation of Big Grain1 significantly improves grain size by regulating auxin transport in rice. Proceedings of the National Academy of Sciences. 2015;112(35):11102–7. doi: 10.1073/pnas.1512748112 26283354

15. Kato T, Shinmura D, Taniguchi A. Activities of enzymes for sucrose–starch conversion in developing endosperm of rice and their association with grain filling in extra-heavy panicle types. Plant Production Science. 2007;10(4):442–50. doi: 10.1626/pps.10.442

16. Xu G, Zhang J, Lam HM, Wang Z, Yang J. Hormonal changes are related to the poor grain filling in the inferior spikelets of rice cultivated under non-flooded and mulched condition. Field Crops Research. 2007;101(1):53–61.

17. Yang J, Zhang J, Wang Z, Liu K, Wang P. Post-anthesis development of inferior and superior spikelets in rice in relation to abscisic acid and ethylene. Journal of Experimental Botany. 2006;57(1):149–60. doi: 10.1093/jxb/erj018 16330527

18. Zhang H, Tan G, Yang L, Yang J, Zhang J, Zhao B. Hormones in the grains and roots in relation to post-anthesis development of inferior and superior spikelets in japonica/indica hybrid rice. Plant Physiology & Biochemistry. 2009;47(3):195–204. doi: 10.1016/j.plaphy.2008.11.012 19117763

19. Chen T, Xu G, Wang Z, Zhang H, Yang J, Zhang J. Expiression of proteins in superior and inferior spikelets of rice during grain filling under different irrigation regimes. Proteomics. 2016;16(1):102–21. doi: 10.1002/pmic.201500070 26442785

20. Zhu G, Ye N, Yang J, Peng X, Zhang J. Regulation of expression of starch synthesis genes by ethylene and ABA in relation to the development of rice inferior and superior spikelets. Journal of Experimental Botany. 2011;62(11):3907–16. doi: 10.1093/jxb/err088 21459770

21. Peng T, Lv Q, Zhang J, Li J, Du Y, Zhao Q. Differential expression of the microRNAs in superior and inferior spikelets in rice (Oryza sativa). Journal of Experimental Botany. 2011;62(14):4943–54. doi: 10.1093/jxb/err205 21791435

22. Peng T, Sun H, Qiao M, Zhao Y, Du Y, Zhang J, et al. Differentially expressed microRNA cohorts in seed development may contribute to poor grain filling of inferior spikelets in rice. BMC Plant Biology. 2014;14(1):1–17. doi: 10.1186/s12870-014-0196-4 25052585

23. Ishimaru T, Matsuda T, Ohsugi R, Yamagishi T. Morphological development of rice caryopses located at the different positions in a panicle from early to middle stage of grain filling. Functional Plant Biology. 2003;30(11):1139–49. doi: 10.1071/FP03122 32689096

24. Bi X, Khush GS, Bennett J. The rice nucellin gene ortholog OsAsp1 encodes an active aspartic protease without a plant-specific insert and is strongly expressed in early embryo. Plant and Cell Physiology. 2005;46(1):87–98. doi: 10.1093/pcp/pci002 15659452

25. Guo J, Wang F, Song J, Sun W, Zhang X-S. The expression of Orysa;CycB1;1 is essential for endosperm formation and causes embryo enlargement in rice. Planta. 2010;231(2):293–303. doi: 10.1007/s00425-009-1051-y 19921249

26. Barrôco RM, Peres A, Droual A-M, De Veylder L, Nguyen LSL, De Wolf J, et al. The cyclin-dependent kinase inhibitor Orysa;KRP1 plays an important role in seed development of rice. Plant Physiology. 2006;142(3):1053–64. doi: 10.1104/pp.106.087056 17012406

27. Su’Udi M, Cha JY, Min HJ, Ermawati N, Han CD, Min GK, et al. Potential role of the rice OsCCS52A gene in endoreduplication. Planta. 2012;235(2):387–97. doi: 10.1007/s00425-011-1515-8 21927949

28. Figueiredo DD, Batista RA, Roszak PJ, Köhler C. Auxin production couples endosperm development to fertilization. Nature Plants. 2015;1(12):1–6. doi: 10.1038/NPLANTS.2015.184 27251719

29. Yoshikawa T, Ito M, Sumikura T, Nakayama A, Nishimura T, Kitano H, et al. The rice FISH BONE gene encodes a tryptophan aminotransferase, which affects pleiotropic auxin-related processes. Plant Journal. 2014;78(6):927–36. doi: 10.1111/tpj.12517 24654985

30. Yamamoto Y, Kamiya N, Morinaka Y, Matsuoka M, Sazuka T. Auxin biosynthesis by the YUCCA genes in rice. Plant Physiology. 2007;143(3):1362–71. doi: 10.1104/pp.106.091561 17220367

31. Aoki Y, Okamura Y, Tadaka S, Kinoshita K, Obayashi T. ATTED-II in 2016: A plant coexpression database towards lineage-specific coexpression. Plant and Cell Physiology. 2016;57(1):e5. doi: 10.1093/pcp/pcv165 26546318

32. Franco-Zorrilla JM, López-Vidriero I, Carrasco JL, Godoy M, Vera P, Solano R. DNA-binding specificities of plant transcription factors and their potential to define target genes. Proceedings of the National Academy of Sciences. 2014;111(6):2367–72. doi: 10.1073/pnas.1316278111 24477691

33. Cao P, Jung K-H, Choi D, Hwang D, Zhu J, Ronald PC. The Rice Oligonucleotide Array Database: An atlas of rice gene expression. Rice. 2012;5(1):17. doi: 10.1186/1939-8433-5-17 24279809

34. Yin L-L, Xue H-W. The MADS29 transcription factor regulates the degradation of the nucellus and the nucellar projection during rice seed development. Plant Cell. 2012;24(3):1049–65. doi: 10.1105/tpc.111.094854 22408076

35. Mathan J, Bhattacharya J, Ranjan A. Enhancing crop yield by optimizing plant developmental features. Development. 2016;143(18):3283–94. doi: 10.1242/dev.134072 27624833

36. Wang B, Smith SM, Li J. Genetic regulation of shoot architecture. Annual Review of Plant Biology. 2018;69(1):437–68. doi: 10.1146/annurev-arplant-042817-040422 29553800

37. Davies DR. The structure and function of the aspartic proteinases. Annual Review of Biophysics and Biophysical Chemistry. 1990;19(1):189–215. doi: 10.1146/annurev.bb.19.060190.001201 2194475

38. Rawlings ND, Barrett AJ. Families of aspartic peptidases, and those of unknown catalytic mechanism. In Methods in Enzymology: Academic Press; 1995. p. 105–20.

39. Chen J, Ouyang Y, Wang L, Xie W, Zhang Q. Aspartic proteases gene family in rice: Gene structure and expression, predicted protein features and phylogenetic relation. Gene. 2009;442(1–2):108–18. doi: 10.1016/j.gene.2009.04.021 19409457

40. Soares A, Ribeiro Carlton SM, Simões I. Atypical and nucellin-like aspartic proteases: Emerging players in plant developmental processes and stress responses. Journal of Experimental Botany. 2019;70(7):2059–76. doi: 10.1093/jxb/erz034 30715463

41. Niu N, Liang W, Yang X, Jin W, Wilson ZA, Hu J, et al. EAT1 promotes tapetal cell death by regulating aspartic proteases during male reproductive development in rice. Nature Communications. 2013;4:1445. http://www.nature.com/ncomms/journal/v4/n2/suppinfo/ncomms2396_S1.html. doi: 10.1038/ncomms2396 23385589

42. Woodward AW, Bartel B. Auxin: Regulation, action, and interaction. Annals of Botany. 2005;95(5):707–35. doi: 10.1093/aob/mci083 15749753

43. Leyser O. Dynamic integration of auxin transport and signalling. Current Biology. 2006;16(11):R424–R33. doi: 10.1016/j.cub.2006.05.014 16753558

44. Casanova-Sáez R, Voß U. Auxin metabolism controls developmental decisions in land plants. Trends in Plant Science. 2019;24(8):741–54. doi: 10.1016/j.tplants.2019.05.006 31230894

45. Gallavotti A, Barazesh S, Malcomber S, Hall D, Jackson D, Schmidt RJ, et al. sparse inflorescence1 encodes a monocot-specific YUCCA-like gene required for vegetative and reproductive development in maize. Proceedings of the National Academy of Sciences. 2008;105(39):15196. doi: 10.1073/pnas.0805596105 18799737

46. Mashiguchi K, Tanaka K, Sakai T, Sugawara S, Kawaide H, Natsume M, et al. The main auxin biosynthesis pathway in Arabidopsis. Proceedings of the National Academy of Sciences. 2011;108(45):18512. doi: 10.1073/pnas.1108434108 22025724

47. Phillips KA, Skirpan AL, Liu X, Christensen A, Slewinski TL, Hudson C, et al. vanishing tassel2 encodes a grass-specific tryptophan aminotransferase required for vegetative and reproductive development in maize. The Plant Cell. 2011;23(2):550–66. Epub 2011/02/18. doi: 10.1105/tpc.110.075267 21335375

48. Abu-Zaitoon YM, Bennett K, Normanly J, Nonhebel HM. A large increase in IAA during development of rice grains correlates with the expression of tryptophan aminotransferase OsTAR1 and a grain-specific YUCCA. Physiologia Plantarum. 2012;146(4):487–99. doi: 10.1111/j.1399-3054.2012.01649.x 22582989

49. Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, et al. The TIGR Rice Genome Annotation Resource: Improvements and new features. Nucleic Acids Research. 2007;35(Suppl. 1):D883–D7. doi: 10.1093/nar/gkl976 17145706

50. Zhao Z, Zhang Y, Liu X, Zhang X, Liu S, Xiaowen YU, et al. A role for a dioxygenase in auxin metabolism and reproductive development in rice. Developmental Cell. 2013;27(1):113–22. doi: 10.1016/j.devcel.2013.09.005 24094741

51. Xu Y, Jin W, Li N, Zhang W, Liu C, Li C, et al. UBIQUITIN-SPECIFIC PROTEASE14 interacts with ULTRAVIOLET-B INSENSITIVE4 to regulate endoreduplication and cell and organ growth in Arabidopsis. Plant Cell. 2016;28(5):1200–14. doi: 10.1105/tpc.16.00007 27099260

52. Pan X, Welti R, Wang X. Quantitative analysis of major plant hormones in crude plant extracts by high-performance liquid chromatography–mass spectrometry. Nature Protocols. 2010;5(6):986–92. doi: 10.1038/nprot.2010.37 20448544

53. Sivamani E, DeLong RK, Qu R. Protamine-mediated DNA coating remarkably improves bombardment transformation efficiency in plant cells. Plant Cell Reports. 2009;28(2):213–21. doi: 10.1007/s00299-008-0636-4 19015859

54. Bowler C, Benvenuto G, PL, Molino D, Probst AV MT, et al. Chromatin techniques for plant cells. Plant Journal. 2004;39(5):776–89. doi: 10.1111/j.1365-313X.2004.02169.x 15315638

55. Cui X, Lu F, Qiu Q, Zhou B, Gu L, Zhang S, et al. REF6 recognizes a specific DNA sequence to demethylate H3K27me3 and regulate organ boundary formation in Arabidopsis. Nature Genetics. 2016;48(6):694. doi: 10.1038/ng.3556 27111035

56. Yoo S-D, Cho Y-H, Sheen J. Arabidopsis mesophyll protoplasts: A versatile cell system for transient gene expression analysis. Nature Protocols. 2007;2(7):1565–72. doi: 10.1038/nprot.2007.199 17585298

57. Wu F-H, Shen S-C, Lee L-Y, Lee S-H, Chan M-T, Lin C-S. Tape-Arabidopsis Sandwich—a simpler Arabidopsis protoplast isolation method. Plant Methods. 2009;5(1):1–10. doi: 10.1186/1746-4811-5-16 19930690

58. Bart R, Chern M, Park C-J, Bartley L, Ronald PC. A novel system for gene silencing using siRNAs in rice leaf and stem-derived protoplasts. Plant Methods. 2006;2(1):1–9. doi: 10.1186/1746-4811-2-13 16808845

59. Ding Y, Li H, Zhang X, Xie Q, Gong Z, Yang S. OST1 kinase modulates freezing tolerance by enhancing ICE1 stability in Arabidopsis. Developmental Cell. 2015;32(3):278–89. doi: 10.1016/j.devcel.2014.12.023 25669882

60. Uchiumi T, Okamoto T. Rice fruit development is associated with an increased IAA content in pollinated ovaries. Planta. 2010;232(3):579–92. doi: 10.1007/s00425-010-1197-7 20512651

61. Ulmasov T, Murfett J, Hagen G, Guilfoyle TJ. Aux/IAA proteins repress expression of reporter genes containing natural and highly active synthetic auxin response elements. Plant Cell. 1997;9(11):1963–71. doi: 10.1105/tpc.9.11.1963 9401121

62. Clough SJ, Bent AF. Floral dip: A simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant Journal. 1998;16(6):735–43. doi: 10.1046/j.1365-313x.1998.00343.x 10069079


Článek vyšel v časopise

PLOS Genetics


2020 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#