Identification zmzfp2 mutant form small kernels
In our previous study, we isolated a candidate gene, RING Zinc Finger Protein (ZmZFP2, Zm00001d033780), from two crucial developmental stages of endosperm in a dent maize inbred Dan232 with large-sized kernels, and a popcorn inbred N04 with small-sized kernels (Li, et al. 2012; Liu, et al. 2010). ZmZFP2 comprises six exons and a C4HC3-type RING zinc finger domain (Fig. S1A-B). We identified a mutant for ZmZFP2 alleles in the Maize EMS induced Mutant Database (MEMD) (http://www.elabcaas.cn/memd/) (Lu, et al. 2018). The zmzfp2-ems allele is caused by a typical EMS-induced C to G mutation, resulting in a codon change of CAG (Arginine) to TAG (stop codon) within the open reading frame (Fig. S1A and C), leading to premature termination of the mutant allele translation (Fig. S1B). Subsequently, homozygous zmzfp2-ems mutant and corresponding wild-type (WT) were selected from self-crossed heterozygous mutants. The zmzfp2-ems mutant plants exhibited a significant decrease in plant height compared to the WT (Fig. 1A and 1F), while the phenotypic difference in ear appearance was not apparent. Quantitative analysis revealed that the 100-kernel weight of zmzfp2-ems decreased by 34.9%, and the length and width of the kernels were reduced by 9.6% and 17.0% respectively, compared to the WT (Fig. 1B-D and 1G-I). Quantitative RT-PCR analysis showed that the transcript levels of ZmZFP2 in the zmzfp2-ems mutant were lower than in its corresponding WT (Fig. 1E). These results indicate that the EMS-induced zmzfp2-ems mutant affects kernel weight and size.
The function of ZmZFP2 in maize transgenic lines
We employed the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated protein 9 (Cas9) technology strategy to target two sites in the second exon of ZmZFP2 and performed genetic transformation in the maize inbred line B104 (Fig. 2A-B). Two CRISPR editing events, zmzfp2-cr1 and zmzfp2-cr2, were generated, with zmzfp2-cr1 showing a 1-bp insertion within the open reading frame and zmzfp2-cr2 showing a 2-bp deletion (Fig. 2A). Phenotypic analysis of the kernels revealed that both zmzfp2-cr1 and zmzfp2-cr2 exhibited significant reductions in kernel 100-kernel weight, length and width compared to their respective WT controls, with reductions of 26.0% and 21.6% in 100-kernel weight, 17.4% and 7.7% in length, and 17.2% and 12.9% in width, respectively (Fig. 2C-D and E-G). These phenotypes were similar to those observed in the zmzfp2-ems mutant.
In Arabidopsis, transgenic plants over-expressing ZmZFP2 under the control of the CaMV35S promoter exhibited a significant increase in leaf area and plant height compared to WT plants (Fig. S2A-B and D-E). Additionally, the 1000-seed weight, seed width, and seed length of Arabidopsis were significantly enhanced (Fig. S2C and F-G). Next, we over-expressed ZmZFP2 in the maize inbred line B104 under the control of a maize ubiquitin promoter (Ubi::ZmZFP2). We assessed the agronomic traits of transgenic over-expression ZmZFP2 plants in the T2 generation, and the ear morphology was slightly larger than that of the WT (Fig. 3A). The transcription levels of ZmZFP2 in the over-expressing cell lines were 3-fold, 7-fold, and 2-fold higher than the WT (Fig. 3D). In comparison to the WT, over-expression of ZmZFP2 in maize increased the 100-kernel weight and kernel length, with little change in the kernel width (Fig. 3B-C and E-G). Statistical analysis showed that the 100-kernel weight in over-expression lines, ZmZFP2-OE1, OE2, and OE3, was enhanced by 5.9%, 9.8%, and 6.6%, respectively (Fig. 3E). Agronomic traits of the three OE lines collectively indicate that over-expressing ZmZFP2 significantly increases kernel weight and enlarges kernel size. These results, in combination with the results in the zmzfp2 mutants, suggest that ZmZFP2 positively regulates kernel weight and size.
The effect of ZmZFP2 on endosperm cell phenotype, starch, and protein accumulation
Given that ZmZFP2 affects the size and weight of maize kernels, the development process of the kernels is determined by the division and expansion of endosperm cells (Dante, et al. 2014; Ji, et al. 2022). We analyzed the longitudinal sections of the endosperm of kernels at 6, 8 and 16 DAP (days after pollination) (Fig. 4A-B and 4G-H). The morphology of the endosperm of the zmzfp2-ems mutant (Fig. 4A-B), zmzfp2-cr1 mutant (Fig. 4G), and ZmZFP2-OE2 (Fig. 4H) showed significant differences compared to their corresponding WT, with Fig. 4A and Fig. 4B being the magnified views of Fig. 4D and Fig. 4E, respectively. From the figures, it can be observed that the grain-filling of the endosperm in the zmzfp2-ems and zmzfp2-cr1 mutants were lower than that of the WT, particularly at the top region, in which the difference between the kernels at 8 DAP and their corresponding WT was more significant compared to that of the 16 DAP kernels and their corresponding WT. We found that the size of the endosperm in the zmzfp2-ems mutant and zmzfp2-cr1 mutant were smaller than that of the WT at both 6, 8 and 16 DAP. Statistical analysis revealed a significant reduction in the number of endosperm cells in the zmzfp2-ems mutant and zmzfp2-cr1 mutant compared to the WT, while the cell size difference was not significant (Fig. 4C-F and 4I-J).
By observing the top region of the endosperm at 16 DAP through transmission electron microscopy, the starch granules of the zmzfp2-ems mutant were smaller and more irregular compared to the WT (Fig. 5A), and the protein bodies were also smaller (Fig. 5B). In the mature dry kernels of the zmzfp2-ems mutant and zmzfp2-cr1 mutants, the starch contents were reduced by 5.1% and 6.6% respectively compared to the WT, and the total protein content decreased by 8.9% and 6.2% respectively (Fig. 5C-D and 5E-F). The content of zein in the zmzfp2-ems and zmzfp2-cr1 mutants significantly decreased, while the difference in non-zein content was not significant (Fig. S3A-D). The starch, total protein, and zein content in the mature kernels of ZmZFP2-OE2 increased by 4.8% and 7.6% respectively. These results indicate that the decrease in kernel size and weight caused by ZmZFP2 may be affected by the number of endosperm cells, as well as related to the accumulation of starch and protein.
Expression profile and subcellular localization of ZmZFP2
We employed quantitative RT-PCR analysis to measure the expression levels of ZmZFP2 in the inbred line B73, allowing to analyze the spatio-temporal expression patterns of ZmZFP2 in different tissues and developmental stages of maize. The highest expression of ZmZFP2 was observed in the shoot apical meristem (SAM) and the ear, with detectable expression in kernels at different developmental stages (Fig. 6A). Additionally, we utilized RNA in situ hybridization to explore the expression of ZmZFP2 in kernels at 10 and 15 DAP. Figure 6B showed abundant expressed of ZmZFP2 in the apical region of the embryo in the kernel, with the highest expression observed in the 15 DAP embryo (Fig. 6B). These in situ hybridization results were consistent with the developmental morphology of ZmZFP2 in the endosperm during grain-filling (Fig. 6A, Fig. 4A-B).
We conducted a search of the National Center for Biotechnology Information (NCBI) database and identified orthologs of ZmZFP2 in eudicots including soybean (Glycine max) and cotton (Gossypium hirsutum), and monocots including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), millet (Setaria italica) and sorghum (Sorghum bicolor) (Fig. S4). The multiple sequence alignment results showed high similarities between the ZmZFP2 protein and its orthologs, with amino acid sequence similarities of 70%, 88%, 94%, 94% and 84% with sorghum and Arabidopsis, rice, millet, sorghum and wheat, respectively (Fig. S5A). Additionally, an analysis of the conserved domain demonstrated that the ZmZFP2 protein contained a single RING domain at the N-terminal and three putative transmembrane domains in its C-terminal region, which were conserved (Fig. S5A). The consensus sequence of RING zinc finger proteins Cys-X2-Cys-X(9−39)-Cys-X1-Cys/His-X2-Cys-X(4−48)-Cys-X2-Cys was highly conserved across different species (Fig. S5A).
To investigate the subcellular localization of ZmZFP2, the full-length ZmZFP2 protein were fused to the N terminus of green fluorescent protein (GFP). The resulting constructs driven by the constitutive 35S promoter were first transiently expressed in tobacco (Nicotiana benthamiana) leaves. ZmZFP2-GFP signals appears to be closely around the DAPI (4',6-diamidino-2-phenylindole) in the nucleus (Fig. 6C). To further determine ZmZFP2 localization, similar around the nucleus localization of ZmZFP2-GFP was observed in maize leaf mesophyll protoplasts (Fig. 6D). Using TMHMM 2.0 analysis (https://services.healthtech.dtu.dk/services/TMHMM-2.0/), we identified three transmembrane domains in ZmZFP2 (Fig. S5A). It revealed that ZmZFP2-GFP containing the entire ZmZFP2 protein localizes to the nuclear outer membrane. Using the yeast GAL4 system, the result showed ZmZFP2 has no intrinsic activation properties (Fig. S5B). Signal peptide was not predicted in ZmZFP2 protein sequence using SignalP 3.0 (https://services.healthtech.dtu.dk/services/SignalP-3.0/). These results showed that ZmZFP2-GFP distinctly localized to the nuclear outer membrane without transcriptional activation.
Transcriptome Analysis of Genes Regulated by ZmZFP2
To access the changes in gene expression associated with ZmZFP2, we performed a transcriptome analysis using RNA-seq on 8-DAP kernels of zmzfp2-ems and WT. A total of 266 upregulated and 1156 downregulated transcripts were identified in the zmzfp2-ems relative to the WT, indicating that ZmZFP2 has a more positive effect on gene expression than a negative effect (Fig. 7A). Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of all the differentially expressed genes revealed that enrichment of genes was mainly involved in organic metabolism, compound biosynthesis, and signaling pathways (Fig. 7B; Supplemental Data 1). For instance, four genes involved in the MAPK signaling pathway in plant were also strongly enriched in plant hormone signaling transduction (Fig. 7C; Supplementary Data 2). In addition, four genes, including Zm00001d005823 (flavonoid 3'-monooxygenase), Zm00001d020628 (trimethyltridecatetraene synthase-like), Zm00001d022351, Zm00001d053617 (brassinosteroid catabolism 1) were significantly downregulated and enriched in brassinosteroid biosynthesis (Fig. 7B-C). Zeatin is a common plant cytokinin and plays a crucial role in cell division. We found that the differentially expressed genes were enriched in zeatin biosynthesis (Fig. 7B-C), including Zm00001d000237 (cis-zeatin O-glucosyltransferase 1), Zm00001d012335 (cis-zeatin O-glucosyltransferase 2), Zm00001d012641 (cytokinin oxidase4b), Zm00001d016417 (UDP-glycosyltransferase), Zm00001d032664 (cytokinin oxidase6), and Zm00001d039643 (cytokinin-O-glucosyltransferase 3), indicating that ZmZFP2 regulates kernel width via cell division. Starch and sucrose metabolism were also significantly enriched, including genes Zm00001d014083 (beta amylase), Zm00001d018082 (trehalose-phosphate phosphatase 4), Zm00001d025354 (invertase cell wall5), Zm00001d032118 (trehalose-6-phosphate synthase3), and Zm00001d035037 (fructokinase2). All of these genes were found to be downregulated in zmzfp2-cr1 compared with WT (Fig. S6A-D). Endoplasmic reticulum associated degradation (ERAD) is a primary method to degrade proteins to maintain cell homeostasis, which requires the ubiquitination of substrate proteins by E1 (ubiquitin-activating enzyme), E2s (ubiquitin-conjugating enzymes) and E3s (ubiquitin E3 ligases) (Chen, et al. 2020). We found that eleven differentially expressed genes were significantly enriched in protein processing in the endoplasmic reticulum, which is consistent with its ubiquitin ligase function (Fig. 7B).