Evolution of RGF/GLV/CLEL Peptide Hormones and Their Roles in Land Plant Growth and Regulation

Rooting is a key innovation during plant terrestrialization. RGFs/GLVs/CLELs are a family of secreted peptides, playing key roles in root stem cell niche maintenance and pattern formation. The origin of this peptide family is not well characterized. RGFs and their receptor genes, RGIs, were investigated comprehensively using phylogenetic and genetic analyses. We identified 203 RGF genes from 24 plant species, representing a variety of land plant lineages. We found that the RGF genes originate from land plants and expand via multiple duplication events. The lineage-specific RGF duplicates are retained due to their regulatory divergence, while a majority of RGFs experienced strong purifying selection in most land plants. Functional analysis indicated that RGFs and their receptor genes, RGIs, isolated from liverwort, tomato, and maize possess similar biological functions with their counterparts from Arabidopsis in root development. RGFs and RGIs are likely coevolved in land plants. Our studies shed light on the origin and functional conservation of this important peptide family in plant root development.     

A Novel DUF569 Gene Is a Positive Regulator of the Drought Stress Response in Arabidopsis

In the last two decades, global environmental change has increased abiotic stress on plants and severely affected crops. For example, drought stress is a serious abiotic stress that rapidly and substantially alters the morphological, physiological, and molecular responses of plants. In Arabidopsis, several drought-responsive genes have been identified; however, the underlying molecular mechanism of drought tolerance in plants remains largely unclear. Here, we report that the “domain of unknown function” novel gene DUF569 (AT1G69890) positively regulates drought stress in Arabidopsis. The Arabidopsis loss-of-function mutant atduf569 showed significant sensitivity to drought stress, i.e., severe wilting at the rosette-leaf stage after water was withheld for 3 days. Importantly, the mutant plant did not recover after rewatering, unlike wild-type (WT) plants. In addition, atduf569 plants showed significantly lower abscisic acid accumulation under optimal and drought-stress conditions, as well as significantly higher electrolyte leakage when compared with WT Col-0 plants. Spectrophotometric analyses also indicated a significantly lower accumulation of polyphenols, flavonoids, carotenoids, and chlorophylls in atduf569 mutant plants. Overall, our results suggest that novel DUF569 is a positive regulator of the response to drought in Arabidopsis.     

Functional Specialization within the EXO70 Gene Family in Arabidopsis

Localized delivery of plasma-membrane and cell-wall components is a crucial process for plant cell growth. One of the regulators of secretory-vesicle targeting is the exocyst tethering complex. The exocyst mediates first interaction between transport vesicles and the target membrane before their fusion is performed by SNARE proteins. In land plants, genes encoding the EXO70 exocyst subunit underwent an extreme proliferation with 23 paralogs present in the Arabidopsis (Arabidopsis thaliana) genome. These paralogs often acquired specialized functions during evolution. Here, we analyzed functional divergence of selected EXO70 paralogs in Arabidopsis. Performing a systematic cross-complementation analysis of exo70a1 and exo70b1 mutants, we found that EXO70A1 was functionally substituted only by its closest paralog, EXO70A2. In contrast, none of the EXO70 isoforms tested were able to substitute EXO70B1, including its closest relative, EXO70B2, pointing to a unique function of this isoform. The presented results document a high degree of functional specialization within the EXO70 gene family in land plants.     

Genome-Wide Analysis of the Apple CBL Family Reveals That Mdcbl10.1 Functions Positively in Modulating Apple Salt Tolerance

Abiotic stresses are increasingly harmful to crop yield and quality. Calcium and its signaling pathway play an important role in modulating plant stress tolerance. As specific Ca²⁺ sensors, calcineurin B-like (CBL) proteins play vital roles in plant stress response and calcium signaling. The CBL family has been identified in many plant species; however, the characterization of the CBL family and the functional study of apple MdCBL proteins in salt response have yet to be conducted in apple. In this study, 11 MdCBL genes were identified from the apple genome. The coding sequences of these MdCBL genes were cloned, and the gene structure and conserved motifs were analyzed in detail. The phylogenetic analysis indicated that these MdCBL proteins could be divided into four groups. The functional identification in Na⁺-sensitive yeast mutant showed that the overexpression of seven MdCBL genes could confer enhanced salt stress resistance in transgenic yeast. The function of MdCBL10.1 in regulating salt tolerance was also verified in cisgenic apple calli and apple plants. These results provided valuable insights for future research examining the function and mechanism of CBL proteins in regulating apple salt tolerance.     

Genome-Wide Identification and Characterization of the Soybean DEAD-Box Gene Family and Expression Response to Rhizobia

DEAD-box proteins are a large family of RNA helicases that play important roles in almost all cellular RNA processes in model plants. However, little is known about this family of proteins in crops such as soybean. Here, we identified 80 DEAD-box family genes in the Glycine max (soybean) genome. These DEAD-box genes were distributed on 19 chromosomes, and some genes were clustered together. The majority of DEAD-box family proteins were highly conserved in Arabidopsis and soybean, but Glyma.08G231300 and Glyma.14G115100 were specific to soybean. The promoters of these DEAD-box genes share cis-acting elements involved in plant responses to MeJA, salicylic acid (SA), low temperature and biotic as well as abiotic stresses; interestingly, half of the genes contain nodulation-related cis elements in their promoters. Microarray data analysis revealed that the DEAD-box genes were differentially expressed in the root and nodule. Notably, 31 genes were induced by rhizobia and/or were highly expressed in the nodule. Real-time quantitative PCR analysis validated the expression patterns of some DEAD-box genes, and among them, Glyma.08G231300 and Glyma.14G115100 were induced by rhizobia in root hair. Thus, we provide a comprehensive view of the DEAD-box family genes in soybean and highlight the crucial role of these genes in symbiotic nodulation.     

Genome-wide Identification and Characterization of FCS-Like Zinc Finger (FLZ) Family Genes in Maize (Zea mays) and Functional Analysis of ZmFLZ25 in Plant Abscisic Acid Response

FCS-like zinc finger family proteins (FLZs), a class of plant-specific scaffold of SnRK1 complex, are involved in the regulation of various aspects of plant growth and stress responses. Most information of FLZ family genes was obtained from the studies in Arabidopsis thaliana, whereas little is known about the potential functions of FLZs in crop plants. In this study, 37 maize FLZ (ZmFLZ) genes were identified to be asymmetrically distributed on 10 chromosomes and can be divided into three subfamilies. Protein interaction and subcellular localization assays demonstrated that eight typical ZmFLZs interacted and partially co-localized with ZmKIN10, the catalytic α-subunit of the SnRK1 complex in maize leaf mesophyll cells. Expression profile analysis revealed that several ZmFLZs were differentially expressed across various tissues and actively responded to diverse abiotic stresses. In addition, ectopic overexpression of ZmFLZ25 in Arabidopsis conferred hypersensitivity to exogenous abscisic acid (ABA) and triggered higher expression of ABA-induced genes, pointing to the positive regulatory role of ZmFLZ25 in plant ABA signaling, a scenario further evidenced by the interactions between ZmFLZ25 and ABA receptors. In summary, these data provide the most comprehensive information on FLZ family genes in maize, and shed light on the biological function of ZmFLZ25 in plant ABA signaling.     

Phylogenomic and Microsynteny Analysis Provides Evidence of Genome Arrangements of High-Affinity Nitrate Transporter Gene Families of Plants

Nitrate transporter 2 (NRT2) and NRT3 or nitrate-assimilation-related 2 (NAR2) proteins families form a two-component, high-affinity nitrate transport system, which is essential for the acquisition of nitrate from soils with low N availability. An extensive phylogenomic analysis across land plants for these families has not been performed. In this study, we performed a microsynteny and orthology analysis on the NRT2 and NRT3 genes families across 132 plants (Sensu lato) to decipher their evolutionary history. We identified significant differences in the number of sequences per taxonomic group and different genomic contexts within the NRT2 family that might have contributed to N acquisition by the plants. We hypothesized that the greater losses of NRT2 sequences correlate with specialized ecological adaptations, such as aquatic, epiphytic, and carnivory lifestyles. We also detected expansion on the NRT2 family in specific lineages that could be a source of key innovations for colonizing contrasting niches in N availability. Microsyntenic analysis on NRT3 family showed a deep conservation on land plants, suggesting a high evolutionary constraint to preserve their function. Our study provides novel information that could be used as guide for functional characterization of these gene families across plant lineages.     

Genome-Wide Investigation and Functional Verification of the ZIP Family Transporters in Wild Emmer Wheat

The zinc/iron-regulated transporter-like protein (ZIP) family has a crucial role in Zn homeostasis of plants. Although the ZIP genes have been systematically studied in many plant species, the significance of this family in wild emmer wheat (Triticum turgidum ssp. dicoccoides) is not yet well understood. In this study, a genome-wide investigation of ZIPs genes based on the wild emmer reference genome was conducted, and 33 TdZIP genes were identified. Protein structure analysis revealed that TdZIP proteins had 1 to 13 transmembrane (TM) domains and most of them were predicted to be located on the plasma membrane. These TdZIPs can be classified into three clades in a phylogenetic tree. They were annotated as being involved in inorganic ion transport and metabolism. Cis-acting analysis showed that several elements were involved in hormone, stresses, grain-filling, and plant development. Expression pattern analysis indicated that TdZIP genes were highly expressed in different tissues. TdZIP genes showed different expression patterns in response to Zn deficiency and that 11 genes were significantly induced in either roots or both roots and shoots of Zn-deficient plants. Yeast complementation analysis showed that TdZIP1A-3, TdZIP6B-1, TdZIP6B-2, TdZIP7A-3, and TdZIP7B-2 have the capacity to transport Zn. Overexpression of TdZIP6B-1 in rice showed increased Zn concentration in roots compared with wild-type plants. The expression levels of TdZIP6B-1 in transgenic rice were upregulated in normal Zn concentration compared to that of no Zn. This work provides a comprehensive understanding of the ZIP gene family in wild emmer wheat and paves the way for future functional analysis and genetic improvement of Zn deficiency tolerance in wheat.     

Diverse Evolution in 111 Plant Genomes Reveals Purifying and Dosage Balancing Selection Models for F-Box Genes

The F-box proteins function as substrate receptors to determine the specificity of Skp1-Cul1-F-box ubiquitin ligases. Genomic studies revealed large and diverse sizes of the F-box gene superfamily across plant species. Our previous studies suggested that the plant F-box gene superfamily is under genomic drift evolution promoted by epigenomic programming. However, how the size of the superfamily drifts across plant genomes is currently unknown. Through a large-scale genomic and phylogenetic comparison of the F-box gene superfamily covering 110 green plants and one red algal species, I discovered four distinct groups of plant F-box genes with diverse evolutionary processes. While the members in Clusters 1 and 2 are species/lineage-specific, those in Clusters 3 and 4 are present in over 46 plant genomes. Statistical modeling suggests that F-box genes from the former two groups are skewed toward fewer species and more paralogs compared to those of the latter two groups whose presence frequency and sizes in plant genomes follow a random statistical model. The enrichment of known Arabidopsis F-box genes in Clusters 3 and 4, along with comprehensive biochemical evidence showing that Arabidopsis members in Cluster 4 interact with the Arabidopsis Skp1-like 1 (ASK1), demonstrates over-representation of active F-box genes in these two groups. Collectively, I propose purifying and dosage balancing selection models to explain the lineage/species-specific duplications and expansions of F-box genes in plant genomes. The purifying selection model suggests that most, if not all, lineage/species-specific F-box genes are detrimental and are thus kept at low frequencies in plant genomes.     

Distinct Evolutionary Profiles and Functions of microRNA156 and microRNA529 in Land Plants

MicroRNA156 (miR156) and miR529 have high sequence similarity and recognize overlapping sites in the same target genes, SQUAMOSA promoter binding protein-like (SPL or SBP box) genes, making it difficult to accurately distinguish their roles in regulatory networks that affect numerous biological functions. Here, we collected data about miR156 and miR529 family members from representative land plants and performed sequence comparisons, phylogenetic analysis, small RNA sequencing, and parallel analysis of RNA ends (PARE) analysis to dissect their evolutionary and functional differences. Although miR156 and miR529 are highly similar, there are differences in their mismatch-sensitive regions, which are essential for target recognition. In land plants, miR156 precursors are conserved mainly within the hairpin region, whereas miR529 precursors are conserved outside the hairpin region, including both the 5’ and 3’ arms. Phylogenetic analysis showed that MIR156 and MIR529 evolved independently, through divergent evolutionary patterns. The two genes also exhibit different expression patterns, with MIR529 preferentially expressed in reproductive tissues and MIR156 in other tissues. PARE analysis revealed that miR156 and miR529 possess specific targets in addition to common targets in maize, pointing to functional differences between them. Based on our findings, we developed a method for the rapid identification of miR529 and miR156 family members and uncovered the evolutionary divergence of these families, providing insights into their different regulatory roles in plant growth and development.     

Genome-Wide Identification,Characterization and Expression Analysis of the Chalcone Synthase Family in Maize

Members of the chalcone synthase (CHS) family participate in the synthesis of a series of secondary metabolites in plants, fungi and bacteria. The metabolites play important roles in protecting land plants against various environmental stresses during the evolutionary process. Our research was conducted on comprehensive investigation of CHS genes in maize (Zea mays L.), including their phylogenetic relationships, gene structures, chromosomal locations and expression analysis. Fourteen CHS genes (ZmCHS01–14) were identified in the genome of maize, representing one of the largest numbers of CHS family members identified in one organism to date. The gene family was classified into four major classes (classes I–IV) based on their phylogenetic relationships. Most of them contained two exons and one intron. The 14 genes were unevenly located on six chromosomes. Two segmental duplication events were identified, which might contribute to the expansion of the maize CHS gene family to some extent. In addition, quantitative real-time PCR and microarray data analyses suggested that ZmCHS genes exhibited various expression patterns, indicating functional diversification of the ZmCHS genes. Our results will contribute to future studies of the complexity of the CHS gene family in maize and provide valuable information for the systematic analysis of the functions of the CHS gene family.