Genome-Wide Analysis of Soybean JmjC Domain-Containing Proteins Suggests Evolutionary Conservation Following Whole-Genome Duplication

A Novel CFEM Effector in Fusarium verticillioides Required for Virulence Involved in Plant Immunity Suppression and Fungal Cell Wall Integrity

Common in Fungal Extracellular Membrane (CFEM) effectors, a unique class of fungal-specific proteins, play critical roles in host-pathogen interactions. While CFEM proteins have been extensively characterized in phytopathogens, their presence and functions in Fusarium verticillioides remained unexplored. Here, we systematically identified 19 CFEM-containing proteins in F. verticillioides, among which FvCFEM12 exhibited secretory activity and plant infection-induced expression. Functional characterization revealed that FvCFEM12 suppressed Bax- and INF1-triggered cell death in Nicotiana benthamiana leaves. Furthermore, heterologous expression of FvCFEM12 in maize leaves using P. syringae strain D36E can compromise immune responses against bacterial pathogens. Deletion of FvCFEM12 impaired fungal virulence, altered hyphal morphology, and reduced cell wall stress tolerance. Interestingly, FvCFEM12 physically interacted with the maize wall-associated receptor kinase ZmWAK17ET, and targeted silencing of ZmWAK17 in maize enhanced susceptibility to F. verticillioides. Our findings revealed that FvCFEM12 is a dual-function effector that suppresses plant immunity and maintains fungal cell wall integrity, thereby orchestrating fungal pathogenicity at the host-pathogen interface.

Systematic identification and functional characterization of the CFEM proteins in fishscale bamboo rhombic-spot pathogen Neostagonosporella sichuanensis

Fungal effectors play a crucial role in the interaction between pathogenic fungi and their hosts. These interactions directly influence the invasion and spread of pathogens, and the development of diseases. Common in fungal extracellular membrane (CFEM) effectors are closely associated with the pathogenicity, cell wall stability, and pathogenic processes of pathogenic fungi. The aim of this study was to investigate the role of CFEM proteins in Neostagonosporella sichuanensis in pathogen-host interactions. We retrieved 19 proteins containing CFEM structural domains from the genome of N. sichuanensis. By systematic analysis, five NsCFEM proteins had signal peptides but lacked transmembrane structural domains, and thus were considered as potential effectors. Among them, NsCFEM1 and NsCFEM2 were successfully cloned and their functions were further investigated. The validation results show that NsCFEM1 was localized in the cell membrane and nucleus, whereas NsCFEM2 was exclusively observed in the cell membrane. Both were identified as secreted proteins. Additionally, NsCFEM1 inhibited Bax-induced programmed cell death in Nicotiana benthamiana, whereas NsCFEM2 did not induce or inhibit this response. NsCFEM1 was implicated as a virulence factor that contributes to fungal growth, development, stress response, and pathogenicity. NsCFEM2 was implicated in maintenance of cell wall stability. This study lays a foundation for elucidating the role of CFEM proteins in the pathogen of fishscale bamboo rhombic-spot caused by N. sichuanensis. In particular, the functional studies of NsCFEM1 and NsCFEM2 revealed their potential roles in the interaction between N. sichuanensis and the host Phyllostachys heteroclada.

Genome-wide identification and expression analysis reveal the role of histone methyltransferase and demethylase genes in heat stress response in potato (Solanum tuberosum L.)

BACKGROUND

Potato (Solanum tuberosum L.), the third most important non-cereal crop, is sensitive to high temperature. Histone modifications have been known to regulate various abiotic stress responses. However, the role of histone methyltransferases and demethylases remain unexplored in potato under heat stress.

METHODS

Potato genome database was used for genome-wide analysis of StPRMT and StHDMA gene families, which were further characterized by analyzing gene structure, conserved motif, domain organization, sub-cellular localization, promoter region and phylogenetic relationships. Additionally, expression profiling under high-temperature stress in leaf and stolon tissue of heat contrasting potato genotypes was done to study their role in response to high temperature stress.

RESULTS

The genome-wide analysis led to identification of nine StPRMT and eleven StHDMA genes. Structural analysis, including conserved motifs, exon/intron structure and phylogenetic relationships classified StPRMT and StHDMA gene families into two classes viz. Class I and Class II. A variety of cis-regulatory elements were explored in the promoter region associated with light, developmental, hormonal and stress responses. Prediction of sub-cellular localization of StPRMT proteins revealed their occurrence in nucleus and cytoplasm, whereas StHDMA proteins were observed in different sub-cellular compartments. Furthermore, expression profiling of StPRMT and StHDMA gene family members revealed genes responding to heat stress. Heat-inducible expression of StPRMT1, StPRMT3, StPRMT4 and StPRMT5 in leaf and stolon tissues of HS and HT cultivar indicated them as probable candidates for enhancing thermotolerance in potato. However, StHDMAs responded dynamically in leaf and stolon tissue of heat contrasting genotypes under high temperature.

CONCLUSION

The current study presents a detailed analysis of histone modifiers in potato and indicates their role as an important epigenetic regulators modulating heat tolerance.

GENERAL SIGNIFICANCE

Understanding epigenetic mechanisms underlying heat tolerance in potato will contribute towards breeding of thermotolerant potato varieties.

Identification of Quantitative Trait Nucleotides and Development of Diagnostic Markers for Nine Fatty Acids in the Peanut

The cultivated peanut (Arachis hypogaea L.) is an important oilseed crop worldwide, and fatty acid composition is a major determinant of peanut oil quality. In the present study, we conducted a genome-wide association study (GWAS) for nine fatty acid traits using the whole genome sequences of 160 representative Chinese peanut landraces and identified 6-1195 significant SNPs for different fatty acid contents. Particularly for oleic acid and linoleic acid, two peak SNP clusters on Arahy.09 and Arahy.19 were found to contain the majority of the significant SNPs associated with these two fatty acids. Additionally, a significant proportion of the candidate genes identified on Arahy.09 overlap with those identified in early studies, among which three candidate genes are of special interest. One possesses a significant missense SNP and encodes a known candidate gene FAD2A. The second gene is the gene closest to the most significant SNP for linoleic acid. It codes for an MYB protein that has been demonstrated to impact fatty acid biosynthesis in Arabidopsis. The third gene harbors a missense SNP and encodes a JmjC domain-containing protein. The significant phenotypic difference in the oleic acid/linoleic acid between the genotypes at the first and third candidate genes was further confirmed with PARMS analysis. In addition, we have also identified different candidate genes (i.e., Arahy.ZV39IJ, Arahy.F9E3EA, Arahy.X9ZZC1, and Arahy.Z0ELT9) for the remaining fatty acids. Our findings can help us gain a better understanding of the genetic foundation of peanut fatty acid contents and may hold great potential for enhancing peanut quality in the future.

The Bax inhibitor GmBI-1α interacts with a Nod factor receptor and plays a dual role in the legume-rhizobia symbiosis

The gene networks surrounding Nod factor receptors that govern the symbiotic process between legumes and rhizobia remain largely unexplored. Here, we identify 13 novel GmNFR1α-associated proteins by yeast two-hybrid screening, and describe a potential interacting protein, GmBI-1α. GmBI-1α had the highest positive correlation with GmNFR1α in a co-expression network analysis, and its expression at the mRNA level in roots was enhanced by rhizobial infection. Moreover, GmBI-1α-GmNFR1α interaction was shown to occur in vitro and in vivo. The GmBI-1α protein was localized to multiple subcellular locations, including the endoplasmic reticulum and plasma membrane. Overexpression of GmBI-1α increased the nodule number in transgenic hairy roots or transgenic soybean, whereas down-regulation of GmBI-1α transcripts by RNA interference reduced the nodule number. In addition, the nodules in GmBI-1α-overexpressing plants became smaller in size and infected area with reduced nitrogenase activity. In GmBI-1α-overexpressing transgenic soybean, the elevated GmBI-1α also promoted plant growth and suppressed the expression of defense signaling-related genes. Infection thread analysis of GmBI-1α-overexpressing plants showed that GmBI-1α promoted rhizobial infection. Collectively, our findings support a GmNFR1α-associated protein in the Nod factor signaling pathway and shed new light on the regulatory mechanism of GmNFR1α in rhizobial symbiosis.

Flagellin and mannitol modulate callose biosynthesis and deposition in soybean seedlings

Callose is a polymer deposited on the cell wall and is necessary for plant growth and development. Callose is synthesized by genes from the glucan synthase-like family (GSL) and dynamically responds to various types of stress. Callose can inhibit pathogenic infection, in the case of biotic stresses, and maintain cell turgor and stiffen the plant cell wall in abiotic stresses. Here, we report the identification of 23 GSL genes (GmGSL) in the soybean genome. We performed phylogenetic analyses, gene structure prediction, duplication patterns, and expression profiles on several RNA-Seq libraries. Our analyses show that WGD/Segmental duplication contributed to expanding this gene family in soybean. Next, we analyzed the callose responses in soybean under abiotic and biotic stresses. The data show that callose is induced by both osmotic stress and flagellin 22 (flg22) and is related to the activity of β-1,3-glucanases. By using RT-qPCR, we evaluated the expression of GSL genes during the treatment of soybean roots with mannitol and flg22. The GmGSL23 gene was upregulated in seedlings treated with osmotic stress or flg22, showing the essential role of this gene in the soybean defense response to pathogenic organisms and osmotic stress. Our results provide an important understanding of the role of callose deposition and regulation of GSL genes in response to osmotic stress and flg22 infection in soybean seedlings.

Genome-wide identification and expression analysis of the cotton patatin-related phospholipase A genes and response to stress tolerance

Patatin-related phospholipase A genes were involved in the response of Gossypium hirsutum to drought and salt tolerance. pPLA (patatin-related phospholipase A) is a key enzyme that catalyzes the initial step of lipid hydrolysis, which is involved in biological processes, such as drought, salt stress, and freezing injury. However, a comprehensive analysis of the pPLA gene family in cotton, especially the role of pPLA in the response to drought and salt tolerance, has not been reported so far. A total of 33 pPLA genes were identified in this study using a genome-wide search approach, and phylogenetic analysis classified these genes into three groups. These genes are unevenly distributed on the 26 chromosomes of cotton, and most of them contain a few introns. The results of the collinear analysis showed that G. hirsutum contained 1-5 copies of each pPLA gene found in G. arboreum and G. raimondii. The subcellular localization analysis of Gh_D08G061200 showed that the protein was localized in the nucleus. In addition, analysis of published upland cotton transcriptome data revealed that six GhPLA genes were expressed in various tissues and organs. Two genes (Gh_A04G142100.1 and Gh_D04G181000.1) were highly expressed in all tissues under normal conditions, showing the expression characteristics of housekeeping genes. Under different drought and salt tolerance stresses, we detected four genes with different expression levels. This study helps to clarify the role of pPLA in the response to drought and salt tolerance.

Identification and analysis of differentially expressed trihelix genes in maize (Zea mays) under abiotic stresses

Background

Trihelix transcription factors play important roles in triggering plant growth and imparting tolerance against biotic and abiotic stresses. However, a systematical analysis of the trihelix transcription factor family under heat and drought stresses in maize has not been reported.

Methods

PlantTFDB and TBtools were employed to identify the trihelix domain-containing genes in the maize genome. The heat-regulated transcriptome data for maize were obtained from NCBI to screen differentially expressed ZmTHs genes through statistical analysis. The basic protein sequences, chromosomal localization, and subcellular localization were analyzed using Maize GDB, Expasy, SOMPA, TBtools, and Plant-mPLoc. The conserved motifs, evolutionary relationships, and cis-elements, were analyzed by MEME, MEGA7.0 and PlantCARE software, respectively. The tissue expression patterns of ZmTHs and their expression profiles under heat and drought stress were detected using quantitative real-time PCR (qRT-PCR).

Results

A total of 44 trihelix family members were discovered, and members were distributed over 10 chromosomes in the maize genome. A total of 11 genes were identified that were regulated by heat stress; these were unevenly distributed on chromosomes 1, 2, 4, 5, and 10. ZmTHs encoded a total of 16 proteins, all of which were located in the nucleus; however, ZmTH04.1 was also distributed in the chloroplast. The protein length varied from 206 to 725 amino acids; the molecular weight ranged from 22.63 to 76.40 kD; and the theoretical isoelectric point (pI) ranged from 5.24 to 11.2. The protein's secondary structures were mainly found to be random coils and α-helices, with fewer instances of elongation chains and β-rotations. Phylogenetic relationship analysis showed that these can be divided into five sub-groups. The conserved domain of ZmTHs was GT1 or MyB_DNA-Bind_4. The protein and gene structure of ZmTHs differed greatly among the subfamilies, while the structures within the subfamilies were similar. The promoter of ZmTHs contained abundant tissue-specific expression cis-acting elements and abiotic stress response elements. qRT-PCR analysis showed that ZmTHs expression levels were significantly different in different tissues. Furthermore, the expression of ZmTH08 was dramatically up-regulated by heat stress, while the expression of ZmTH03, ZmTH04, ZmTH05, ZmTH06, ZmTH07, ZmTH09, ZmTH10, and ZmTH11 were down-regulated by heat stress. Upon PEG-simulated drought stress, ZmTH06 was significantly up-regulated, while ZmTH01 and ZmTH07 were down-regulated.

Conclusions

We performed a genome-wide, systematic identification and analysis of differentially expressed trihelix genes under heat and drought stresses in maize.

Systematic identification and functional characterization of the CFEM proteins in poplar fungus Marssonina brunnea

Proteins containing Common in Fungal Extracellular Membrane (CFEM) domains uniquely exist in fungi and play significant roles in their whole life history. In this study, a total of 11 MbCFEM proteins were identified from Marssonina brunnea f. sp. multigermtubi (MULT), a hemibiotrophic pathogenic fungus on poplars that causes severe leaf diseases. Phylogenic analysis showed that the 11 proteins (MbCFEM1-11) were divided into three clades based on the trans-membrane domain and the CFEM domain. Sequence alignment and WebLogo analysis of CFEM domains verified the amino acids conservatism therein. All of them possess eight cysteines except MbCFEM4 and MbCFEM11, which lack two cysteines each. Six MbCFEM proteins with a signal peptide and without trans-membrane domain were considered as candidate effectors for further functional analysis. Three-dimensional (3D) models of their CFEM domains presented a helical-basket structure homologous to the crucial virulence factor Csa2 of Candida albicans. Afterward, four (MbCFEM1, 6, 8, and 9) out of six candidate effectors were successfully cloned and a yeast signal sequence trap (YSST) assay confirmed their secretion activity. Pathogen challenge assays demonstrated that the transient expression of four candidate MbCFEM effectors in Nicotiana benthamiana promoted Fusarium proliferatum infection, respectively. In an N. benthamiana heterogeneous expression system, MbCFEM1, MbCFEM6, and MbCFEM9 appeared to suppress both BAX/INF1-triggered PCD, whereas MbCFEM8 could only defeat BAX-triggered PCD. Additionally, subcellular localization analysis indicated that the four candidate MbCFEM effectors accumulate in the cell membrane, nucleus, chloroplast, and cytosolic bodies. These results demonstrate that MbCFEM1, MbCFEM6, MbCFEM8, and MbCFEM9 are effectors of M. brunnea and provide valuable targets for further dissection of the molecular mechanisms underlying the poplar-M. brunnea interaction.

Genome-wide identification and expression analysis of JmjC domain-containing genes in grape under MTA treatment

Histone demethylases containing the JmjC domain play an extremely important role in maintaining the homeostasis of histone methylation and are closely related to plant growth and development. Currently, the JmjC domain-containing proteins have been reported in many species; however, they have not been systematically studied in grapes. In this paper, 21 VviJMJ gene family members were identified from the whole grape genome, and the VviJMJ genes were classified into five subfamilies: KDM3, KDM4, KDM5, JMJD6, and JMJ-only based on the phylogenetic relationship and structural features of Arabidopsis and grape. After that, the conserved sites of VviJMJ genes were revealed by protein sequence analysis. In addition, chromosomal localization and gene structure analysis revealed the heterogeneous distribution of VviJMJ genes on grape chromosomes and the structural features of VviJMJ genes, respectively. Analysis of promoter cis-acting elements demonstrated numerous hormone, light, and stress response elements in the promoter region of the VviJMJ genes. Subsequently, the grape fruit was treated with MTA (an H3K4 methylation inhibitor), which significantly resulted in the early ripening of grape fruits. The qRT-PCR analysis showed that VviJMJ genes (except VviJMJ13c) had different expression patterns during grape fruit development. The expression of VviJMJ genes in the treatment group was significantly higher than that in the control group. The results indicate that VviJMJ genes are closely related to grape fruit ripening.

Insights to the superoxide dismutase genes and its roles in Hevea brasiliensis under abiotic stress

The superoxide dismutase (SOD) protein significantly influences the development and growth of plants and their reaction to abiotic stresses. However, little is known about the characteristics of rubber tree SOD genes and their expression changes under abiotic stresses. The present study recognized 11 SOD genes in the rubber tree genome, including 7 Cu/ZnSODs, 2 MnSODs, and 2 FeSODs. Except for HbFSD1, SODs were scattered on five chromosomes. The phylogenetic analysis of SOD proteins in rubber trees and a few other plants demonstrated that the SOD proteins contained three major subgroups. Moreover, the genes belonging to the same clade contained similar gene structures, which confirmed their classification further. The extension of the SOD gene family in the rubber tree was mainly induced by the segmental duplication events. The cis-acting components analysis showed that HbSODs were utilized in many biological procedures. The transcriptomics data indicated that the phosphorylation of the C-terminal domain of RNA polymerase II might control the cold response genes through the CBF pathway and activate the SOD system to respond to cold stress. The qRT-PCR results showed that the expression of HbCSD1 was significantly downregulated under drought and salt stresses, which might dominate the adaption capability to different stresses. Additionally, salt promoted the expression levels of HbMSD1 and HbMSD2, exhibiting their indispensable role in the salinity reaction. The study results will provide a theoretical basis for deep research on HbSODs in rubber trees.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03328-7.

Genome-wide identification of the jumonji C domain- containing histone demethylase gene family in wheat and their expression analysis under drought stress

Methylation and demethylation of histone play a crucial role in regulating chromatin formation and gene expression. The jumonji C (JmjC) domain-containing proteins are demethylases that are involved in regulating epigenetic modification in plants. In our study, the JmjC genes in Triticum aestivum L., Triticum turgidum L., Triticum dicoccoides L., Triticum urartu L., and Aegilops tauschii L. were identified. Phylogenetic relationship and colinearity analysis revealed that the wheat JmjC genes were conserved in A, B, and D subgenomes during evolution. Cis-acting elements analysis showed that elements related to stress response, hormone response, and light response were found in wheat JmjC genes. The expression of JmjC genes was affected by tissue types and developmental stages, and members of the same subfamily tended to have similar expression patterns in wheat. They also showed a unique expression pattern in root during PEG (Polyethylene glycol) treatment. In conclusion, comprehensive analysis indicated that three members (Tr-1A-JMJ2, Tr-1B-JMJ2, and Tr-1D-JMJ2) might be regulated by several hormones and function in the early stages of drought stress, while eight members (Tr-1B-JMJ3, Tr-4B-JMJ1, Tr-7A-JMJ1, etc.) displayed a significantly high expression after 24 h of PEG treatment, indicating a role in the later stages of drought stress. This research presents the first genome-wide study of the JmjC family in wheat, and lays the foundation for promoting the study of their functional characterization in wheat drought resistance.

Evolutionary History and Functional Diversification of the JmjC Domain-Containing Histone Demethylase Gene Family in Plants

Histone demethylases containing JumonjiC (JmjC) domains regulate gene transcription and chromatin structure by changing the methylation status of lysine residues and play an important role in plant growth and development. In this study, a total of 332 JmjC family genes were identified from 21 different plant species. The evolutionary analysis results showed that the JmjC gene was detected in each species, that is, the gene has already appeared in algae. The phylogenetic analysis showed that the KDM3/JHDM2 subfamily genes may have appeared when plants transitioned from water to land, but were lost in lycophytes (Selaginella moellendorffii). During the evolutionary process, some subfamily genes may have been lost in individual species. According to the analysis of the conserved domains, all of the plant JmjC genes contained a typical JmjC domain, which was highly conserved during plant evolution. The analysis of cis-acting elements showed that the promoter region of the JmjC gene was rich in phytohormones and biotic and abiotic stress-related elements. The transcriptome data analysis and protein interaction analyses showed that JmjC genes play an important role in plant growth and development. The results clarified the evolutionary history of JmjC family genes in plants and lay the foundation for the analysis of the biological functions of JmjC family genes.

Genome-wide identification of nitrate transporter 2 (NRT2) gene family and functional analysis of MeNRT2.2 in cassava (Manihot esculenta Crantz)

Nitrate transporter 2 (NRT2) proteins play an important role in nitrate uptake and utilization in plants. The NRT2 family has been identified and functionally characterized in many plants. However, no systematic identification of NRT2 family members has been reported in cassava (Manihot esculenta Crantz). In this study, six MeNRT2 genes were identified from cassava genome and named as MeNRT2.1-2.6 according to their chromosomal locations. Phylogenetic tree showed that NRT2 proteins were divided into four main subgroups, which was further supported by their gene structure and conserved motifs. All six MeNRT2 genes are randomly distributed on 4 chromosomes (LG8, LG11, LG13, and LG17), two tandem duplicated genes (MeNRT2.3/MeNRT2.4) and a pair of segmental duplicated gene (MeNRT2.1/MeNRT2.2) was detected. Subsequently, expression profiles of MeNRT2 genes in eight different tissues and in response to nitrate deficient treatment were analyzed. The results showed that the MeNRT2 genes had differential expression patterns. All of MeNRT2 genes induced by nitrate deficiency, of them the MeNRT2.2 had the highest expression level after treatment. Arabidopis transformed with MeNRT2.2 gene showed higher fresh weight than wild type plants in response to N starvation, suggesting that MeNRT2.2 play important role in adapting to low nitrogen. Taken together, our results provide the reference for further analyses of the molecular functions of the MeNRT2 gene family, but also some candidate genes for developing nitrogen efficient crops.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid Evolution

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

Comprehensive Analysis of Jumonji Domain C Family from Citrus grandis and Expression Profilings in the Exocarps of “Huajuhong” (Citrus grandis “Tomentosa”) during Various Development Stages

Citrus grandis “Tomentosa” (“Huajuhong”) is a famous Traditional Chinese Medicine. In this study, a total of 18 jumonji C (JMJC) domain-containing proteins were identified from C. grandis. The 18 CgJMJCs were unevenly located on six chromosomes of C. grandis. Phylogenetic analysis revealed that they could be classified into five groups, namely KDM3, KDM4, KDM5, JMJC, and JMJD6. The domain structures and motif architectures in the five groups were diversified. Cis-acting elements on the promoters of 18 CgJMJC genes were also investigated, and the abscisic acid-responsive element (ABRE) was distributed on 15 CgJMJC genes. Furthermore, the expression profiles of 18 CgJMJCs members in the exocarps of three varieties of “Huajuhong”, for different developmental stages, were examined. The results were validated by quantitative real-time PCR (qRT-PCR). The present study provides a comprehensive characterization of JMJC domain-containing proteins in C. grandis and their expression patterns in the exocarps of C. grandis “Tomentosa” for three varieties with various development stages.

Genome-Wide Analysis of WRKY Gene Family and the Dynamic Responses of Key WRKY Genes Involved in Ostrinia furnacalis Attack in Zea mays

WRKY transcription factors comprise one of the largest gene families and serve as key regulators of plant defenses against herbivore attack. However, studies related to the roles of WRKY genes in response to herbivory are limited in maize. In this study, a total of 128 putative maize WRKY genes (ZmWRKYs) were identified from the new maize genome (v4). These genes were divided into seven subgroups (groups I, IIa–e, and III) based on phylogenomic analysis, with distinct motif compositions in each subgroup. Syntenic analysis revealed that 72 (56.3%) of the genes were derived from either segmental or tandem duplication events (69 and 3, respectively), suggesting a pivotal role of segmental duplication in the expansion of the ZmWRKY family. Importantly, transcriptional regulation prediction showed that six key WRKY genes contribute to four major defense-related pathways: L-phenylalanine biosynthesis II and flavonoid, benzoxazinoid, and jasmonic acid (JA) biosynthesis. These key WRKY genes were strongly induced in commercial maize (Jingke968) infested with the Asian corn borer, Ostrinia furnacalis, for 0, 2, 4, 12 and 24 h in the field, and their expression levels were highly correlated with predicted target genes, suggesting that these genes have important functions in the response to O. furnacalis. Our results provide a comprehensive understanding of the WRKY gene family based on the new assembly of the maize genome and lay the foundation for further studies into functional characteristics of ZmWRKY genes in commercial maize defenses against O. furnacalis in the field.

"The PLCP gene family of grapevine (Vitis vinifera L.): characterization and differential expression in response to Plasmopara Viticola"

BACKGROUND

Papain-like cysteine proteases (PLCPs), a large group of cysteine proteases, are structurally related to papain. The members belonging to PLCPs family contribute to plant immunity, senescence, and defense responses in plants. The PLCP gene family has been identified in Arabidopsis, rice, soybean, and cotton. However, no systematic analysis of PLCP genes has been undertaken in grapevine. Since Plasmopara viticola as a destructive pathogen could affect immunity of grapes in the field, we considered that the members belonged to PLCPs family could play a crucial role in defensive mechanisms or programmed cell death. We aimed to evaluate the role of PLCPs in 2 different varieties of grapevines and compared the changes of their expressions with the transcriptional data in response to P. viticola.

RESULTS

In this study, 23 grapevine PLCP (VvPLCP) genes were identified by comprehensive bioinformatics analysis. Subsequently, the chromosomal localizations, gene structure, conserved domains, phylogenetic relationship, gene duplication, and cis-acting elements were analyzed. Numerous cis-acting elements related to plant development, hormone, and stress responses were identified in the promoter of the VvPLCP genes. Phylogenetic analysis grouped the VvPLCP genes into nine subgroups. The transcription of VvPLCP in different inoculation time points and varieties indicated that VvPLCP may have vital functions in grapevine defense against Plasmopara viticola. According to transcriptome data and qPCR analysis, we observed the increasing expression levels of VvRD21-1 at 72 h after inoculation in resistant variety, inferring that it was related to grape downy mildew resistance. Meanwhile, 3 genes including VvXBCP1, VvSAG12-1, and VvALP1 showed higher expression at 24 h after pathogen inoculation in the susceptible variety and might be related to the downy mildew phenotype. We nominated these four genes to function during hypersensitive response (HR) process, inferring that these genes could be associated with downy mildew resistance in grapes.

CONCLUSIONS

Our results provide the reference for functional studies of PLCP gene family, and highlight its functions in grapevine defense against P. viticola. The results help us to better understand the complexity of the PLCP gene family in plant immunity and provide valuable information for future functional characterization of specific genes in grapevine.

Genome-wide identification, classification, and expression analysis of the JmjC domain-containing histone demethylase gene family in birch

BACKGROUND

Histone methylation occurs primarily on lysine residues and requires a set of enzymes capable of reading, writing, and erasing to control its establishment and deletion, which is essential for maintaining chromatin structure and gene expression. Histone methylation and demethylation are contributed to plant growth and development, and are involved in adapting to environmental stresses. The JmjC domain-containing proteins are extensively studied for their function in histone lysine demethylation in plants, and play a critical role in sustaining histone methylation homeostasis.

RESULTS

In this study, a total of 21 JmjC domain-containing histone demethylase proteins (JHDMs) in birch were identified and classified into five subfamilies based on structural characteristics and phylogenetic relationships among Arabidopsis, rice, maize, and birch. Although the BpJMJ genes displayed significant schematic variation, their distribution on the chromosomes is relatively uniform. Additionally, the BpJMJ genes in birch have never experienced a tandem-duplication event proved by WGD analysis and were remaining underwent purifying selection (Ka/Ks < < 1). A typical JmjC domain was found in all BpJMJ genes, some of which have other essential domains for their functions. In the promoter regions of BpJMJ genes, cis-acting elements associated with hormone and abiotic stress responses were overrepresented. Under abiotic stresses, the transcriptome profile reveals two contrasting expression patterns within 21 BpJMJ genes. Furthermore, it was established that most BpJMJ genes had higher expression in young tissues under normal conditions, with BpJMJ06/16 having the highest expression in germinating seeds and participating in the regulation of BpGA3ox1/2 gene expression. Eventually, BpJMJ genes were found to directly interact with genes involved in the "intracellular membrane" in respond to cold stress.

CONCLUSIONS

The present study will provide a foundation for future experiments on histone demethylases in birch and a theoretical basis for epigenetic research on growth and development in response to abiotic stresses.

Linkage mapping and QTL analysis of flowering time in faba bean

Flowering time marks the transition from vegetative to reproductive growth and is key for optimal yield in any crop. The molecular mechanisms controlling this trait have been extensively studied in model plants such as Arabidopsis thaliana and rice. While knowledge on the molecular regulation of this trait is rapidly increasing in sequenced galegoid legume crops, understanding in faba bean remains limited. Here we exploited translational genomics from model legume crops to identify and fine map QTLs linked to flowering time in faba bean. Among the 31 candidate genes relevant for flowering control in A. thaliana and Cicer arietinum assayed, 25 could be mapped in a segregating faba bean RIL population. While most of the genes showed conserved synteny among related legume species, none of them co-localized with the 9 significant QTL regions identified. The FT gene, previously implicated in the control of flowering time in numerous members of the temperate legume clade, mapped close to the most relevant stable and conserved QTL in chromosome V. Interestingly, QTL analysis suggests an important role of epigenetic modifications in faba bean flowering control. The new QTLs and candidate genes assayed here provide a robust framework for further genetic studies and will contribute to the elucidation of the molecular mechanisms controlling this trait.

Linkage mapping and QTL analysis of flowering time in faba bean

Flowering time marks the transition from vegetative to reproductive growth and is key for optimal yield in any crop. The molecular mechanisms controlling this trait have been extensively studied in model plants such as Arabidopsis thaliana and rice. While knowledge on the molecular regulation of this trait is rapidly increasing in sequenced galegoid legume crops, understanding in faba bean remains limited. Here we exploited translational genomics from model legume crops to identify and fine map QTLs linked to flowering time in faba bean. Among the 31 candidate genes relevant for flowering control in A. thaliana and Cicer arietinum assayed, 25 could be mapped in a segregating faba bean RIL population. While most of the genes showed conserved synteny among related legume species, none of them co-localized with the 9 significant QTL regions identified. The FT gene, previously implicated in the control of flowering time in numerous members of the temperate legume clade, mapped close to the most relevant stable and conserved QTL in chromosome V. Interestingly, QTL analysis suggests an important role of epigenetic modifications in faba bean flowering control. The new QTLs and candidate genes assayed here provide a robust framework for further genetic studies and will contribute to the elucidation of the molecular mechanisms controlling this trait.

Genome-wide analysis of JMJ-C histone demethylase family involved in salt-tolerance in Gossypium hirsutum L

The jumonji C (JMJ-C) domain-containing protein is a histone demethylase and is involved in plant stress. However, the function of the JMJ-C gene family in cotton is still not confirmed. Herein, 25, 26, 52, and 53 members belonging to the JMJ-C gene family were identified in Gossypium raimondii, Gossypium arboreum, Gossypium hirsutum, and Gossypium barbadense, respectively. Based on phylogenetic relationships and conserved domains, the JMJ-C genes were categorized into five subfamilies, KDM3, KDM4, KDM5, JMJC, and JMJD6. The chromosomal location, gene structure, motif compositions, and cis-elements have been displayed. The collinear investigation showed that whole-genome duplication event is the mainly power to drive JMJ-C gene family expansion. Transcriptome and qRT-PCR analysis revealed that eight GhJMJs were induced by salt and PEG treatment. Further assays confirmed that GhJMJ34/40 greatly improved salt and osmotic tolerance in Saccharomyces cerevisiae. These results help clarify JMJ-C protein functions in preparation for further study.

Genome-wide identification, classification and expression analysis of the Hsf and Hsp70 gene families in maize

Heat shock factors (Hsfs) and heat shock proteins (Hsps) play a critical role in the molecular mechanisms such as plant development and defense against abiotic. As an important food crop, maize is vulnerable to adverse environment such as heat stress and water logging, which leads to a decline in yield and quality. To date, very little is known regarding the structure and function of Hsf and Hsp genes in maize. Although some Hsf and Hsp genes have been characterized in maize, analysis of the entire Hsf and Hsp70 gene families were not completed following Maize (B73) Genome Sequencing Project. Therefore, studying their molecular mechanism and revealing their biological function in plant stress resistance process will contribute to reveal important theoretical significance and application value for improving corn yield and quality. In this study, we have identified 25 ZmHsf and 22 ZmHsp70 genes in maize. The structural characteristics and phylogenetic relationships of the Hsf and Hsp70 gene families of Arabidopsis thaliana, rice and maize were compared. The final 25 ZmHsf proteins and 22 ZmHsp70 proteins were divided into three and four subfamilies, respectively. In addition, chromosomal localization indicated that the ZmHsf and ZmHsp70 genes were unevenly distributed on the chromosome, and the gene structure map revealed the characteristics of their structures. Finally, transcriptome analysis indicated that most of the ZmHsf and ZmHsp70 genes showed different expression patterns at different developmental stages of maize. Further, by semi-quantitative RT-PCR and quantitative real-time PCR analysis, all 25 ZmHsf and 22 ZmHsp70 genes were confirmed to respond to heat stress treatment, indicating that they have potential effects in heat stress response. The analyses performed by combining co-expression network with protein-protein interaction network among the members of the Hsf and Hsp70 gene families in maize further enabled us to recognize components involved in the regulatory network associated with hsfs and hsp70s complex. The predicted subcellular location revealed that maize Hsp70 proteins exhibited a various subcellular distribution, which may be associated with functional diversification in heat stress response. Taken together, our study provides comprehensive information on the members of Hsf and Hsp70 gene families and will help in elucidating their exact function in maize.

Characterization and Stress Response of the JmjC Domain-Containing Histone Demethylase Gene Family in the Allotetraploid Cotton Species Gossypium hirsutum

Histone modification is an important epigenetic modification that controls gene transcriptional regulation in eukaryotes. Histone methylation is accomplished by histone methyltransferase and can occur on two amino acid residues, arginine and lysine. JumonjiC (JmjC) domain-containing histone demethylase regulates gene transcription and chromatin structure by changing the methylation state of the lysine residue site and plays an important role in plant growth and development. In this study, we carried out genome-wide identification and comprehensive analysis of JmjC genes in the allotetraploid cotton species Gossypium hirsutum. In total, 50 JmjC genes were identified and in G. hirsutum, and 25 JmjC genes were identified in its two diploid progenitors, G. arboreum and G. raimondii, respectively. Phylogenetic analysis divided these JmjC genes into five subfamilies. A collinearity analysis of the two subgenomes of G. hirsutum and the genomes of G. arboreum and G. raimondii uncovered a one-to-one relationship between homologous genes of the JmjC gene family. Most homologs in the JmjC gene family between A and D subgenomes of G. hirsutum have similar exon-intron structures, which indicated that JmjC family genes were conserved after the polyploidization. All G. hirsutumJmjC genes were found to have a typical JmjC domain, and some genes also possess other special domains important for their function. Analysis of promoter regions revealed that cis-acting elements, such as those related to hormone and abiotic stress response, were enriched in G. hirsutum JmjC genes. According to a reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis, most G. hirsutumJmjC genes had high abundance expression at developmental stages of fibers, suggesting that they might participate in cotton fiber development. In addition, some G. hirsutumJmjC genes were found to have different degrees of response to cold or osmotic stress, thus indicating their potential role in these types of abiotic stress response. Our results provide useful information for understanding the evolutionary history and biological function of JmjC genes in cotton.

Genome-wide survey of soybean papain-like cysteine proteases and their expression analysis in root nodule symbiosis

BACKGROUND

Plant papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes and play important roles in root nodule symbiosis (RNS), while the whole-genome studies of PLCP family genes in legume are quite limited, and the roles of Glycine max PLCPs (GmPLCPs) in nodulation, nodule development and senescence are not fully understood.

RESULTS

In the present study, we identified 97 GmPLCPs and performed a genome-wide survey to explore the expansion of soybean PLCP family genes and their relationships to RNS. Nineteen paralogous pairs of genomic segments, consisting of 77 GmPLCPs, formed by whole-genome duplication (WGD) events were identified, showing a high degree of complexity in duplication. Phylogenetic analysis among different species showed that the lineage differentiation of GmPLCPs occurred after family expansion, and large tandem repeat segment were specifically in soybean. The expression patterns of GmPLCPs in symbiosis-related tissues and nodules identified RNS-related GmPLCPs and provided insights into their putative symbiotic functions in soybean. The symbiotic function analyses showed that a RNS-related GmPLCP gene (Glyma.04G190700) really participate in nodulation and nodule development.

CONCLUSIONS

Our findings improved our understanding of the functional diversity of legume PLCP family genes, and provided insights into the putative roles of the legume PLCPs in nodulation, nodule development and senescence.

Genome-wide identification and functional analysis of JmjC domain-containing genes in flower development of Rosa chinensis

We genome-wide identified 28 JmjC domain-containing genes, further spatio-temporal expression profiling and genetic analysis defined them as epigenetic regulators in flowering initiation of Rosa chinensis. The JmjC domain-containing histone demethylases play critical roles in maintaining homeostasis of histone methylations, thus are vital for plant growth and development. Genome-wide identification of the JmjC domain-containing genes have been reported in several species, however, no systematic study has been performed in rose plants. In this paper, we identified 28 JmjC domain-containing genes from the newly published genome database of Rosa chinensis. The JmjC domain-containing proteins in R. chinensis were divided into seven groups, KDM3 was the largest group with 13 members, and JmjC domain-only A and KDM5B were the smallest clades both with only one member. Although all the JmjC domain proteins having a conserved JmjC domain, the gene and protein structure experienced differentiation and specification during the evolution, especially in KDM3 clade, one gene (RcJMJ40) was found carrying site deletions for cofactors Fe (II) and α-KG binding which were crucial for demethylase activities, three genes (RcJMJ41, RcJMJ43 and RcJMJ44) had no intron while two of them had tandem JmjC domains. Spatial expression pattern analysis of these JmjC domain-containing genes in different tissues showed most of them were highly expressed in reproductive tissues such as floral meristem and closed flowers than vegetative tissues, demonstrating their important functions in developmental switch from vegetative to reproductive growth of roses. Temporal expression profiling indicated majority of JmjC domain-containing genes from R. chinensis fluctuated along with floral bud differentiation and development, further proving their essential roles in flower organogenesis. VIGS induced silencing of RcJMJ12 led to delayed flowering time, and decreased the expression levels of flowering integrator such as RcFT, RcSOC1, RcFUL, RcLFY and RcAP1, therefore providing the genetic evidence of RcJMJ12 in flowering initiation. Collectively, spatio-temporal expression profiling and genetic analysis defined the JmjC domain-containing genes as important epigenetic regulators in flower development of R. chinensis.

Catalpol Promotes the Proliferation and Differentiation of Osteoblasts Induced by High Glucose by Inhibiting KDM7A

INTRODUCTION

The protective effect of catalpol on diabetic osteoporosis (DOP) and its mechanism remain unclear. This study aimed to explore whether catalpol enhanced the proliferation and differentiation of MC3T3 cells induced by high glucose by inhibiting the expression of KDM7A.

METHODS

MC3T3 cells were induced by high glucose (HG) and treated with different concentrations of catalpol. The proliferation and mineralization abilities of MC3T3-E1 cells were determined by CCK-8 assay and Alizarin Red Staining, respectively. The expression of differentiation-related osteogenic proteins, KDM7A and related proteins of Wnt/β-catenin signaling pathway was analyzed by Western blot analysis. The alkaline phosphatase (ALP) activity was detected by ALP assay kits.

RESULTS

MC3T3-E1 cells induced by high glucose showed decreased proliferation and mineralization abilities and decreased ALP activity, which were all reversed by the treatment of catalpol. High glucose induction inhibited the expression of KDM7A, Total-β-catenin, Nuclear-β-catenin and p-GSK3β, which was reversed by the treatment of catalpol. And KDM7A interference up-regulated the expression of Total-β-catenin, Nuclear-β-catenin and p-GSK3β, which was down-regulated by KDM7A overexpression. Furthermore, the proliferation and mineralization abilities and ALP activity were improved when treated with KDM7A interference and decreased when treated with KDM7A overexpression. However, SKL2001 could improve the proliferation and mineralization abilities and ALP activity of MC3T3-E1 cells.

DISCUSSION

Catalpol promotes the proliferation and differentiation of osteoblasts induced by high glucose by regulating the Wnt/β-catenin signaling pathway through KDM7A.

Identification, Evolution, and Expression Profiling of Histone Lysine Methylation Moderators in Brassica rapa

Histone modifications, such as methylation and demethylation, are vital for regulating chromatin structure, thus affecting its expression patterns. The objective of this study is to understand the phylogenetic relationships, genomic organization, diversification of motif modules, gene duplications, co-regulatory network analysis, and expression dynamics of histone lysine methyltransferases and histone demethylase in Brassica rapa. We identified 60 SET (HKMTases), 53 JmjC, and 4 LSD (HDMases) genes in B. rapa. The domain composition analysis subcategorized them into seven and nine subgroups, respectively. Duplication analysis for paralogous pairs of SET and JmjC (eight and nine pairs, respectively) exhibited variation. Interestingly, three pairs of SET exhibited Ka/Ks > 1.00 values, signifying positive selection, whereas the remaining underwent purifying selection with values less than 1.00. Furthermore, RT-PCR validation analysis and RNA-sequence data acquired on six different tissues (i.e., leaf, stem, callus, silique, flower, and root) revealed dynamic expression patterns. This comprehensive study on the abundance, classification, co-regulatory network analysis, gene duplication, and responses to heat and cold stress of SET and JmjC provides insights into the structure and diversification of these family members in B. rapa. This study will be helpful to reveal functions of these putative SET and JmjC genes in B. rapa.

Genome-Wide Identification, Expression Pattern Analysis and Evolution of the Ces/Csl Gene Superfamily in Pineapple (Ananas comosus)

The cellulose synthase (Ces) and cellulose synthase-like (Csl) gene families belonging to the cellulose synthase gene superfamily, are responsible for the biosynthesis of cellulose and hemicellulose of the plant cell wall, and play critical roles in plant development, growth and evolution. However, the Ces/Csl gene family remains to be characterized in pineapple, a highly valued and delicious tropical fruit. Here, we carried out genome-wide study and identified a total of seven Ces genes and 25 Csl genes in pineapple. Genomic features and phylogeny analysis of Ces/Csl genes were carried out, including phylogenetic tree, chromosomal locations, gene structures, and conserved motifs identification. In addition, we identified 32 pineapple AcoCes/Csl genes with 31 Arabidopsis AtCes/Csl genes as orthologs by the syntenic and phylogenetic approaches. Furthermore, a RNA-seq investigation exhibited the expression profile of several AcoCes/Csl genes in various tissues and multiple developmental stages. Collectively, we provided comprehensive information of the evolution and function of pineapple Ces/Csl gene superfamily, which would be useful for screening out and characterization of the putative genes responsible for tissue development in pineapple. The present study laid the foundation for future functional characterization of Ces/Csl genes in pineapple.

Genome-wide identification, classification and expression analysis of the JmjC domain-containing histone demethylase gene family in maize

Background

Histone methylation mainly occurs on the lysine residues and plays a crucial role during flowering and stress responses of plants, through changing the methylation status or ratio of lysine residues. Histone lysine residues of plants can arise in three forms of methylation (single, double and triple) and the corresponding demethylation can also ensue on certain occasions, by which the plants can accommodate the homeostasis of histone methylation by means of lysine methyltransferase and demethylase. The JmjC domain-containing proteins, an important family of histone lysine demethylases, play a vital role in maintaining homeostasis of histone methylation in vivo.

Results

In this study, we have identified 19 JmjC domain-containing histone demethylase (JHDM) proteins in maize. Based on structural characteristics and a comparison of phylogenetic relationships of JHDM gene families from Arabidopsis, rice and maize, all 19 JHDM proteins in maize were categorized into three different subfamilies. Furthermore, chromosome location and schematic structure revealed an unevenly distribution on chromosomes and structure features of ZmJMJ genes in maize, respectively. Eventually, the 19 ZmJMJ genes displayed different expression patterns at diverse developmental stages of maize based on transcriptome analysis. Further, quantitative real-time PCR analysis showed that all 19 ZmJMJ genes were responsive to heat stress treatment, suggesting their potential roles in heat stress response.

Conclusions

Overall, our study will serve to present an important theoretical basis for future functional verification of JHDM genes to further unravel the mechanisms of epigenetic regulation in plants.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5633-1) contains supplementary material, which is available to authorized users.

Histone H3K9 demethylase JMJ25 epigenetically modulates anthocyanin biosynthesis in poplar

Anthocyanins are involved in several aspects of development and defence in poplar (Populus spp.). Although, over the past decades, significant progress has been made in uncovering these anthocyanin biosynthetic and regulatory mechanisms, the fundamental understanding of the epigenetic regulation in this pathway is still largely unclear. Here, we isolated a histone H3K9 demethylase gene JMJ25 from Populus and characterized its role in anthocyanin biosynthesis by genetic and biochemical approaches. JMJ25 was induced by continuous dark treatment. Overexpression of JMJ25 led to downregulated expression of anthocyanin biosynthetic genes in transgenic poplar, resulting in a significant reduction in anthocyanin content. ChIP-qPCR assays showed that JMJ25 could directly associate with MYB182 chromatin and dynamically demethylate at H3K9me2. Furthermore, JMJ25 also affected the DNA methylation levels of MYB182. By contrast, knockout of JMJ25 by CRISPR/Cas9 resulted in ectopic anthocyanin accumulation under dark condition and increased expression of anthocyanin biosynthetic genes. Our results support a model in which JMJ25 directly affects MYB182 expression by altering the histone methylation status of its chromatin and DNA methylation, resulting in repression of anthocyanin accumulation. This study uncovered an epigenetic mechanism that modulates anthocyanin biosynthesis in poplar.

Identification of jumonjiC domain containing gene family among the Oryza species and their expression analysis in FL478, a salt tolerant rice genotype

The jumonji (JMJ)-C domain containing proteins belong to histone demethylases family with the ability to demethylate the tri-methylated histone residues. They act as chromatin regulators to regulate many physiological functions in plants. The present study deals with the characterization of JMJ-C gene family members in wild as well as cultivated rice species and their expression analysis in salt tolerant rice genotype, FL478. The genome wide study identified 151 members belonging to JMJ-C gene family in 11 different Oryza species. We also studied their structure, genomic location, gene duplication events, phylogenetic relationship, in silico expression analysis and identified cis elements in their promoters. We also found a few JMJ-C gene family members in rice which underwent duplication before the whole genome duplication event of the rice. The qRT-PCR based expression profiling revealed that out of the total 15 rice JMJ-C members, two were highly expressed in the flag leaf stage of FL478 under salt treatment. These two candidate JMJ-C members were also found to render salinity tolerance when over-expressed in yeast cells. Thus, the present study helps in further structural as well as functional characterization of JMJ-C genes under salinity stress in Oryza species.

Genome-Wide Identification and Classification of Soybean C2H2 Zinc Finger Proteins and Their Expression Analysis in Legume-Rhizobium Symbiosis

Root nodule symbiosis (RNS) is one of the most productive and economical systems for nitrogen fixation, and previous studies have shown that several nodule-specific C2H2-zinc finger proteins (ZFPs) play important roles in symbiosis establishment and nodule function. However, C2H2-ZFPs are the most widespread ZFPs in eukaryotes, and a great variation of structure and function exist among the family members. It remains largely unclear whether or not special types of C2H2-ZF genes participate in symbiosis, especially in soybean. In the present study, we performed a genome-wide survey of soybean C2H2-ZF genes, and 321 soybean C2H2-ZF genes were identified and classified into 11 clearly distinguishable subsets (Gm-t1-SF, Gm-t2-SF, Gm-1i-Q-SF, Gm-1i-M-SF, Gm-1i-Z-SF, Gm-1i-D-SF, Gm-2i-Q-SF, Gm-2i-M-SF, Gm-2i-Mix-SF, Gm-3i-SF, and Gm-4i-SF) based on the arrangements, numbers, and types of C2H2-ZF domains. Phylogenetic and gene ontology analyses were carried out to assess the conserved sequence and GO function among these subsets, and the results showed that the classification of soybean C2H2-ZFPs was reasonable. The expression profile of soybean C2H2-ZFPs in multiple tissues showed that nearly half of soybean C2H2-ZFPs within different subsets had expressions in nodules, including a clustering branch consisting of 11 Gm-1i-Q-SF genes specifically expressed in symbiotic-relative tissues. RNA-Seq was used to identify symbiosis-related soybean C2H2-ZFPs, and the expression pattern of the soybean C2H2-ZFPs in roots and nodules at different development stages showed that soybean C2H2-ZFPs mainly played roles in nodule development or nodule function rather than nodulation signal transduction, and nearly half of these genes had high expressions and/or different expression patterns during soybean nodule development, especially for the six clustering branches of genes consisting of different subsets of C2H2-ZFPs. Furthermore, the selected symbiosis-related soybean C2H2-ZFPs might function in legume-rhizobium symbiosis through regulating or interacting with other key proteins. Taken together, our findings provided useful information for the study on classification and conservative function of C2H2-ZFPs, and offered solid evidence for investigation of rhizobium symbiosis-related C2H2-ZFPs in soybean or other legumes.