- freely available
Plants 2014, 3(1), 95-112; doi:10.3390/plants3010095
Published: 7 February 2014
Abstract: Indole-3-acetic acid (IAA) synthesis is required for grain-fill in maize and appears to be regulated by cell-wall invertase (CWIN) activity. OsYUC12 is one of three IAA biosynthesis genes we previously reported as expressed during early rice grain development, correlating with a large increase in IAA content of the grain. This work aimed to investigate further the role of OsYUC12 and its relationship to CWIN activity and invertase inhibitors (INVINH). The analysis shows a brief peak of OsYUC12 expression early in endosperm development. Meta-analysis of microarray data, confirmed by quantitative expression analysis, revealed that OsYUC12 is coexpressed with OsIAA29, which encodes an unusual AUX/IAA transcription factor previously reported as poorly expressed. Maximum expression of OsYUC12 and OsIAA29 coincided with maximum CWIN activity, but also with a peak in INVINH expression. Unlike ZmYUC1, OsYUC12 expression is not reduced in the rice CWIN mutant, gif1. Several reports have investigated CWIN expression in rice grains but none has reported on expression of INVINH in this species. We show that rice has 54 genes encoding putative invertase/pectin methylesterase inhibitors, seven of which are expressed exclusively during grain development. Our results suggest a more complex relationship between IAA, CWIN, and INVINH than previously proposed.
Several recent publications indicate key intersecting signaling roles for indole-3-acetic acid (IAA) and cell wall invertases (CWIN) during cereal grain development. The maize defective endosperm18 (de18) phenotype was recently shown to result from loss of expression of the IAA biosynthesis gene ZmYUC1 and related low levels of IAA in the developing grain . The cell wall invertase miniature1 (mn1) mutant also shows poor grain fill, the low levels of cell wall invertase resulting in a defective basal endosperm transfer layer (BETL) with poorly developed wall in-growths . Maize mn1 mutants have low levels of IAA and low expression of ZmYUC1; glucose was able to increase ZmYUC1 transcript levels in cultured kernels . Forestan et al.  showed that IAA accumulates in the BETL, aleurone and embryo surrounding region (ESR) just before the endosperm starts to accumulate starch. In addition, both BETL and ESR showed a high level of auxin transporter ZmPIN1 transcript and protein. The rice ortholog of Mn1, GIF1/OsCIN2, also appears to be important for grain development, with gif1 mutants (that have much lower CWIN activity in developing grains) showing poor grain-fill; expression of GIF1 during grain development is localized to the vascular trace . The relationship between IAA and invertase in rice has not been investigated.
In a previous paper on IAA synthesis in developing rice grains we showed that IAA accumulates more than 50-fold in developing kernels during endosperm cellularisation and early starch deposition . The rise in IAA content was correlated with a large increase in expression of IAA biosynthesis genes OsYUC9, OsYUC11, and OsTAR1. We also obtained some rather inconsistent data on a third endosperm-specific YUCCA gene, OsYUC12. Preliminary results suggested that OsYUC12 was expressed only briefly during grain development and may have a distinct role from the more highly expressed OsYUC9 and OsYUC11. The aim of the work reported here was to clarify the expression of OsYUC12 in developing grains and to investigate its relationship to CWIN, as well as the unexpected observation that OsYUC12 is coexpressed with a putative invertase inhibitor (INVINH).
A number of reports, e.g., Jin et al. , indicate that invertases in dicots are regulated by INVINH. In addition, ZM-INVINH1 inhibited maize CWIN activity in vitro, bound to a glycoprotein fraction including CWIN, and was localized to the ESR of young developing maize kernels . Several papers report on expression of CWIN in rice, e.g., , however the role of invertase inhibitors in rice has been neglected. Furthermore, extraction of invertase for assay of its activity has normally been carried out at pH 7.5, e.g., , conditions under which any bound inhibitor would be expected to dissociate from the enzyme . In this work, we carried out a comprehensive phylogenetic analysis of invertase inhibitor homologs in rice; we investigated the expression of Os04g49720, a co-ortholog of ZM-INVINH1 previously designated as OsINVINH3 ; we extracted and assayed CWIN activity under conditions that have been shown to preserve an enzyme/inhibitor complex . Our data on expression of OsYUC12 and OsINVINH3, as well as CWIN activity suggest that the relationship between IAA and CWIN is more complex than previously proposed.
In addition, preliminary work suggested a second auxin related gene, OsIAA29, encoding a putative AUX/IAA transcription factor is coexpressed with OsYUC12 and OsINVINH3. Coexpression analysis has been used productively to predict networks of genes that may act together in a common pathway, e.g., . To explore other genes that may act in the same signaling pathway as OsYUC12, we carried out a comprehensive coexpression meta-analysis using OsYUC12, OsIAA29, and OsINVINH3 as bait genes. The search for coexpressed genes revealed a number of sequences encoding uncharacterized proteins, some of which are homologs of proteins with an identified signaling role earlier in plant reproductive development. In addition, a number of the coexpressed genes encode homologs of molecules specifically found in endosperm transfer cells of other cereals. We suggest that the coexpressed genes may encode novel signaling factors.
2.1. OsYUC12 and Other Endosperm YUCCA Genes Are Conserved in Cereal Species
In preparation for investigating further the role of OsYUC12 in grain-fill we asked the question whether this gene had conserved orthologs in other cereal species. A phylogenetic analysis of OsYUC9, OsYUC11, and OsYUC12 in maize, sorghum, and rice (Supplemental Figure S1) showed that each of these YUCCA genes has clear orthologs in each species. Secondly, we explored the expression of GRMZM2G107761_T01, the maize ortholog of OsYUC12, using microarray data from Sekhon et al.  accessed via PLEXdb  as experiment ZM37. This revealed that GRMZM2G107761_T01 was also expressed primarily during grain development with a peak at eight days after pollination (DAP) (Supplemental Figure S2).
2.2. Expression of OsYUCCA Genes in the gif1 CWIN Mutant
In light of the low expression and activity of ZmYUC1 in the maize CWIN mutant mn1, we investigated expression of OsYUC9, OsYUC11, and OsYUC12 in the equivalent gif1 mutant of rice using data from Wang et al.  “Global gene expression profiles of Oryza sativa wild type Zhonghua11 and mutant gif1 in filling stage” accessible as accession OS17 in PLEXdb . The data shown in Figure 1 indicates that expression of OsYUC9 only is reduced in the gif1 mutant.
2.3. Invertase Inhibitor Phylogeny
A BlastP search of the rice proteome using query sequences ZM-INVINH1, ZM-INVINH2, and ZM-INVINH3 initially identified 20 homologous proteins. Additional rice homologs, as well as sequences from sorghum, were obtained using the homologous protein function in Phytozome. The phylogenetic tree shown in Figure 2 resulted from analysis of 54 protein sequences from rice, 34 sequences from sorghum, as well as experimentally characterized pectin methylesterase inhibitors (PMEIs) from wheat  and Arabidopsis , INVINHs from tomato , potato , and Arabidopsis  and F2DCS9 HORVD, an invertase inhibitor homolog from barley identified as up-regulated under heat stress by Mangelsen et al. . The comprehensive set of homologs from the sorghum proteome was included to identify rice sequences that were conserved in more than one cereal species.
A number of bootstrap values for larger clades are low due to the highly diverse nature of the sequences. However, an examination of the multiple sequence alignment indicated that the four conserved cysteine residues previously shown to form disulfide bridges in PMEIs and INVINHs  were present in most sequences and were correctly aligned (Supplemental Figure S3). Branches with high bootstrap values indicate reliable prediction of phylogenetic relationships. These show that sorghum has unambiguous orthologs of ZM-INVINH1, ZM-INVINH2, and ZM-INVINH3. Rice also has clear orthologs of both ZM-INVINH1 (Os04g49730) and ZM-INVINH3 (Os01g20970). Further analysis of sequences most similar to ZM-INVINH1 indicated that the tandem repeats Os04g49730 and Os04g49720, as well as Os02g46360 are co-orthologous to ZM-INVINH1, the only characterized invertase inhibitor from a cereal . Interestingly, characterized invertase inhibitors from dicot species are in separate clades from ZM-INVINH1. Arabidopsis invertase inhibitor, C/VIF1, as well as potato and tomato invertase inhibitors do not have any unambiguous cereal orthologs. On the other hand, C/VIF2 falls into a well-supported clade containing ZM-INVINH2, two rice sequences and four sorghum sequences. ZM-INVINH3 and Os01g20970 are orthologs of TaPMEI and are, therefore, likely to be PMEIs rather than invertase inhibitors. The Arabidopsis PMEIs appear in a separate branch from cereal PMEIs. The tree contains many orthologous pairs from rice and sorghum; however rice has many more sequences which appear to be due to recent gene duplication as evidenced by groups of tandem repeats, particularly on chromosome 8. The barley protein F2DCS9_HORVD, which has been given the name INVINH1 , is only distantly related to other characterized invertase inhibitors but has unambiguous orthologs in both rice and sorghum.
2.4. Analysis of OsCIN2/GIF1 and Putative INVINH Expression in Developing Caryopses Using Published Microarray Data
The extracellular invertase GIF1/OsCIN2 has been shown to be responsible for most of the CWIN activity in developing rice grains . We mined microarray data for information on expression of GIF1 as well as the 54 invertase inhibitor homologs during grain development. Data was accessed via PLEXdb  from experiment accessions OS5 “Expression data for reproductive development in rice” , OS8 “Expression data from rice embryo, endosperm, root, leaf, and seedling” , OS16 “Genome-wide gene expression profiling of rice stigma” , OS44 “Rice expression atlas (3): Early embryogenesis”  and OS89 “Expression data from rice embryo and endosperm development” .
Seven INVINH/PMEI homologs (Table 1) were found to be expressed exclusively during caryopsis development with expression highest in samples three to four days after pollination (DAP); expression appeared specific to the endosperm. Caryopsis-expressed INVINH/PMEI homologs included Os04g49720 which is co-orthologous with Os04g49730 and Os02g46360 to ZM-INVINH1. Os04g49730 and Os02g46360 appeared to be only weakly expressed; Os04g49720 has previously been designated as OsINVINH3 . In contrast to the large changes in expression of OsINVINH3, there was little change in the expression of GIF1 from prior to anthesis until 20 DAP (Figure 3). Interestingly examination of endosperm-specific genes in Figure 2 showed them to be present in diverse clades. We also noted that the expression profiles of genes in Table 1 matched that of IAA biosynthesis gene OsYUC12 and the gene encoding putative AUX/IAA transcriptional regulator OsIAA29.
|Table 1. Invertase inhibitor/pectin methylesterase inhibitor (INVINH/PMEI) homologs from rice identified as expressed during caryopsis development using microarray data accessed via PLEXdb  from experiment accessions OS5 “Expression data for reproductive development in rice” , OS8 “Expression data from rice embryo, endosperm, root, leaf, and seedling” , OS16 “Genome-wide gene expression profiling of rice stigma” , OS44 “Rice expression atlas (3): Early embryogenesis” , and OS89 “Expression data from rice embryo and endosperm development” .|
|Os02g46290.1||2 to 4 DAP, endosperm|
|Os03g43820.1||3-4 DAP, 5-10 DAP, endosperm|
|Os04g49720.1||3-4 DAP, 5-10 DAP, endosperm|
|Os08g04670.1||3-4 DAP, endosperm|
|Os08g04710.1||3-4 DAP, endosperm|
|Os08g04740.1||3-4 DAP, 5-10 DAP, endosperm|
|Os08g42890.1||3-4 DAP, 5-10 DAP, endosperm|
2.5. Coexpression Meta-Analysis Using Bait Genes OsIAA29, OsYUC12 and OsINVINH3
Several on-line platforms are available for coexpression analysis including RiceNet , GeneCAT , and Rice Oligonucleotide Array Database . Initial investigations indicated that the limited data set used in RiceNet led to the identification of some genes that had poorly correlated expression with our genes of interest in data sets, such as OS44 , which includes samples collected daily during the early stages after pollination. Rice Oligonucleotide Array Database contains more data sets but as most of these did not contain caryopsis tissue and the site does not allow for condition dependent selection of data sets, results were similarly unsatisfactory. We therefore downloaded RMA normalized datasets via PLEXdb for previously mentioned accessions OS5, OS8, OS16, OS44, and OS89 and calculated pairwise mutual coexpression ranks  for three bait genes: OsIAA29, OsYUC12, and OsINVINH3. Examination of the expression profiles of genes with a mutual co-expression rank <50 in individual experiments indicated very close correspondence with expression of bait genes. All appear to be expressed exclusively in the endosperm of developing grains with maximal expression at 3–4 DAP.
Genes with a mutual co-expression rank less than 20 (MR < 20) are shown in Table 2, with a larger set (MR < 50) shown in Supplemental Table S1. All database annotations have been manually checked against the literature; where the annotation differs from the literature, the published protein designation has been used. Table 2 includes two lipid transfer protein-like genes (LTPL) and four defensin/defensin like genes (DEF/DEFL). This group of genes encoding cysteine-rich peptides (CRPs) was also highly represented in the larger MR < 50 group, with a total of 13 genes present. Also present was a ZmEBE-like (embryo sac/basal endosperm transfer layer/embryo surrounding region) gene; a further four homologs of this gene are found in the MR <50 group. A third group of homologous proteins found in the MR <50 group included two additional INVINH/PMEI homologs. Table 2 also includes genes encoding OsHAP3D/OsLEC1A (Leafy cotyledon1)  and a leucine-rich repeat receptor-like kinase. Genes with a likely signaling role in the MR <50 group also include OsRR33, a type-B response regulator  an ortholog of OBERON1/2  and at least two other proteins with a predicted nuclear location.
|Table 2. Genes with a pairwise mutual coexpression rank < 20 with each of the bait genes, INVINH3, OsIAA29, and OsYUC12. Genes listed in order of increasing mutual co-expression rank.|
|Os12g13960.1||Lipid transfer protein-like LTPL33 |
|Os02g07628.1||Defensin DEF5 |
|Os11g45360.1||Defensin-like DEFL15 |
|Os01g70680.1||Defensin DEF1 |
|Os02g07624.1||Defensin DEF4 |
|Os11g18140.1||homologous to ZmEBE 1 (embryo sac/basal endosperm transfer layer/embryo surrounding region |
|Os06g17480.1||Nuclear factor YC OsHAP3D /OsLEC1A |
|Os05g07850.1||Leucine-rich repeat receptor-like kinase |
|Os11g14880.1||Lipid transfer protein-like LTPL32 |
2.6. Quantitative Expression Analysis of Genes Coexpressed with OsIAA29, OsYUC12, and OsINVINH3
To confirm the expression profile of bait genes OsYUC12, OsIAA29, and INVINH3 along with coexpressed genes LTPL32, DEFL13, ZmEBE-like Os11g18140, and OsHAP3D in plants of the same variety and grown under the same conditions as our previously published IAA measurements, we carried out quantitative RT-PCR analysis. LTPL32, DEFL13, ZmEBE-like Os11g18140 were chosen as representatives from the LTPL, DEFL, and ZmEBE-like gene families, as these genes contained introns allowing for primers to be designed across exon boundaries. Expression data are presented in Figure 4. All genes showed a short-lived increase in expression during early grain development; under our plant growth conditions an increase in expression was seen at 4 DAP (apart from OsINVINH3) with highest expression apparent in 7 DAP samples. Large error bars for some 4 DAP samples resulted from one of the biological replicates in which the major increase in expression had occurred earlier than the other two samples. The early increase in expression was seen in most genes when using the same RNA template.
2.7. Invertase Activity in Developing Rice Grains
Cell wall invertase activity during grain development is shown in Figure 5. Total cell wall invertase shown on a per grain basis increased to a maximum at 7 DAP; when calculated on a per g fresh weight basis, no significant difference was seen between samples from 1–7 DAP with a substantial decrease at 10 DAP.
3.1. Expression of IAA Related Genes, OsYUC12 and OsIAA29 Coincides with Endosperm Cellularisation
Previous research using maize and rice has suggested a key role for IAA and CWIN during early grain-fill however the exact roles of IAA and CWIN remain to be established. Our results indicated that the relationship between IAA and CWIN is complex, with differences in expression of three IAA biosynthesis genes OsYUC9, OsYUC11, and OsYUC12 during grain development. In this paper we explore the role of one IAA biosynthesis gene, OsYUC12. We chose to investigate this gene as to our knowledge there have been no other studies of it or its orthologs. We have shown that OsYUC12 has conserved orthologs in maize and sorghum, it is specifically expressed in endosperm for a short time during early grain development and this expression pattern also appears in maize. We, therefore, suggest that this gene has a specific and conserved signaling role during grain development. The aim of this work was to identify other possible components of the signaling network that could be investigated.
An examination of on-line accessible microarray data indicated that AUX/IAA gene OsIAA29 has a very similar expression pattern to OsYUC12; this was confirmed by qRT-PCR. Our quantitative expression data confirm that expression of OsYUC12 and OsIAA29 is limited to a short period during early grain-fill. The increase in expression coincided with our previously reported large increase in the IAA content of developing grains at 7 DAP in plants grown under the same conditions . Examination of material under the light microscope indicated that, at 4 DAP, the endosperm syncytium contained many nuclei and a few small irregularly shaped cells. At day seven there were clearly defined endosperm cells containing amyloplasts, as well as differentiated outer cellular layers. The coexpression of OsYUC12 and OsIAA29 in the endosperm over a very limited period during grain development, suggests a specific role for IAA during this period. However, OsIAA29 is an unusual AUX/IAA protein that lacks domains I and II . As domain II is responsible for interacting with the F-box protein, IAA co-receptor  it is unclear how OsIAA29 may interact with IAA. Dimerization domains III and IV, which mediate interactions between AUX/IAA proteins and ARFs are still present in OsIAA29 . In work with seedling tissue, Jain et al.  reported very low expression of OsIAA29 and no expression of its tandem repeat OsIAA28. Our observations showing high expression at a very specific stage of development suggest that OsIAA29 may in fact play a key but as yet unidentified role. This is further supported by the existence of OsIAA29 orthologs in both sorghum and maize and our recent unpublished observation that the sorghum ortholog is also expressed during grain-fill.
3.2. Expression of OsINVINH3 Coincides with Maximum Extracellular Invertase Activity
LeClere et al.  have provided evidence that ZmYUC1 expression may be regulated by glucose levels produced as a result of CWIN activity. We therefore investigated the expression of GIF1 during grain development. Examination of microarray data indicated GIF1 expression appeared to change little from prior to pollination until 20 DAP. However, we were surprised to find that expression of Os04g49720 (OsINVINH3), a co-ortholog of ZM-INVINH1  changed dramatically during early grain-development. Invertase inhibitor ZM-INVINH1 has been reported to bind CWIN during grain development in maize ; both ZM-INVINH1 and OsINVINH3 have predicted signal peptides. Quantitative expression analysis confirmed that Os04g49720 showed a large, short-lived increase in expression during early grain development. In plants grown under our conditions maximum expression was found in 7-DAP samples, also coinciding with the peak in OsYUC12 and OsIAA29 expression and rapidly increasing IAA levels . We were surprised to also observe maximum invertase activity on a per grain basis at 7 DAP. The unexpected result was confirmed by repetition of the experiment by a second person in the laboratory. Clearly, further work is required to clarify the relationship between CWIN and OsINVINH3, however, our data support the suggestion by Wang et al.  that CWIN in developing grains may be tightly regulated.
3.3. Rice Has 54 INVINH Homologs Several of Which Are Expressed Specifically in Developing Endosperm
Invertase inhibitors have similar structure and are homologous to pectin methylesterase inhibitors . A comprehensive search for members of this family revealed 54 diverse homologous proteins in rice. Sorghum homologs were also included to distinguish between species-specific proteins and those that have been conserved. Although the high sequence diversity resulted in some ambiguous relationships, both sorghum and rice have clear orthologs of ZM-INVINH1. In addition, a small well-supported clade containing the experimentally characterized wheat PMEI suggested that ZM-INVINH3, Sb03g012790, Os01g20970 and Os05g29740 are likely to be PMEIs. OsINVINH3 is not the only member from the family to be expressed exclusively during grain development. Os02g46290, Os03g43820, Os08g42890, Os08g04670, Os08g04710, and Os08g04740 are also grain specific and expressed maximally at 3–4 DAP. Exploration of data from OS8 , OS16 , and OS89  indicated that expression is restricted to the endosperm. The seven endosperm-specific proteins are found in four clades; Os03g43820 is present in a large clade quite separate from the others, but including the barley protein F2DCS9 HORVD. Os08g04670, Os08g04710, and Os08g04740 are part of a cluster of nine genes found in tandem; Os08g42890 appears to be part of the same clade though this is not strongly supported. Os02g46290 is a possible ortholog of C/VIF1, an experimentally characterized vacuolar invertase inhibitor from Arabidopsis. The high diversity of these proteins in separate clades of a tree that is known to contain proteins of different biological functions raises the possibility that not all are invertase inhibitors or PMEIs. We suggest that some members of the INVINH/PMEI protein family may regulate other carbohydrate metabolizing enzymes. An investigation of the predicted sub-cellular locations of the gene products, using CELLO  and TargetP , suggests that all proteins have a signal peptide, however an examination of extracellular and vacuolar invertase inhibitors from potato  shows that domain prediction programs cannot reliably distinguish between these two sub-cellular destinations.
3.4. Putative Signaling Proteins Coexpressed with OsYUC12, OsIAA29 and OsINVINH3
The observation that OsYUC12, OsIAA29, and OsINVINH3 are expressed exclusively in endosperm samples from developing grains 3–10 DAP suggest they may be co-regulated and may be part of a common signaling pathway. A meta-analysis of microarray data indicated genes from three families, CRPs, ZmEBE-like, and INVINH/PMEI-like, all encoding putative extracellular proteins, are prominent among genes coexpressed with OsYUC12, OsIAA29, and OsINVINH1. The CRPs come from two distinct groups; the so-called lipid transfer protein-like group and the defensin/defensin-like group. We note that despite their reported antimicrobial activity, defensins have recently been shown to have a signaling role earlier in maize reproductive development . In addition, the CRP MEG1 is essential for the differentiation of ETCs in maize . OsLTPL29/OsPR602 and DEFL OsPR9a, as well as ZmEBE1, ZmEBE2, and ZM-INVINH1, are also expressed primarily in ETCs and ESR [38,47]. Differentiating transfer cells in maize kernels accumulate high levels of auxin transporter OsPIN1 transcripts ; the BETL and ESR also show high levels of PIN proteins as well as an accumulation of IAA itself. The specific localization of rice proteins found in Table 2 should be investigated to determine whether they are also restricted to the ETCs and/or ESR.
Genes with a mutual co-expression rank of <20 included also OsHAP3D  and Os05g07850, a putative leucine-rich repeat receptor-like kinase. OsHAP3D, also referred to as OsLEC1a, encodes a nuclear factor YC protein, homologous to Arabidopsis LEC1/L1L. These Arabidopsis proteins have a crucial role in embryogenesis  including the up-regulation of YUCCA10 . Rice has two homologs of LEC1, OsHAP3D, and OsHAP3E. An examination of their expression profile via PLEXdb indicates that both are expressed in developing grains with a peak at 3–4 DAP. However, whereas OsHAP3D is expressed in endosperm, OsHAP3E expression appears to be restricted to the embryo. Os05g07850 encodes a leucine-rich repeat receptor-like kinase which has a predicted trans-membrane domain and serine/threonine kinase domain. It is homologous but not orthologous to the brassinolide receptor BRI1 and aligns in the same clade as Arabidopsis NSP interacting kinases. A partially homologous receptor-like kinase has been shown to form part of a signaling network that also includes a DEFL during fertilization in maize . The identification of these coexpressed genes provides a framework for investigating protein-protein interactions that could form part of a signaling network.
4.1. Bioinformatic Analysis
Rice protein sequences were downloaded via Phytozome 8.0  from MSU Release 7.0 of the Rice Genome Annotation; protein sequences for ZM-INVINH1, ZM-INVINH2, ZM-INVINH3, and experimentally characterized pectin methyl esterase inhibitors (PMEI) and invertase inhibitors (INVINH) from wheat, potato, tomato, and Arabidopsis were downloaded from Uniprot. The phylogenetic analysis of protein sequences was carried out in MEGA5.2  using the Neighbor-Joining method . MUSCLE  was used for multiple sequence alignment. The bootstrap consensus tree was inferred from 500 replicates . Evolutionary distances were computed using the Poisson correction method . All ambiguous positions were removed for each sequence pair. Normalized microarray data from accessions OS5 , OS8 , OS16 , OS17 , OS44 , and OS89  was downloaded via PLEXdb  and analyzed in Excel. Pairwise mutual coexpression ranks were calculated as per Obayashi et al. .
4.2. Plant Material
Rice grains (ssp. Japonica cv. Jarrah) were direct sown into flooded soil-filled cylindrical plastic pots and thinned out to three plants per pot when the seedlings had reached the 2–3 leaf growth stage. Plants were grown in a glasshouse under a natural light with day/night temperatures 28 °C/18 °C. Plants were fertilized fortnightly using Aquasol® (8 g/5 L). When spikelets in the top half of the panicle had reached anthesis, panicles were tagged to record the date; panicles were harvested at 1, 4, 7, 10, and 15 days after pollination (DAP). Grain samples were placed in microfuge tubes, the weight and number of grains recorded before freezing in liquid nitrogen and storage at −80 °C.
4.3. Quantitative Reverse Transcriptase PCR
Total RNA was extracted from 60–80 mg grain samples following instructions for the Bioline ISOLATE Plant RNA Mini Kit. RNA concentration and 260/280 ratio were determined using a NanoDrop ND-1000 Spectrophotometer (NanoDrop Technologies, Inc., Wilmington, DE, USA). Only samples with A260/A280 of 1.8–2.0 were used for further analysis. The RNA quality was also checked using agarose gel electrophoresis for two clear bands of 18s and 28s rRNA .
Gene-specific primers (See Supplemental Table S2) were designed with melting points in the range 59–61 °C and product sizes between 100–150 bp. Primers were tested for amplification of a single product of the expected size using a QIAGEN OneStep RT-PCR kit (QIAGEN, Limburg, The Netherlands), with agarose gel analysis of products. To avoid the possibility of amplifying genomic DNA, either left or right primers were designed to span an exon-exon boundary. Controls (no reverse transcriptase) were included for INVINH1 and OsHAP3D genes, which had no introns; amplification was not detected in the controls. Quantitative RT-PCR was carried out following manufacturer’s instructions using 50 ng RNA per 20 μL reaction, a final primer concentration of 0.5 μM and reagents from a Bioline SensiFAST™SYBR No-ROX RT-PCR kit (Bioline, London, UK). Reaction tubes were prepared by using a QIAgility robotic liquid handling system (QIAGEN). Reactions were monitored in a Rotor-Gene Q (QIAGEN) thermocycler fitted with a 100-tube gene disc rotor. The amplification program included 48 °C min, 95 °C 2 min and 45 cycles of 95 °C 5 s, 58 °C 10 s, and 72 °C 5 s. Expression was calculated relative to reference gene for ubiquitin-conjugating enzyme E2 (OsUBC)  using ROTORGENE software (QIAGEN). Melt curve analysis confirmed that a single product was obtained for each gene amplified.
4.4. Invertase Assay
Rice grain samples were extracted using a procedure modified from Tomlinson et al. . Samples (200 mg) were ground to a fine powder in liquid N2 in a mortar and pestle then transferred to a close fitting glass homogenizer and ground further with 1 mL of extraction buffer (pH 4.8, 100 mM acetic acetate buffer, containing 8 mM MgCl2, 2 mM EDTA and 12.5% (w/v) glycerol, 10 mM dithiothreitol and 1 mM phenylmethanesulfonyl fluoride). The homogenate plus a 400 μL wash were transferred to a 2 mL microfuge tube and centrifuged at 16,000 ×g for 10 min. The pellet was washed and re-centrifuged twice with 500 μL of extraction buffer before resuspending in two volumes of salt extraction buffer (pH 4.8 100 mM acetate buffer, 1 M NaCl). Samples were extracted for 1 h on ice with gentle agitation, as well as 3× sonication for 60 s. Following salt extraction, samples were centrifuged at 16,000 ×g for 10 min and the supernatant retained for assay of CWIN activity.
Reaction mixtures containing 80 μL of enzyme preparation, 10 mM sucrose and pH 4.8, 100 mM acetate buffer to a total volume of 440 μL were incubated at 37 °C for 30 min with gentle shaking. Three 40 μL aliquots were taken for glucose determination at 0 and 30 min. Glucose in the reaction mixture was measured using the Somogyi-Nelson method ; and calculated in comparison with a set of standards containing 5 to 100 μg glucose.
In conclusion, we have shown that the endosperm-specific IAA synthesis gene OsYUC12 and the AUX/IAA gene OsIAA29 are coexpressed, with expression limited to a short period during endosperm cellularisation. Extractable CWIN activity, as well as expression of putative invertase inhibitor OsINVINH3, also reached a maximum at 7 DAP, coinciding with maximal expression of OsYUC12 and OsIAA29. The possible regulation of CWIN activity by OsINVINH3 needs to be considered as part of any investigation of the role of CWIN during grain development. We have shown that rice has at least 54 genes encoding diverse proteins homologous to invertase inhibitors and PMEIs. In addition to OsINVINH3, a rice ortholog of ZM-INVINH1, six other INVINH/PMEI genes from widely divergent clades were specifically expressed in endosperm during early grain-fill. We explored other genes co-expressed with OsYUC12, OsIAA29, and OsINVINH3. Seven out of the nine genes with MR<20 appear to encode extracellular proteins homologous to molecules expressed exclusively in ETCs/ESR of other plants; these include six cysteine-rich proteins, some with similarity to BETL1 from maize. In addition, the presence of an endosperm-specific LEC1 ortholog, OsHADP3D and a receptor–like kinase gene suggests that these genes and their encoded proteins are a priority for experimental investigation in regard to a signaling network regulating differentiation of ETCs in rice. The expression of OsYUC12 is distinct from OsYUC9 and OsYUC11. It has orthologs in sorghum and maize; the maize ortholog is also expressed for a short time during grain development, suggesting conserved subfunctionalization. Expression of OsYUC9 but not OsYUC12 is reduced in the gif1 CWIN mutant. We suggest that further investigation of the specific roles of three IAA biosynthesis genes will be necessary to reveal the full extent of the relationship between CWIN and IAA.
This work was supported by funding from the University of New England.
Sarah Russell French carried out quantitative gene expression analysis, Yousef Abu Zaitoon completed the coexpression analysis, Myn Uddin carried out the search for invertase inhibitor/PMEI sequences, invertase enzyme assays were done by Karina Bennett. Heather Nonhebel supervised all work, carried out data analysis and wrote the paper.
Conflicts of Interest
The authors declare no conflict of interest.
- Bernardi, J.; Lanubile, A.; Li, Q.-B.; Kumar, D.; Kladnik, A.; Cook, S.D.; Ross, J.J.; Marocco, A.; Chourey, P.S. Impaired auxin biosynthesis in the defective endosperm18 mutant is due to mutational loss of expression in the ZmYUC1 gene encoding endosperm-specific YUCCA1 protein in maize. Plant Physiol. 2012, 160, 1318–1328. [Google Scholar] [CrossRef]
- Kang, B.-H.; Xiong, Y.; Williams, D.S.; Pozueta-Romero, D.; Chourey, P.S. Miniature1-encoded cell wall invertase is essential for assembly and function of wall-in-growth in the maize endosperm transfer cell. Plant Physiol. 2009, 151, 1366–1376. [Google Scholar] [CrossRef]
- LeClere, S.; Schmelz, E.A.; Chourey, P.S. Sugar levels regulate tryptophan-dependent auxin biosynthesis in developing maize kernels. Plant Physiol. 2010, 153, 306–318. [Google Scholar] [CrossRef]
- Forestan, C.; Meda, S.; Varotto, S. ZmPIN1-mediated auxin transport is related to cellular differentiation during maize embryogenesis and endosperm development. Plant Physiol. 2010, 152, 1373–1390. [Google Scholar] [CrossRef]
- Wang, E.; Wang, J.; Zhu, X.D.; Hao, W.; Wang, L.Y.; Li, Q.; Zhang, L.X.; He, W.; Lu, B.R.; Lin, H.X.; et al. Control of rice grain-filling and yield by a gene with a potential signature of domestication. Nat. Genet. 2008, 40, 1370–1374. [Google Scholar] [CrossRef]
- Abu-Zaitoon, Y.M.; Bennett, K.; Normanly, J.; Nonhebel, H.M. A large increase in IAA during development of rice grains correlates with the expression of tryptophan aminotransferase OsTAR1 and a grain-specific yucca. Physiol. Plant. 2012, 146, 487–499. [Google Scholar] [CrossRef]
- Jin, Y.; Ni, D.-A.; Ruan, Y.-L. Posttranslational elevation of cell wall invertase activity by silencing its inhibitor in tomato delays leaf senescence and increases seed weight and fruit hexose level. Plant Cell Online 2009, 21, 2072–2089. [Google Scholar] [CrossRef]
- Bate, N.J.; Niu, X.P.; Wang, Y.W.; Reimann, K.S.; Helentjaris, T.G. An invertase inhibitor from maize localizes to the embryo surrounding region during early kernel development. Plant Physiol. 2004, 134, 246–254. [Google Scholar] [CrossRef]
- Cho, J.I.; Lee, S.K.; Ko, S.H.; Kim, H.K.; Jun, S.H.; Lee, Y.H.; Bhoo, S.H.; Lee, K.W.; An, G.H.; Hahn, T.R.; et al. Molecular cloning and expression analysis of the cell-wall invertase gene family in rice (Oryza sativa L.). Plant Cell Rep. 2005, 24, 225–236. [Google Scholar] [CrossRef]
- Ishimaru, T.; Hirose, T.; Matsuda, T.; Goto, A.; Takahashi, K.; Sasaki, H.; Terao, T.; Ishii, R.; Ohsugi, R.; Yamagishi, T. Expression patterns of genes encoding carbohydrate-metabolizing enzymes and their relationship to grain filling in rice (Oryza sativa L.): Comparison of caryopses located at different positions in a panicle. Plant Cell Physiol. 2005, 46, 620–628. [Google Scholar] [CrossRef]
- Hothorn, M.; van den Ende, W.; Lammens, W.; Rybin, V.; Scheffzek, K. Structural insights into the pH-controlled targeting of plant cell-wall invertase by a specific inhibitor protein. Proc. Natl. Acad. Sci. USA 2010, 107, 17427–17432. [Google Scholar] [CrossRef]
- Fu, F.-F.; Xue, H.-W. Coexpression analysis identifies rice starch regulator1, a rice AP2/EREBP family transcription factor, as a novel rice starch biosynthesis regulator. Plant Physiol. 2010, 154, 927–938. [Google Scholar] [CrossRef]
- Sekhon, R.S.; Lin, H.N.; Childs, K.L.; Hansey, C.N.; Buell, C.R.; de Leon, N.; Kaeppler, S.M. Genome-wide atlas of transcription during maize development. Plant J. 2011, 66, 553–563. [Google Scholar] [CrossRef]
- Dash, S.; van Hemert, J.; Hong, L.; Wise, R.P.; Dickerson, J.A. PLEXdb: Gene expression resources for plants and plant pathogens. Nucleic Acids Res. 2011, 40, D1194–D1201. [Google Scholar]
- Hong, M.J.; Kim, D.Y.; Lee, T.G.; Jeon, W.B.; Seo, Y.W. Functional characterization of pectin methylesterase inhibitor (PMEI) in wheat. Genes Genet. Syst. 2010, 85, 97–106. [Google Scholar] [CrossRef]
- Rausch, T.; Greiner, S. Plant protein inhibitors of invertases. Biochim. Biophys. Acta 1696, 1696, 253–261. [Google Scholar] [CrossRef]
- Liu, X.; Song, B.T.; Zhang, H.L.; Li, X.Q.; Xie, C.H.; Liu, J. Cloning and molecular characterization of putative invertase inhibitor genes and their possible contributions to cold-induced sweetening of potato tubers. Mol. Genet. Genomics 2010, 284, 147–159. [Google Scholar] [CrossRef]
- Mangelsen, E.; Kilian, J.; Harter, K.; Jansson, C.; Wanke, D.; Sundberg, E. Transcriptome analysis of high-temperature stress in developing barley caryopses: Early stress responses and effects on storage compound biosynthesis. Mol. Plant 2011, 4, 97–115. [Google Scholar] [CrossRef]
- Scognamiglio, M.A.; Ciardiello, M.A.; Tamburrini, M.; Carratore, V.; Rausch, T.; Camardella, L. The plant invertase inhibitor shares structural properties and disulfide bridges arrangement with the pectin methylesterase inhibitor. J. Protein Chem. 2003, 22, 363–369. [Google Scholar] [CrossRef]
- Jain, M.; Nijhawan, A.; Arora, R.; Agarwal, P.; Ray, S.; Sharma, P.; Kapoor, S.; Tyagi, A.K.; Khurana, J.P. F-box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress. Plant Physiol. 2007, 143, 1467–1483. [Google Scholar] [CrossRef]
- Xue, H.W. Institute of Plant Physiology and Ecology, Shanghai Institutes of Biological Sciences, Shanghai, China, Unpublished data. 2009.
- Li, M.N.; Xu, W.Y.; Yang, W.Q.; Kong, Z.S.; Xue, Y.B. Genome-wide gene expression profiling reveals conserved and novel molecular functions of the stigma in rice. Plant Physiol. 2007, 144, 1797–1812. [Google Scholar] [CrossRef]
- Fujita, M.; Horiuchi, Y.; Ueda, Y.; Mizuta, Y.; Kubo, T.; Yano, K.; Yamaki, S.; Tsuda, K.; Nagata, T.; Niihama, M.; et al. Rice expression atlas in reproductive development. Plant Cell Physiol. 2010, 51, 2060–2081. [Google Scholar] [CrossRef]
- Gao, L.-L.; Xue, H.-W. Global analysis of expression profiles of rice receptor-like kinase genes. Mol. Plant 2012, 5, 143–153. [Google Scholar] [CrossRef]
- Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. Mega5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef]
- Saitou, N.; Nei, M. The neighbor-joining method—A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar]
- Edgar, R.C. Muscle: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef]
- Felsenstein, J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef]
- Zuckerkandl, E.; Pauling, L. Evolutionary divergence and convergence in proteins. In Evolving Genes and Proteins; Bryson, V., Vogel, H.J., Eds.; Academic Press: New York, NY, USA, 1965; pp. 97–166. [Google Scholar]
- Mutwil, M.; Klie, S.; Tohge, T.; Giorgi, F.M.; Wilkins, O.; Campbell, M.M.; Fernie, A.R.; Usadel, B.; Nikoloski, Z.; Persson, S. PlaNet: Combined sequence and expression comparisons across plant networks derived from seven species. Plant Cell 2011, 23, 895–910. [Google Scholar] [CrossRef]
- Mutwil, M.; Obro, J.; Willats, W.G.T.; Persson, S. GeneCAT—Novel webtools that combine blast and co-expression analyses. Nucleic Acids Res. 2008, 36, W320–W326. [Google Scholar] [CrossRef]
- Cao, P.; Jung, K.-H.; Choi, D.; Hwang, D.; Zhu, J.; Ronald, P. The rice oligonucleotide array database: An atlas of rice gene expression. Rice 2012, 2. [Google Scholar] [CrossRef]
- Obayashi, T.; Hayashi, S.; Saeki, M.; Ohta, H.; Kinoshita, K. ATTED-II provides coexpressed gene networks for arabidopsis. Nucleic Acids Res. 2009, 37, D987–D991. [Google Scholar] [CrossRef]
- Thirumurugan, T.; Ito, Y.; Kubo, T.; Serizawa, A.; Kurata, N. Identification, characterization and interaction of hap family genes in rice. Mol. Genet. Genomics 2008, 279, 279–289. [Google Scholar] [CrossRef]
- Tsai, Y.-C.; Weir, N.R.; Hill, K.; Zhang, W.; Kim, H.J.; Shiu, S.-H.; Schaller, G.E.; Kieber, J.J. Characterization of genes involved in cytokinin signaling and metabolism from rice. Plant Physiol. 2012, 158, 1666–1684. [Google Scholar] [CrossRef]
- Saiga, S.; Möller, B.; Watanabe-Taneda, A.; Abe, M.; Weijers, D.; Komeda, Y. Control of embryonic meristem initiation in arabidopsis by phd-finger protein complexes. Development 2012, 139, 1391–1398. [Google Scholar] [CrossRef]
- Silverstein, K.A.T.; Moskal, W.A.; Wu, H.C.; Underwood, B.A.; Graham, M.A.; Town, C.D.; VandenBosch, K.A. Small cysteine-rich peptides resembling antimicrobial peptides have been under-predicted in plants. Plant J. 2007, 51, 262–280. [Google Scholar] [CrossRef]
- Magnard, J.L.; Lehouque, G.; Massonneau, A.E.; Frangne, N.; Heckel, T.; Gutierrez-Marcos, J.F.; Perez, P.; Dumas, C.; Rogowsky, P.M. ZmEBE genes show a novel, continuous expression pattern in the central cell before fertilization and in specific domains of the resulting endosperm after fertilization. Plant Mol. Biol. 2003, 53, 821–836. [Google Scholar] [CrossRef]
- Hwang, S.G.; Kim, D.S.; Jang, C.S. Comparative analysis of evolutionary dynamics of genes encoding leucine-rich repeat receptor-like kinase between rice and arabidopsis. Genetica 2011, 139, 1023–1032. [Google Scholar] [CrossRef]
- Jain, M.; Kaur, N.; Garg, R.; Thakur, J.K.; Tyagi, A.K.; Khurana, J.P. Structure and expression analysis of early auxin-responsive AUX/IAA gene family in rice (Oryza sativa). Funct. Integr. Genomics 2006, 6, 47–59. [Google Scholar] [CrossRef]
- Gray, W.M.; Kepinski, S.; Rouse, D.; Leyser, O.; Estelle, M. Auxin regulates SCFTIR1-dependent degradation of AUX/IAA proteins. Nature 2001, 414, 271–276. [Google Scholar] [CrossRef]
- LeClere, S.; Schmelz, E.A.; Chourey, P.S. Cell wall invertase-deficient miniature1 kernels have altered phytohormone levels. Phytochemistry 2008, 69, 692–699. [Google Scholar] [CrossRef]
- Yu, C.S.; Chen, Y.C.; Lu, C.H.; Hwang, J.K. Prediction of protein subcellular localization. Proteins 2006, 64, 643–651. [Google Scholar] [CrossRef]
- Emanuelsson, O.; Nielsen, H.; Brunak, S.; von Heijne, G. Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J. Mol. Biol. 2000, 300, 1005–1016. [Google Scholar] [CrossRef]
- Amien, S.; Kliwer, I.; Marton, M.L.; Debener, T.; Geiger, D.; Becker, D.; Dresselhaus, T. Defensin-like ZmES4 mediates pollen tube burst in maize via opening of the potassium channel kzm1. PLoS Biol. 2010, 8. [Google Scholar] [CrossRef]
- Costa, L.M.; Yuan, J.; Rouster, J.; Paul, W.; Dickinson, H.; Gutierrez-Marcos, J.F. Maternal control of nutrient allocation in plant seeds by genomic imprinting. Curr. Biol. 2012, 22, 160–165. [Google Scholar]
- Li, M.; Singh, R.; Bazanova, N.; Milligan, A.S.; Shirley, N.; Langridge, P.; Lopato, S. Spatial and temporal expression of endosperm transfer cell-specific promoters in transgenic rice and barley. Plant Biotechnol. J. 2008, 6, 465–476. [Google Scholar] [CrossRef]
- Kwong, R.W.; Bui, A.Q.; Lee, H.; Kwong, L.W.; Fischer, R.L.; Goldberg, R.B.; Harada, J.J. Leafy cotyledon1-like defines a class of regulators essential for embryo development. Plant Cell 2003, 15, 5–18. [Google Scholar] [CrossRef]
- Junker, A.; Mönke, G.; Rutten, T.; Keilwagen, J.; Seifert, M.; Thi, T.M.N.; Renou, J.-P.; Balzergue, S.; Viehöver, P.; Hähnel, U.; et al. Elongation-related functions of leafy cotyledon1 during the development of Arabidopsis thaliana. Plant J. 2012, 71, 427–442. [Google Scholar]
- Boisson-Dernier, A.; Kessler, S.A.; Grossniklaus, U. The walls have ears: The role of plant CrRLK1Ls in sensing and transducing extracellular signals. J. Exp. Botany 2011, 62, 1581–1591. [Google Scholar] [CrossRef]
- Goodstein, D.M.; Shu, S.; Howson, R.; Neupane, R.; Hayes, R.D.; Fazo, J.; Mitros, T.; Dirks, W.; Hellsten, U.; Putnam, N.; et al. Phytozome: A comparative platform for green plant genomics. Nucleic Acids Res. 2012, 40, D1178–D1186. [Google Scholar] [CrossRef]
- Nolan, T.; Hands, R.E.; Bustin, S.A. Quantification of mrna using real-time RT-PCR. Nat. Protoc. 2006, 1, 1559–1582. [Google Scholar] [CrossRef]
- Li, Q.F.; Sun, S.S.M.; Yuan, D.Y.; Yu, H.X.; Gu, M.H.; Liu, Q.Q. Validation of candidate reference genes for the accurate normalization of real-time quantitative RT-PCR data in rice during seed development. Plant Mol. Biol. Rep. 2010, 28, 49–57. [Google Scholar]
- Tomlinson, K.L.; McHugh, S.; Labbe, H.; Grainger, J.L.; James, L.E.; Pomeroy, K.M.; Mullin, J.W.; Miller, S.S.; Dennis, D.T.; Miki, B.L.A. Evidence that the hexose-to-sucrose ratio does not control the switch to storage product accumulation in oilseeds: Analysis of tobacco seed development and effects of overexpressing apoplastic invertase. J. Exp. Botany 2004, 55, 2291–2303. [Google Scholar] [CrossRef]
- Fournier, E. Colorimetric quantification of carbohydrates. In Current Protocols in Food Analytical Chemistry; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2001. [Google Scholar]
© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).