Genome-Wide Identification and Characterization of Polygalacturonase Gene Family in Maize (Zea mays L.)

Polygalacturonase (PG, EC 3.2.1.15) is a crucial enzyme for pectin degradation and is involved in various developmental processes such as fruit ripening, pollen development, cell expansion, and organ abscission. However, information on the PG gene family in the maize (Zea mays L.) genome and the specific members involved in maize anther development are still lacking. In this study, we identified 55 PG family genes from the maize genome and further characterized their evolutionary relationship and expression patterns. Phylogenetic analysis revealed that ZmPGs are grouped into six Clades, and gene structures of the same Clade are highly conserved, suggesting their functional conservation. The ZmPGs are randomly distributed across maize chromosomes, and collinearity analysis showed that many ZmPGs might be derived from tandem duplications and segmental duplications, and these genes are under purifying selection. Furthermore, gene expression analysis provided insights into possible functional divergence among ZmPGs. Based on the RNA-seq data analysis, we found that many ZmPGs are expressed in various tissues while 18 ZmPGs are highly expressed in maize anther, and their detailed expression profiles in different anther developmental stages were further investigated by using RT-qPCR analysis. These results provide valuable information for further functional characterization and application of the ZmPGs in maize.     

Genome-Wide Identification and Characterization of Soybean GmLOR Gene Family and Expression Analysis in Response to Abiotic Stresses

The LOR (LURP-one related) family genes encode proteins containing a conserved LOR domain. Several members of the LOR family genes are required for defense against Hyaloperonospora parasitica (Hpa) in Arabidopsis. However, there are few reports of LOR genes in response to abiotic stresses in plants. In this study, a genome-wide survey and expression levels in response to abiotic stresses of 36 LOR genes from Glycine max were conducted. The results indicated that the GmLOR gene family was divided into eight subgroups, distributed on 14 chromosomes. A majority of members contained three extremely conservative motifs. There were four pairs of tandem duplicated GmLORs and nineteen pairs of segmental duplicated genes identified, which led to the expansion of the number of GmLOR genes. The expansion patterns of the GmLOR family were mainly segmental duplication. A heatmap of soybean LOR family genes showed that 36 GmLOR genes exhibited various expression patterns in different tissues. The cis-acting elements in promoter regions of GmLORs include abiotic stress-responsive elements, such as dehydration-responsive elements and drought-inducible elements. Real-time quantitative PCR was used to detect the expression level of GmLOR genes, and most of them were expressed in the leaf or root except that GmLOR6 was induced by osmotic and salt stresses. Moreover, GmLOR4/10/14/19 were significantly upregulated after PEG and salt treatments, indicating important roles in the improvement of plant tolerance to abiotic stress. Overall, our study provides a foundation for future investigations of GmLOR gene functions in soybean.     

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.     

Genome-Wide Identification,Characterization and Expression Analysis of the JAZ Gene Family in Resistance to Gray Leaf Spots in Tomato

The plant disease resistance system involves a very complex regulatory network in which jasmonates play a key role in response to external biotic or abiotic stresses. As inhibitors of the jasmonic acid (JA) signaling pathway, JASMONATE ZIM domain (JAZ) proteins have been identified in many plant species, and their functions are gradually being clarified. In this study, 26 JAZ genes were identified in tomato. The physical and chemical properties, predicted subcellular localization, gene structure, cis-acting elements, and interspecies collinearity of 26 SlJAZ genes were subsequently analyzed. RNA-seq data combined with qRT-PCR analysis data showed that the expression of most SlJAZ genes were induced in response to Stemphylium lycopersici, methyl jasmonate (MeJA) and salicylic acid (SA). Tobacco rattle virus RNA2-based VIGS vector (TRV2)-SlJAZ25 plants were more resistant to tomato gray leaf spots than TRV2-00 plants. Therefore, we speculated that SlJAZ25 played a negative regulatory role in tomato resistance to gray leaf spots. Based on combining the results of previous studies and those of our experiments, we speculated that SlJAZ25 might be closely related to JA and SA hormone regulation. SlJAZ25 interacted with SlJAR1, SlCOI1, SlMYC2, and other resistance-related genes to form a regulatory network, and these genes played an important role in the regulation of tomato gray leaf spots. The subcellular localization results showed that the SlJAZ25 gene was located in the nucleus. Overall, this study is the first to identify and analyze JAZ family genes in tomato via bioinformatics approaches, clarifying the regulatory role of SlJAZ25 genes in tomato resistance to gray leaf spots and providing new ideas for improving plant disease resistance.     

Gibberellin Oxidase Gene Family in L. chinense: Genome-Wide Identification and Gene Expression Analysis

GAox is a key enzyme for the transformation of gibberellins, and belongs to the 2-ketoglutarate dependent dioxygenase gene family (2ODD). However, a systematic analysis of GAox in the angiosperm L. chinense has not yet been reported. Here, we identified all LcGAox gene family members in L. chinense, which were classified into the three subgroups of GA20ox, C19GA2ox, and C20GA2ox. Comparison of the gene structure, conserve motifs, phylogenetic relationships, and syntenic relationships of gibberellin oxidase gene families in different species indicated that the gene functional differences may be due to the partial deletion of their domains during evolution. Furthermore, evidence for purifying selection was detected between orthologous GAox genes in rice, grape, Arabidopsis, and L. chinense. Analysis of the codon usage patterns showed that mutation pressure and natural selection might have induced codon usage bias in angiosperms; however, the LcGAox genes in mosses, lycophytes, and ambarella plants exhibited no obvious codon usage preference. These results suggested that the gibberellin oxidase genes were more primitive. The gene expression pattern was analyzed in different organs subjected to multiple abiotic stresses, including GA, abscisic acid (ABA), and chlormequat (CCC) treatment, by RNA-seq and qRT-PCR, and the stress- and phytohormone-responsive cis-elements were counted. The results showed that the synthesis and decomposition of GA were regulated by different LcGAox genes in the vegetative and reproductive organs of L. chinense, and only LcGA2ox1,4, and 7 responded to the NaCl, polyethylene glycol, 4 °C, GA, ABA, and CCC treatment in the roots, stems, and leaves of seedlings at different time periods, revealing the potential role of LcGAox in stress resistance.     

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.     

Genome-Wide Identification of GDSL-Type Esterase/Lipase Gene Family in Dasypyrum villosum L. Reveals That DvGELP53 Is Related to BSMV Infection

GDSL-type esterase/lipase proteins (GELPs) characterized by a conserved GDSL motif at their N-terminus belong to the lipid hydrolysis enzyme superfamily. In plants, GELPs play an important role in plant growth, development and stress response. The studies of the identification and characterization of the GELP gene family in Triticeae have not been reported. In this study, 193 DvGELPs were identified in Dasypyrum villosum and classified into 11 groups (clade A–K) by means of phylogenetic analysis. Most DvGELPs contain only one GDSL domain, only four DvGELPs contain other domains besides the GDSL domain. Gene structure analysis indicated 35.2% DvGELP genes have four introns and five exons. In the promoter regions of the identified DvGELPs, we detected 4502 putative cis-elements, which were associated with plant hormones, plant growth, environmental stress and light responsiveness. Expression profiling revealed 36, 44 and 17 DvGELPs were highly expressed in the spike, the root and the grain, respectively. Further investigation of a root-specific expressing GELP, DvGELP53, indicated it was induced by a variety of biotic and abiotic stresses. The knockdown of DvGELP53 inhibited long-distance movement of BSMV in the tissue of D. villosum. This research provides a genome-wide glimpse of the D. villosum GELP genes and hints at the participation of DvGELP53 in the interaction between virus and plants.     

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.     

Identification and Characterization of Salt- and Drought-Responsive AQP Family Genes in Medicago sativa L.

Aquaporins (AQP) are distributed ubiquitously in plants, and they play important roles in multiple aspects of plant growth and development, as well as in plant resistance to various environmental stresses. In this study, 43 MsAQP genes were identified in the forage crop Medicago sativa. All the MsAQP proteins were clustered into four subfamilies based on sequence similarity and phylogenetic relationship, including 17 TIPs, 14 NIPs, 9 PIPs and 3 SIPs. Analyses of gene structure and conserved domains indicated that the majority of the deduced MsAQP proteins contained the signature transmembrane domains and the NPA motifs. Analyses on cis-acting elements in the promoter region of MsAQP genes revealed the presence of multiple and diverse stress-responsive and hormone-responsive cis-acting elements. In addition, by analyzing the available and comprehensive gene expression data of M. truncatula, we screened ten representative MtAQP genes that were responsive to NaCl or drought stress. By analyzing the sequence similarity and phylogenetic relationship, we finally identified the corresponding ten salt- or drought-responsive AQP genes in M. sativa, including three MsTIPs, three MsPIPs and four MsNIPs. The qPCRs showed that the relative expression levels of these ten selected MsAQP genes responded differently to NaCl or drought treatment in M. sativa. Gene expression patterns showed that most MsAQP genes were preferentially expressed in roots or in leaves, which may reflect their tissue-specific functions associated with development. Our results lay an important foundation for the future characterization of the functions of MsAQP genes, and provide candidate genes for stress resistance improvement through genetic breeding in M. sativa.     

Transcriptome-Wide Identification and Functional Characterization of CIPK Gene Family Members in Actinidia valvata under Salt Stress

Fruit plants are severely constrained by salt stress in the soil due to their sessile nature. Ca²⁺ sensors, which are known as CBL-interacting protein kinases (CIPKs), transmit abiotic stress signals to plants. Therefore, it is imperative to investigate the molecular regulatory role of CIPKs underlying salt stress tolerance in kiwifruit. In the current study, we have identified 42 CIPK genes from Actinidia. valvata (A.valvata). All the AvCIPKs were divided into four different phylogenetic groups. Moreover, these genes showed different conserved motifs. The expression pattern analysis showed that AvCIPK11 was specifically highly expressed under salt stress. The overexpression of AvCIPK11 in ‘Hongyang’ (a salt sensitive commercial cultivar from Actinidia chinensis) enhanced salt tolerance by maintaining K⁺/Na⁺ homeostasis in the leaf and positively improving the activity of POD. In addition, the salt-related genes AcCBL1 and AcNHX1 had higher expression in overexpression lines. Collectively, our study suggested that AvCIPK11 is involved in the positive regulation of salt tolerance in kiwifruit.     

Genome-Wide Identification and Expression Profile of OSCA Gene Family Members in Triticum aestivum L.

Hyperosmolality and various other stimuli can trigger an increase in cytoplasmic-free calcium concentration ([Ca²⁺]_cyt). Members of the Arabidopsis thaliana (L.) reduced hyperosmolality-gated calcium-permeable channels (OSCA) gene family are reported to be involved in sensing extracellular changes to trigger hyperosmolality-induced [Ca²⁺]_cyt increases and controlling stomatal closure during immune signaling. Wheat (Triticum aestivum L.) is a very important food crop, but there are few studies of its OSCA gene family members. In this study, 42 OSCA members were identified in the wheat genome, and phylogenetic analysis can divide them into four clades. The members of each clade have similar gene structures, conserved motifs, and domains. TaOSCA genes were predicted to be regulated by cis-acting elements such as STRE, MBS, DRE1, ABRE, etc. Quantitative PCR results showed that they have different expression patterns in different tissues. The expression profiles of 15 selected TaOSCAs were examined after PEG (polyethylene glycol), NaCl, and ABA (abscisic acid) treatment. All 15 TaOSCA members responded to PEG treatment, while TaOSCA12/-39 responded simultaneously to PEG and ABA. This study informs research into the biological function and evolution of TaOSCA and lays the foundation for the breeding and genetic improvement of wheat.