Yanping Wang
Qiang Xu1- 1National Engineering Laboratory of Crop Stress Resistance Breeding, Anhui Agricultural University, Hefei, China
- 2Yunnan Key Laboratory of Plant Reproductive Adaptation and Evolutionary Ecology and Institute of Biodiversity, School of Ecology and Environmental Science, Yunnan University, Kunming, China
14-3-3 proteins (regulatory protein family) are phosphate serine-binding proteins. A number of transcription factors and signaling proteins have been shown to bind to the 14-3-3 protein in plants, which plays a role in regulating their growth (seed dormancy, cell elongation and division, vegetative and reproduction growth and stress response (salt stress, drought stress, cold stress). Therefore, the 14-3-3 genes are crucial in controlling how plants respond to stress and develop. However, little is known about the function of 14-3-3 gene families in gramineae. In this study, 49 14-3-3 genes were identified from four gramineae, including maize, rice, sorghum and brachypodium, and their phylogeny, structure, collinearity and expression patterns of these genes were systematically analyzed. Genome synchronization analysis showed large-scale replication events of 14-3-3 genes in these gramineae plants. Moreover, gene expression revealed that the 14-3-3 genes respond to biotic and abiotic stresses differently in different tissues. Upon arbuscular mycorrhizal (AM) symbiosis, the expression level of 14-3-3 genes in maize significantly increased, suggesting the important role of 14-3-3 genes in maize-AM symbiosis. Our results provide a better understanding on the occurrence of 14-3-3 genes in Gramineae plants, and several important candidate genes were found for futher study on AMF symbiotic regulation in maize.
Introduction
The 14-3-3 protein was originally isolated from bovine brain tissues by Moore and Perez (Carlson and Hubbard, 1968). The 14-3-3 family is composed of several subunits, exists in all eukaryotes, such as yeast (Heusden et al., 2010), human (Boston and Jackson, 1980) and Arabidopsis thaliana (Rosenquist et al., 2001). In Arabidopsis thaliana, they major components of the protein/G-box complex, and they were named “G-box factor 14-3-3” or “GF14” (Lu et al., 1992). 14-3-3 proteins (regulatory protein family) are phosphate serine-binding proteins that regulate the activity of a series of targets through direct protein-protein interactions (Bridges and Moorhead, 2005). A number of transcription factors and signaling proteins have been shown to bind to the 14-3-3 protein in plants, which plays a role in regulating their growth [seed dormancy (Darling and Yingling, 2005), cell elongation and division (Zhang et al., 2010), vegetative and reproduction growth (Asih et al., 2009; John, 2012)] and stress response [salt stress (Wei et al., 2009), drought stress (Chevalier et al., 2009), cold stress (Yao et al., 2007)].
Arbuscular mycorrhizal (AM) fungi belong to the ancient spherical phylum and are soil-borne microorganisms. They can establish the most extensive symbiotic relationship with around 80% of terrestrial flowering plant species, i.e. AM symbiosis (Simon et al., 1993; Remy et al., 1994). AM fungi consume plant photosynthates (Bago et al., 2000) and lipids for the duration of their life cycle (Bravo et al., 2017; Jiang et al., 2017). In return, AM fungi not only enhance the absorption of mineral nutrients and water for host plants, but also protect them from fungal pathogens (Jung et al., 2012; Chitarra et al., 2016) and all kinds of abiotic stresses (for example, heavy metals, low available nutrients, drought, extreme temperatures, acidic soils, salinity, aluminum toxins and contaminants) (Augé, 2001; Schützendübel and Polle, 2002; Lenoir et al., 2016). Therefore, AM fungi are the key endophytes of plant symbiosis. Symbiosis has an important impact on plant productivity and ecosystem functions (Van der Heijden and Klironomos, 1998), and is critical to agricultural sustainability (Gianinazzi et al., 2010).
At present, the AM fungal induced functions of 14-3-3 gene are mainly related to drought resistance. For example, Gi14-3-3 gene is involved in the protection of AM symbiosis against drought stress in host plants (Porcel et al., 2006). AM symbiosis decreases the expression of two plasma membrane aquaporin genes (PIP) activated by ABA signal, which reduces cell membrane permeability and improves cell water retention (Porcel et al., 2006). D-inositol-3-phosphate synthase (IPS) and 14-3-3 protein GF14 (14-3GF) are also regulated by AM fungi to activate cross-talk between symbiotic partners and increase the drought resilience of maize plants (Li et al., 2016). AM symbiosis may regulate stomatal behavior and maintain water use efficiency by regulating 14-3-3 genes in ABA signal pathway, thus improving drought resistance of plants (Xu et al., 2018). Moreover, it has also been reported that 14-3-3 genes can promote symbiosis. For example, three 14-3-3 protein genes have been identified in AMF: Fm201, Ri14-3-3, and RiBMH2. their transcriptional levels are highly induced in the early stage of symbiosis and symbiotic stage, including germinating spores, stem hyphae and roots (Sun et al., 2018).
Therefore, understanding the evolutionary relationship of 14-3-3 of Gramineae and the induction and regulation of 14-3-3 genes by AM symbiosis is of great significance to improve the response of crops to biotic and abiotic stresses
Materials and methods
Identification of the 14-3-3 proteins in maize, sorghum, rice, and brachypodium
The genome and predicted protein sequences for maize (Zea mays), rice (Oryza sativa-IRGSP-1.0), sorghum (Sorghum bicolor), and brachypodium (Brachypodium distachyon-v3.0) were obtained from the respective databases: MaizeGDB (http://www.maizegdb.org/, accessed on 20 October 2022), Phytozome (http://www.phytozome.net/ accessed on 20 October 2022), and Ensembl Plants (http://plants.ensembl.org, accessed on 20 October 2022). The version number for four plant databases are Zea mays-v4, Oryza sativa Japonica Group-IRGSP-1.0, Sorghum bicolor-v3, Brachypodium distachyon-v3. Two techniques were used to identify the gene encoding 14-3-3 from four plants at the whole genome level. Initially, functionally known 14-3-3 protein sequences from Arabidopsis (AtGRF1) and Oryza sativa (OsGF14a) were employed to query homologous proteins ofmaize, rice, sorghum, and brachypodium using BLASTp with E-values of less than le-5 (Huang et al., 2015). Then, Pfam scan (http://www.ebi.ac.uk/Tools/pfa/pfamscan/, accessed on 20 October 2022) was used to validate the existence of a 14-3-3-related domain (PF00244.23) in all 14-3-3 protein sequences.
Protein sequence analysis, chromosomal locations, gene structure and protein motif identification analysis
The theoretical isoelectric point (pI) and molecular weight (MW, kDa) were calculated using the ProtParam tool in ExPASy (http://web.expasy.org/protparam/, accessed on 22 October 2022) (Elisabeth et al., 2003)). Each protein’s subcellular localization was predicted using the Cell-PLoc 2.0 server (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/, accessed on 22 October 2022). The ensemblplant Database (http://plants.ensembl.org/index.html, accessed on 22 October 2022) was used to find the chromosomal locations of the 14-3-3 genes. The exon/intron structures of the 14-3-3 genes were determined using the Gene Structure Display Server (GSDS 2.0: http://gsds.cbi.pku.edu.cn/ (Hu et al., 2014), accessed on 22 October 2022) The Multiple Expectation Maximization for Motif Elicitation program (MEME, http://alternate.meme-suite.org, accessed on 22 October 2022) was used to identify the conserved motif of 14-3-3 proteins in four gramineae, the parameter was set as any number of repetitions and the maximum value of motif was 15 to identify the conserved motif of 14-3-3 protein in four gramineae plants (Bailey et al., 2006).
Duplication analysis
Evolutionary selection is measured by the ratio of non-synonym to synonymous substitution (Ka/Ks). In the TBtools software (v1.098726) (Chen et al., 2020), the genome databases and gene annotation files of four species were analyzed by using the MCScanX function, and the homologous gene pairs were brought into the Ka/Ks function for calculation, thus the evolutionary relationship was obtained. Homologous gene pairs can also be used for collinear analysis to evaluate the assembly effect of genomes and the retention and loss of homologous genes by homologous comparison, and to study the evolutionary relationship of materials. We integrated GFF file format files, gene link format files and Chr Layout format files using File Merge for MCScanX package on TBtools (v1.098726). Then the synteny relationships of 14-3-3 genes in four species were constructed.
Gene expression profile
Gramene database (http://www.gramene.org/, accessed on 1 November 2022) was utilized to gather tissue expression data. We selected the expression data of the required species from the website and downloaded it to search for the required 14-3-3 gene family genes. Then we selected the heatmap function in Tbtools (v1.098726) and imported the organized expression pattern data into it and click star. The 14-3-3 gene expression data of four gramineous plants were investigated and heat maps were drawn.
Plant materials, growth conditions, and treatments
Maize B73 seeds were germinated in hydroponic culture under sterile conditions. After germination, we moved seedlings of maize B73 to a sand-based pot culture to inoculate them with arbuscular mycorrhizal fungus. The control group was not inoculated with arbuscular mycorrhizal fungi. Culture medium was mixed with vermiculite:perlite at a 4:1 ratio and sterilised by 40 min high-pressure steam at 121℃. The AMF species was Glomus intraradices (Gi, provided by Sun Yat-Sen University). We grew maize at 28°C with 16 hours of light and 8 hours of darkness in a greenhouse. After 60 days post-treatment, the plants were collected. For RNA isolation, samples were stored at 80°C.
TransZol Up Plus RNA Kit was used to extract total RNA (TransGen Biotech, Beijing, China). cDNAs were obtained using reverse transcriptase (Vazyme, Nanjing, China). ZmTubulin (Zm00001d009780) and ZmGAPDH (Zm00001d049641) were used as the endogenous reference genes to standardize the relative expression levels of ZmGRF genes. Our previous published work provided extensive knowledge regarding qRT-PCR reactions (Wang et al., 2022). The 10-(ΔCt/3)method was used to determine the relative gene expression levels (Li et al., 2005; Liao et al., 2014). The qPCR assays were performed with three biological replicates.
Construction of regulatory network of symbiosis-related ZmGRFs
According to the qRT-PCR results, the co-expression network of 14-3-3 genes induced by AMF with significant differences was constructed. First, the ZmGRFs co-expression data was downloaded from PlantRegMap (http://plantregmap.gao-lab.org/, accessed on 2 November 2022), a plant transcriptional regulation map database. Secondly, we measured the RNAseq data of maize symbiosis [PRJNA918441(SAMN32395817-SAMN32395900)]. Co-expression networks were created using the WGCNA package (v1.63) in R with a soft threshold of 18 (Langfelder and Horvath, 2008). Expression modules were obtained using the automatic network construction function block wise Modules with default settings. The modules were grouped using a stringency threshold of 0.75. The two regulatory networks were intersected to obtain the primary expression regulatory network. The co-expression network was visualized using the Cytoscape program (v3.6.1) (Shannon et al., 2003).
Results
Identification of 14-3-3 gene family in maize, sorghum, rice, and brachypodium
In the current study, 49 genes were identified as members of 14-3-3 family in the genome of four plants. Maize 14-3-3 gene family was named according to the reported naming methods of Arabidopsis 14-3-3 gene family. This study found 28 ZmGRF genes in the maize genome, 6 in the sorghum genome, 7 in the brachypodium genome, and 8 in the rice genome. The number of 14-3-3 genes in maize is about four times that of other three species. It is speculated that corn needs 14-3-3 to play different functions to adapt to the environment. These 14-3-3s were designated as ZmGRF, SbGF14, OsGF14, and BdGF14 followed by number or alphabet. The detailed information of 14-3-3 genes were listed in Table 1, Tables S1–S3, including chromosome distribution, gene length, isoelectric point, and molecular weight. The number of amino acid residues in the four plants 14-3-3 proteins ranged from 94 (ZmGRF27) to 501 (ZmGRF3), and their relative molecular weights (MWs) ranged from 10.87 kDa (ZmGRF27) to 53.42 kDa (ZmGRF3). Theoretical pI (isoelectric point) values ranged from 4.71 (OsGF14e 、 OsGF14g 、 BdGF14b) to 9.69 (ZmGRF26) and rice, sorghum and brachypodium were all less than 7. Subcellular localization prediction of 14-3-3 gene family showed that only ZmGRF4 was located in Cytoplasm and ZmGRF27、ZmGRF28 、OsGF14g、 OsGF14h、 SbGF14f、 BdGF14g were located in cytoplasm and Nucleus, while the other 42 14-3-3 proteins were all located in Nucleus. The chromosomal location maps revealed that ZmGRFs were found on nine of ten chromosomes, whereas SbGF14s were found on three of ten chromosomes (Chr. 5,6, and 7), and OsGF14s were mainly existed on six of twelve chromosomes (Chr.1, 2, 3, 4, 8, and 11), and BdGF14s were dispersed on four of five chromosomes (Chr.1, 3, 4, and 5).
Phylogenetic and structural analysis of 14-3-3s in maize, sorghum, rice, and Brachypodium
14-3-3 genes from maize, sorghum, rice, and brachypodium were studied for evolutionary relationships. Firstly, the 14-3-3 proteins of four species and the reported 14-3-3 proteins of Arabidopsis thaliana were analyzed by phylogenetic tree. According to the group of Arabidopsis thaliana, the 14-3-3s of five plants were divided into ϵ and non-ϵ group (Figure 1). The phylogenetic differentiation patterns of the five plants were similar, but the four gramineae were more closely related than arabidopsis, reflecting different 14-3-3 conserved evolution and function. Among the four gramineae species, ϵ-group contains 8 members, accounting for 16.3% of the total 14-3-3s. Non ϵ-group contains 41 members, accounting for 81.7% of the total 14-3-3s. In ϵ group, ZmGRF25 has a distant genetic relationship with other genes, which may be caused by evolutionary selection to adapt to the environment. In non- ϵ classes, it can be seen from the evolutionary tree that there are three branches. Each branch can be regarded as an ancestor’s 14-3-3 gene. Through replication and species differentiation, each branch represents gene pedigree-specific replication and loss. This is of instructive significance for us to judge the evolution and origin of 14-3-3 gene family in Gramineae.
Figure 1 Phylogenetic relationship among 14-3-3 proteins in plants. The 14-3-3 protein sequences of Arabidopsis thaliana, maize, rice, sorghum and Brachypodium were mapped using Multiple Sequence Alignment by CLUSTALW. The 14-3-3 proteins in the phylogenetic tree were divided into two groups (ϵ group and non-ϵ classes). At, Arabidopsis thaliana; Zm, Zea mays; Os, Oryza sativa; Sb, Sorghum bicolor; Bd, Brachypodium distachyon.
To further understand the structural diversity of 14-3-3 genes and the diversity of protein structure, The exon/intron structure of 14-3-3 genes and the conserved motif of protein sequences were compared (Figure 2). Highly conserved distribution of CDS regions (Almost all the 14-3-3s contained 3 to 5 in numbers, except ZmGRF7, ZmGRF14, ZmGRF15, ZmGRF27, ZmGRF28, OsGF14h, SbGF14f and BdGF14g) was found among 14-3-3 genes. To further understand the structural variety of 14-3-3 proteins, online MEME/MAST tools were used to analyze the conserved motifs. The MEME program was used to predict fifteen putative protein motifs. In non - ϵ group and ϵ group, the motifs were not regular, but they all contained at least one of the first 8 motifs. In ϵ group, except ZmGRF26, other seven genes all contained motif7 and motif8. In non - ϵ group, all 14-3-3 proteins contained motif1, motif2, Motif3 and Motif5 except ZmGRF2, ZmGRF3, ZmGRF4, ZmGRF8, ZmGRF15 and ZmGRF16. Motif sequences were shown in Table S4.
Figure 2 Phylogenetic tree, motif, and gene structure of maize, sorghum, brachypodium, and rice. (A). The phylogenetic tree of four plants. non - ϵ group and ϵ group indicate the two 14-3-3 subfamilies. 14-3-3s from maize (ZmGRF1 to ZmGRF28), sorghum (SbGF14a to SbGF14f), brachypodium (BdGF14a to Bd GF14g), and rice (OsGF14a to OsGF14h). (B). Gene structure of four plants. Exons are indicated by pink boxes, introns are indicated by single lines, and gray boxes represent untranslated regions (UTR); (C). Motifs’ analysis. Different colors represent different motifs.
Duplication events of 14-3-3 genes among maize, sorghum, rice, and brachypodium
According to the whole-genome analysis in the four gramineae plants, maize had the most putative duplicated gene pairs (ZmGRF17/ZmGRF18, ZmGRF21/ZmGRF22, and ZmGRF23/ZmGRF24), while rice contained two pairs (OsGF14b/OsGF14e, OsGF14c/OsGF14e), brachypodium contained one pair (BdGF14b/BdGF14e), and sorghum had no gene pairs (Figure 3).
Figure 3 Synteny of four plants 14-3-3 genes. Zm1–10 represented maize chromosome 1–10, Sb1–10 represented sorghum chromosome, Bd1–5 represented brachypodium chromosome and Os1–12 represented rice chromosome, all chromosome indicated by colored boxes. Chromosome box numbers represent sequence lengths in megabases. All the syntenic genes were located in the map.
For better comprehending of these genes, we explored the gene duplication events of them in four Gramineae plants. To better understand the evolution history of the 14-3-3 gene family, putative orthologous links between all the 14-3-3 genes in maize and other three plants were established (Figure 4, Table 2).
Figure 4 Circle plot showing segmentally duplicated 14-3-3 genes on maize, sorghum, brachypodium, and rice chromosomes. Zm1–10 represented maize chromosome 1–10, Sb1–10 represented sorghum chromosome, Bd1–5 represented brachypodium chromosome and Os1–12 represented rice chromosome, all chromosome indicated by colored boxes. Chromosome box numbers represent sequence lengths in megabases. lines indicated duplicated 14-3-3 gene pairs.
We identified 6 homologous gene pairs between maize and rice (ZmGRF5 and OsGF14d;ZmGRF21 and OsGF14b; ZmGRF22 and OsGF14c; ZmGRF23 and OsGF14b; ZmGRF24 and OsGF14e; ZmGRF28 and OsGF14h), and 8 homologous gene pairs between maize and sorghum (ZmGRF1 and SbGF14a; ZmGRF5 and SbGF14d ; ZmGRF17 and SbGF14c; ZmGRF18 and SbGF14c; ZmGRF22 and SbGF14e; ZmGRF23 and SbGF14b; ZmGRF24 and SbGF14b; ZmGRF28 and SbGF14f), 8 homologous gene pairs between maize and brachypodium (ZmGRF1 and BdGF14d ; ZmGRF5 and BdGF14f ; ZmGRF13 and BdGF14f; ZmGRF21 and BdGF14c;ZmGRF22 and BdGF14c; ZmGRF23 and BdGF14b ;ZmGRF24 and BdGF14e; ZmGRF28 and BdGF14g) (Figure 4). The difference number of orthologous gene pairs between maize and the other three plants may be due to different loss rates of duplicated genes during the evolutionary processes. Except for 10 genes in the above gene pairs, other ZmGRFs did not appear in any duplication block in maize. OsGF14g of rice also did not appear in any duplication block.
The ratio of nonsynonymous-to-synonymous substitutions (Ka/Ks) was used to calculate evolutionary selection. When the ratio of a pair of sequences Ka/Ks is less than one, purifying selection is inferred; when it is equal to one, neutral drift is inferred; and when it is greater than one, positive or Darwinian selection is inferred (Juretic, 2005; Chen et al., 2014). In four gramineae plants, the selection analysis revealed that duplicated gene pairs were mostly subjected to purifying selection (Ka/Ks < 1.0). (Table 2).
Tissue specificity expression of 14-3-3s in maize sorghum, rice, and brachypodium
With reported data, the expression levels of the 14-3-3 genes in four gramineae plants from different tissues at different times were used to create a heat map (Figure 5) (Davidson et al., 2012; Sekhon et al., 2013). According to the heat map, ZmGRF1、ZmGRF5、ZmGRF13、ZmGRF21、ZmGRF22、ZmGRF23、ZmGRF24 and ZmGRF28 were expressed in all tissues. ZmGRF5 and ZmGRF13 were highly expressed in pollen, which might be closely related to reproductive development; ZmGRF22、ZmGRF23 and ZmGRF28 were highly expressed in internodes, which might be related to water and nutrient transport; In addition, ZmGRF28 was also highly expressed in 12 days immature embryos and roots at all stages, which was speculated to be related to reproduction and nutrient transport. Other ZmGRFs were low or not expressed in all tissues. OsGF14d, OsGF14f and BdGF14f were highly expressed in anthers in rice and brachypodium, which might be closely related to reproductive development; Other genes were moderately or low expressed in the tissues shown by the heat map. In sorghum, all genes were expressed in all tissues, with SbGF14b being significantly expressed in flowers, stems, and meristems; SbGF14d was highly expressed in embryos. But the repeat pairs of maize, rice and Brachypodium showed different expression patterns, suggesting that 14-3-3s have different functions in maize, rice and Brachypodium.
Figure 5 Heat maps of maize, sorghum, brachypodium, and rice 14-3-3 genes. Color scale at the right represents log2-transformed foldchange values. Red indicates high expression level; yellow indicates medium expression level; and white indicates low expression level.
Expression patterns of ZmGRF genes in mycorrhizal roots
In order to understand the expression of 14-3-3 genes induced by AM fungi in maize, RT-qPCR was performed using cDNA of maize roots inoculated with and without AM fungi as templates. We totally found 28 sequences of ZmGRF genes, but only 21 of genes were identified by qRT-PCR, while ZmGRF2, ZmGRF5, ZmGRF7, ZmGRF15, ZmGRF20, ZmGRF21 and ZmGRF23 could not be identified in maize. The expression level of the 21 ZmGRF genes were analyzed by RT-PCR (Figure 6). The primer pairs are listed in Table S5. Compared with non-AM symbiosis, there was no significant difference in the expression of ZmGRF1, ZmGRF10, ZmGRF11, ZmGRF14, ZmGRF17, ZmGRF19 and ZmGRF22 genes after AM symbiosis, indicating that these six genes were not induced by AM symbiosis. The expression of ZmGRF12 gene was significantly up-regulated (P < 0.05); The expression of ZmGRF3, ZmGRF6, ZmGRF16, ZmGRF24, ZmGRF25 and ZmGRF26 genes were significantly up-regulated (P < 0.01), The expression of ZmGRF4, ZmGRF8, ZmGRF9, ZmGRF13, ZmGRF18, ZmGRF27 and ZmGRF28 genes were significantly up-regulated (P < 0.001), indicating that these 14 genes were up-regulated by AM symbiosis. It was speculated that these genes might had a significant impact in the symbiotic pathway, and further research is needed to determine the specific mechanism.
Figure 6 Relative expressions of 14-3-3 genes in roots of wild-type. Maize was inoculated with AM symbiosis and sampled 60 days after inoculation. (p*<0.05、p**<0.01、p***<0.001).
To further verify that the above genes could be induced by AM fungi, the promoter cis-acting elements of the above 13 14-3-3 genes were analyzed (Figure 7). W-BOX (TTGACC), OSEROOTNODULE (AAAGAT) (JLQ et al., 2008), CTTC element (TCTTGT) (Lota et al., 2013), were cis-acting elements associated with AM symbiosis, and NODCON2GM (CTCTT) was a cis-acting element associated with root nodule symbiosis. The cis-acting element map was drawn for the 2000 bp upstream promoter through the RSAT website (http://rsat.eead.csic.es/plants/dna-pattern_form.cgi, accessed on 25 October 2022). It can be seen from the figure that the promoters all contain the above four cis-acting elements.
Figure 7 The transcriptional regulatory networks of 13 symbiosis-related maize 14-3-3s. These transcriptional regulatory networks were constructed by Cytoscape software. Red represents regulatory genes, and green represents symbiosis-related genes in the 14-3-3 gene family.
Interaction and regulatory networks prediction of symbiosis-related ZmGRFs
Transcriptional regulatory networks of TFs are particularly important for identifying target regulatory genes and interacting proteins of TFs (Barabási and Oltvai, 2004; Jang et al., 2015). It is still largely unknown how 14-3-3s regulate the graminaceous - AM symbiosis, particularly in maize. In this study, we selected 13 symbiotic-related maize 14-3-3 genes based on the results of RT-qPCR and constructed their co-expression networks through the plant transcriptional regulation map database and symbiotic-related RNA-seq data. According to the predicted regulatory network, ZmGRF4, ZmGRF6, ZmGRF25 and ZmGRF26 had relatively more regulatory genes, about 10-15. ZmGRF 3, ZmGRF 8, ZmGRF 9, ZmGRF 16, ZmGRF 18 and ZmGRF 28 were followed by about 5-9 regulated genes. ZmGRF13, ZmGRF 24 and ZmGRF 27 were the least, with only 1-2 regulated genes (Table S6 and Figure 7). Among the predicted regulatory genes, the following five transcription factor gene families account for most of them, namely Ethylene-responsive transcription factor, AP2/EREBP transcription factor, MADS transcription factor, MYB transcription factor and NAC transcription factor. These transcription factor families had been reported in plant-AM symbiosis. For example, the Medicago AP2/EREBP transcription factor WRI5a plays a role in controlling bidirectional nutrient exchange and periarbuscular membrane formation by regulating genes involved in fatty acid biosynthesis and phosphate uptake in arbuscular cells (Jiang et al., 2018).MYB1 regulated the transcriptional program of arbuscular degeneration during AM symbiosis (Floss et al., 2017).Expression of MtNAC969 in symbiotic root nodules was induced by nitrate treatment and antagonized by salt in roots and nodules, similar to markers of senescence (Zélicourt et al., 2012).The above results indicated that AM symbiosis was a very complex process involving multiple regulatory pathways and exerting different functions during the symbiosis process.
Discussion
14-3-3 protein have been identified and studied in many species (Cheng et al., 2018; Wang et al., 2019; Liu et al., 2020), but its evolutionary origin and phylogenetic relationship have not been fully clarified. In this study, there are 49 genes in the 14-3-3 gene family found in four plants. The number of 14-3-3 genes in maize is about 4 times that of the other three species, which may be in the process of evolution, and the function of the 14-3-3 genes is very different, which is particularly important for plant environment adaptation. Moreover, we further analyze the properties (amino acid length, molecular weight and isoelectric point) of 14-3-3 proteins (Table 1, Tables S1-S3). Physicochemical properties of 14-3-3 protein shows that the proteins of the other three species were acidic and stable, consistent with previous reports. But ZmGRFs proteins are partially acidic and partially alkaline, which is similar to Pe14-3-3s (Guo et al., 2022).Similarly, the diversity of gene structure can provide more in-depth insights for the study of structure, evolution and functional relationships (Xu et al., 2012). Except maize, the other three species all contain more than 4 introns. This indicates that evolution may lead to the diversity of this structure. Phylogenetic analysis shows that 14-3-3 proteins in four gramineae are divided into ϵ and non-ϵ groups, which is consistent with previous reports in other species (Figure 1 and Figure 2) (Rosenquist et al., 2001; Wei, 2006; Wang et al., 2019).14-3-3 proteins are highly conserved in eukaryotes (Aitken et al., 1993) and can form homologous and heterodimers (Jones et al., 1995). They can interact with two different target proteins to form complexes at the same time (Yaffe et al., 1997; Ottmann et al., 2007; Taoka et al., 2011). The analysis of conserved motifs show that most of the 14-3-3 gene families of four gramineae have 8 conserved motifs, which are the core structures that bind with many ligands (Yang et al., 2006).
Table 1 Information of 14-3-3s in maize.
The expansion of plant gene family has experienced whole-genome replication and tandem replication, which has led to plant evolution and created diversity of gene functions (Bowers et al., 2003; Doyle et al., 2008). The difference between the four gramineae plants may be caused by the whole-genome replication event (Sun et al., 2017). It can be seen from the results, a total of 6 fragment repetition events were identified in the other three plants except sorghum (Figure 3), supporting this hypothesis. We also performed a collinear analysis of the four species (Figure 4), guessing that they originated from the same ancestor. Therefore, genes may have similar functions. The Ka/Ks ratio can be used to indicate the direction of gene selection (Cannon et al., 2004; Shiu, 2004; Sun et al., 2017). In this study, 14-3-3 genes are subjected to purifying selection (Table 2), and it is speculated that there was little difference in function (Hurst, 2002). These results provide reference for the study of 14-3-3 genes.
Table 2 Evolutionary selection between 14-3-3s in maize and other three plants.
The expression patterns of 14-3-3 genes are different in different tissues of different plants (Li et al., 2016). In this study, 14-3-3 genes of the four plant species show different expressions in tissues (Figure 5), indicating that they play different roles in plant responses to different biological processes. The same has been reported in other plants (Cheng et al., 2018; Xu et al., 2021; Xia et al., 2022). For example, ZmGRF5, ZmGRF13, OsGF14d, OsGF14f and BdGF14f are highly expressed in anthers in rice and brachypodium, which may be closely related to reproductive development. It has been previously reported that AtGRF12 is highly expressed in flowers and floral organs (Tian et al., 2015). These results indicate that tissue expression pattern provides a reference direction for studying gene function.
Gramineae is an important food crop. People are very interested in AM fungal symbiosis to improve crop resistance. It has previously been reported that the expression of 14-3-3 genes induced by AM fungi can improve drought tolerance in plants (Porcel et al., 2006; Porcel et al., 2006; Li et al., 2016). In addition, the overexpression of Fm201, Ri14-3-3, and RiBMH2 genes can promote the symbiosis between plants and AM fungi (Xu et al., 2018). The results show that 21 ZmGRFs genes were significantly up-regulated after induction by AM fungi, while the remaining 7 ZmGRFs are not expressed or had no significant differences (Figure 6). Next, the cis-acting elements of the 13 14-3-3 gene promoters are analyzed, and it is found that they all containe elements related to AM fungal induction (Figure 8) (JLQ et al., 2008; Lota et al., 2013). These results further confirm that maize 14-3-3 gene family plays an important role in AM fungal symbiosis, but the detailed mechanism remains to be studied.
Figure 8 Promoter cis-acting element analysis. The promoter fragments are 2.0 kb in length. The yellow block represents the NODCON2GM motif. The green block represents the W-box motif. The red block represents the CTTC motif. The purple block represents the OSEROOTNODULE motif.
Different proteins interact to form complexes to perform various biological functions (Morsy et al., 2008; Braun et al., 2013). A co-expression network is constructed for 13 AM fungi-induced 14-3-3 genes. Three of the maize 14-3-3 genes have interacting proteins that differ from the other ten members, according to the predicted interacting proteins. The three genes are ZmGRF13, ZmGRF16 and ZmGRF24 respectively (Figure 7). It is speculated that the difference may be caused by the need to participate in other response pathways in order to adapt to environmental changes during evolution. Among the regulated transcription factors, there are 30 AP2/EREBP transcription factor family genes, 10 MYB transcription factor family genes, 7 MADS transcription factor family genes, 5 NAC transcription factor family genes, and the number of other transcription factor family genes is less than 5 (Table S6). These transcription factor families have been reported in other species to play different roles in the symbiotic process of AM fungi. This co expression network has certain guiding significance for studying the mechanism of transcription factor regulating 14-3-3 genes in the symbiotic process.
In conclusion, Plant 14-3-3 proteins can interact with proteins of various biological processes to participate in many key processes. For example, related proteins such as growth and development, material transport and biological and microbial stress. At present, the related factors have been isolated, and the molecular mechanism of 14-3-3 protein has been deeply studied (Aitken, 2006; Toker et al., 2010; Kumar et al., 2015; Camoni et al., 2018). For example, 14-3-3 proteins play a vital function in regulating stomatal movement by regulating binding partners in guard cells (Cotelle and Leonhardt, 2016). By binding to other proteins,14-3-3 proteins regulate the plant’s internal physiological processes in response to changing external environmental conditions (Denison et al., 2011). For example, overexpression of 14-3-3 protein in Arabidopsis thaliana enhances cold tolerance (Visconti et al., 2019). The 14-3-3 proteins are curcial for the response to biological stress (Yang et al., 2009). We aim to study the mechanism of improving plant tolerance to biotic and abiotic stresses by altering the expression pattern of 14-3-3 genes through AM symbiosis. This study identified 14 maize 14-3-3 genes that were up-regulated by AM symbiosis. Maize is an important food crop and cash crop. This study is of great significance to the breeding of excellent maize varieties.
Conclusion
14-3-3 proteins are involved in a variety of metabolic activities of plants. Although great progresses have been made in the function analysis of 14-3-3 genes in some plants, there are few studies on the 14-3-3 gene family in gramineae. According to the study, we identified 49 14-3-3 genes in maize, rice, sorghum, and brachypodium and demonstrated that the Gramineae 14-3-3 genes underwent large-scale replication and purifying selection. The expression analysis of 14-3-3 genes revelated that 14-3-3 genes in maize, sorghum, rice and short brachypodium may have different functions. Other species have reported that 14-3-3 genes are induced by AM fungi in response to variety of abiotic stresses. These studies provide basic information for 14-3-3 genes.In this study, the 14-3-3 gene in maize was induced by AM fungi, and 13 genes up-regulated by AM fungi were identified. Promoter analysis and co-expression network analysis were performed for them, which has certain guiding significance for studying the mechanism of transcription factors regulating 14-3-3 gene in response to various biotic and abiotic stresses during symbiosis.
Data availability statement
The data presented in the study are deposited in the NCBI repository, accession number PRJNA9184 (SAMN32395817-SAMN32395900).
Author contributions
The authors have made the following declaration about their contributions. Conceived and designed the experiments: YW, QX and HS. Performed the experiments: YW and QX. Analyzed and interpreted the data: HS, YX and YN. Performed statistical analysis: MX. Drafted the manuscript: YW, QX, HS, YN and MX. Edited and finalized the manuscript: BC and YX. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by Key research and development Program of Anhui Province (No. 202204c06020021) and National Natural Science Foundation of China (NSFC, No. U21A20235 and No. 31870415).
Acknowledgments
We thank the Key Research and Development Program of Anhui Province (No. 202204c06020021) and National Natural Science Foundation of China (NSFC, No. U21A20235 and No. 31870415) for their financial support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1117879/full#supplementary-material
References
Aitken, A. (2006). 14-3-3 proteins: a historic overview. Semin. Cancer Biol. 16 (3), 162–172. doi: 10.1016/j.semcancer.2006.03.005
Aitken, A., Collinge, D. B., Van Heusden, B. P. H., Isobe, T., Roseboom, P. H., Rosenfeld, G., et al. (1993). 14-3-3 proteins: a highly conserved, widespread family of eukaryotic proteins. Trends Biochem. Sci. 17 (12), 498–501. doi: 10.1016/0968-0004(92)90339-b
Asih, P. Y., Yuka, O., Shojiro, T., Hiroyuki, T., Ko, S. (2009). The 14-3-3 protein GF14c acts as a negative regulator of flowering in rice by interacting with the florigen Hd3a. Plant Cell Physiol. 3), 429. doi: 10.1093/pcp/pcp012
Augé, R. (2001). Water relations, drought and vesicular-arbuscular mycorrhizal symbiosis. Mycorrhiza 11 (1), 3–42. doi: 10.1007/s005720100097
Bago, B., Pfeffer, P. E., Shachar-Hill, Y. (2000). Carbon metabolism and transport in arbuscular mycorrhizas. Plant Physiol. 124 (3), 949–958. doi: 10.1104/pp.124.3.949
Bailey, T. L., Williams, N., Misleh, C., Li, W. W. (2006). MEME: Discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 1 (34), 369–373. doi: 10.1093/nar/gkl198
Barabási, Albert-László, Oltvai, Zoltán (2004). Network biology: Understanding the cell's functional organization. Nat. Rev. Genet. 5), 101–113. doi: 10.1038/nrg1272
Boston, P., Jackson, P. (1980). Purification and properties of a brain-specific protein, human 14-3-3 protein. Biochem. Soc. Trans. 8 (5), 617–618. doi: 10.1042/bst0080617
Bowers, J. E., Chapman, B. A., Rong, J., Paterson, A. H. (2003). Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature 422 (6930), 433–438. doi: 10.1038/nature01521
Braun, P., Aubourg, S., Van Leene, J., De Jaeger, G., Lurin, C. (2013). Plant protein interactomes. Annu. Rev. Plant Biol. 64, 161–187. doi: 10.1146/annurev-arplant-050312-120140
Bravo, A., Brands, M., Wewer, V., Dörmann, P., Harrison, M. J. (2017). Arbuscular mycorrhiza-specific enzymes FatM and RAM2 fine-tune lipid biosynthesis to promote development of arbuscular mycorrhiza. New Phytol. 214 (4), 1631–1645. doi: 10.1111/nph.14533
Bridges, D., Moorhead, G. (2005). 14-3-3 proteins: A number of functions for a numbered protein. Sci. s STKE 2005 (296), re10. doi: 10.1126/stke.2962005re10
Camoni, L., Visconti, S., Aducci, P., Marra, M. (2018). 14-3-3 proteins in plant hormone signaling: Doing several things at once. Front. Plant Sci. 9, 297. doi: 10.3389/fpls.2018.00297
Cannon, S. B., Mitra, A., Baumgarten, A., Young, N. D., May, G. (2004). The roles of segmental and tandem gene duplication in the evolution of large gene families in arabidopsis thaliana. BMC Plant Biol. 4 (10), 1471–2229. doi: 10.1186/1471-2229-4-10
Carlson, F. D., Hubbard, J. (1968). Physiological and biochemical aspects of nervous integration. Physiol. Zoology 42 (3), 348–349. doi: 10.1086/physzool.42.3.30155498
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13 (8), 1194–1202. doi: 10.1016/j.molp.2020.06.009
Chen, X., Chen, Z., Zhao, H., Zhao, Y., Cheng, B., Xiang, Y. (2014). Genome-wide analysis of soybean HD-zip gene family and expression profiling under salinity and drought treatments. PLos One 9 (2), e87156. doi: 10.1371/journal.pone.0087156
Cheng, C., Wang, Y., Chai, F., Li, S., Xin, H., Liang, Z. (2018). Genome-wide identification and characterization of the 14–3-3 family in vitis vinifera l. during berry Dev. cold- heat-stress response. BMC Genomics 19 (1), 579. doi: 10.1186/s12864-018-4955-8
Chevalier, D., Morris, E. R., Walker, J. C. (2009). 14-3-3 and FHA domains mediate phosphoprotein interactions. Annu. Rev. Plant Biol. 60 (1), 67–91. doi: 10.1146/annurev.arplant.59.032607.092844
Chitarra, W., Pagliarani, C., Maserti, B., Lumini, E., Siciliano, I. (2016). Insights on the impact of arbuscular mycorrhizal symbiosis on tomato tolerance to water stress. Plant Physiol. 171 (2), 1009–1023. doi: 10.1104/pp.16.00307
Cotelle, V., Leonhardt, N. (2016). 14-3-3 proteins in guard cell signaling. Front. Plant Sci. 6 (e90734), 1210. doi: 10.3389/fpls.2015.01210
Darling, D. L., Yingling, J. (2005). Wynshaw〣oris a: Role of 14–3–3 proteins in eukaryotic signaling and development. Curr. Topics Dev. Biol. 68 (4), 281–315. doi: 10.1016/S0070-2153(05)68010-6
Davidson, R. M., Gowda, M., Moghe, G., Lin, H., Vaillancourt, B., Shiu, S. H., et al. (2012). Comparative transcriptomics of three poaceae species reveals patterns of gene expression evolution. Plant J. 71 (3), 492–502. doi: 10.1111/j.1365-313X.2012.05005.x
Denison, F. C., Paul, A. L., Zupanska, A. K., Ferl, R. J. (2011). 14-3-3 proteins in plant physiology. Semin. Cell Dev. Biol. 22 (7), 720–727. doi: 10.1016/j.semcdb.2011.08.006
Doyle, J. J., Paterson, A. H., Soltis, D. E., Wendel, J. F. (2008). Evolutionary genetics of genome merger and doubling in plants. Annu. Rev. Genet. 42 (0), 443–461. doi: 10.1146/annurev.genet.42.110807.091524
Elisabeth, G., Alexandre, G., Christine, H., Ivan, I., Appel, R. D., Amos, B. (2003). ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 31 (13), 3784–3788. doi: 10.1093/nar/gkg563
Floss, D. S., Gomez, S. K., Park, H. J., MacLean, A. M., Müller, L. M., Bhattarai, K. K., et al. (2017). A transcriptional program for arbuscule degeneration during AM symbiosis is regulated by MYB1. Curr. Biol. 27 (8), 1206–1212. doi: 10.1016/j.cub.2017.03.003
Gianinazzi, S., Gollotte, A., Binet, M. N., Tuinen, D. V., Redecker, D., Wipf, D. (2010). Agroecology: The key role of arbuscular mycorrhizas in ecosystem services. Mycorrhiza 20 (8), 519–530. doi: 10.1007/s00572-010-0333-3
Guo, D., Zhu, C., Yang, K., Liu, Y., Xiao, X., Li, Z., et al. (2022). Identification of the 14-3-3 gene family in bamboo and characterization of Pe14-3-3b reveals its potential role in promoting growth. Int. J. Mol. Sci. 23 (19), 11221. doi: 10.3390/ijms231911221
Heusden, G., Griffiths, D., Ford, J. C., Chin-A-Woeng, T., Schrader, P., Carr, A. M., et al. (2010). The 14-3-3 proteins encoded by the BMH1 and BMH2 genes are essential in the yeast saccharomyces cerevisiae and can be replaced by a plant homologue. FEBS J. 229 (1), 45–53. doi: 10.1111/j.1432-1033.1995.0045l.x
Hu, B., Jin, J., Guo, A. Y., Zhang, H., Gao, G. (2014). GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics 31 (8), 1296. doi: 10.1093/bioinformatics/btu817
Huang, W., Xian, Z., Kang, X., Tang, N., Li, Z. (2015). Genome-wide identification, phylogeny and expression analysis of GRAS gene family in tomato. BMC Plant Biology,15,1(2015-08-25) 15 (1), 209. doi: 10.1186/s12870-015-0590-6
Hurst, L. D. (2002). The Ka/Ks ratio: Diagnosing the form of sequence evolution. Trends Genet. 18 (9), 486. doi: 10.1016/S0168-9525(02)02722-1
Jang, S. M., Kim, C. H., Kim, J. W., Choi, K. H. (2015). Transcriptional regulatory network of SOX4 during myoblast differentiation. Biochem. Biophys. Res. Commun. 462 (4), 365–370. doi: 10.1016/j.bbrc.2015.04.142
Jiang, Y., Wang, W., Xie, Q., Liu, N, Lixia, D. (2017). Plants transfer lipids to sustain colonization by mutualistic mycorrhizal and parasitic fungi. Science 356 (6343), 1172–1175. doi: 10.1126/science.aam9970
Jiang, Y., Xie, Q., Wang, W., Yang, J., Zhang, X., Yu, N., et al. (2018). Medicago AP2-domain transcription factor WRI5a is a master regulator of lipid biosynthesis and transfer during mycorrhizal symbiosis. Molucelar Plant. 11 (11), 16. doi: 10.1016/j.molp.2018.09.006
JLQ, A., WJZ, A., SHL, A., HW, A., DYS, A., GHX, A., et al. (2008). A HMZ: Expression analysis of light-regulated genes isolated from a full-length-enriched cDNA library of onosma paniculatum cell cultures. J. Plant Physiol. 165 (14), 1474–1482. doi: 10.1016/j.jplph.2007.11.003
John, D. (2012). Mayfield* a-LP, ferl RJ: The 14-3-3 proteins of arabidopsis regulate root growth and chloroplast development as components of the photosensory system. J. Exp. Bot. 63 (8), 3061–3070. doi: 10.1093/jxb/ers022
Jones, D. H., Ley, S., Aitken, A. (1995). Isoforms of 14-3-3 protein can form homo- and heterodimers in vivo and in vitro: implications for function as adapter proteins. FEBS Lett. 368 (1), 55–58. doi: 10.1016/0014-5793(95)00598-4
Jung, S. C., Martinez-Medina, A., Lopez-Raez, J. A., Pozo, M. J. (2012). Mycorrhiza-induced resistance and priming of plant defenses. J. Chem. Ecol. 38, 651–664. doi: 10.1007/s10886-012-0134-6
Juretic, N. (2005). The evolutionary fate of MULE-mediated duplications of host gene fragments in rice. Genome Res. 15 (9), 1292–1297. doi: 10.1101/gr.4064205
Kumar, K., Muthamilarasan, M., Bonthala, V. S., Roy, R., Prasad, M. (2015). Unraveling 14-3-3 proteins in C4 panicoids with emphasis on model plant setaria italica reveals phosphorylation-dependent subcellular localization of RS splicing factor. PLos One 10 (4), e0123236. doi: 10.1371/journal.pone.0123236
Langfelder, P., Horvath, S. (2008). WGCNA: an r package for weighted correlation network analysis. BMC Bioinf. 9, 559. doi: 10.1186/1471-2105-9-559
Lenoir, I., Fontaine, J., Sahraoui, L. H. (2016). Arbuscular mycorrhizal fungal responses to abiotic stresses: A review. Phytochemistry 123, 4–15. doi: 10.1016/j.phytochem.2016.01.002
Li, X. B., Fan, X. P., Wang, X. L., Cai, L., Yang, W. C. (2005). The cotton ACTIN1 gene is functionally expressed in fibers and participates in fiber elongation. Plant Cell 17 (3), 859–875. doi: 10.1105/tpc.104.029629
Li, M., Ren, L., Xu, B., Yang, X., Xia, Q., He, P., et al. (2016). Genome-wide identification, phylogeny, and expression analyses of the 14-3-3 family reveal their involvement in the development, ripening, and abiotic stress response in banana. Front. Plant Sci. 7 (22), 1442. doi: 10.3389/fpls.2016.01442
Li, T., Sun, Y., Ruan, Y., Xu, L., Hu, Y., Hao, Z., et al. (2016). Potential role of d-myo-inositol-3-phosphate synthase and 14-3-3 genes in the crosstalk between zea mays and rhizophagus intraradices under drought stress. Mycorrhiza 26 (8), 879–893. doi: 10.1007/s00572-016-0723-2
Liao, D., Chen, X., Chen, A., Wang, H., Liu, J., Liu, J., et al. (2014). The characterization of six auxin-induced tomato GH3 genes uncovers a member, SlGH3.4, strongly responsive to arbuscular mycorrhizal symbiosis. Plant Cell Physiol. 56 (4), 674–687. doi: 10.1093/pcp/pcu212
Liu, J., Jiang, C., Kang, L., Zhang, H., Zheng, W. (2020). Over-expression of a 14-3-3 protein from foxtail millet improves plant tolerance to salinity stress in arabidopsis thaliana. Front. Plant Sci. 11, 449. doi: 10.3389/fpls.2020.00449
Lota, F., Wegmüller, S., Buer, B., Sato, S., Brutigam, A., Hanf, B., et al. (2013). The cis-acting CTTC-P1BS module is indicative for gene function of LjVTI12, a qb-SNARE protein gene that is required for arbuscule formation in lotus japonicus. Plant J. 74 (2), 280–293. doi: 10.1111/tpj.12120
Lu, G., Delisle, A. J., Vetten, N. D., Ferl, R. J. (1992). Brain proteins in plants: an arabidopsis homolog to neurotransmitter pathway activators is part of a DNA binding complex. Proc. Natl. Acad. Sci. U.S.A. 89 (23), 11490–11494. doi: 10.1073/pnas.89.23.11490
Morsy, M., Gouthu, S., Orchard, S., Thorneycroft, D., Harper, J. F., Mittler, R., et al. (2008). Charting plant interactomes: possibilities and challenges. Trends Plant Sci. 13 (4), 183–191. doi: 10.1016/j.tplants.2008.01.006
Ottmann, Marco, Jaspert, Marcon, Schauer, Weyand, et al. (2007). : Structure of a 14-3-3 coordinated hexamer of the plant plasma membrane h+-ATPase by combining X-ray crystallography and electron cryomicroscopy. Mol. Cell 25 (3), 427–440. doi: 10.1016/j.molcel.2006.12.017
Porcel, R., Aroca, R., Azcón, R., Ruiz-Lozano, J. M. (2006). PIP aquaporin gene expression in arbuscular mycorrhizal glycinemax and lactuca sativa plants in relation to drought stress tolerance. Plant Mol. Biol. 60 (3), 389–404. doi: 10.1007/s11103-005-4210-y
Porcel, R., Aroca, R., Cano, C., Bago, A., Ruiz-Lozano, J. M. (2006). Identification of a gene from the arbuscular mycorrhizal fungus glomus intraradices encoding for a 14-3-3 protein that is up-regulated by drought stress during the AM symbiosis. Microb. Ecol. 52 (3), 575–582. doi: 10.1007/s00248-006-9015-2
Remy, W., Taylor, T. N., Hass, H., Kerp, H. (1994). Four hundred-million-year-old vesicular arbuscular mycorrhizae. Proc. Natl. Acad. Sci. 91 (25), 11841–11843. doi: 10.1073/pnas.91.25.11841
Rosenquist, M., Alsterfjord, M., Larsson, C., Sommarin, M. (2001). Data mining the arabidopsis genome reveals fifteen 14-3-3 genes. Expression is demonstrated two out five novel genes. Plant Physiol. 127 (1), 142–149. doi: 10.1104/pp.127.1.142
Schützendübel, A., Polle, A. (2002). Plant responses to abiotic stresses: heavy metal-induced oxidative stress and protection by mycorrhization. J. Exp. Bot. 53 (372), 1351–1365. doi: 10.1093/jexbot/53.372.1351
Sekhon, R. S., Roman, B., Hirsch, C. N., Myers, C. L., Springer, N. M., Robin, B. C., et al. (2013). Maize gene atlas developed by RNA sequencing and comparative evaluation of transcriptomes based on RNA sequencing and microarrays. PLos One 8 (4), e61005. doi: 10.1371/journal.pone.0061005
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498–2504. doi: 10.1101/gr.1239303
Shiu, S. H. (2004). Comparative analysis of the receptor-like kinase family in arabidopsis and rice. THE Plant Cell Online 16 (5), 1220–1234. doi: 10.1105/tpc.020834
Simon, L., Bousquet, J., Lévesque, R., Lalonde, M. (1993). Origin and diversification of endomycorrhizal fungi and coincidence with vascular land plants. Nature 363 (6424), 67–69. doi: 10.1038/363067a0
Sun, Z., Song, J., Xin, X., Xie, X., Zhao, B. (2018). Arbuscular mycorrhizal fungal 14-3-3 proteins are involved in arbuscule formation and responses to abiotic stresses during AM symbiosis. Front. Microbiol. 9, 91. doi: 10.3389/fmicb.2018.00091
Sun, H., Wu, S., Zhang, G., Jiao, C., Guo, S. (2017). Karyotype stability and unbiased fractionation in the paleo-allotetraploid cucurbita genomes. Mol. Plant (Cell Press) 10 (10), 1293–1306. doi: 10.1016/j.molp.2017.09.003
Taoka, K., Ohki, I., Tsuji, H., Furuita, K., Hayashi, K., Yanase, T., et al. (2011). et al: 14-3-3 proteins act as intracellular receptors for rice Hd3a florigen. Nature 476 (7360), 332–335. doi: 10.1038/nature10272
Tian, F., Wang, T., Xie, Y., Zhang, J., Hu, J. (2015). Genome-wide identification, classification, and expression analysis of 14-3-3 gene family in populus. PLos One 10 (4), e0123225. doi: 10.1371/journal.pone.0123225
Toker, A., Ellis, C. A., Sellers, L. A., Aitken, A. (2010). Protein kinase c inhibitor proteins. Purification sheep Brain sequence similarity to lipocortins 14-3-3 protein. FEBS J. 191 (2), 421–429. doi: 10.1111/j.1432-1033.1990.tb19138.x
Van der Heijden, M. G.A., Klironomos, J. N. (1998). Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396 (6706), 69–72. doi: 10.1038/23932
Visconti, S., D'Ambrosio, C., Fiorillo, A., Arena, S., Camoni, L. (2019). Overexpression of 14-3-3 proteins enhances cold tolerance and increases levels of stress-responsive proteins of arabidopsis plants. Plant Sci. 289 (340), 110215. doi: 10.1016/j.plantsci.2019.110215
Wang, Y., Ling, L., Jiang, Z., Tan, W., Li, W. (2019). Genome-wide identification and expression analysis of the 14-3-3 gene family in soybean ( glycine max ). PeerJ 7 (6), e7950. doi: 10.7717/peerj.7950
Wang, Y., Zhou, W., Wu, J., Xie, K., Li, X. (2022). LjAMT2;2 promotes ammonium nitrogen transport during arbuscular mycorrhizal fungi symbiosis in lotus japonicus. Int. J. Mol. Sci. 23 (17), 9522. doi: 10.3390/ijms23179522
Wei, F. X. (2006). Shi, Ming W: Expression profiling of the 14-3-3 gene family in response to salt stress and potassium and iron deficiencies in young tomato (Solanum lycopersicum) roots: Analysis by real-time RT–PCR. Ann. Bot. 98 (5), 965–974. doi: 10.1093/aob/mcl189
Wei, X., Zhang, Z., Li, Y., Wang, X., Shao, S., Chen, L., et al. (2009). Expression analysis of two novel cotton 14-3-3 genes in root development and in response to salt stress. Prog. Natural Sci. 19 (2), 173–178. doi: 10.1016/j.pnsc.2008.06.016
Xia, L., He, X., Huang, X., Yu, H., Luo, C. (2022). Genome-wide identification and expression analysis of the 14-3-3 gene family in mango (Mangifera indica l.). Int. J. Mol. Sci. 23 (23), 1593. doi: 10.3390/ijms23031593
Xu, G., Guo, C., Shan, H., Kong, H. (2012). Divergence of duplicate genes in exon-intron structure. Proc. Natl. Acad. Sci. United States America 109 (4), 1187–1192. doi: 10.1073/pnas.1109047109
Xu, M., Hu, Z., Lai, W., Liu, S., Wu, H., Zhou, Y. (2021). Comprehensive analysis of 14-3-3 family genes and their responses to cold and drought stress in cucumber. Funct. Plant Biol. 48 (12), 1264–1276. doi: 10.1071/FP21022
Xu, L., Li, T., Wu, Z., Feng, H., Yu, M., Zhang, X., et al. (2018). Arbuscular mycorrhiza enhances drought tolerance of tomato plants by regulating the 14-3-3 genes in the ABA signaling pathway. Appl. Soil Ecol. 125, 213–221. doi: 10.1016/j.apsoil.2018.01.012
Yaffe, M. B., Rittinger, K., Volinia, S., Caron, P. R., Cantley, L. C. (1997). The structural basis for 14-3-3:phosphopeptide binding specificity. Cell 91 (7), 961–971. doi: 10.1016/S0092-8674(00)80487-0
Yang, X., Wang, W., Coleman, M., Orgil, U., Feng, J., Ma, X., et al. (2009). Arabidopsis 14-3-3 lambda is a positive regulator of RPW8-mediated disease resistance. Plant J. 60 (3), 539–550. doi: 10.1111/j.1365-313X.2009.03978.x
Yang, X., Wen, H. L., Sobott, F., Papagrigoriou, E., Robinson, C. V., Grossmann, J. G., et al. (2006). Structural basis for protein–protein interactions in the 14-3-3 protein family. Proc. Natl. Acad. Sci. United States America 103 (46), 17237–17242. doi: 10.1073/pnas.0605779103
Yao, Y., Du, Y., Jiang, L., Liu, J. Y. (2007). Molecular analysis and expression patterns of the 14-3-3 gene family from oryza sativa. J. Biochem. Mol. Biol. 40 (3), 349–357. doi: 10.5483/BMBRep.2007.40.3.349
Zélicourt, A., Diet, A., Marion, J., Laffont, C., Ariel, F., Moison, M., et al. (2012). Dual involvement of a medicago truncatula NAC transcription factor in root abiotic stress response and symbiotic nodule senescence. The Plant Journal. 70, 2, 220–230. doi: 10.1111/j.1365-313X.2011.04859.x
Keywords: gramineous, 14-3-3, maize, evolution, AM symbiosis
Citation: Wang Y, Xu Q, Shan H, Ni Y, Xu M, Xu Y, Cheng B and Li X (2023) Genome-wide analysis of 14-3-3 gene family in four gramineae and its response to mycorrhizal symbiosis in maize. Front. Plant Sci. 14:1117879. doi: 10.3389/fpls.2023.1117879
Received: 07 December 2022; Accepted: 31 January 2023;
Published: 17 February 2023.
Edited by:
Renesh Bedre, Texas A&M AgriLife Research & Extension Center, Weslaco, United StatesReviewed by:
Mohit Verma, Mississippi State University, United StatesHimanshu Kumar, University of Texas Health Science Center at Houston, United States
Raviraj M. Kalunke, Donald Danforth Plant Science Center, United States
Copyright © 2023 Wang, Xu, Shan, Ni, Xu, Xu, Cheng and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Beijiu Cheng, cbj@ahau.edu.cn; Xiaoyu Li, lixiaoyu@ahau.edu.cn