Proteomic identification of genes associated with maize grain-filling rate

Stimulus-responsive proteins involved in multi-process regulation of storage substance accumulation during rice grain filling under elevated temperature

BACKGROUND

The intensified global warming during grain filling deteriorated rice quality, in particular increasing the frequency of chalky grains which markedly impact market value. The formation of rice quality is a complex process influenced by multiple genes, proteins and physiological metabolic processes. Proteins responsive to stimulus can adjust the ability of plants to respond to unfavorable environments, which may be an important protein involved in the regulation of quality formation under elevated temperature. However, relatively few studies have hindered our further understanding of rice quality formation under elevated temperature.

RESULTS

We conducted the actual field elevated temperature experiment and performed proteomic analysis of rice grains at the early stage of grain filling. Starting with the response to stimulus in GO annotation, 22 key proteins responsive to stimulus were identified in the regulation of grain filling and response to elevated temperature. Among the proteins responsive to stimulus, during grain filling, an increased abundance of signal transduction and other stress response proteins, a decreased abundance of reactive oxygen species-related proteins, and an increased accumulation of storage substance metabolism proteins consistently contributed to grain filling. However, the abundance of probable indole-3-acetic acid-amido synthetase GH3.4, probable indole-3-acetic acid-amido synthetase GH3.8 and CBL-interacting protein kinase 9 belonged to signal transduction were inhibited under elevated temperature. In the reactive oxygen species-related protein, elevated temperature increased the accumulation of cationic peroxidase SPC4 and persulfide dioxygenase ETHE1 homolog to maintain normal physiological homeostasis. The increased abundance of alpha-amylase isozyme 3E and seed allergy protein RA5 was related to the storage substance metabolism, which regulated starch and protein accumulation under elevated temperature.

CONCLUSION

Auxin synthesis and calcium signal associated with signal transduction, other stress responses, protein transport and modification, and reactive oxygen species-related proteins may be key proteins responsive to stimulus in response to elevated temperature. Alpha-amylase isozyme 3E and seed allergy protein RA5 may be the key proteins to regulate grain storage substance accumulation and further influence quality under elevated temperature. This study enriched the regulatory factors involved in the response to elevated temperature and provided a new idea for a better understanding of grain response to temperature.

Physiological and omics analysis of maize inbred lines during late grain development

BACKGROUND

There were significant differences in the change of moisture content and grain composition at the late stage of grain development among different maize varieties, but the regulation mechanism is not clear.

OBJECTIVE

To explore the key genes causing the variation in physiological traits of two typical maize inbred lines in late grain development.

METHODS

The grains at different development stages were selected as materials to determine the content of water, sucrose, starch and ABA. Transcriptomic and proteomic analysis of the materials were performed to screen relevant genes.

RESULTS

The grain dehydration rate and the content of sucrose, starch and ABA were showed significant differences between two varieties in the late stage of grain development. The enrichment analysis of common differentially expressed genes (proteins) showed that most of the genes (proteins) were enriched in the extracellular region. The downregulated genes were mainly concentrated in carbohydrate metabolism and lipid metabolism, while the upregulated genes were mainly in response to stress. Furthermore, this study also identified many key candidate genes (dehydrin genes, pathogenesis-related genes, sucrose synthase and secondary metabolites related genes) related to late grain development of maize.

CONCLUSIONS

The suggested genes related to late grain development of maize can be candidates for further functional study.

Comparative Transcriptome and Proteome Analysis Provides New Insights Into the Mechanism of Protein Synthesis in Kenaf (Hibiscus cannabinus L.) Leaves

Given the rising domestic demand and increasing global prices of corn and soybean, China is looking for alternatives for these imports to produce animal fodder. Kenaf (Hibiscus cannabinus L.) has great potential as a new forage source, due to abundant proteins, phenols and flavonoids in its leaves. However, few studies have evaluated the mechanism of protein synthesis in kenaf leaves. In the current work, compared with kenaf material "L332," the percentage of crude protein content in leaves of material "Q303" increased by 6.13%; combined with transcriptome and proteome data, the kenaf samples were systematically studied to obtain mRNA-protein correlation. Then, the genes/proteins related to protein synthesis in the kenaf leaves were obtained. Moreover, this work detected mRNA expression of 20 differentially expressed genes (DEGs). Meanwhile, 20 differentially expressed proteins (DEPs) related to protein synthesis were performed parallel reaction monitoring. Fructose-1,6-bisphosphatase (FBP), nitrite reductase (NirA), prolyl tRNA synthase (PARS) and glycine dehydrogenase (GLDC) presented increased mRNA and protein levels within kenaf leaves with high protein content. Based on the obtained findings, FBP, NirA, PARS, and GLDC genes may exert a vital function in the protein synthesis of kenaf leaves. The results provide a new idea for further studying the potential genes affecting the quality trait of protein content in kenaf leaves and provide gene resources and a theoretical foundation for further cultivating high protein kenaf varieties.

Comparative transcriptome analysis provides insights into grain filling commonalities and differences between foxtail millet [Setaria italica (L.) P. Beauv.] varieties with different panicle types

Grain filling affects grain weight and quality and is among the most critical factors in determining the yield and quality of cereal crops. Though hybrids have larger panicles and numerous spikelets with a larger sink capacity than conventional varieties, data on the grain filling commonalities and differences between foxtail millet varieties with different panicle types remain sparse. In this study, we found that "Zhang Gu 13" (ZG, large panicle) exhibits a significantly higher panicle weight than "Yu Gu 18" (YG, conventional panicle) at the early stage of grain filling, but the weight of YG increased rapidly and gradually overtook ZG during the middle stages. A temporal expression pattern analysis demonstrated that the genes involved in photosynthesis, metabolic pathways, and phenylpropanoid biosynthesis were downregulated, while those related to peroxisome function, purine metabolism, and zeatin biosynthesis were upregulated during grain filling in both varieties. A total of 6,832 differentially expressed genes (DEGs) were identified in both varieties, with the majority identified at the early and late stages. A Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis further revealed that the upregulated DEGs in YG were associated with gibberellin (GA) biosynthesis, ATP-binding cassette (ABC) transporters, and plant hormone signal transduction. Photosynthesis-related DEGs, such as photosystem and antenna proteins, were significantly upregulated in ZG. This study provides preliminary insights into the differences in gene expression and molecular mechanisms of grain filling between ZG and YG in the North China summer-sowing region.

Identification of a new mutant allele of ZmYSL2 that regulates kernel development and nutritional quality in maize

The discovery and characterization of the opaque endosperm gene provide ideas and resources for the production and application of maize. We found an o213 mutant whose phenotype was opaque and shrunken endosperm with semi-dwarf plant height. The protein, lipid, and starch contents in the o213 endosperm were significantly decreased, while the free amino acid content in the o213 endosperm significantly increased. The aspartic acid, asparagine, and lysine contents were raised in the o213 endosperm by 6.5-, 8.5-, and 1.7-fold, respectively. Genetic analysis showed that this o213 mutant is a recessive single-gene mutation. The position mapping indicated that o213 is located in a 468-kb region that contains 11 protein-encoding genes on the long arm of chromosome 5. The coding sequence analysis of candidate genes between the WT and o213 showed that ZmYSL2 had only a single-base substitution (A-G) in the fifth exon, which caused methionine substitution to valine. Sequence analysis and the allelic test showed that o213 is a new mutant allele of ZmYSL2. The qRT-PCR results indicated that o213 is highly expressed in the stalks and anthers. Subcellular localization studies showed that o213 is a membrane transporter. In the variation analysis of o213, the amplification of 65 inbred lines in GWAS showed that this 3-bp deletion of the first exon of o213 was found only in temperate inbred lines, implying that the gene was artificially affected in the selection process. Our results suggest that o213 is an important endosperm development gene and may serve as a genetic resource.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01278-9.

Spatial distribution of proteins and metabolites in developing wheat grain and their differential regulatory response during the grain filling process

Grain filling and grain development are essential biological processes in the plant's life cycle, eventually contributing to the final seed yield and quality in all cereal crops. Studies of how the different wheat (Triticum aestivum L.) grain components contribute to the overall development of the seed are very scarce. We performed a proteomics and metabolomics analysis in four different developing components of the wheat grain (seed coat, embryo, endosperm, and cavity fluid) to characterize molecular processes during early and late grain development. In-gel shotgun proteomics analysis at 12, 15, 20, and 26 days after anthesis (DAA) revealed 15 484 identified and quantified proteins, out of which 410 differentially expressed proteins were identified in the seed coat, 815 in the embryo, 372 in the endosperm, and 492 in the cavity fluid. The abundance of selected protein candidates revealed spatially and temporally resolved protein functions associated with development and grain filling. Multiple wheat protein isoforms involved in starch synthesis such as sucrose synthases, starch phosphorylase, granule-bound and soluble starch synthase, pyruvate phosphate dikinase, 14-3-3 proteins as well as sugar precursors undergo a major tissue-dependent change in abundance during wheat grain development suggesting an intimate interplay of starch biosynthesis control. Different isoforms of the protein disulfide isomerase family as well as glutamine levels, both involved in the glutenin macropolymer pattern, showed distinct spatial and temporal abundance, revealing their specific role as indicators of wheat gluten quality. Proteins binned into the functional category of cell growth/division and protein synthesis/degradation were more abundant in the early stages (12 and 15 DAA). At the metabolome level all tissues and especially the cavity fluid showed highly distinct metabolite profiles. The tissue-specific data are integrated with biochemical networks to generate a comprehensive map of molecular processes during grain filling and developmental processes.

Combined linkage mapping and association analysis reveals genetic control of maize kernel moisture content

The free moisture in crop kernels after being naturally dried is referred to as kernel moisture content (KMC). Maize KMC reflects grain quality and influences transportation and storage of seeds. We used an IBM Syn10 DH maize population consisting of 249 lines and an association panel comprising 310 maize inbred lines to identify the genetic loci affecting maize KMC in three environments. Using the IBM population detected 13 QTL on seven chromosomes, which were clustered into nine common QTL. Genome-wide association analysis (GWAS) identified 16 significant SNPs across the 3 environments, which were linked to 158 genes across the three environments. Combined QTL mapping and GWAS found two SNPs that were located in two of the mapped QTL, respectively. Twenty-three genes were linked with the loci co-localized in both populations. Of these 181 genes, five have previously been reported to be associated with KMC or to regulate seed development. These associations were verified by candidate gene association analysis. Two superior alleles and one favorable haplotype for Zm00001d007774 and Zm00001d047868 were found to influence KMC. These findings provide insights into molecular mechanisms underlying maize KMC and contribute to the use of marker-assisted selection for breeding low-KMC maize.

Label-Free Comparative Proteomic Analysis Combined with Laser-Capture Microdissection Suggests Important Roles of Stress Responses in the Black Layer of Maize Kernels

The black layer (BL) is traditionally used as an indicator for kernel harvesting in maize, as it turns visibly dark when the kernel reaches physiological maturity. However, the molecular roles of BL in kernel development have not been fully elucidated. In this work, microscopy images showed that BL began to appear at a growth stage earlier than 10 days after pollination (DAP), and its color gradually deepened to become dark as the development period progressed. Scanning electron microscopy observations revealed that BL is a tissue structure composed of several layers of cells that are gradually squeezed and compressed during kernel development. Laser-capture microdissection (LCM) was used to sample BL and its neighboring inner tissue, basal endosperm transfer layer (BETL), and outer tissue, inner epidermis (IEP), from 20 DAP of kernels. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry profiling (MALDI-TOF MS profiling) detected 41, 104, and 120 proteins from LCM-sampled BL, BETL, and IEP, respectively. Gene ontology (GO) analysis indicated that the 41 BL proteins were primarily involved in the response to stress and stimuli. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis found that the BL proteins were enriched in several defense pathways, such as the ascorbate and aldarate metabolic pathways. Among the 41 BL proteins, six were BL-specific proteins that were only detected from BL. Annotations of five BL-specific proteins were related to stress responses. During kernel development, transcriptional expression of most BL proteins showed an increase, followed by a decrease, and reached a maximum zero to 20 DAP. These results suggest a role for BL in stress responses for protecting filial tissue against threats from maternal sides, which helps to elucidate the biological functions of BL.

Global Dynamics in Protein Disorder during Maize Seed Development

Intrinsic protein disorder is a physicochemical attribute of some proteins lacking tridimensional structure and is collectively known as intrinsically disordered proteins (IDPs). Interestingly, several IDPs have been associated with protective functions in plants and with their response to external stimuli. To correlate the modulation of the IDPs content with the developmental progression in seed, we describe the expression of transcripts according to the disorder content of the proteins that they codify during seed development, from the early embryogenesis to the beginning of the desiccation tolerance acquisition stage. We found that the total expression profile of transcripts encoding for structured proteins is highly increased during middle phase. However, the relative content of protein disorder is increased as seed development progresses. We identified several intrinsically disordered transcription factors that seem to play important roles throughout seed development. On the other hand, we detected a gene cluster encoding for IDPs at the end of the late phase, which coincides with the beginning of the acquisition of desiccation tolerance. In conclusion, the expression pattern of IDPs is highly dependent on the developmental stage, and there is a general reduction in the expression of transcripts encoding for structured proteins as seed development progresses. We proposed maize seeds as a model to study the regulation of protein disorder in plant development and its involvement in the acquisition of desiccation tolerance in plants.

Proteomic analysis reveals that auxin homeostasis influences the eighth internode length heterosis in maize (Zea mays)

Ear height is an important maize morphological trait that influences plant lodging resistance in the field, and is based on the number and length of internodes under the ear. To explore the effect of internodes on ear height, the internodes under the ear were analysed in four commercial hybrids (Jinsai6850, Zhengdan958, Xundan20, and Yuyu22) from different heterotic groups in China. The eighth internode, which is the third aboveground extended internode, exhibited high-parent or over high-parent heterosis and contributed considerably to ear height. Thus, the proteome of the eighth internode was examined. Sixty-six protein spots with >1.5-fold differences in accumulation (P < 0.05) among the four hybrids were identified by mass spectrometry and data analyses. Most of the differentially accumulated proteins exhibited additive accumulation patterns, but with epistatic effects on heterosis performance. Proteins involved in phenylpropanoid and benzoxazinoid metabolic pathways were observed to influence indole-3-acetic acid biosynthesis and polar auxin transport during internode development. Moreover, indole-3-acetic acid content was positively correlated with the eighth internode length, but negatively correlated with the extent of the heterosis of the eighth internode length.

Identification of genes alternatively spliced in developing maize endosperm

The process of alternative splicing is critical for the regulation of growth and development of plants. Thus far, little is known about the role of alternative splicing in the regulation of maize (Zea mays L.) endosperm development. RNA sequencing (RNA-seq) data of endosperms from two maize inbred lines, Mo17 and Ji419, at 15 and 25 days after pollination (DAP), respectively, were used to identify genes that were alternatively spliced during endosperm development. Intron retention (IR) in GRMZM2G005887 was further validated using PCR and re-sequencing technologies. In total, 49,000 alternatively spliced events and ca. 20,000 alternatively spliced genes were identified in the two maize inbred lines. Of these, 30 genes involved in amino acid biosynthesis and starch biosynthesis were identified, with IR occurring only in a specific sample, and were significantly co-expressed with ten well-known genes related to maize endosperm development. Moreover, IR in GRMZM2G005887, which encodes a cysteine synthase, was confirmed to occur only in the endosperm of Mo17 at 15 DAP, resulting in the retention of a 121-bp fragment in its 5' untranslated region. Two cis-acting regulatory elements, CAAT-box and TATA-box were observed in the retained fragment in Mo17 at 15 DAP; this could regulate the expression of this gene and influence endosperm development. The results suggest that the 30 genes with IR identified herein might be associated with maize endosperm development, and are likely to play important roles in the developing maize endosperm.

Proteomic analysis of maize grain development using iTRAQ reveals temporal programs of diverse metabolic processes

BACKGROUND

Grain development in maize is an essential process in the plant's life cycle and is vital for use of the plant as a crop for animals and humans. However, little is known regarding the protein regulatory networks that control grain development. Here, isobaric tag for relative and absolute quantification (iTRAQ) technology was used to analyze temporal changes in protein expression during maize grain development.

RESULTS

Maize grain proteins and changes in protein expression at eight developmental stages from 3 to 50 d after pollination (DAP) were performed using iTRAQ-based proteomics. Overall, 4751 proteins were identified; 2639 of these were quantified and 1235 showed at least 1.5-fold changes in expression levels at different developmental stages and were identified as differentially expressed proteins (DEPs). The DEPs were involved in different cellular and metabolic processes with a preferential distribution to protein synthesis/destination and metabolism categories. A K-means clustering analysis revealed coordinated protein expression associated with different functional categories/subcategories at different development stages.

CONCLUSIONS

Our results revealed developing maize grain display different proteomic characteristics at distinct stages, such as numerous DEPs for cell growth/division were highly expressed during early stages, whereas those for starch biosynthesis and defense/stress accumulated in middle and late stages, respectively. We also observed coordinated expression of multiple proteins of the antioxidant system, which are essential for the maintenance of reactive oxygen species (ROS) homeostasis during grain development. Particularly, some DEPs, such as zinc metallothionein class II, pyruvate orthophosphate dikinase (PPDK) and 14-3-3 proteins, undergo major changes in expression at specific developmental stages, suggesting their roles in maize grain development. These results provide a valuable resource for analyzing protein function on a global scale and also provide new insights into the potential protein regulatory networks that control grain yield and quality.

Genome-Wide Identification and Characterization of microRNAs in Developing Grains of Zea mays L

The development and maturation of maize kernel involves meticulous and fine gene regulation at transcriptional and post-transcriptional levels, and miRNAs play important roles during this process. Although a number of miRNAs have been identified in maize seed, the ones involved in the early development of grains and in different lines of maize have not been well studied. Here, we profiled four small RNA libraries, each constructed from groups of immature grains of Zea mays inbred line Chang 7–2 collected 4–6, 7–9, 12–14, and 18–23 days after pollination (DAP). A total of 40 known (containing 111 unique miRNAs) and 162 novel (containing 196 unique miRNA candidates) miRNA families were identified. For conserved and novel miRNAs with over 100 total reads, 44% had higher accumulation before the 9^th DAP, especially miR166 family members. 42% of miRNAs had highest accumulation during 12–14 DAP (which is the transition stage from embryogenesis to nutrient storage). Only 14% of miRNAs had higher expression 18–23 DAP. Prediction of potential targets of all miRNAs showed that 165 miRNA families had 377 target genes. For miR164 and miR166, we showed that the transcriptional levels of their target genes were significantly decreased when co-expressed with their cognate miRNA precursors in vivo. Further analysis shows miR159, miR164, miR166, miR171, miR390, miR399, and miR529 families have putative roles in the embryogenesis of maize grain development by participating in transcriptional regulation and morphogenesis, while miR167 and miR528 families participate in metabolism process and stress response during nutrient storage. Our study is the first to present an integrated dynamic expression pattern of miRNAs during maize kernel formation and maturation.

Identification and Characterization of microRNAs during Maize Grain Filling

The grain filling rate is closely associated with final grain yield of maize during the period of maize grain filling. To identify the key microRNAs (miRNAs) and miRNA-dependent gene regulation networks of grain filling in maize, a deep-sequencing technique was used to research the dynamic expression patternsof miRNAs at four distinct developmental grain filling stages in Zhengdan 958, which is an elite hybrid and cultivated widely in China. The sequencing result showed that the expression amount of almost all miRNAs was changing with the development of the grain filling and formed in seven groups. After normalization, 77 conserved miRNAs and 74 novel miRNAs were co-detected in these four samples. Eighty-one out of 162 targets of the conserved miRNAs belonged to transcriptional regulation (81, 50%), followed by oxidoreductase activity (18, 11%), signal transduction (16, 10%) and development (15, 9%). The result showed that miRNA 156, 393, 396 and 397, with their respective targets, might play key roles in the grain filling rate by regulating maize growth, development and environment stress response. The result also offered novel insights into the dynamic change of miRNAs during the developing process of maize kernels and assistedin the understanding of how miRNAs are functioning about the grain filling rate.

Proteomics of seed development, desiccation tolerance, germination and vigor

Proteomics, the large-scale study of the total complement of proteins in a given sample, has been applied to all aspects of seed biology mainly using model species such as Arabidopsis or important agricultural crops such as corn and rice. Proteins extracted from the sample have typically been separated and quantified by 2-dimensional polyacrylamide gel electrophoresis followed by liquid chromatography and mass spectrometry to identify the proteins in the gel spots. In this way, qualitative and quantitative changes in the proteome during seed development, desiccation tolerance, germination, dormancy release, vigor alteration and responses to environmental factors have all been studied. Many proteins or biological processes potentially important for each seed process have been highlighted by these studies, which greatly expands our knowledge of seed biology. Proteins that have been identified to be particularly important for at least two of the seed processes are involved in detoxification of reactive oxygen species, the cytoskeleton, glycolysis, protein biosynthesis, post-translational modifications, methionine metabolism, and late embryogenesis-abundant (LEA) proteins. It will be useful for molecular biologists and molecular plant breeders to identify and study genes encoding particularly interesting target proteins with the aim to improve the yield, stress tolerance or other critical properties of our crop species.

Genome-wide identification and characterization of maize expansin genes expressed in endosperm

By promoting cell wall loosening, expansins contribute to cell enlargement during various developmental processes. Nevertheless, the role of expansins in the expansion and development of endosperm--a major seed component whose cell size is significantly associated with grain yield--is poorly understood. To explore associated biological processes and the evolution of expansins in maize, we performed a systematic analysis of the expansin gene family encompassing gene structure, phylogeny, chromosomal location, gene duplication, and gene ontology. A total of 88 maize expansin genes (ZmEXPs) were identified and categorized into three subfamilies according to their phylogenetic relationships. Expression patterns of ZmEXPs were also investigated in nine different tissues by semi-quantitative RT-PCR. The expression of eight ZmEXPs was detected in endosperm, with five showing endosperm-specific expression. Quantitative RT-PCR was used to analyze expression patterns of the eight ZmEXPs in endosperm (10 days after pollination) under abscisic acid (ABA) and gibberellic acid (GA3) treatments. All eight ZmEXPs were found to be significantly regulated by ABA and GA3 in endosperm, suggesting important roles for these hormones in the regulation of ZmEXPs during endosperm development. Our results provide essential information for ZmEXPs cloning and functional exploration, which will assist research on expansin-related mechanisms and contribute to future enhancement of maize grain yield.

Maize yield response to water supply and fertilizer input in a semi-arid environment of Northeast China

Maize grain yield varies highly with water availability as well as with fertilization and relevant agricultural management practices. With a 311-A optimized saturation design, field experiments were conducted between 2006 and 2009 to examine the yield response of spring maize (Zhengdan 958, Zea mays L) to irrigation (I), nitrogen fertilization (total nitrogen, urea-46% nitrogen,) and phosphorus fertilization (P₂O₅, calcium superphosphate-13% P₂O₅) in a semi-arid area environment of Northeast China. According to our estimated yield function, the results showed that N is the dominant factor in determining maize grain yield followed by I, while P plays a relatively minor role. The strength of interaction effects among I, N and P on maize grain yield follows the sequence N+I >P+I>N+P. Individually, the interaction effects of N+I and N+P on maize grain yield are positive, whereas that of P+I is negative. To achieve maximum grain yield (10506.0 kg·ha⁻¹) for spring maize in the study area, the optimum application rates of I, N and P are 930.4 m³·ha⁻¹, 304.9 kg·ha⁻¹ and 133.2 kg·ha⁻¹ respectively that leads to a possible economic profit (EP) of 10548.4 CNY·ha⁻¹ (CNY, Chinese Yuan). Alternately, to obtain the best EP (10827.3 CNY·ha⁻¹), the optimum application rates of I, N and P are 682.4 m³·ha⁻¹, 241.0 kg·ha⁻¹ and 111.7 kg·ha⁻¹ respectively that produces a potential grain yield of 10289.5 kg·ha⁻¹.