Next Article in Journal
Developmental Instability and Gene Dysregulation in an Extracted Tetraploid from Hexaploid Wheat
Next Article in Special Issue
Transcriptomic and Metabolomic Analyses Reveal the Key Genes Related to Shade Tolerance in Soybean
Previous Article in Journal
GMM-Based Expanded Feature Space as a Way to Extract Useful Information for Rare Cell Subtypes Identification in Single-Cell Mass Cytometry
Previous Article in Special Issue
Genomic Prediction of Resistance to Tan Spot, Spot Blotch and Septoria Nodorum Blotch in Synthetic Hexaploid Wheat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Identification of the Pum Gene Family and Its Involvement in Kernel Development in Maize

1
Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Maize Research Institute, Sichuan Agricultural University, Chengdu 611130, China
2
National Research Centre of Intercropping, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(18), 14036; https://doi.org/10.3390/ijms241814036
Submission received: 3 August 2023 / Revised: 7 September 2023 / Accepted: 12 September 2023 / Published: 13 September 2023

Abstract

:
The Pumilio (Pum) RNA-binding protein family regulates post-transcription and plays crucial roles in stress response and growth. However, little is known about Pum in plants. In this study, a total of 19 ZmPum genes were identified and classified into two groups in maize. Although each ZmPum contains the conserved Pum domain, the ZmPum members show diversity in the gene and protein architectures, physicochemical properties, chromosomal location, collinearity, cis-elements, and expression patterns. The typical ZmPum proteins have eight α-helices repeats, except for ZmPum2, 3, 5, 7, and 14, which have fewer α-helices. Moreover, we examined the expression profiles of ZmPum genes and found their involvement in kernel development. Except for ZmPum2, ZmPum genes are expressed in maize embryos, endosperms, or whole seeds. Notably, ZmPum4, 7, and 13 exhibited dramatically high expression levels during seed development. The study not only contributes valuable information for further validating the functions of ZmPum genes but also provides insights for improvement and enhancing maize yield.

1. Introduction

Post-transcriptional regulation of gene expression employs a wide range of RNA-binding proteins (RBPs) and plays crucial roles in finely controlling protein synthesis in a spatial and temporal manner. RBPs contribute to regulating RNA processing, mRNA transport, stability, and translation by targeting specific 3′-untranslated regions (UTRs) of target mRNA [1,2,3]. Additionally, RBPs can collaborate with ribosomal protein binding sites within the 5′-UTR or microRNAs (miRNAs) to regulate mRNA metabolism [4,5]. The Pumilio (Pum) RNA-binding proteins, known as Puf proteins, are a kind of RBP and are well characterized in animals and fungi [6,7], but little is known about the Pum family in plants.
The Pum proteins exhibit high conservation of the Pumilio Homology Domain (Pum-HD) in various organisms [8,9,10]. The Pum-HD possesses a unique crescent-shaped structure and is necessary for RNA binding [11,12,13]. Typically, this domain consists of imperfect tandem Puf repeats forming as α-helices each containing approximately 36 amino acids, and allows Pum proteins specifically to interact with mRNA to regulate post-transcription processes [11,14,15,16]. Within each Puf repeat, the second α-helix serves as the primary binding interface between the Pum and the target RNA [17]. The Puf repeat binds to a single RNA base through hydrogen bonds, van der Waals interactions, and base stacking. The binding of Pum and RNA can be facilitated by three conserved amino acid side chains within each repeat, which allows Pum proteins to selectively bind to specific mRNA sequences [8]. Generally, Pum proteins are recognized by the conserved UGUA core motif situated within the 3′-UTR of target mRNA [18,19]. Interestingly, it was also found that Pum proteins interact with other proteins to inhibit translation or the initiation of mRNA decay processes [17]. The human Pum interacts with NORAD (non-coding RNA activated by DNA damage) to preserve genomic stability [20]. The PUF5p forms a complex with Pop2p, a component of the deadenylase complex, to regulate mRNA decay in yeast [21]. Specifically, Pum proteins stimulate deadenylation and decapping to accelerate mRNA turnover and reduce translation efficiency [22].
Although Pum proteins show high conservation in the Pum-HD sequence, there is a high diversity of Pum gene members in plants. For example, there are 2, 2, 2, 6, 10, and 11 Pum genes in the genome of Drosophila melanogaster, human, mouse, Saccharomyces cerevisiae, Trypanosoma cruzi, and Caenorhabditis elegans, respectively [8,11,23,24]. On the contrary, there are 31, 26, 28, 26, 22, and 20 Pum genes in the genome of Arabidopsis lyrata, Arabidopsis thaliana, Malus domestica, Glycine max, Oryza sativa ssp. Indica, and Oryza sativa ssp. Japonica, respectively [17,24,25]. This suggests that Pum proteins are involved in a wide range of post-transcriptional/translational regulations to control growth and development and cope with environmental stresses in plants. However, little is known about the roles of Pum genes in plants. In Arabidopsis, AtPum5 is involved in cucumber mosaic virus (CMV) and abiotic stress response and negatively regulates salt and drought tolerance by binding to 3′-UTR of the abiotic stress-responsive genes containing the Pum RNA-binding motifs at the 3′-UTR [26,27]. AtPum9 binds to target transcripts to trigger mRNA degradation via Pum-HD at the C-terminal and interacts with DCP2 (the catalytic subunit of the decapping complex) to positively regulate heat stress and seed dormancy mediated by REDUCED DORMANCY5 encoding a PP2C phosphatase [28,29]. AtPum23, a nuclear-localized protein, is required for normal plant growth including leaf development and organ polarity, as well as being involved in salt response mediated by ABA signaling via regulating rRNA processing [30,31,32]. AtPum24 is an atypical Pum protein and reduces mRNA stability of the BTB/POZMATH (BPM) gene family by directly binding to their 3′-UTR to regulate plant development, seed maturation, and starch, protein, and oil biosynthesis [33,34]. To date, the Pum gene family has only been identified genome-wide in Arabidopsis and rice [17]. In addition, the roles of Pum genes in plants are largely unknown.
Maize is one of the most important crops and plays a crucial role in ensuring food and economic security [35]. The primary components including starch, protein, and oil are stored within maize seeds, which account for approximately 90% of the total dry seed weight. The content and composition of these components in maize kernels have a significant impact on their quality [36]. Hence, the maize kernel is a valuable resource for human consumption, animal feed, and bioenergy applications. In this study, we focused on comprehensively exploring the maize Pum gene family. We identified 19 ZmPum genes and investigated their physicochemical properties, phylogenetic relationships, chromosome localization, gene and protein-conserved domain structure, cis-acting elements, and expression patterns in kernel development. The main objective is to provide valuable insights into the underlying role of ZmPum genes in regulating seed development in maize and contribute to the analysis of the Pum family in plants.

2. Results

2.1. The ZmPum Family in Maize

In the maize genome, a total of 19 candidate genes encoding Pum proteins were identified and defined as ZmPum1 to ZmPum19. The coding sequence of the ZmPum gene was 1191 to 3009 bp in length, encoding 396 to 1002 amino acids (aa), with a molecular weight (MW) of 42.49 to 109.02 kDa. The isoelectric points (PIs) of ZmPum1, 2, 3, 5, 7, 13, 14, and 18 proteins were more than 7.00, the other eleven ZmPum proteins had PIs ranging from 5.70 to 6.73. The instable indices of ZmPum1, 6, 17, 18, and 19 were less than 40.00, while others were more than 40.00. All ZmPum proteins were hydrophilic proteins with a grand average of hydropathicity (GRAVY) < 0. Fourteen ZmPum proteins were predicted to show nuclear localization, only ZmPum1, 2, and 13 showed cytoplasm localization, ZmPum6 showed vacuole localization, and ZmPum18 showed chloroplast localization (Table 1). The diversity of properties of ZmPum proteins may imply their different roles.
The phylogenetic analysis showed that maize ZmPum members were clustered into four subclades within the phylogenetic tree based on the similarity of their amino acid sequences with AtPum proteins. ZmPum3, 4, 6, 8, 9, 11, 12, 13, 15, and 19 were grouped in subclade I. ZmPum1, 10, 16, 17, and 18 were clustered into clade II. Two (ZmPum7 and 14) and three (ZmPum2, 5, and 24) of them were grouped in clades III and IV, respectively (Figure 1).

2.2. Protein Architectures of ZmPum

As shown in Figure 2, six conserved motifs were discovered in ZmPum proteins and named motifs 1–6. The majority of ZmPum proteins contained motifs 1, 2, and 3, excluding ZmPum3 and ZmPum5. Among them, motifs 1, 2, 3, 4, and 5 contribute to the composition of Pum domains. Conserved domain analysis showed that the ZmPum proteins could be divided into two groups: typical Pum and atypical Pum. ZmPum3, 5, 7, and 14 belonged to one subgroup and were atypical Pum families because they possessed few Pum domains. The other 14 ZmPum members were grouped into another clade and were typical Pum proteins.
Three-dimensional structures of 19 ZmPum proteins were predicated using structure-based analysis in the Ensembl database. All ZmPum proteins contained the conserved Pum domain with a different number of α-helix repeats. The ZmPum proteins from the same subgroup exhibited a similar three-dimensional structure. However, the members from subclade I were atypical Pum proteins containing fewer imperfect Pum domains. The ZmPum members from group II, except ZmPum16, had eight Pum domains at the C-terminal region (Figure 2 and Figure 3).

2.3. Chromosomal Location and Gene Duplication of ZmPum

The information on chromosome location of ZmPum genes was obtained from the maizeGDB database and used for visualization of 19 ZmPum mapping to the maize genome (Figure 4). In detail, there was no ZmPum gene on chromosome 3. Other ZmPum genes were unevenly distributed on the other 9 maize chromosomes. There were five, three, three, two, and two ZmPum genes on chromosomes 1 (ZmPum1, 2, 3, 4, and 5), 4 (ZmPum8, 9, and 10), 10 (ZmPum17, 18, and 19), 2 (ZmPum6 and 7), and 7 (ZmPum13 and 14), respectively. The ZmPum11, 12, 15, and 16 genes were mapped on chromosomes 5, 6, 8, and 9, respectively. The results of gene duplication analysis showed that seven segmental duplication events were detected among 19 ZmPum genes, and each gene pair was located on a distinct chromosome, including pairs of ZmPum1 and 16, ZmPum2 and 16, ZmPum4 and 8, ZmPum6 and 19, ZmPum7 and 14, ZmPum9 and 11, as well as ZmPum12 and 15 (Figure 4).
Additionally, the synteny between the Pum gene families in the maize and rice, as well as maize and Arabidopsis genomes, was also examined. It was revealed that there were nineteen pairs of Pum orthologous genes in maize and rice, and five Pum gene pairs in maize and Arabidopsis (Figure 5; Table S1).

2.4. Gene Structure and Cis-Elements of ZmPum

To further examine the organization of the exons and introns of ZmPum genes, the CDS and the corresponding gDNA sequence of each ZmPum gene were analyzed using GSDS. It showed that the numbers of exons and introns varied greatly among different ZmPum members (Figure 6), which ranged from 5 to 10 exons unevenly. For example, ZmPum3, 5, 9, and 11 had 10 exons. ZmPum6, 12, and 15 had 9 exons, ZmPum4 and 19 both possessed 8 exons, and the other ZmPum genes had fewer exons.
Cis-elements analysis showed that abundant elements involved in hormone response were identified in promoter sequences of ZmPum genes, such as ABRE, AuxRRcore, P- box/GARE, TCA, and TGACG/CGTCA motif elements, which were responsive to ABA, auxin, gibberellin, SA, and MeJA, respectively (Figure 7; Table S2). In general, 13 ZmPum genes (68.4%) possessed auxin-responsive elements (AuxRRcore). Thirteen ZmPum genes had gibberellin-responsive elements (P-box or GARE). Moreover, ZmPum12 carried RY-elements involved in seed-specific regulation. Meanwhile, MBS, MBSI, and LTR elements involved in drought, flavonoid biosynthetic regulation, and low-temperature response, respectively, were found in their promoters. In total, 11 ZmPum genes had MBS elements, 10 ZmPum genes (52.6%) had LTR elements, and 2 ZmPum genes had MBSI elements.

2.5. Tissue-Specific Expression Patterns of ZmPum

The results of tissue-specific expression analysis showed that the ZmPum genes were classified into three groups in terms of different expression patterns (Figure 8). ZmPum3, 4, 5, 9, 11, 15, and 18 could be clustered into one group and highly expressed in pollinated internodes, embryos, endosperm, and whole seeds. ZmPum1, 2, 6, 7, 8, 17, and 19 were clustered into one group and expressed in bicellular male gametophytes, microspores, and sperm cells. Interestingly, ZmPum6 and 17 were also slightly expressed during seed development. In addition, ZmPum10, 12, 13, 14, and 16 were grouped into one branch and dominantly expressed in endosperm and whole seeds. Meanwhile, ZmPum12 and 14 exhibit high expression in pollinated internodes and slight expression in seminal, silk, roots, and bicellular male gametophytes. The results suggest that the ZmPum genes may play essential roles in regulating maize growth and development, particularly in the formation of seeds.

2.6. ZmPum Regulates Kernel Development

High-resolution transcriptome data in maize endosperms ranging from 48 to 144 h (En48-144) after pollination (HAP) with a time interval of 24 h were recently reported by Fu et al. [37]. Interestingly, except for the ZmPum2 gene, the expression of the other 18 ZmPum genes was detected during En48-144 HAP (Figure 9). ZmPum1, 4, 8, 9, 11, 16, 17, and 18 exhibited a high expression level in En48 HAP and a low transcript level in En72, En96, En120, and En144 HAP. While ZmPum3, 5, 6, 7, 10, 12, 13, 14, 15, and 19 showed high expression levels in En72, En96, En120, and En144 HAP, but low expression in En48 HAP.
Subsequently, qRT-PCR was performed and used to confirm the expression of ZmPum genes in maize kernel development at 15, 20, and 25 days after pollination (DAP). As shown in Figure 10, the expression of ZmPum2 was not detected in any samples. However, the ZmPum4, 7, and 13 genes showed extremely high expression levels in the kernel of 15, 20, and 25 DAP. The ZmPum3, 8, 9, 10, 11, 12, 14, 16, and 17 genes exhibited a high transcript level. Inversely, the expression of the ZmPum1, 5, 6, 18, and 19 genes in the kernel of 15, 20, and 25 DAP was lower than others.
The above results suggest that the ZmPum genes play a key role in kernel development and may have distinct roles and functions during endosperm development.

3. Discussion

To date, Pum proteins have been identified as a kind of RBP to regulate gene post-transcription via conserved Pum-HD [11,13,14,29]. In the present study, a total of 19 ZmPum members were identified in the maize genome (Table 1; Figure 1). Meanwhile, ZmPum were grouped into typical and atypical Pum with different numbers of Pum domains and showed the diversity of gene structures (Figure 2, Figure 3 and Figure 6), which was similar to the AtPum family [17,30]. However, the number of Pums in different organisms is variable and showed higher diversity in plants, which could be explained by whole-genome duplications [17]. It was also found that there were some paralogous and orthologous Pum gene pairs in the maize genome and between the rice and Arabidopsis genomes (Figure 4 and Figure 5). Tandem duplication and segmental duplication have played essential roles in expanding gene families during the species’ evolutionary history [38,39].
Plants evolved various mechanisms, including physiological, biochemical, and molecular changes, for survival under adverse conditions [40]. Herein, improvement of crop performance and yield under environmental stimuli is a crucial goal during sustainable agriculture production to ensure food security. Although it is not well known for Pum in plants, few available reports show that Pum can respond to stress such as heat, drought, salt, osmotic, ABA, dark, light, brassinolide, and glucose, as well as regulate development [26,27,28,32,34,41]. In Arabidopsis, AtPum1 to 6 are specifically associated with genes related to shoot stem cell maintenance genes [25]. AtPum5 regulates CMV infection and salt tolerance [26,27,42]. Similarly, rice Pum genes respond to biotic and abiotic stress including Magnaporthe oryzae and Nilaparvata lugens infections, cold, drought, auxin, and cytokinin [30]. However, the atpum23 mutant exhibited delayed germination rates compared to wild-type plants [32]. Reducing AtPum24 expression resulted in abnormal seed maturation, wrinkled seeds, and lower seed oil contents, but higher starch and sugar contents. Inversely, overexpression of AtPum24 increased seed fatty acid, size, and weight [33]. The findings suggest that Pum regulates seed development in plants.
In maize, the endosperm is the main nutritive tissue and accounts for approximately 90% of the total dry seed weight. Starch, protein, and oil are the primary storage components within maize endosperm [43]. The improvement of kernel traits holds significant importance for cultivating new maize germplasm with superior quality and high yield. Statistical analysis of transcriptome data from MaizeGDB showed that 84% (16/19) of ZmPum genes are highly expressed in seeds after pollination, except ZmPum2, 8, and 19 (Figure 8). In maize, the early endosperm development phase plays a key role in kernel development and comes to an end at 144 h HAP. Afterward, the endosperm shifts to rapid cell proliferation and differentiation, and enters the filling stage [38,44,45,46]. Here, it was also found that 18 ZmPum genes exhibited high transcription activity in En48, En72, En96, En120, or En144 HAP (Figure 9). Most ZmPum genes maintained high expression levels in maize seeds (Figure 10). These results indicate that the ZmPum family is involved in the regulation of early endosperm development.
Insights into the molecular interactions between Pum proteins and RNA bases have been well revealed in some eukaryotes but are still urgently needed to be explored in plants. Plants possess highly complex genomes with a high number of Pum members, which implies that Pum has specific target mRNA and function in plants [25]. In summary, we identified 19 ZmPum genes and found their involvement in kernel development in maize. In a further study, the function and molecular mechanism of ZmPum genes in regulating seed traits will be revealed. Overall, the study provides a valuable reference to improve crops through genetic engineering approaches.

4. Materials and Methods

4.1. Identification of ZmPum and Phylogenetic Analysis

To identify the ZmPum in maize, the database of amino acid of Zm-B73 V5.0 and the AtPum protein sequences were retrieved from MaizeGDB (https://maizegdb.org/, accessed on 5 April 2023) and TAIR (https://www.arabidopsis.org/accessed on 5 April 2023), respectively. The local BLASTP was conducted for ZmPum searching using the AtPum sequences as a reference. Additionally, the hidden Markov model (HMM) files of Pum-HD were acquired from the Pfam database (http://pfam.xfam.org/, accessed on 7 April 2023) and used to search for ZmPUM protein sequences [47]. The properties of the ZmPUM proteins such as molecular weight, isoelectric point, hydrophilic index, stability coefficient, and grand average of hydropathicity (GRAVY) were analyzed using the ExPASy tool (www.expasy.org/tools/, accessed on 7 April 2023) [48]. The subcellular localization was predicted using the BUSCA tool (http://busca.biocomp.unibo.it/, accessed on 7 April 2023) [49]. For the phylogenetic analysis, the amino acid sequences of Pum from maize and Arabidopsis were aligned and used for constructing a phylogenetic tree by Mega 7.0, employing 1000 bootstrap replicates to assess the reliability of the tree topology.

4.2. Conserved Motif, Domain, and Structures Analysis

The conserved motifs of ZmPum proteins were identified using the MEME online program (http://meme.sdsc.edu/meme/intro.html, accessed on 5 April 2023) [50]. To verify the conserved Pum-HD of ZmPum proteins, the protein sequences of ZmPum were analyzed using an online tool in the SMART database (https://smart.embl.de, accessed on 12 April 2023) and NCBI-CDD database (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 14 April 2023) [51,52]. Meanwhile, the 3D structures of ZmPum proteins were predicted by the Swiss-Model tool (https://swissmodel.expasy.org/interactive/, accessed on 14 April 2023) [53]. Subsequently, the quality of the predicted protein structures was evaluated using the SAVES server (http://nihserver.mbi.ucla.edu/SAVES/, accessed on 15 April 2023) [54].

4.3. Gene Structure and Duplication Analysis

For the examination of exon/intron structures, the coding sequences (CDS) and genomic DNA (gDNA) sequences of ZmPum genes were analyzed using the Gene Structure Display Server 2.0 (GSDS) (http://gsds.cbi.pku.edu.cn/, accessed on 18 April 2023) [55]. Chromosome localization of the ZmPum gene was obtained from maizeGDB and visualized using TBtools [56]. Gene duplications of the ZmPum genes were identified and examined by TBtools. The collinearity of the orthologous Pum genes between maize, rice, and Arabidopsis was determined and plotted using MCScanX Circos within TBtools, respectively.

4.4. Cis-Acting Elements and Expression Analysis of ZmPum

To determine the cis-elements in ZmPum gene promoters, the 2000 bp promoter sequences of ZmPum genes were acquired from the maizeGDB and analyzed using the PlantCARE software (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 20 April 2023) [57]. Furthermore, the TBtools software was employed to visualize the composition of cis-elements in promoters.
To examine the specific expression patterns of ZmPum genes in maize, expression data for maize’s different developmental stages and tissues and high-resolution transcriptome data from time points ranging from 48 to 144 h after pollination (HAP) were obtained from qTeller in MaizeGDB (https://qteller.maizegdb.org/, accessed on 22 April 2023) and Fu et al. [37], respectively. Then, the expression of ZmPum genes was analyzed and used to create a heatmap using TBtools.

4.5. RNA Extraction and Quantitative Real-Time PCR (qRT-PCR)

The maize seeds were sampled at 15, 20, and 25 DAP and used to extract total RNA using RNAiso plus kit (TaKaRa, Dalian, China) according to the manufacturer’s instructions. The concentration and purity of the RNA samples were determined using NanoDropTM OneC (ThermoScientific, Waltham, USA). The qualified RNA samples were then reverse-transcribed into cDNA using the PrimeScriptTM reagent kit (TaKaRa, Dalian, China) and used for qRT-PCR. The qRT-PCR was conducted using the TransScript® II Two-Step RT-PCR SuperMix (Transgen, Beijing, China) in the Bio-Rad CFX96TM Real-Time PCR system. The specific primers of each ZmPum gene were designed using Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_LOC=BlastHome, accessed on 25 May 2023), synthesized at Sangon Biotech (Chengdu, China), and listed in Table S3. In addition, a 139 bp fragment of maize ZmTUB gene was amplified using primers T-F/T-R (Table S3) and used as an internal reference for normalization. To determine the relative expression levels of the ZmPum genes, the 2−ΔΔCT method was employed [58]. This assay was conducted with three biological and technological replicates.

5. Conclusions

In summary, we identified 19 ZmPum genes and found their involvement in kernel development in maize. In the next study, the function and molecular mechanism of ZmPum genes in regulating seed traits will be revealed. Overall, the study provides a valuable reference to improve crops through genetic engineering approaches.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241814036/s1.

Author Contributions

Conceptualization, H.Y. and W.F.; methodology, W.F. and H.Z.; software, Y.C. and C.Y.; resources, M.H.B.K. and Q.Y.; writing—original draft preparation, W.F.; writing—review and editing, H.Y.; supervision, W.L., F.F. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (32102226), and the Sichuan Science and Technology Program (2022YFH0067).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are included in the article and Supplementary Files.

Acknowledgments

The authors are thankful for the technical support from the Key Laboratory of Biology and Genetic Improvement of Maize in the Southwest Region.

Conflicts of Interest

The authors declare that there are no conflict of interest.

References

  1. Lee, M.H.; Wu, X.; Zhu, Y. RNA-Binding Protein PUM2 Regulates Mesenchymal Stem Cell Fate via Repression of JAK2 and RUNX2 MRNAs. J. Cell. Physiol. 2020, 235, 3874–3885. [Google Scholar] [CrossRef]
  2. Yang, H.; Deng, Z.; Pan, X.; Shen, H.-B.; Choi, K.S.; Wang, L.; Wang, S.; Wu, J. RNA-Binding Protein Recognition Based on Multi-View Deep Feature and Multi-Label Learning. Brief. Bioinform. 2021, 22, bbaa174. [Google Scholar] [CrossRef]
  3. Goldstrohm, A.C.; Hall, T.M.T.; McKenney, K.M. Post-Transcriptional Regulatory Functions of Mammalian Pumilio Proteins. Trends Genet. 2018, 34, 972–990. [Google Scholar] [CrossRef]
  4. Qiu, C.; Dutcher, R.C.; Porter, D.F.; Arava, Y.; Wickens, M.; Hall, T.M.T. Distinct RNA-Binding Modules in a Single PUF Protein Cooperate to Determine RNA Specificity. Nucleic Acids Res. 2019, 47, 8770–8784. [Google Scholar] [CrossRef]
  5. Yan, Y.; Ham, B.K.; Chong, Y.H.; Yeh, S.D.; Lucas, W.J. A Plant SMALL RNA-BINDING PROTEIN 1 Family Mediates Cell-to-Cell Trafficking of RNAi Signals. Mol. Plant 2020, 13, 321–335. [Google Scholar] [CrossRef]
  6. Nishanth, M.J.; Simon, B. Functions, Mechanisms and Regulation of Pumilio/Puf Family RNA Binding Proteins: A Comprehensive Review. Mol. Biol. Rep. 2020, 47, 785–807. [Google Scholar] [CrossRef]
  7. Galgano, A.; Forrer, M.; Jaskiewicz, L.; Kanitz, A.; Zavolan, M.; Gerber, A.P. Comparative Analysis of MRNA Targets for Human PUF-Family Proteins Suggests Extensive Interaction with the MiRNA Regulatory System. PLoS ONE 2008, 3, e3164. [Google Scholar] [CrossRef]
  8. Wang, M.; Ogé, L.; Perez-Garcia, M.D.; Hamama, L.; Sakr, S. The PUF Protein Family: Overview on PUF RNA Targets, Biological Functions, and Post Transcriptional Regulation. Int. J. Mol. Sci. 2018, 19, 410. [Google Scholar] [CrossRef]
  9. Joshna, C.R.; Saha, P.; Atugala, D.; Chua, G.; Muench, D.G. Plant PUF RNA-Binding Proteins: A Wealth of Diversity for Post-Transcriptional Gene Regulation. Plant Sci. 2020, 297, 110505. [Google Scholar] [CrossRef]
  10. Spassov, D.S.; Jurecic, R. The PUF Family of RNA-Binding Proteins: Does Evolutionarily Conserved Structure Equal Conserved Function? IUBMB Life 2003, 55, 359–366. [Google Scholar] [CrossRef]
  11. Wang, X.; Mclachlan, J.; Zamore, P.D.; Tanaka Hall, T.M. Modular Recognition of RNA by a Human Pumilio-Homology Domain. Cell 2002, 110, 501–512. [Google Scholar] [CrossRef] [PubMed]
  12. Wharton, R.P.; Sonoda, J.; Lee, T.; Patterson, M.; Murata, Y.; Carolina, N. The Pumilio RNA-Binding Domain Is Also a Translational Regulator. Mol. Cell 1998, 1, 863–872. [Google Scholar] [CrossRef] [PubMed]
  13. Najdrová, V.; Stairs, C.W.; Vinopalová, M.; Voleman, L.; Doležal, P. The Evolution of the Puf Superfamily of Proteins across the Tree of Eukaryotes. BMC Biol. 2020, 18, 77. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, B.; Gallegos, M.; Puoti, A.; Durkin, E.; Fields, S.; Kimble, J.; Wickens, M.P. A Conserved RNA-Binding Protein That Regulates Sexual Fates in the C. Elegans Hermaphrodite Germ Line. Nature 1997, 390, 477–484. [Google Scholar] [CrossRef] [PubMed]
  15. Edwards, T.A.; Pyle, S.E.; Wharton, R.P.; Aggarwal, A.K. Structure of Pumilio Reveals Similarity between RNA and Peptide Binding Motifs. Cell 2001, 105, 281–289. [Google Scholar] [CrossRef]
  16. Wang, X.; Zamore, P.D.; Tanaka Hall, T.M. Crystal Structure of a Pumilio Homology Domain. Mol. Cell 2001, 7, 855–865. [Google Scholar] [CrossRef]
  17. Tam, P.P.C.; Barrette-Ng, I.H.; Simon, D.M.; Tam, M.W.C.; Ang, A.L.; Muench, D.G. The Puf Family of RNA-Binding Proteins in Plants: Phylogeny, Structural Modeling, Activity and Subcellular Localization. BMC Plant Biol. 2010, 10, 44. [Google Scholar] [CrossRef]
  18. Zamore, P.D.; Williamson, J.R.; Lehmann, R. The Pumilio Protein Binds RNA through a Conserved Domain That Defines a New Class of RNA-Binding Proteins. RNA 1997, 3, 1421–1433. [Google Scholar] [CrossRef]
  19. Gerber, A.P.; Herschlag, D.; Brown, P.O. Extensive Association of Functionally and Cytotopically Related MRNAs with Puf Family RNA-Binding Proteins in Yeast. PLoS Biol. 2004, 2, e79. [Google Scholar] [CrossRef]
  20. Lee, S.; Kopp, F.; Chang, T.C.; Sataluri, A.; Chen, B.; Sivakumar, S.; Yu, H.; Xie, Y.; Mendell, J.T. Noncoding RNA NORAD Regulates Genomic Stability by Sequestering PUMILIO Proteins. Cell 2016, 164, 69–80. [Google Scholar] [CrossRef]
  21. Goldstrohm, A.C.; Hook, B.A.; Seay, D.J.; Wickens, M. PUF Proteins Bind Pop2p to Regulate Messenger RNAs. Nat. Struct. Mol. Biol. 2006, 13, 533–539. [Google Scholar] [CrossRef] [PubMed]
  22. Miller, M.A.; Olivas, W.M. Roles of Puf Proteins in MRNA Degradation and Translation. Wiley Interdiscip. Rev. RNA 2011, 2, 471–492. [Google Scholar] [CrossRef] [PubMed]
  23. Miller, M.T.; Higgin, J.J.; Tanaka Hall, T.M. Basis of Altered RNA-Binding Specificity by PUF Proteins Revealed by Crystal Structures of Yeast Puf4p. Nat. Struct. Mol. Biol. 2008, 15, 397–402. [Google Scholar] [CrossRef] [PubMed]
  24. Huh, S.U. The Role of Pumilio Rna Binding Protein in Plants. Biomolecules 2021, 11, 1851. [Google Scholar] [CrossRef] [PubMed]
  25. Francischini, C.W.; Quaggio, R.B. Molecular Characterization of Arabidopsis Thaliana PUF Proteins-Binding Specificity and Target Candidates. FEBS J. 2009, 276, 5456–5470. [Google Scholar] [CrossRef]
  26. Un Huh, S.; Paek, K.-H. Role of Arabidopsis Pumilio RNA Binding Protein 5 in Virus Infection. Plant Signal. Behav. 2013, 8, e23975. [Google Scholar] [CrossRef]
  27. Huh, S.U.; Paek, K.H. APUM5, Encoding a Pumilio RNA Binding Protein, Negatively Regulates Abiotic Stress Responsive Gene Expression. BMC Plant Biol. 2014, 14, 75. [Google Scholar] [CrossRef]
  28. Nyikó, T.; Auber, A.; Bucher, E. Functional and Molecular Characterization of the Conserved Arabidopsis PUMILIO Protein, APUM9. Plant Mol. Biol. 2019, 100, 199–214. [Google Scholar] [CrossRef] [PubMed]
  29. Xiang, Y.; Nakabayashi, K.; Ding, J.; He, F.; Bentsink, L.; Soppe, W.J.J. REDUCED DORMANCY5 Encodes a Protein Phosphatase 2C That Is Required for Seed Dormancy in Arabidopsis. Plant Cell 2014, 26, 4362–4375. [Google Scholar] [CrossRef]
  30. Abbasi, N.; Kim, H.B.; Park, N., II; Kim, H.S.; Kim, Y.K.; Park, Y., II; Choi, S.B. APUM23, a Nucleolar Puf Domain Protein, Is Involved in Pre-Ribosomal RNA Processing and Normal Growth Patterning in Arabidopsis. Plant J. 2010, 64, 960–976. [Google Scholar] [CrossRef]
  31. Huang, T.; Kerstetter, R.A.; Irish, V.F. APUM23, a PUF Family Protein, Functions in Leaf Development and Organ Polarity in Arabidopsis. J. Exp. Bot. 2014, 65, 1181–1191. [Google Scholar] [CrossRef]
  32. Huang, K.C.; Lin, W.C.; Cheng, W.H. Salt Hypersensitive Mutant 9, a Nucleolar APUM23 Protein, Is Essential for Salt Sensitivity in Association with the ABA Signaling Pathway in Arabidopsis. BMC Plant Biol. 2018, 18, 40. [Google Scholar] [CrossRef]
  33. Huang, R.; Liu, M.; Gong, G.; Wu, P.; Patra, B.; Yuan, L.; Qin, H.; Wang, X.; Wang, G.; Liao, H.; et al. The Pumilio RNA-Binding Protein APUM24 Regulates Seed Maturation by Fine-Tuning the BPM-WRI1 Module in Arabidopsis. J. Integr. Plant Biol. 2021, 63, 1240–1259. [Google Scholar] [CrossRef] [PubMed]
  34. Shanmugam, T.; Abbasi, N.; Kim, H.S.; Kim, H.B.; Park, N., II; Park, G.T.; Oh, S.A.; Park, S.K.; Muench, D.G.; Choi, Y.; et al. An Arabidopsis Divergent Pumilio Protein, APUM24, Is Essential for Embryogenesis and Required for Faithful Pre-RNA Processing. Plant J. 2017, 92, 1092–1105. [Google Scholar] [CrossRef] [PubMed]
  35. Klopfenstein, T.J.; Erickson, G.E.; Berger, L.L. Maize Is a Critically Important Source of Food, Feed, Energy and Forage in the USA. Field Crop. Res. 2013, 153, 5–11. [Google Scholar] [CrossRef]
  36. García-Lara, S.; Serna-Saldivar, S.O. Corn History and Culture. Corn Chem. Technol. 2018, 1–18. [Google Scholar] [CrossRef]
  37. Fu, Y.; Li, S.; Xu, L.; Ji, C.; Xiao, Q.; Shi, D.; Wang, G.; Wang, W.; Wang, J.; Wang, J.; et al. RNA Sequencing of Cleanly Isolated Early Endosperms Reveals Coenocyte-to-Cellularization Transition Features in Maize. Seed Biol. 2023, 2, 8. [Google Scholar] [CrossRef]
  38. Lu, Y.; Liu, L.; Wang, X.; Han, Z.; Ouyang, B.; Zhang, J.; Li, H. Genome-Wide Identification and Expression Analysis of the Expansin Gene Family in Tomato. Mol. Genet. Genomics 2016, 291, 597–608. [Google Scholar] [CrossRef]
  39. Zhang, C.; Liu, S.; Liu, D.; Guo, F.; Yang, Y.; Dong, T.; Zhang, Y.; Ma, C.; Tang, Z.; Li, F.; et al. Genome-Wide Survey and Expression Analysis of GRAS Transcription Factor Family in Sweetpotato Provides Insights into Their Potential Roles in Stress Response. BMC Plant Biol. 2022, 22, 232. [Google Scholar] [CrossRef]
  40. Zhu, J.K. Abiotic Stress Signaling and Responses in Plants. Cell 2016, 167, 313–324. [Google Scholar] [CrossRef]
  41. Park, S.H.; Kim, H.S.; Kalita, P.J.; Choi, S.B. Structural and Functional Similarities and Differences in Nucleolar Pumilio RNA-Binding Proteins between Arabidopsis and the Charophyte Chara Corallina. BMC Plant Biol. 2020, 20, 230. [Google Scholar] [CrossRef] [PubMed]
  42. Huh, S.U.; Kim, M.J.; Paek, K.H. Arabidopsis Pumilio Protein APUM5 Suppresses Cucumber Mosaic Virus Infection via Direct Binding of Viral RNAs. Proc. Natl. Acad. Sci. USA 2013, 110, 779–784. [Google Scholar] [CrossRef] [PubMed]
  43. Yang, T.; Guo, L.; Ji, C.; Wang, H.; Wang, J.; Zheng, X.; Xiao, Q.; Wu, Y. The B3 Domain-Containing Transcription Factor ZmABI19 Coordinates Expression of Key Factors Required for Maize Seed Development and Grain Filling. Plant Cell 2021, 33, 104–128. [Google Scholar] [CrossRef] [PubMed]
  44. Lopes, M.A.; Larkins, B.A. Endosperm origin, development, and function. Plant Cell 1993, 5, 1383–1399. [Google Scholar] [CrossRef]
  45. Olsen, O.A. Endosperm Development: Cellularization and Cell Fate Specification. Annu. Rev. Plant Biol. 2001, 52, 233–267. [Google Scholar] [CrossRef]
  46. Sabelli, P.A.; Larkins, B.A. The Development of Endosperm in Grasses. Plant Physiol. 2009, 149, 14–26. [Google Scholar] [CrossRef]
  47. Mistry, J.; Chuguransky, S.; Williams, L.; Qureshi, M.; Salazar, G.A.; Sonnhammer, E.L.L.; Tosatto, S.C.E.; Paladin, L.; Raj, S.; Richardson, L.J.; et al. Pfam: The Protein Families Database in 2021. Nucleic Acids Res. 2021, 49, D412–D419. [Google Scholar] [CrossRef]
  48. Artimo, P.; Jonnalagedda, M.; Arnold, K.; Baratin, D.; Csardi, G.; De Castro, E.; Duvaud, S.; Flegel, V.; Fortier, A.; Gasteiger, E.; et al. ExPASy: SIB Bioinformatics Resource Portal. Nucleic Acids Res. 2012, 40, 597–603. [Google Scholar] [CrossRef]
  49. Savojardo, C.; Martelli, P.L.; Fariselli, P.; Profiti, G.; Casadio, R. BUSCA: An Integrative Web Server to Predict Subcellular Localization of Proteins. Nucleic Acids Res. 2018, 46, W459–W466. [Google Scholar] [CrossRef]
  50. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME Suite: Tools for Motif Discovery and Searching. Nucleic Acids Res. 2009, 37, 202–208. [Google Scholar] [CrossRef]
  51. Letunic, I.; Khedkar, S.; Bork, P. SMART: Recent Updates, New Developments and Status in 2020. Nucleic Acids Res. 2021, 49, D458–D460. [Google Scholar] [CrossRef] [PubMed]
  52. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional Classification of Proteins via Subfamily Domain Architectures. Nucleic Acids Res. 2017, 45, D200–D203. [Google Scholar] [CrossRef] [PubMed]
  53. Arnold, K.; Bordoli, L.; Kopp, J.; Schwede, T. The SWISS-MODEL Workspace: A Web-Based Environment for Protein Structure Homology Modelling. Bioinformatics 2006, 22, 195–201. [Google Scholar] [CrossRef]
  54. Luthy, R.; Bowei, J.; Einsenberg, D. Verify3D: Assessment of Protein Models with Three-Dimensional Profiles. Methods Enzymol. 1997, 277, 396–404. [Google Scholar] [CrossRef]
  55. Hu, B.; Jin, J.; Guo, A.Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An Upgraded Gene Feature Visualization Server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef]
  56. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  57. Rombauts, S.; Déhais, P.; Van Montagu, M.; Rouzé, P. PlantCARE, a Plant Cis-Acting Regulatory Element Database. Nucleic Acids Res. 1999, 27, 295–296. [Google Scholar] [CrossRef]
  58. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic tree of ZmPum and AtPum proteins. The phylogenetic tree was constructed using MEGA7 software. The full-length amino acid sequences of 19 ZmPum and 26 AtPum proteins were aligned and used to construct a tree with 1000 bootstrap replicates to support the branching patterns.
Figure 1. Phylogenetic tree of ZmPum and AtPum proteins. The phylogenetic tree was constructed using MEGA7 software. The full-length amino acid sequences of 19 ZmPum and 26 AtPum proteins were aligned and used to construct a tree with 1000 bootstrap replicates to support the branching patterns.
Ijms 24 14036 g001
Figure 2. Diagram of motifs and domains of the ZmPum proteins.
Figure 2. Diagram of motifs and domains of the ZmPum proteins.
Ijms 24 14036 g002
Figure 3. The structure modeling of ZmPum proteins.
Figure 3. The structure modeling of ZmPum proteins.
Ijms 24 14036 g003
Figure 4. The location and segmental duplication of ZmPum members. The outer circle with chromosome numbers illustrates different maize chromosomes. The maize gene density is displayed by short-green lines in the inner circle. The gray lines indicate all segmental duplications in the maize genome, and the colored lines indicate segmentally duplicated ZmPum gene pairs.
Figure 4. The location and segmental duplication of ZmPum members. The outer circle with chromosome numbers illustrates different maize chromosomes. The maize gene density is displayed by short-green lines in the inner circle. The gray lines indicate all segmental duplications in the maize genome, and the colored lines indicate segmentally duplicated ZmPum gene pairs.
Ijms 24 14036 g004
Figure 5. Synteny analysis of Pum genes. Gray lines: all collinear blocks within maize, rice, and Arabidopsis genomes. Red lines: the synteny of Pum gene pairs. The species names with the prefixes Zm, Os, and At indicate maize, rice, and Arabidopsis, respectively.
Figure 5. Synteny analysis of Pum genes. Gray lines: all collinear blocks within maize, rice, and Arabidopsis genomes. Red lines: the synteny of Pum gene pairs. The species names with the prefixes Zm, Os, and At indicate maize, rice, and Arabidopsis, respectively.
Ijms 24 14036 g005
Figure 6. Exon and intron composition of ZmPum genes. The CDS and gDNA sequences of ZmPum genes were retrieved from maizeGDB and analyzed using the Gene Structure Display Server 2.0 (GSDS).
Figure 6. Exon and intron composition of ZmPum genes. The CDS and gDNA sequences of ZmPum genes were retrieved from maizeGDB and analyzed using the Gene Structure Display Server 2.0 (GSDS).
Ijms 24 14036 g006
Figure 7. Predicted cis-regulatory elements in ZmPum promoters. The promoter sequence of each ZmPum gene was analyzed by PlantCARE.
Figure 7. Predicted cis-regulatory elements in ZmPum promoters. The promoter sequence of each ZmPum gene was analyzed by PlantCARE.
Ijms 24 14036 g007
Figure 8. Expression profiles of the ZmPum genes in different tissues. The colored scale represents expression data.
Figure 8. Expression profiles of the ZmPum genes in different tissues. The colored scale represents expression data.
Ijms 24 14036 g008
Figure 9. Heatmap of ZmPum expression after pollination in the endosperm. En48, En72, En96, En120, and En144 indicate the endosperm at 48, 72, 96, 120, and 144 h after pollination, respectively.
Figure 9. Heatmap of ZmPum expression after pollination in the endosperm. En48, En72, En96, En120, and En144 indicate the endosperm at 48, 72, 96, 120, and 144 h after pollination, respectively.
Ijms 24 14036 g009
Figure 10. Relative expression level of ZmPum genes in maize kernel development. Here, 15, 20, and 25 DAP indicate the seed at 15, 20, and 25 days after pollination, respectively.
Figure 10. Relative expression level of ZmPum genes in maize kernel development. Here, 15, 20, and 25 DAP indicate the seed at 15, 20, and 25 days after pollination, respectively.
Ijms 24 14036 g010
Table 1. The ZmPum members in maize.
Table 1. The ZmPum members in maize.
Gene NameIDLength (bp) aLength (aa) bMW (kDa)PIIIGRAVYLocal c
ZmPum1Zm00001eb006950208569476.228.1239.11−0.219C
ZmPum2Zm00001eb008320119139644.958.0453.99−0.002C
ZmPum3Zm00001eb028000251783892.829.3242.50−0.548N
ZmPum4Zm00001eb0357002985994108.976.0247.08−0.470N
ZmPum5Zm00001eb063600114938242.499.3742.96−0.245N
ZmPum6Zm00001eb087930232877586.485.7039.96−0.420V
ZmPum7Zm00001eb098710141947253.147.6366.34−0.202N
ZmPum8Zm00001eb1746302961986108.056.0049.57−0.476N
ZmPum9Zm00001eb19129030091002109.026.1650.92−0.450N
ZmPum10Zm00001eb200460226575483.936.4340.51−0.479N
ZmPum11Zm00001eb25905030031000108.946.3650.17−0.459N
ZmPum12Zm00001eb266980261687194.926.1254.14−0.381N
ZmPum13Zm00001eb301130161453758.778.5849.65−0.165C
ZmPum14Zm00001eb311340144348054.238.1662.40−0.248N
ZmPum15Zm00001eb355970262287395.186.0852.47−0.390N
ZmPum16Zm00001eb401310203167674.146.7346.54−0.182N
ZmPum17Zm00001eb408690205868576.046.1536.44−0.336N
ZmPum18Zm00001eb409750144047952.638.3331.04−0.037Ch
ZmPum19Zm00001eb419690234378087.275.9239.25−0.403N
a The CDS length of each ZmPum gene. b The amino acid length of each ZmPum protein. c Subcellular localization of each ZmPum protein. MW: molecular weight. PI: isoelectric point. II: instable index. GRAVY: grand average of hydropathicity. C: cytoplast. N: nuclear. V: vacuole. Ch: chloroplast.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Feng, W.; Zhang, H.; Cao, Y.; Yang, C.; Khalid, M.H.B.; Yang, Q.; Li, W.; Wang, Y.; Fu, F.; Yu, H. Comprehensive Identification of the Pum Gene Family and Its Involvement in Kernel Development in Maize. Int. J. Mol. Sci. 2023, 24, 14036. https://doi.org/10.3390/ijms241814036

AMA Style

Feng W, Zhang H, Cao Y, Yang C, Khalid MHB, Yang Q, Li W, Wang Y, Fu F, Yu H. Comprehensive Identification of the Pum Gene Family and Its Involvement in Kernel Development in Maize. International Journal of Molecular Sciences. 2023; 24(18):14036. https://doi.org/10.3390/ijms241814036

Chicago/Turabian Style

Feng, Wenqi, Hongwanjun Zhang, Yang Cao, Cheng Yang, Muhammad Hayder Bin Khalid, Qingqing Yang, Wanchen Li, Yingge Wang, Fengling Fu, and Haoqiang Yu. 2023. "Comprehensive Identification of the Pum Gene Family and Its Involvement in Kernel Development in Maize" International Journal of Molecular Sciences 24, no. 18: 14036. https://doi.org/10.3390/ijms241814036

APA Style

Feng, W., Zhang, H., Cao, Y., Yang, C., Khalid, M. H. B., Yang, Q., Li, W., Wang, Y., Fu, F., & Yu, H. (2023). Comprehensive Identification of the Pum Gene Family and Its Involvement in Kernel Development in Maize. International Journal of Molecular Sciences, 24(18), 14036. https://doi.org/10.3390/ijms241814036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop